The metagenomic datasets of the microbial DNA from tropical bed bugs (Cimex hemipterus) after feeding on human blood were presented. Next-generation sequencing of the community DNA was carried out on an Illumina Miseq platform and the raw fastq files were analyzed using QIIME (version 1.9.1). The metagenome of three samples comprised of 108,198 sequences representing 44,646,263 bps with a mean length of 412.63 bps. The sequence data is accessible at the NCBI SRA under the bioproject number PRJNA600667. Community analysis showed Proteobacteria was the most abundance (more than 99%) microbial community that present in the guts of fully fed tropical bed bugs.
Here, we present the first complete mitochondrial genome of Malayan Gaur (Bos gaurus hubbacki) inferred using next-generation sequencing. The mitogenome is 16,367 bp in length with the structural organization of a typical bovine mitochondrial arrangement comprising 13 protein-coding genes, 21 tRNAs, and 2 rRNAs. No internal stop codon was found in the protein-coding genes. Phylogenetic tree analysis revealed that Malayan gaur is more closely related to Burmese banteng instead of gaur.
Known for its valuable agarwood, Aquilaria malaccensis Lam. is an evergreen tropical forest tree species endemic to the Indo-malesian region. Indiscriminate damaging and harvesting of the trees in the wild have resulted in it being listed in the CITES Appendix II for controlled trade and in the IUCN Red List as 'Vulnerable (VU)'. In this study, the complete chloroplast genome of A. malaccensis was assembled using data from high-throughput Illumina sequencing. The chloroplast genome was 174,832 bp in size, which included two inverted repeat regions of 42,091 bp each, separated by a large single copy region of 87,302 bp and a small single copy region of 3,348 bp. A total of 139 genes were predicted, including 39 tRNA, 8 rRNA, and 92 protein-coding genes. Phylogenetic analysis placed A. malaccensis within the family Thymelaeaceae. The chloroplast genome sequence of A. malaccensis offers a useful resource for future studies on the taxonomy and conservation of the threatened Aquilaria trees.
The complete mitochondrial genome of the cavity-nesting honeybee Apis cerana from Sabah on Borneo Island was analyzed using next-generation sequencing. The mitochondrial genome of A. cerana was a circular molecule of 15,884 bp and was similar to that of the other cavity-nesting honeybee species. The average AT content in the A. cerana mitochondrial genome was 84.4%. It was predicted to contain 13 protein-coding, 22 tRNA, and two rRNA genes, along with one A + T-rich control region.
Analysis of repetitive DNA sequence content and divergence among the repetitive functional classes is a well-accepted approach for estimation of inter- and intra-generic differences in plant genomes. Among these elements, microsatellites, or Simple Sequence Repeats (SSRs), have been widely demonstrated as powerful genetic markers for species and varieties discrimination. We present PlantFuncSSRs platform having more than 364 plant species with more than 2 million functional SSRs. They are provided with detailed annotations for easy functional browsing of SSRs and with information on primer pairs and associated functional domains. PlantFuncSSRs can be leveraged to identify functional-based genic variability among the species of interest, which might be of particular interest in developing functional markers in plants. This comprehensive on-line portal unifies mining of SSRs from first and next generation sequencing datasets, corresponding primer pairs and associated in-depth functional annotation such as gene ontology annotation, gene interactions and its identification from reference protein databases. PlantFuncSSRs is freely accessible at: http://www.bioinfocabd.upo.es/plantssr.
MAIN CONCLUSION: Karyotyping using high-density genome-wide SNP markers identified various chromosomal aberrations in oil palm (Elaeis guineensis Jacq.) with supporting evidence from the 2C DNA content measurements (determined using FCM) and chromosome counts. Oil palm produces a quarter of the world's total vegetable oil. In line with its global importance, an initiative to sequence the oil palm genome was carried out successfully, producing huge amounts of sequence information, allowing SNP discovery. High-capacity SNP genotyping platforms have been widely used for marker-trait association studies in oil palm. Besides genotyping, a SNP array is also an attractive tool for understanding aberrations in chromosome inheritance. Exploiting this, the present study utilized chromosome-wide SNP allelic distributions to determine the ploidy composition of over 1,000 oil palms from a commercial F1 family, including 197 derived from twin-embryo seeds. Our method consisted of an inspection of the allelic intensity ratio using SNP markers. For palms with a shifted or abnormal distribution ratio, the SNP allelic frequencies were plotted along the pseudo-chromosomes. This method proved to be efficient in identifying whole genome duplication (triploids) and aneuploidy. We also detected several loss of heterozygosity regions which may indicate small chromosomal deletions and/or inheritance of identical by descent regions from both parents. The SNP analysis was validated by flow cytometry and chromosome counts. The triploids were all derived from twin-embryo seeds. This is the first report on the efficiency and reliability of SNP array data for karyotyping oil palm chromosomes, as an alternative to the conventional cytogenetic technique. Information on the ploidy composition and chromosomal structural variation can help to better understand the genetic makeup of samples and lead to a more robust interpretation of the genomic data in marker-trait association analyses.
Salmonella Typhi is a human restricted pathogen with a significant number of individuals as asymptomatic carriers of the bacterium. Salmonella infection can be effectively controlled if a reliable method for identification of these carriers is developed. In this context, the availability of whole genomes of carrier strains through high- throughput sequencing and further downstream analysis by comparative genomics approaches is very promising. Herein we describe the genome sequence of a Salmonella Typhi isolate representing an asymptomatic carrier individual during a prolonged outbreak of typhoid fever in Kelantan, Malaysia. Putative genomic coordinates relevant in pathogenesis and persistence of this carrier strain are identified and discussed.
The Tapirus indicus, also known as Malayan tapir, has been listed as a rapidly declining animal species in the past decades, along with being declared and categorized as an endangered species by the International Union for Conservation of Nature (IUCN) 2016. This tapir species is geographically distributed across several countries in Southeast Asia such as Peninsular Malaysia, Indonesia (Sumatra), South Thailand, and Myanmar. Amongst these countries, the Peninsula Malaysia forest is recorded to contain the highest number of Malayan tapir population. Unfortunately, in the past decades, the population of Malayan tapirs has declined swiftly due to serious deforestation, habitat fragmentation, and heavy vehicle accidents during road crossings at forest routes. Concerned by this predicament, the Department of Wildlife and National Parks (DWNP) Peninsular Malaysia collaborated with a few local universities to conduct various studies aimed at increasing the population number of tapirs in Malaysia. Several studies were conducted with the aim of enhancing the well-being of tapirs in captivity. Veterinarians face problems when it comes to selecting healthy and suitable tapirs for breeding programs at conservation centers. Conventional molecular methods using high-throughput sequencing provides a solution in determining the health condition of Malayan tapirs using the Next-Generation Sequencing (NGS) technology. Unaware by most, gut microbiome plays an important role in determining the health condition of an organism by various aspects: (1) digestion control; (2) benefiting the immune system; and (3) playing a role as a "second brain." Commensal gut bacterial communities (microbiomes) are predicted to influence organism health and disease. Imbalance of unhealthy and healthy microbes in the gut may contribute to weight gain, high blood sugar, high cholesterol, and other disorders. In infancy, neonatal gut microbiomes are colonized with maternal and environmental flora, and mature toward a stable composition in two to three years. Interactions between the microorganism communities and the host allow for the establishment of microbiological roles. Identifying the core microbiome(s) are essential in the prediction of diseases and changes in environmental behavior of microorganisms. The dataset of 16S rRNA amplicon sequencing of Malayan tapir was deposited in the MG-RAST portal. Parameters such as quality control, taxonomic prediction (unknown and predicted), diversity (rarefaction), and diversity (alpha) were analyzed using sequencing approaches (Amplicon sequencing). Comparisons of parameters, according to the type of sequencing, showed significant differences, except for the prediction variable. In the Amplicon sequencing datasets, the parameters Rarefaction and Unknown had the highest correlation, while Alpha and Predicted had the lowest. Firmicutes, Bacteroidetes, Proteobacteria, Bacilli, and Bacteroidia were the most representative genera in Malayan tapir amplicon sequences, which indicated that most of the tapirs were healthy. However, continuous assessment to maintain the well-being of tapir for long term is still required. This chapter focuses on the introduction of 16S rRNA amplicon metagenomics in analyzing Malayan tapir gut microbiome dataset.
The Malaysian giant freshwater prawn, Macrobrachium rosenbergii, is an economically important crustacean worldwide. However, production of this prawn is facing a serious threat from Vibriosis disease caused by Vibrio species such as Vibrio parahaemolyticus. Unfortunately, the mechanisms involved in the immune response of this species to bacterial infection are not fully understood. We therefore used a high-throughput deep sequencing technology to investigate the transcriptome and comparative expression profiles of the hepatopancreas from this freshwater prawn infected with V. parahaemolyticus to gain an increased understanding of the molecular mechanisms underlying the species' immune response to this pathogenic bacteria.
Mycobacterium brisbanense is a member of Mycobacterium fortuitum third biovariant complex, which includes rapidly growing Mycobacterium spp. that normally inhabit soil, dust and water, and can sometimes cause respiratory tract infections in humans. We present the first whole-genome analysis of M. brisbanense UM_WWY which was isolated from a 70-year-old Malaysian patient. Molecular phylogenetic analyses confirmed the identification of this strain as M. brisbanense and showed that it has an unusually large genome compared with related mycobacteria. The large genome size of M. brisbanense UM_WWY (~7.7Mbp) is consistent with further findings that this strain has a highly variable genome structure that contains many putative horizontally transferred genomic islands and prophage. Comparative analysis showed that M. brisbanense UM_WWY is the only Mycobacterium species that possesses a complete set of genes encoding enzymes involved in the urea cycle, suggesting that this soil bacterium is able to synthesize urea for use as plant fertilizers. It is likely that M. brisbanense UM_WWY is adapted to live in soil as its primary habitat since the genome contains many genes associated with nitrogen metabolism. Nevertheless, a large number of predicted virulence genes were identified in M. brisbanense UM_WWY that are mostly shared with well-studied mycobacterial pathogens such as Mycobacterium tuberculosis and Mycobacterium abscessus. These findings are consistent with the role of M. brisbanense as an opportunistic pathogen of humans. The whole-genome study of UM_WWY has provided the basis for future work of M. brisbanense.
The marine clam Lutraria rhynchaena is gaining popularity as an aquaculture species in Asia. Lutraria populations are present in the wild throughout Vietnam and several stocks have been established and translocated for breeding and aquaculture grow-out purposes. In this study, we demonstrate the feasibility of utilising Illumina next-generation sequencing technology to streamline the identification and genotyping of microsatellite loci from this clam species. Based on an initial partial genome scan, 48 microsatellite markers with similar melting temperatures were identified and characterised. The 12 most suitable polymorphic loci were then genotyped using 51 individuals from a population in Quang Ninh Province, North Vietnam. Genetic variation was low (mean number of alleles per locus = 2.6; mean expected heterozygosity = 0.41). Two loci showed significant deviation from Hardy-Weinberg equilibrium (HWE) and the presence of null alleles, but there was no evidence of linkage disequilibrium among loci. Three additional populations were screened (n = 7-36) to test the geographic utility of the 12 loci, which revealed 100 % successful genotyping in two populations from central Vietnam (Nha Trang). However, a second population from north Vietnam (Co To) could not be successfully genotyped and morphological evidence and mitochondrial variation suggests that this population represents a cryptic species of Lutraria. Comparisons of the Qang Ninh and Nha Trang populations, excluding the 2 loci out of HWE, revealed statistically significant allelic variation at 4 loci. We reported the first microsatellite loci set for the marine clam Lutraria rhynchaena and demonstrated its potential in differentiating clam populations. Additionally, a cryptic species population of Lutraria rhynchaena was identified during initial loci development, underscoring the overlooked diversity of marine clam species in Vietnam and the need to genetically characterise population representatives prior to microsatellite development. The rapid identification and validation of microsatellite loci using next-generation sequencing technology warrant its integration into future microsatellite loci development for key aquaculture species in Vietnam and more generally, aquaculture countries in the South East Asia region.
Acute myeloid leukaemia with normal karyotype (AML-NK) is a cytogenetically cryptic yet heterogeneous clinical subgroup that lacks structural chromosomal abnormalities. Although classified as an intermediate-risk group, AML-NK exhibits variability in patient outcomes, underscoring the need for refined molecular stratification. In this study, high-throughput deep sequencing was performed on 51 AML-NK patients, with a comprehensive analysis of a Southeast Asian cohort to identify cryptic and novel fusion genes potentially implicated in disease pathogenesis and prognosis. Two pipelines (Arriba and STAR-fusion) were utilised in this study, which identified 68 findings involving 27 fusion genes, of which approximately 70% (n = 19/27) were novel. The majority (85%) were intrachromosomal fusion events, and 26% involved non-coding RNAs. The most recurrent and novel fusions were LATS2::SAP18 (17.6%) and HOXA3::HOXA9 (15.7%). In addition, known prognostic fusion genes, KMT2A-PTD, NUP98::NSD1, and NPM1::MLF1, were detected. The integration of these prognostic fusion genes into the ELN 2022 criteria reclassified 73% (n = 8/11) of intermediate-risk patients as poor-risk. Notably, the presence of KMT2A-PTD and/or LATS2::SAP18 was associated with reduced overall survival (p = 0.039, Log-Rank), although this significance was not retained in multivariate analysis. This study's findings revealed the complexity of the AML-NK transcriptome, highlighting the utility of RNA sequencing for prognostic risk assessment and therapeutic strategy development. KEY MESSAGE: RNA-seq uncovered 27 fusion genes in AML-NK, 70% of which were novel. Recurrent LATS2::SAP18 and HOXA3::HOXA9 fusions were identified. Fusion-based ELN 2022 reclassified most intermediate-risk patients. RNA-seq enhances prognostic assessment in AML-NK.
The MiSeq FGx Forensic Genomics System (Illumina) enables amplification and massively parallel sequencing of 59 STRs, 94 identity informative SNPs, 54 ancestry informative SNPs, and 24 phenotypic informative SNPs. Allele frequency and population statistics data were generated for the 172 SNP loci included in this panel on four major population groups (Chinese, African Americans, US Caucasians, and Southwest Hispanics). Single-locus and combined random match probability values were generated for the identity informative SNPs. The average combined STR and identity informative SNP random match probabilities (assuming independence) across all four populations were 1.75E-67 and 2.30E-71 with length-based and sequence-based STR alleles, respectively. Ancestry and phenotype predictions were obtained using the ForenSeq™ Universal Analysis System (UAS; Illumina) based on the ancestry informative and phenotype informative SNP profiles generated for each sample. Additionally, performance metrics, including profile completeness, read depth, relative locus performance, and allele coverage ratios, were evaluated and detailed for the 725 samples included in this study. While some genetic markers included in this panel performed notably better than others, performance across populations was generally consistent. The performance and population data included in this study support that accurate and reliable profiles were generated and provide valuable background information for laboratories considering internal validation studies and implementation.
BACKGROUND: Patient-Derived Tumour Xenografts (PDTXs) have emerged as the pre-clinical models that best represent clinical tumour diversity and intra-tumour heterogeneity. The molecular characterization of PDTXs using High-Throughput Sequencing (HTS) is essential; however, the presence of mouse stroma is challenging for HTS data analysis. Indeed, the high homology between the two genomes results in a proportion of mouse reads being mapped as human.
RESULTS: In this study we generated Whole Exome Sequencing (WES), Reduced Representation Bisulfite Sequencing (RRBS) and RNA sequencing (RNA-seq) data from samples with known mixtures of mouse and human DNA or RNA and from a cohort of human breast cancers and their derived PDTXs. We show that using an In silico Combined human-mouse Reference Genome (ICRG) for alignment discriminates between human and mouse reads with up to 99.9% accuracy and decreases the number of false positive somatic mutations caused by misalignment by >99.9%. We also derived a model to estimate the human DNA content in independent PDTX samples. For RNA-seq and RRBS data analysis, the use of the ICRG allows dissecting computationally the transcriptome and methylome of human tumour cells and mouse stroma. In a direct comparison with previously reported approaches, our method showed similar or higher accuracy while requiring significantly less computing time.
CONCLUSIONS: The computational pipeline we describe here is a valuable tool for the molecular analysis of PDTXs as well as any other mixture of DNA or RNA species.
The crucifix crab, Charybdis feriatus, which mainly inhabits Indo-Pacific region, is regarded as one of the most high-potential species for domestication and incorporation into the aquaculture sector. However, the regulatory mechanisms of sex determination and differentiation of this species remain unclear. To identify candidate genes involved in sex determination and differentiation, high throughput sequencing of transcriptome from the testis and ovary of C. feriatus was performed by the Illumina platform. After removing adaptor primers, low-quality sequences and very short (<50 nt) reads, we obtained 80.9 million and 66.2 million clean reads from testis and ovary, respectively. A total of 86,433 unigenes were assembled, and ~43% (37,500 unigenes) were successfully annotated to the NR, NT, Swiss-Prot, KEGG, COG, GO databases. By comparing the testis and ovary libraries, we obtained 27,636 differentially expressed genes. Some candidate genes involved in the sex determination and differentiation of C. feriatus were identified, such as vasa, pgds, vgr, hsp90, dsx-f, fem-1, and gpr. In addition, 88,608 simple sequence repeats were obtained, and 61,929 and 77,473 single nucleotide polymorphisms from testis and ovary were detected, respectively. The transcriptome profiling was validated by quantitative real-time PCR in 30 selected genes, which showed a good consistency. The present study is the first high-throughput transcriptome sequencing of C. feriatus. These findings will be useful for future functional analysis of sex-associated genes and molecular marker-assisted selections in C. feriatus.
High-depth next generation sequencing data provide valuable insights into the number and distribution of RNA editing events. Here, we report the RNA editing events at cellular level of human primary monocyte using high-depth whole genomic and transcriptomic sequencing data. We identified over a ten thousand putative RNA editing sites and 69% of the sites were A-to-I editing sites. The sites enriched in repetitive sequences and intronic regions. High-depth sequencing datasets revealed that 90% of the canonical sites were edited at lower frequencies (<0.7). Single and multiple human monocytes and brain tissues samples were analyzed through genome sequence independent approach. The later approach was observed to identify more editing sites. Monocytes was observed to contain more C-to-U editing sites compared to brain tissues. Our results establish comparable pipeline that can address current limitations as well as demonstrate the potential for highly sensitive detection of RNA editing events in single cell type.
Understanding the relationship between genotypes and phenotypes is essential to disentangle biological mechanisms and to unravel the molecular basis of diseases. Genes and proteins are closely linked in biological systems. However, genomics and proteomics have developed separately into two distinct disciplines whereby crosstalk among scientists from the two domains is limited and this constrains the integration of both fields into a single data modality of useful information. The emerging field of proteogenomics attempts to address this by building bridges between the two disciplines. In this review, how genomics and transcriptomics data in different formats can be utilized to assist proteogenomics application is briefly discussed. Subsequently, a much larger part of this review focuses on proteogenomics research articles that are published in the last five years that answer two important questions. First, how proteogenomics can be applied to tackle biological problems is discussed, covering genome annotation and precision medicine. Second, the latest developments in analytical technologies for data acquisition and the bioinformatics tools to interpret and visualize proteogenomics data are covered.
The incomplete identification of structural variants (SVs) from whole-genome sequencing data limits studies of human genetic diversity and disease association. Here, we apply a suite of long-read, short-read, strand-specific sequencing technologies, optical mapping, and variant discovery algorithms to comprehensively analyze three trios to define the full spectrum of human genetic variation in a haplotype-resolved manner. We identify 818,054 indel variants (<50 bp) and 27,622 SVs (≥50 bp) per genome. We also discover 156 inversions per genome and 58 of the inversions intersect with the critical regions of recurrent microdeletion and microduplication syndromes. Taken together, our SV callsets represent a three to sevenfold increase in SV detection compared to most standard high-throughput sequencing studies, including those from the 1000 Genomes Project. The methods and the dataset presented serve as a gold standard for the scientific community allowing us to make recommendations for maximizing structural variation sensitivity for future genome sequencing studies.
Antibodies leverage on their unique architecture to bind with an array of antigens. The strength of interaction has a direct relation to the affinity of the antibodies towards the antigen. In vivo affinity maturation is performed through multiple rounds of somatic hypermutation and selection in the germinal centre. This unique process involves intricate sequence rearrangements at the gene level via molecular mechanisms. The emergence of in vitro display technologies, mainly phage display and recombinant DNA technology, has helped revolutionize the way antibody improvements are being carried out in the laboratory. The adaptation of molecular approaches in vitro to replicate the in vivo processes has allowed for improvements in the way recombinant antibodies are designed and tuned. Combinatorial libraries, consisting of a myriad of possible antibodies, are capable of replicating the diversity of the natural human antibody repertoire. The isolation of target-specific antibodies with specific affinity characteristics can also be accomplished through modification of stringent protocols. Despite the ability to screen and select for high-affinity binders, some 'fine tuning' may be required to enhance antibody binding in terms of its affinity. This review will provide a brief account of phage display technology used for antibody generation followed by a summary of different combinatorial library characteristics. The review will focus on available strategies, which include molecular approaches, next generation sequencing, and in silico approaches used for antibody affinity maturation in both therapeutic and diagnostic applications.
The emergence of Massively Parallel Sequencing technologies enabled the analysis of full mitochondrial (mt)DNA sequences from forensically relevant samples that have, so far, only been typed in the control region or its hypervariable segments. In this study, we evaluated the performance of a commercially available multiplex-PCR-based assay, the Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific), for the amplification and sequencing of the entire mitochondrial genome (mitogenome) from even degraded forensic specimens. For this purpose, more than 500 samples from 24 different populations were selected to cover the vast majority of established superhaplogroups. These are known to harbor different signature sequence motifs corresponding to their phylogenetic background that could have an effect on primer binding and, thus, could limit a broad application of this molecular genetic tool. The selected samples derived from various forensically relevant tissue sources and were DNA extracted using different methods. We evaluated sequence concordance and heteroplasmy detection and compared the findings to conventional Sanger sequencing as well as an orthogonal MPS platform. We discuss advantages and limitations of this approach with respect to forensic genetic workflow and analytical requirements.