Illumina Next-generation Sequencing Research Articles

Genomic patterns of selection in morphometric traits across diverse Indian cattle breeds.

This study seeks a comprehensive exploration of genome-wide selective processes impacting morphometric traits across diverse cattle breeds, utilizing an array of statistical methods. Morphometric traits, encompassing both qualitative and quantitative variables, play a pivotal role in characterizing and selecting livestock breeds based on their external appearance, size, and physical attributes. While qualitative traits, such as color, horn structure, and coat type, contribute to adaptive features and breed identification, quantitative traits like body weight and conformation measurements bear a closer correlation with production characteristics. This study employs advanced genotyping technologies, including the Illumina BovineSNP50 Bead Chip and next-generation sequencing methods like Reduced Representation sequencing, to identify genomic signatures associated with these traits. We applied four intra-population methods to find evidence of selection, such as Tajima's D, CLR, iHS, and ROH. We found a total of 40 genes under the selection signature, that were associated with morphometric traits in five cattle breeds (Kankrej, Tharparkar, Nelore, Sahiwal, and Gir). Crucial genes such as ADIPDQ, DPP6, INSIG1, SLC35D2 in Kankrej, LPL, ATP6V1B2, CDC14B in Tharparkar, HPSE2, PLAG1 in Nelore, PCSK1, PRKD1 in Sahiwal, and GNAQ, HPCAL1 in Gir were identified in our study. This approach provides valuable insights into the genetic basis of variations in body weight and conformation traits, facilitating informed selection processes and offering a deeper understanding of the evolutionary and domestication processes in diverse cattle breeds.

Fluorogenic Aptamer Optimizations on a Massively Parallel Sequencing Platform.

F luorogenic ap tamers (FAPs) have become an increasingly important tool in cellular sensing and pathogen diagnostics. However, fine-tuning FAPs for enhanced performance remains challenging even with the structural details provided by X-ray crystallography. Here we present a novel approach to optimize a DNA-based FAP (D-FAP), Lettuce, on repurposed Illumina next-generation sequencing (NGS) chips. When substituting its cognate chromophore, DFHBI-1T, with TO1-biotin, Lettuce not only shows a red-shifted emission peak by 53 nm (from 505 to 558 nm), but also a 4-fold bulk fluorescence enhancement. After screening 8,821 Lettuce variants complexed with TO1-biotin, the C14T mutation is found to exhibit an improved apparent dissociated constant ( vs. 0.82 µM), an increased quantum yield (QY: 0.62 vs. 0.59) and an elongated fluorescence lifetime (τ: 6.00 vs. 5.77 ns), giving 45% more ensemble fluorescence than the canonical Lettuce/TO1-biotin complex. Molecular dynamic simulations further indicate that the π-π stacking interaction is key to determining the coordination structure of TO1-biotin in Lettuce. Our screening-and-simulation pipeline can effectively optimize FAPs without any prior structural knowledge of the canonical FAP/chromophore complexes, providing not only improved molecular probes for fluorescence sensing but also insights into aptamer-chromophore interactions.

Sequence comparison of the mitochondrial genomes of Plesionika species (Caridea: Pandalidae), gene rearrangement and phylogenetic relationships of Caridea.

Despite the Caridean shrimps' vast species richness and ecological diversity, controversies persist in their molecular classification. Within Caridea, the Pandalidae family exemplifies significant taxonomic diversity. As of June 25, 2023, GenBank hosts only nine complete mitochondrial genomes (mitogenomes) for this family. The Plesionika genus within Pandalidae is recognized as polyphyletic. To improve our understanding of the mitogenome evolution and phylogenetic relationships of Caridea, this study introduces three novel mitogenome sequences from the Plesionika genus: P. ortmanni, P. izumiae and P. lophotes. The complete mitochondrial genomes of three Plesionika species were sequenced utilizing Illumina's next-generation sequencing (NGS) technology. After assembling and annotating the mitogenomes, we conducted structural analyses to examine circular maps, sequence structure characteristics, base composition, amino acid content, and synonymous codon usage frequency. Additionally, phylogenetic analysis was performed by integrating existing mitogenome sequences of true shrimp available in GenBank. The complete mitogenomes of the three Plesionika species encompass 37 canonical genes, comprising 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region (CR). The lengths of these mitogenomes are as follows: 15,908 bp for P. ortmanni, 16,074 bp for P. izumiae and 15,933 bp for P. lophotes. Our analyses extended to their genomic features and structural functions, detailing base composition, gene arrangement, and codon usage. Additionally, we performed selection pressure analysis on the PCGs of all Pandalidae species available in Genbank, indicating evolutionary purification selection acted on the PCGs across Pandalidae species. Compared with the ancestral Caridea, translocation of two tRNA genes, i.e., trnP or trnT, were found in the two newly sequenced Plesionika species-P. izumiae and P. lophotes. We constructed a phylogenetic tree of Caridea using the sequences of 13 PCGs in mitogenomes. The results revealed that family Pandalidae exhibited robust monophyly, while genus Plesionika appeared to be a polyphyletic group. Gene rearrangements within the Pandalidae family were observed for the first time. Furthermore, a significant correlation was discovered between phylogenetics of the Caridea clade and arrangement of mitochondrial genes. Our findings offer a detailed exploration of Plesionika mitogenomes, laying a crucial groundwork for subsequent investigations into genetic diversity, phylogenetic evolution, and selective breeding within this genus.

Profiling of microRNAs by next-generation sequencing: Potential biomarkers for diffuse large B-cell lymphoma

BackgroundLymphoma ranks fifth in prevalence among common cancer types worldwide. This lymphatic system cancer arises from T or B cells. Diffuse large B cell lymphomas (DLBCLs) are associated with most non-Hodgkin lymphomas. Non-coding microRNAs (miRNAs) greatly affect gene expression. A single miRNA can target numerous genes, thus largely influencing gene expression networks. MiRNAs can act as oncogenes or tumor suppressors in controlling DLBCL progression. This study investigated the roles of miRNAs in patients with DLBCL through next-generation sequencing, which was found to be sensitive, accurate, and robust. MethodsThe study involved seven patients with DLBCLs and three controls at a hematology-oncology clinic. MiRNA was extracted from existing formalin-fixed, paraffin-embedded (FFPE) tissue specimens. Illumina next-generation sequencing was used to sequence samples for miRNA profiling. ResultsSamples from patients showed expression of various hsa-mir miRNAs (1248, 3607, 21, 142, 1244, 182, 6516, 766, 1291, 4449, and 181a), whereas those from healthy individuals showed expression of hsa-mir 1248, 3607, 21, 142, and 877. Hsa-mir-877-3p is known to target multiple genes, and miRNAs such as hsa-mir-877-3p, hsa-mir-1291, and hsa-mir-181a-5p interact primarily with target genes. ConclusionsMiRNA profiling in FFPE tissues from patients with DLBCL suggested that miRNA levels can distinguish patients with DLBCL from controls, and therefore may provide prognostic or diagnostic biomarkers for DLBCL. Altered genes and miRNAs may also be potential therapeutic targets.

Sequencing of the Complete Mitochondrial Genome of the Big Brown Mactra Clam, Mactra grandis (Venerida: Mactridae).

Mitochondrial genomes are playing an increasingly important role in molluscan taxonomy, germplasm, and evolution studies. The first complete mitochondrial genome of the commercial big brown mactra clam, Mactra grandis, was characterized using Illumina next-generation sequencing in this study. The 17,289 bp circular genome has a typical gene organization of 13 protein-coding genes (PCGs), 2 rRNAs, and 22 tRNAs, with an obvious (A + T)-bias of 64.54%. All PCGs exhibited a homogeneous bias in nucleotide composition with a (A + T)-bias, a positive GC skew, and a negative AT skew. Results of phylogenetic analysis showed that Mactra grandis was most closely related to Mactra cygnus. The functional gene arrangement of the two species was identical but different from other Mactra species. The congeneric relationships among Mactra species were demonstrated by genetic distance analysis. Additionally, the selective pressure analysis suggested that cox1 was highly efficient for discriminating closely related species in genus Mactra, while nad2 was the most appropriate marker for population genetic analysis.

Insights into soil nematode diversity and bacterial community of Thai jasmine rice rhizosphere from different paddy fields in Thailand.

Globally, phytonematodes cause significant crop losses. Understanding the functions played by the plant rhizosphere soil microbiome during phytonematodes infection is crucial. This study examined the distribution of phytonematodes in the paddy fields of five provinces in Thailand, as well as determining the keystone microbial taxa in response to environmental factors that could be considered in the development of efficient biocontrol tactics in agriculture. The results demonstrated that Meloidogyne graminicola and Hirschmanniella spp. were the major and dominant phytonematodes distributed across the paddy fields of Thailand. Soil parameters (total P, Cu, Mg, and Zn) were the important factors affecting the abundance of both nematodes. Illumina next-generation sequencing demonstrated that the levels of bacterial diversity among all locations were not significantly different. The Acidobacteriota, Proteobacteria, Firmicutes, Actinobacteriota, Myxococcota, Chloroflexi, Verrucomicrobiota, Bacteroidota, Gemmatimonadota, and Desulfobacterota were the most abundant bacterial phyla observed at all sites. The number of classes of the Acidobacteriae, Clostridia, Bacilli, and Bacteroidia influenced the proportions of Hirschmanniella spp., Tylenchorhynchus spp., and free-living nematodes in the sampling dirt, whereas the number of classes of the Polyangia and Actinobacteria affected the amounts of Pratylenchus spp. in both roots and soils. Soil organic matter, N, and Mn were the main factors that influenced the structure of the bacterial community. Correlations among rhizosphere microbiota, soil nematodes, and soil properties will be informative data in considering phytonematode management in a rice production system.

A novel homozygous missense mutation in the RPE65 gene in an Iranian family with retinitis pigmentosa

Hereditary retinal dystrophies such as retinitis pigmentosa (RP) are genetically and phenotypically heterogeneous and result in total blindness. One of the key genes causing RP is RPE65. In 2017, a gene therapy drug, voretigene neparvovec-rzyl, for RP associated with RPE65 mutations, received Food and Drug Administration approval. So, identification of RPE65 gene mutations in affected patients which may be amenable to genetic treatment is critical. Here, we examine three related Iranian patients with RP which have not been previously described in the population. Whole exome sequencing (WES) using next-generation Illumina sequencing was utilized to identify the disease-causing mutation in the proband. Subsequently, Sanger sequencing was performed to confirm the identified mutation in the patients and their family members. Finally, the possible effect of the reported mutation on the target protein was also evaluated using bioinformatics tools. WES identified a novel homozygous missense mutation c.T170C: p.F57S in exon 3 of the RPE65 gene. Sanger sequencing confirmed the mutation found in the patients and their healthy family members. This mutation segregated with the autosomal recessive inheritance pattern of RP disorder. Furthermore, the pathogenic nature of the identified mutation was confirmed by in silico approaches and genotype-phenotype correlation analysis. Our findings revealed a novel missense mutation c.T170C: p.F57S in the RPE65 gene and could expand the pathogenic mutations spectrum of the RPE65 gene. Additionally, such studies can aid to conduct genetic counseling, prenatal diagnosis, clinical management and even genetic treatment more accurately for all affected patients with RPE65 gene mutations.

De Novo Transcriptome Analysis by PacBio SMRT-Seq and Illumina RNA-Seq Provides New Insights into Polyphenol Biosynthesis in Chinese Olive Fruit

Polyphenols play a crucial role in fruit flavor. To elucidate the mechanism of fruit polyphenol metabolism, we constructed a transcriptome atlas through PacBio single-molecule real-time (SMRT) sequencing and Illumina next-generation sequencing (NGS) using Canarium album (Lour.) Raeusch., which is a fantastic fruit rich in polyphenolic compounds. In this work, PacBio full-length transcriptome assembly generated 135,439 isoforms with an average length of all isoforms of 2687.94 bp and an N50 length of 3224 bp. To gain deeper insights into the molecular mechanisms of polyphenol biosynthesis in C. album, we constructed twelve RNA-Seq libraries from four developmental stages of the fruits. We identified a total of 28,658 differentially expressed genes (DEGs). We found that many DEGs were involved in metabolic pathways, biosynthesis of secondary metabolites, biosynthesis of antibiotics, starch and sucrose metabolism, and plant hormone signal transduction. Here, we report the expression profiles of 215 DEGs encoding 27 enzymes involved in the polyphenol biosynthesis pathway in C. album. In addition, 285 differentially expressed transcription factors (TFs) continuously down-regulated in four developmental periods of C. album fruit, which may indicate their potential role in the response to polyphenol metabolism and phenylpropanoid biosynthesis pathways. This report will help us understand polyphenol biosynthesis’s functions and metabolic mechanism in C. album. The transcriptome data provide a valuable resource for genetic and genomics research. They will facilitate future work exploiting C. album and other fruits used as medicine and food.

Deciphering bacterial community structure in Cola lepidota (K.) Schum fruits using Illumina next-generation sequencing

Cola lepidota cultivation is hampered by pathogenic microorganisms which affect the quality and quantity of produce after harvest. Bacterial organisms cause softening of tissues which may lead to the spoilage of the whole fruit. This study aimed to detect the bacterial organisms and the bacterial community structure associated with monkey cola (Cola lepidota K. Schum) fruits using Illumina next-generation sequencing (NGS). Amplification of the v4 regions of the 16S rDNA gene was performed through Polymerase chain reaction (PCR) using the primers, 515F and 806R. Picking of the Operational Taxonomic Units (OTUs) and assignment of the OTUs into different taxonomic levels were performed on GREENGENES database. The most abundant phyla obtained from Cola lepidota mesocarp and seeds were Proteobacteria (40.83 %), Bacteroidetes (31.01 %) and Firmicutes (24.52 %). The abundance of Bacteroides, which was the predominant genus in the majority of the samples, was between 31 and 72 % among the total bacteria identified in most samples. Other relatively high abundant genera obtained across the total bacterial community of C. lepidota were Pseudomonas (19.20 %), Tepidimonas (13.76 %) and Acinetobacter (13.33 %). We obtained potentially beneficial, human and phytopathogenic microorganisms which play important roles in evolutionary and ecological processes. Among these genera include Acinetobacter and Staphylococcus which are opportunistic human pathogens. This study revealed the bacterial organisms associated with monkey cola fruits and indicates that Illumina next generation sequencing can be used to determine the bacterial community structure and bacterial organisms in fruits.

Phenotypic and genomic characteristics of clinical IMP-producing Klebsiella spp. Isolates in China

BackgroundIMP-producing Klebsiella spp. (IMPKsp) strains have spread globally, including in China. Currently, the prevalence and genomic characterization of IMPKsp is largely unknown nationwide. Here we aimed to provide a general overview of the phenotypic and genomic characteristics of IMPKsp strains.Methods61 IMPKsp strains were obtained from 13 provinces in China during 2016-2021. All strains were tested for their susceptibility to antimicrobial agents by the microdilution broth method and sequenced with Illumina next-generation sequencing. We performed conjugation experiments on thirteen representative strains which were also sequenced by Oxford nanopore sequencing technology to characterize blaIMP-encoding plasmids.ResultsWe find that all IMPKsp strains display multidrug-resistant (MDR) phenotypes. All strains belong to 27 different STs. ST307 emerges as a principal IMP-producing sublineage. blaIMP-4 is found to be the major isoform, followed by blaIMP-38. Seven incompatibility types of blaIMP-encoding plasmids are identified, including IncHI5 (32/61, 52.5%), IncN-IncR (10/61, 16.4%), IncFIB(K)-HI1B (7/61, 11.5%), IncN (5/61, 8.2%), IncN-IncFII (2/61, 3.3%), IncFII (1/61, 1.6%) and IncP (1/61, 1.6%). The strains carrying IncHI5 and IncN plasmids belong to diverse ST types, indicating that these two plasmids may play an important role in the transmission of blaIMP genes among Klebsiella spp. strains.ConclusionsOur results highlight that multi-clonal transmission, multiple genetic environments and plasmid types play a major role in the dissemination process of blaIMP genes among Klebsiella spp. IncHI5 type plasmids have the potential to be the main vectors mediating the spread of the blaIMP genes in Klebsiella spp.

Microbial diversity in dairy manure environment under liquid-solid separation systems

ABSTRACT In dairy manure, a wide array of microorganisms, including many pathogens, survive and grow under suitable conditions. This microbial community offers a tremendous opportunity for studying animal health, the transport of microbes into the soil, air, and water, and consequential impacts on public health. The aim of this study was to assess the impacts of manure management practices on the microbial community of manure. The key novelty of this work is to identify the impacts of various stages of manure management on microbes living in dairy manure. In general, the majority of dairy farms in California use a flush system to manage dairy manure, which involves liquid-solid separations. To separate liquid and solid in manure, Multi-stage Alternate Dairy Effluent Management Systems (ADEMS) that use mechanical separation systems (MSS) or weeping wall separation systems (WWSS) are used. Thus, this study was conducted to understand how these manure management systems affect the microbial community. We studied the microbial communities in the WWSS and MSS separation systems, as well as in the four stages of the ADEMS. The 16S rRNA gene from the extracted genomic DNA of dairy manure was amplified using the NovoSeq Illumina next-generation sequencing platform. The sequencing data were used to perform the analysis of similarity (ANOSIM) and multi-response permutation procedure (MRRP) statistical tests, and the results showed that microbial communities among WWSS and MSS were significantly different (p < 0.05). These findings have significant practical implications for the design and implementation of manure management practices in dairy farms.

Gastric Perineurioma: A Rare Entity with Molecular Analysis and Literature Review.

Perineuriomas of the gastrointestinal tract are benign neoplasms that commonly develop in the distal colon and are identified during screening colonoscopy; however, perineuriomas of the stomach are exceedingly rare and less frequently identified. Differentiating gastric perineuriomas from other more serious gastric neoplasms is critical to avoid unnecessarily aggressive treatments. Thus far, only six patients with gastric perineurioma have been described, and the molecular characterization of this entity is still lacking. We report a 52-year-old woman who presented with abdominal pain and gastric acid reflux and was found to have a 1.5 cm subepithelial gastric neoplasm composed of bland spindle cells displacing the gastric glands with no cytologic atypia or mitotic activity, suggesting a benign spindle cell neoplasm. Immunohistochemical analysis showed reactivity for perineurial markers glucose transporter-1 and epithelial membrane antigen, consistent with benign gastric perineurioma. DNA extracted from the tissue was used for a capture-based target sequence enrichment panel followed by Illumina next-generation sequencing and targeted bioinformatic analysis for oncogenic alterations within defined disease-associated target regions. No sequence variants in the BRAF gene were identified. This rare case of gastric perineurioma helps solidify our understanding of how to discern various types of gastric neoplasms through traditional laboratory analysis alongside genetic sequencing approaches. Although extremely rare, gastric perineurioma should be kept in the differential diagnosis when assessing spindle cell gastric tumors to avoid unnecessary therapies, and physicians should understand the molecular characteristics of benign versus malignant tumors.

Changes in DNA methylation associated with a specific mode of delivery: a pilot study.

The mode of delivery represents an epigenetic factor with potential to affect further development of the individual by multiple mechanisms. DNA methylation may be one of them, representing a major epigenetic mechanism involving direct chemical modification of the individual's DNA. This pilot study aims to examine whether a specific mode of delivery induces changes of DNA methylation by comparing the umbilical cord blood and peripheral blood of the newborns. Blood samples from infants born by vaginal delivery and caesarean section were analysed to prepare the Methylseq library according to NEBNext enzymatic Methyl-seq Methylation Library Preparation Kit with further generation of target-enriched DNA libraries using the Twist Human Methylome Panel. DNA methylation status was determined using Illumina next-generation sequencing (NGS). We identified 168 differentially methylated regions in umbilical cord blood samples and 157 regions in peripheral blood samples. These were associated with 59 common biological, metabolic and signalling pathways for umbilical cord and peripheral blood samples. Caesarean section is likely to represent an important epigenetic factor with the potential to induce changes in the genome that could play an important role in development of a broad spectrum of disorders. Our results could contribute to the elucidation of how epigenetic factors, such as a specific mode of delivery, could have adverse impact on health of an individual later in their life.

Effects of co-inoculation of indole-3-acetic acid- and ammonia-producing bacteria on plant growth and nutrition, soil elements, and the relationships of soil microbiomes with soil physicochemical parameters

Abstract To cope with the problem of the reduction in agricultural productivity, the massive usage of chemical fertilizers for boosting agricultural productivity is popular across the world. In return, it thereby declines soil quality and biodiversity. We therefore determined the potential of co-inoculation of the indole-3-acetic-acid- and ammonia-producing bacteria, Pseudomonas sp. strain NK2/7 and Bacillus sp. strain NK1/19, immobilized in agar to promote mustard greens growth and investigated the effects of co-inoculation of both strains on plant and soil elements as well as soil bacterial community structure, in comparison with a chemical fertilizer and single-strain inoculations. The knowledge on the relationships of both species with plant and soil microbiota in field conditions is limited. The results show that the co-inoculation of both strains immobilized in agar increased plant dry weight up to 62.02% and plant elements, including N (38.46%), P (28.37%), K (203.16%), Ca (10.89%), Mg (25.8%), and Na (70.3%). The co-inoculation also conferred the highest organic matter (OM) (by 429.45%) and increased soil elements including Fe (92.74%), Mn (126.88%), Cu (12.05%), and Zn (107.46%). Illumina next-generation sequencing demonstrated that the establishment of both strains affected the soil bacterial community by reducing bacterial diversity and richness. Bacterial classes Bacteroidia, alpha-Proteobacteria, Vicinamibacteria, Gemmatimonadetes, Verrucomicrobiae, Polyangia, Nitrososphaeria, and Blastocatellia were significantly decreased, whereas the amounts of gamma-Proteobacteria and Bacilli were significantly increased. Soil parameters, including pH, OM, total Fe, Mn, Cu, and Zn, were the major factors influencing the soil bacterial community structure.

Phenotypic and Genotypic Characterizations of Antimicrobial-Resistant Escherichia coli Isolates from Diverse Retail Meat Samples in North Carolina During 2018-2019.

Surveillance of antimicrobial-resistant pathogens in U.S. retail meats is conducted to identify potential risks of foodborne illness. In this study, we conducted a phenotypic and genotypic analysis of Escherichia coli recovered from a diverse range of retail meat types during 2018-2019 in North Carolina. The investigation was conducted as part of the National Antimicrobial Resistance Monitoring System (NARMS). Retail meat sampling and E. coli isolation were performed in accordance with NARMS retail meat isolation protocols. We used the Sensititre™ broth microdilution system to determine phenotypic resistance to 14 antimicrobial agents and the Illumina next-generation sequencing platform for genotypic resistance profiling. The highest prevalence of E. coli isolates was found in ground turkey (n = 57, 42.9%) and chicken (n = 27, 20.3%), followed by ground beef (n = 25, 18.9%) and pork (n = 24, 18%). The isolates were divided into seven different phylogroups using the Clermont typing tool, with B1 (n = 59, 44.4%) and A (n = 39, 29.3%) being the most dominant, followed by B2 (n = 14, 10.5%), D (n = 7, 5.3%), F (n = 6, 4.5%), E (n = 3, 2.3%), and C (n = 2, 1.5%). Using multilocus sequence typing (MLST), 128 Sequence types (STs) were identified indicating high diversity. Phenotypic and genotypic resistance was observed toward aminoglycosides, sulfonamides, beta-lactams, macrolides, tetracyclines, phenicols, and fluoroquinolones. Ground turkey samples were more resistant to the panel of tested antimicrobials than chicken, beef, or pork (p < 0.05). All isolates were found to be susceptible to meropenem. A high percentage of turkey isolates (n = 16, 28%) were multidrug-resistant (MDR) compared with 18.5% of chicken (n = 5), 8.4% of pork (n = 2), and 8% of beef isolates (n = 2). This study highlights the benefit of surveillance to identify MDR E. coli for epidemiologic tracking and is a comprehensive report of the phenotypic and genotypic characterization of E. coli isolated from retail meats in North Carolina.

An overlooked phenomenon: complex interactions of potential error sources on the quality of bacterial de novo genome assemblies

BackgroundParameters adversely affecting the contiguity and accuracy of the assemblies from Illumina next-generation sequencing (NGS) are well described. However, past studies generally focused on their additive effects, overlooking their potential interactions possibly exacerbating one another’s effects in a multiplicative manner. To investigate whether or not they act interactively on de novo genome assembly quality, we simulated sequencing data for 13 bacterial reference genomes, with varying levels of error rate, sequencing depth, PCR and optical duplicate ratios.ResultsWe assessed the quality of assemblies from the simulated sequencing data with a number of contiguity and accuracy metrics, which we used to quantify both additive and multiplicative effects of the four parameters. We found that the tested parameters are engaged in complex interactions, exerting multiplicative, rather than additive, effects on assembly quality. Also, the ratio of non-repeated regions and GC% of the original genomes can shape how the four parameters affect assembly quality.ConclusionsWe provide a framework for consideration in future studies using de novo genome assembly of bacterial genomes, e.g. in choosing the optimal sequencing depth, balancing between its positive effect on contiguity and negative effect on accuracy due to its interaction with error rate. Furthermore, the properties of the genomes to be sequenced also should be taken into account, as they might influence the effects of error sources themselves.

Comparative transcriptomic analysis of Illumina and MGI next-generation sequencing platforms using RUNX3- and ZBTB46-instructed embryonic stem cells.

Introduction: We have previously observed phenotypic and developmental changes upon the ectopic expression of the RUNX3 or the ZBTB46 transcription factors in mouse embryonic stem cell (ESC) derived progenitors. In this study, we evaluated the gene expression profiles of the RUNX3- and the ZBTB46-instructed murine ESCs with RNA-seq testing two next-generation sequencing technologies. Methods: We compared the DNA nanoball-based DNBSEQ G400 sequencer (MGI) with the bridge-PCR-based NextSeq 500 instrument (Illumina) for RNA sequencing. Moreover, we also compared two types of MGI sequencing reagents (Standard versus Hot-massive parallel sequencing (MPS)) with the DNBSEQ G400. Results: We observed that both sequencing platforms showed comparable levels of quality, sequencing uniformity, and gene expression profiles. For example, highly overlapping RUNX3- and ZBTB46-regulated gene lists were obtained from both sequencing datasets. Moreover, we observed that the Standard and the Hot-MPS-derived RUNX3- and ZBTB46-regulated gene lists were also considerably overlapped. This transcriptome analysis also helped us to identify differently expressed genes in the presence of the transgenic RUNX3 or ZBTB46. For example, we found that Gzmb, Gzmd, Gzme, Gdf6, and Ccr7 genes were robustly upregulated upon the forced expression of Runx3; on the other hand, Gpx2, Tdpoz4, and Arg2 were induced alongside the ectopic expression of Zbtb46. Discussion: Similar gene expression profile and greatly overlapping RUNX3- and ZBTB46-regulated gene sets were detected with both DNA sequencing platforms. Our analyses demonstrate that both sequencing technologies are suitable for transcriptome profiling and target gene selection. These findings suggest that DNBSEQ G400 represents a cost-effective alternative sequencing platform for gene expression monitoring. Moreover, this analysis provides a resource for exploration of the RUNX3- and ZBTB46-dependent gene regulatory networks.

Amplicon-Based Bisulfite Conversion-NGS DNA Methylation Analysis Protocol.

DNA methylation is an important epigenetic modification that regulates chromatin structure and the cell-type-specific expression of genes. The association of aberrant DNA methylation with many diseases, as well as the increasing interest in modifying the methylation mark in a directed manner at genomic sites using epigenome editing for research and therapeutic purposes, increases the need for easy and efficient DNA methylation analysis methods. The standard approach to analyze DNA methylation with a single-cytosine resolution is bisulfite conversion of DNA followed by next-generation sequencing (NGS). In this chapter, we describe a robust, powerful, and cost-efficient protocol for the amplification of target regions from bisulfite-converted DNA, followed by a second PCR step to generate libraries for Illumina NGS. In the two consecutive PCR steps, first, barcodes are added to individual amplicons, and in the second PCR, indices and Illumina adapters are added to the samples. Finally, we describe a detailed bioinformatics approach to extract DNA methylation levels of the target regions from the sequencing data. Combining barcodes with indices enables a high level of multiplexing allowing to sequence multiple pooled samples in the same sequencing run. Therefore, this method is a robust, accurate, quantitative, and cheap approach for the readout of DNA methylation patterns at defined genomic regions.

Exploring the Diversity of Plant-Associated Viruses and Related Viruses in Riverine Freshwater Samples Collected in Berlin, Germany.

Plant-infecting RNA viruses from 30 families and floating genera, as well as a great number of uncultured as yet-unclassified plant-associated viruses have been described. Even so, the plant RNA virosphere is still underexplored. RNA extracted from enriched virus particles of 50 L water samples from the Teltow Canal and the Havel River in Berlin, Germany, was sequenced using Illumina next-generation sequencing. Sequences were searched for plant viruses with BLAST and DIAMOND. Phylogenetic analyses were conducted with IQ-TREE 2. Altogether, 647 virus sequences greater than 1 kb were detected and further analyzed. These data revealed the presence of accepted and novel viruses related to Albetovirus, Alphaflexiviridae, Aspiviridae, Bromoviridae, Endornaviridae, Partitiviridae, Potyviridae, Solemoviridae, Tombusviridae and Virgaviridae. The vast majority of the sequences were novel and could not be taxonomically assigned. Several tombus- and endorna-like viruses make use of alternative translation tables that suggest unicellular green algae, ciliates, or diplomonades as their hosts. The identification of 27 albeto-like satellite viruses increases available sequence data five-fold. Sixteen new poty-like viruses align with other poty-like viruses in a link that combines the Astroviridae and Potyviridae families. Further, the identification of viruses with peptidase A6-like and peptidase A21-like capsid proteins suggests horizontal gene transfer in the evolution of these viruses.

Characterization of ‘&lt;em&gt;Ca&lt;/em&gt;. Phytoplasma asteris’, ‘&lt;em&gt;Ca&lt;/em&gt;. Phytoplasma australiense’ and ‘&lt;em&gt;Ca&lt;/em&gt;. Phytoplasma mali’ Associated with Russian Olive Witches’-Broom Disease in Iran Using Next-Generation Sequencing

Russian olive trees showing witches’-broom disease symptoms in urban green spaces and orchards in northwest Iran were sampled for phytoplasma detection. PCR assays and Sanger sequencing of 16S rRNA gene confirmed that ‘&lt;em&gt;Ca&lt;/em&gt;. Phytoplasma asteris’ was associated with Russian olive witches’-broom disease. For genomic characterization of phytoplasma associated with disease, the total DNA of an infected tree was analyzed by Illumina next-generation sequencing (NGS). The NGS analysis generated 46011389 read pairs of raw data with about 6.8 × 10&lt;sup&gt;9&lt;/sup&gt; bp sequences and 31% GC content. Then, MetaPhlAn2 analysis of NGS raw data predicted three phytoplasma species inside the infected Russian olive tree DNA. Finally, the genomes of ‘&lt;em&gt;Ca&lt;/em&gt;. Phytoplasma asteris’ with 833199 bp length, ‘&lt;em&gt;Ca&lt;/em&gt;. Phytoplasma australiense’ with 762261 bp length and ‘&lt;em&gt;Ca&lt;/em&gt;. Phytoplasma mali’ with 518942 bp length were obtained from NGS raw data. This study represents the first report of ‘&lt;em&gt;Ca&lt;/em&gt;. Phytoplasma australiense’ and ‘&lt;em&gt;Ca&lt;/em&gt;. Phytoplasma mali’ association with Russian olive witches’-broom disease using NGS analyses worldwide.