Biostatistical Aspects of Whole Genome Sequencing Studies: Preprocessing and Quality Control.

Rapid advances in high-throughput DNA sequencing technologies have enabled large-scale whole genome sequencing (WGS) studies. Before performing association analysis between phenotypes and genotypes, preprocessing and quality control (QC) of the raw sequence data need to be performed. Because many biostatisticians have not been working with WGS data so far, we first sketch Illumina's short-read sequencing technology. Second, we explain the general preprocessing pipeline for WGS studies. Third, we provide an overview of important QC metrics, which are applied to WGS data: on the raw data, after mapping and alignment, after variant calling, and after multisample variant calling. Fourth, we illustrate the QC with the data from the GENEtic SequencIng Study Hamburg-Davos (GENESIS-HD), a study involving more than 9000 human whole genomes. All samples were sequenced on an Illumina NovaSeq 6000 with an average coverage of 35× using a PCR-free protocol. For QC, one genome in a bottle (GIAB) trio was sequenced in four replicates, and one GIAB sample was successfully sequenced 70 times in different runs. Fifth, we provide empirical data on the compression of raw data using the DRAGEN original read archive (ORA). The most important quality metrics in the application were genetic similarity, sample cross-contamination, deviations from the expected Het/Hom ratio, relatedness, and coverage. The compression ratio of the raw files using DRAGEN ORA was 5.6:1, and compression time was linear by genome coverage. In summary, the preprocessing, joint calling, and QC of large WGS studies are feasible within a reasonable time, and efficient QC procedures are readily available.

Comparison of the performance of MiSeq and NovaSeq in oral microbiome study.

Objective: Next generation sequencing is commonly used to characterize the microbiome structure. MiSeq is most commonly used to analyze the microbiome due to its relatively long read length. Illumina also introduced the 250 × 2 chip for NovaSeq. The purpose of this study was to compare the performance of MiSeq and NovaSeq in the context of oral microbiome study. Methods: Total read count, read quality score, relative bacterial abundance, community diversity, and correlation between two platforms were analyzed. Phylogenetic trees were analyzed for Streptococcus and periodontopathogens. Results: NovaSeq produced significantly more read counts and assigned more operational taxonomic units (OTUs) compared to MiSeq. Community diversity was similar between MiSeq and NovaSeq. NovaSeq were able to detect more unique OTUs compared to MiSeq. When phylogenetic trees were constructed for Streptococcus and periodontopathogens, both platforms detected OTUs for most of the clades. Conclusion: Taken together, while both MiSeq and NovaSeq platforms effectively characterize the oral microbiome, NovaSeq outperformed MiSeq in terms of read counts and detection of unique OTUs, highlighting its potential as a valuable tool for large scale oral microbiome studies.

Transcriptome dataset for comparative analysis of differentially expressed genes between wild type and transgenic potato plants overexpressing Nuclear Factor Y subunit A7.

Nuclear Factor Y (NF-Y) is divided into three different types of subunits, A, B, and C. NF-Ys play crucial roles in plants for controlling gene expression associated with various developmental processes and abiotic stresses, but it is mostly unknown the downstream genes regulated by NF-Ys in plant. One of the potato NF-Y genes, StNF-YA7, increased potato's drought tolerance when overexpressed under the control of constitutive CaMV 35S promoter. Therefore, it was of interest what genes are regulated by the increased expression level of StNF-YA7. To investigate the downstream genes of StNF-YA7, the transcriptome sequencing was carried out for four potato lines, including Solanum tuberosum L 'Superior' as wild type (WT), empty vector control (VC), and two StNF-YA7 overexpressor lines (designated to StNF-YA7 #19 & #26). The RNA sequencing data was produced by the Illumina NovaSeq 6000 sequencing system. The number of total raw reads obtained from the RNA sequencing was 36.7 million for WT, 36.2 for VC, 29.3 for StNF-YA7 #19, and 29.5 million for StNF-YA7 #26, respectively. The length of total raw reads for each sample was between 5.92 Gb (StNF-YA7 #19) and 7.42 Gb (WT), and after filtering raw quality reads, the total length was between 5.81 Gb (StNF-YA7 #19) and 7.29 Gb (WT). Each filtered clear read set of four transcriptomes was mapped on the potato reference genome, SolTub_3.0, and the percentage of mapped reads ranged from 89.8 % (VC) to 90.3 % (WT). GC contents range between 43.01 % (StNF-YA7 #19) and 42.44 % (StNF-YA7 #26). Q20 quality score ranges between 98.63 % (StNF-YA7 #26) and 98.74 % (VC).

Genomic dataset of a multiple-drug resistant Pseudomonas sp. strain RAC1 isolated from a flacherie infected Nistari race of Bombyx mori L.

Species belonging to the genus Pseudomonas is a rod shaped Gram-negative bacteria emerged as an important silkworm pathogen with broad-level multi-drug resistance. The extensive usage of antimicrobials in sericulture farming is gradually leading to the emergence of multi-drug resistance (MDR) strains, posing a significant threat to the well-being of both Bombyx mori L. and serifarmers. Pseudomonas spp. with MDR level may gets transmitted from the infected silkworm to human handlers either via direct contact or through contaminated feces. To understand the emerging concern of antimicrobial resistance (AMR) in Pseudomonas spp. provides insights into their genomic information. Here, we present the draft genome sequence data of Pseudomonas sp. strain RAC1 isolated from a flacherie infected Nistari race of Bombyx mori L. from the silkworm rearing house of Raiganj University, India and sequenced using the Illumina NovaSeq 6000 platform. The estimated genome size of the strain was 4494347 bp with a G + C content of 63.5%. The de novo assembly of the genome generated 38 contigs with an N50 of 200 kb. Our data might help to reveal the genetic diversity, underlying mechanisms of AMR and virulence potential of Pseudomonas spp. This draft-genome shotgun project has been deposited under the NCBI GenBank accession number NZ_JAUTXS000000000.

Comparison of capture-based mtDNA sequencing performance between MGI and illumina sequencing platforms in various sample types.

BACKGROUND: Mitochondrial genome abnormalities can lead to mitochondrial dysfunction, which in turn affects cellular biology and is closely associated with the development of various diseases. The demand for mitochondrial DNA (mtDNA) sequencing has been increasing, and Illumina and MGI are two commonly used sequencing platforms for capture-based mtDNA sequencing. However, there is currently no systematic comparison of mtDNA sequencing performance between these two platforms. To address this gap, we compared the performance of capture-based mtDNA sequencing between Illumina's NovaSeq 6000 and MGI's DNBSEQ-T7 using tissue, peripheral blood mononuclear cell (PBMC), formalin-fixed paraffin-embedded (FFPE) tissue, plasma, and urine samples. RESULTS: Our analysis indicated a high degree of consistency between the two platforms in terms of sequencing quality, GC content, and coverage. In terms of data output, DNBSEQ-T7 showed higher rates of clean data and duplication compared to NovaSeq 6000. Conversely, the amount of mtDNA data obtained by per gigabyte sequencing data was significantly lower in DNBSEQ-T7 compared to NovaSeq 6000. In terms of detection mtDNA copy number, both platforms exhibited good consistency in all sample types. When it comes to detection of mtDNA mutations in tissue, FFPE, and PBMC samples, the two platforms also showed good consistency. However, when detecting mtDNA mutations in plasma and urine samples, significant differenceof themutation number detected was observed between the two platforms. For mtDNA sequencing of plasma and urine samples, a wider range of DNA fragment size distribution was found in NovaSeq 6000 when compared to DNBSEQ-T7. Additionally, two platforms exhibited different characteristics of mtDNA fragment end preference. CONCLUSIONS: In summary, the two platforms generally showed good consistency in capture-based mtDNA sequencing. However, it is necessary to consider the data preferences generated by two sequencing platforms when plasma and urine samples were analyzed.

Different Responses of Bacteria and Archaea to Environmental Variables in Brines of the Mahai Potash Mine, Qinghai-Tibet Plateau.

Salt mines feature both autochthonous and allochthonous microbial communities introduced by industrialization. It is important to generate the information on the diversity of the microbial communities present in the salt mines and how they are shaped by the environment representing ecological diversification. Brine from Mahai potash mine (Qianghai, China), an extreme hypersaline environment, is used to produce potash salts for hundreds of millions of people. However, halophiles preserved in this niche during deposition are still unknown. In this study, using high-throughput 16S rRNA gene amplicon sequencing and estimation of physicochemical variables, we examined brine samples collected from locations with the gradient of industrial activity intensity and discrete hydrochemical compositions in the Mahai potash mine. Our findings revealed a highly diverse bacterial community, mainly composed of Pseudomonadota in the hypersaline brines from the industrial area, whereas in the natural brine collected from the upstream Mahai salt lake, most of the 16S rRNA gene reads were assigned to Bacteroidota. Halobacteria and halophilic methanogens dominated archaeal populations. Furthermore, we discovered that in the Mahai potash mining area, bacterial communities tended to respond to anthropogenic influences. In contrast, archaeal diversity and compositions were primarily shaped by the chemical properties of the hypersaline brines. Conspicuously, distinct methanogenic communities were discovered in sets of samples with varying ionic compositions, indicating their strong sensitivity to the brine hydrochemical alterations. Our findings provide the first taxonomic snapshot of microbial communities from the Mahai potash mine and reveal the different responses of bacteria and archaea to environmental variations in this high-altitude aquatic ecosystem.

Bamboo charcoal affects soil properties and bacterial community in tea plantations.

Bamboo charcoal, a type of manufactured biochar, is produced by pyrolyzing bamboo residue under anoxic conditions. Its beneficial properties in absorption, catalyst support, and agricultural function have attracted significant attention; however, relatively few studies have examined its effects on the soil microbiota. In this study, we analyzed the effects of bamboo charcoal on soil physicochemical properties, enzymes, and microbial community structure in tea plantations and investigated the optimal amount of bamboo charcoal to be added to organic fertilizer. The results show that bamboo charcoal can further increase soil available nitrogen, total and available phosphorus and potassium, organic carbon content, pH, and urease activity. However, only the combined use of bamboo charcoal and organic fertilizer significantly increased total nitrogen, sucrase, and ß-glucosidase activities in the soil. Bamboo charcoal also significantly increased the Chao1 and Shannon indices of microbiota diversity in a concentration-dependent manner. The structure of the bacterial community changed significantly after the bamboo charcoal addition, with Proteobacteria, Actinobacteria, and Firmicutes increasing and Acidobacteria decreasing. This study provides fundamental insights into the suitability of bamboo charcoal application for the ecological remediation of diseased soils.

Systematic and benchmarking studies of pipelines for mammal WGBS data in the novel NGS platform.

BACKGROUND: Whole genome bisulfite sequencing (WGBS), possesses the aptitude to dissect methylation status at the nucleotide-level resolution of 5-methylcytosine (5-mC) on a genome-wide scale. It is a powerful technique for epigenome in various cell types, and tissues. As a recently established next-generation sequencing (NGS) platform, GenoLab M is a promising alternative platform. However, its comprehensive evaluation for WGBS has not been reported. We sequenced two bisulfite-converted mammal DNA in this research using our GenoLab M and NovaSeq 6000, respectively. Then, we systematically compared those data via four widely used WGBS tools (BSMAP, Bismark, BatMeth2, BS-Seeker2) and a new bisulfite-seq tool (BSBolt). We interrogated their computational time, genome depth and coverage, and evaluated their percentage of methylated Cs. RESULT: Here, benchmarking a combination of pre- and post-processing methods, we found that trimming improved the performance of mapping efficiency in eight datasets. The data from two platforms uncovered ~ 80% of CpG sites genome-wide in the human cell line. Those data sequenced by GenoLab M achieved a far lower proportion of duplicates (~ 5.5%). Among pipelines, BSMAP provided an intriguing representation of 5-mC distribution at CpG sites with 5-mC levels > ~ 78% in datasets from human cell lines, especially in the GenoLab M. BSMAP performed more advantages in running time, uniquely mapped reads percentages, genomic coverage, and quantitative accuracy. Finally, compared with the previous methylation pattern of human cell line and mouse tissue, we confirmed that the data from GenoLab M performed similar consistency and accuracy in methylation levels of CpG sites with that from NovaSeq 6000. CONCLUSION: Together we confirmed that GenoLab M was a qualified NGS platform for WGBS with high performance. Our results showed that BSMAP was the suitable pipeline that allowed for WGBS studies on the GenoLab M platform.

Exploring the value of microorganisms in the appendix for inferring postmortem interval in Sprague-Dawley rats using high-throughput sequencing.

Various microorganisms play an important role in daily functions in the body and continue to flourish after death. Our prior investigation using frozen cadavers revealed that the appendix, rather than the transverse colon, was a superior sampling site for intestinal bacteria because the appendiceal flora had higher diversity than that in the transverse colon in the majority of experimental periods after death. We sought to explore out more about whether the appendicular flora is significantly related to postmortem interval (PMI) at natural temperatures following the host's death. In this work, we employed high-throughput sequencing to evaluate the contents of rats' appendices within 2 weeks after death and then utilized the random forest algorithm to build a PMI prediction model after completing basic visual analyses on the sequencing data. The findings revealed that Firmicutes was the absolute dominant species of appendicular flora; alpha-diversity of appendix flora first increased and then decreased, with the highest point appearing at 36 h after death; and the primary metabolic functions were carbohydrate metabolism, amino acid metabolism, as well as cofactors and vitamin metabolism. Finally, a random forest regression model for PMI prediction was built by the training data at the family level, with the mean absolute error of 10.27 h for prediction within 14 days postmortem, and the test set data subsequently proved the model's reliability. Changes in appendicular flora were strongly related to the PMI following rats' deaths, so we have reason to believe that the appendicular flora is valuable in predicting PMI.

A long-read and short-read transcriptomics approach provides the first high-quality reference transcriptome and genome annotation for Pseudotsuga menziesii (Douglas-fir).

Douglas-fir (Pseudotsuga menziesii) is native to western North America. It grows in a wide range of environmental conditions and is an important timber tree. Although there are several studies on the gene expression responses of Douglas-fir to abiotic cues, the absence of high-quality transcriptome and genome data is a barrier to further investigation. Like for most conifers, the available transcriptome and genome reference dataset for Douglas-fir remains fragmented and requires refinement. We aimed to generate a highly accurate, and complete reference transcriptome and genome annotation. We deep-sequenced the transcriptome of Douglas-fir needles from seedlings that were grown under nonstress control conditions or a combination of heat and drought stress conditions using long-read (LR) and short-read (SR) sequencing platforms. We used 2 computational approaches, namely de novo and genome-guided LR transcriptome assembly. Using the LR de novo assembly, we identified 1.3X more high-quality transcripts, 1.85X more "complete" genes, and 2.7X more functionally annotated genes compared to the genome-guided assembly approach. We predicted 666 long noncoding RNAs and 12,778 unique protein-coding transcripts including 2,016 putative transcription factors. We leveraged the LR de novo assembled transcriptome with paired-end SR and a published single-end SR transcriptome to generate an improved genome annotation. This was conducted with BRAKER2 and refined based on functional annotation, repetitive content, and transcriptome alignment. This high-quality genome annotation has 51,419 unique gene models derived from 322,631 initial predictions. Overall, our informatics approach provides a new reference Douglas-fir transcriptome assembly and genome annotation with considerably improved completeness and functional annotation.

Comparison of genotyping by sequencing procedures to determine population genetic structure.

Advancements in technology over the past few decades have resulted in the development of genome sequencing at lower costs. Protocols, costs, and the amount of data produced by different sequencing technologies are highly variable. Ion Torrent and Illumina sequencing instruments are two sequencing technologies which use very similar library preparation procedures. Enzymatic combinations can be changed in genotyping by sequencing (GbS) library protocols without significant adjustments. To compare the outputs from two different GbS procedures, we sequenced samples of two sister species of yellow-nosed albatross collected at multiple geographic locations. The data sets involving different sequencing instruments and enzymatic combinations were analysed using the Stacks pipeline and aligned to the same reference genome. Both procedures identified the same genetic clusters separating Atlantic and Indian yellow-nosed albatross and substructure within Indian yellow-nosed albatross.

The complete chloroplast genome sequencing data of Juniperus sabina L. (Cupressaceae Bartl.) from Kazakhstan.

The genus Juniperus L. (Cupressaceae Bartl.) is consisting of about 75 species that are divided into sections Caryocedrus Endlicher, Sabina (Miller) Spach, and Juniperus (syn: sect. Oxycedrus Spach). Juniperus sabina L. from section Sabina is an important shrub for the maintenance of the ecosystem in mountainous regions and a source of medicinal compounds. The species is monoecious, rarely dioecious, and distributed in Europe, Central Asia, China, and Mongolia. The goal of the present study was to sequence and reconstruct the complete chloroplast genome of J. sabina. De novo chloroplast (cp) genome assembly for J. sabina was conducted using Illumina paired-end reads. The assembled cp genome size is 127,646 bp in length and has a typical circular DNA molecule. The genome encodes 118 genes, including 82 protein-coding genes, 32 tRNA genes, and 4 rRNA genes, the overall GC content is 34,36%. The complete cp genome nucleotide sequence of J. sabina was deposited to the NCBI (National Center for Biotechnology Information) under accession number OL467323. The raw data in fastq format was deposited to the NCBI sequence read archive under accession number SRR21515769.

Effects of Insect-Resistant Maize HGK60 on Community Diversity of Bacteria and Fungi in Rhizosphere Soil.

The influence of biotech crops on microbial communities in rhizosphere soil is an important issue in biosafety assessments. The transgenic maize HGK60 harboring the Bt cry1Ah gene enhanced the resistance to lepidopteran pests, while the ecological risk of HGK60 maize on rhizosphere microorganisms is unclear. In this study, we comprehensively analyzed the diversity and composition of bacterial and fungal communities in the rhizosphere soil around Bt maize HGK60 and the near-isogenic non-Bt maize ZD958 at four growth stages via a high-throughput sequencing technique. The results showed that HGK60 maize unleashed temporary effects on the bacterial and fungal diversity and richness during the study plant's development, which would be restored after one cycle of plant cultivation due to the application of the same agricultural management. The differences of bacterial and fungal communities were marked by seasonality, while the different growth stage was the important factor as opposed to the cultivar contributing to the shifts in the bacterial and fungal communities' structure. This study will provide useful information regarding the impact of Bt transgenic maize on the soil microbiome and a theoretical basis for the development of a safety assessment approach for Bt maize in China.

Whole-Genome Survey and Microsatellite Marker Detection of Antarctic Crocodile Icefish, Chionobathyscus dewitti.

The crocodile icefish, Chionobathyscus dewitti, belonging to the family Channichthyidae, is an endemic species of the Southern Ocean. The study of its biological features and genetics is challenging as the fish inhabits the deep sea around Antarctic waters. The icefish, the sole cryopelagic species, shows unique physiological and genetic features, unlike other teleosts. It lacks hemoglobin and has evolved antifreeze proteins. Here, we report the genome sequencing data of crocodile icefish produced using the Illumina Novaseq 6000 platform. The estimated genome size was 0.88 Gb with a K-value of 19, and the unique sequence, heterozygosity, error, and duplication rates were 57.4%, 0.421%, 0.317%, and 0.738%, respectively. A genome assembly of 880.69 Mb, with an N50 scaffold length of 2401 bp, was conducted. We identified 2,252,265 microsatellite motifs from the genome assembly data, and dinucleotide repeats (1,920,127; 85.25%) had the highest rate. We selected 84 primer pairs from the genome survey assembly and randomly selected 30 primer pairs for validation. As a result, 15 primer pairs were validated as microsatellite markers.

Development of de-novo transcriptome assembly and SSRs in allohexaploid Brassica with functional annotations and identification of heat-shock proteins for thermotolerance.

Crop Brassicas contain monogenomic and digenomic species, with no evidence of a trigenomic Brassica in nature. Through somatic fusion (Sinapis alba + B. juncea), a novel allohexaploid trigenomic Brassica (H1 = AABBSS; 2n = 60) was produced and used for transcriptome analysis to uncover genes for thermotolerance, annotations, and microsatellite markers for future molecular breeding. Illumina Novaseq 6000 generated a total of 76,055,546 paired-end raw reads, which were used for de-novo assembly, resulting in the development of 486,066 transcripts. A total of 133,167 coding sequences (CDSs) were predicted from transcripts with a mean length of 507.12 bp and 46.15% GC content. The BLASTX search of CDSs against public protein databases showed a maximum of 126,131 (94.72%) and a minimum of 29,810 (22.39%) positive hits. Furthermore, 953,773 gene ontology (GO) terms were found in 77,613 (58.28%) CDSs, which were divided into biological processes (49.06%), cellular components (31.67%), and molecular functions (19.27%). CDSs were assigned to 144 pathways by a pathway study using the KEGG database and 1,551 pathways by a similar analysis using the Reactome database. Further investigation led to the discovery of genes encoding over 2,000 heat shock proteins (HSPs). The discovery of a large number of HSPs in allohexaploid Brassica validated our earlier findings for heat tolerance at seed maturity. A total of 15,736 SSRs have been found in 13,595 CDSs, with an average of one SSR per 4.29 kb length and an SSR frequency of 11.82%. The first transcriptome assembly of a meiotically stable allohexaploid Brassica has been given in this article, along with functional annotations and the presence of SSRs, which could aid future genetic and genomic studies.

Analysis of the blood bacterial composition of patients with acute coronary syndrome and chronic coronary syndrome.

Emerging evidence revealed that the blood microbiota plays a role in several non-communicable diseases, including cardiovascular disease. However, the role of circulating microbes in atherosclerosis remains understudied. To test this hypothesis, we performed this study to investigate the microbial profile in the blood of Chines atherosclerosis volunteers. A total of seventy Acute Coronary Syndrome patients, seventy Chronic Coronary Syndrome patients, and seventy healthy individuals were examined using high-throughput Illumina Novaseq targeting the V3-V4 regions of the 16S rRNA gene. The relationship between atherosclerosis and blood microbiome, clinical variables, and their functional pathways were also investigated. Our study observed significantly higher alpha diversity indices (Chao1, p = 0.001, and Shannon, p = 0.004) in the acute coronary syndrome group compared with chronic coronary syndrome and healthy group, although a significantly lower alpha diversity was observed in the chronic coronary syndrome compared to acute coronary syndrome and healthy group. Beta diversity based on principal coordinate analysis demonstrated a major separation among the three groups. In addition, using linear discriminant analysis, a significant distinct taxon such as Actinobacteria _ phylum, and Staphylococcus_ genus in the healthy group; Firmicutes_ phylum, and Lactobacillus_ genus in the chronic coronary syndrome group, and Proteobacteria and Acidobacteriota _ phyla in acute coronary syndrome group were observed among three groups. Clusters of Orthologous Genes grouped and Kyoto Encyclopedia of Genes and Genomes pathways suggested a significant variation among all groups (p < 0.05). The blood microbiota analysis provides potential biomarkers for the detection of coronary syndromes in this population.

Microbial community succession in the fermentation of Qingzhuan tea at various temperatures and their correlations with the quality formation.

With the aim to reveal the microbial community succession at various temperatures in the fermentation of Qingzhuan tea (QZT), the Illumina NovaSeq sequencing was carried out to analyze bacterial and fungal community structure in tea samples collected from the fermentation set at various temperatures, i.e., 25 °C, 30 °C, 37 °C, 45 °C, 55 °C, and room temperature. The results showed that fermentation temperature profoundly affected the microbial community succession in the QZT fermentation. Microbial richness and community diversity decreased along with the increase of fermentation temperature. Despite the differences between microorganisms and their metabolic types among various temperatures, most bacteria and fungi showed positive correlations at the genera level. Klebsiella, Paenibacillus, Cohnella, and Pantoea were confirmed as the main bacterial genera, and Aspergillus and Cyberlindnera were the main fungal genera in QZT fermentation. The microbial genera (i.e. Aspergillus, Rhizomucor, Thermomyces, Ralstonia, Castellaniella, and Vibrio) were positively correlated with fermentation temperature (P < 0.05), while Klebsiella, Paenibacillus, and Aspergillus had good adaptability at different temperatures. Conversely, Pantoea and Cyberlindnera were only suitable for low temperature (≤37 °C) growth, and Thermomyces was only suitable for high temperature (>37 °C) growth. Aspergillus had a significant positive correlation with tea aroma quality (r = 0.64, p < 0.05). This study would help to understand the formation mechanism of QZT from microflora perspective.

Identification and profiling of the community structure and potential function of bacteria from the fruiting bodies of Sanghuangporus vaninii.

Sanghuangporus sp., a medicinal and edible homologous macrofungus known as 'forest gold', which has good effects on antitumor, hypolipidemia and the treatment of gynecological diseases. However, the natural resources of fruiting body are on the verge of depletion due to its long growth cycle and over exploitation. The growth and metabolism of macrofungi are known to depend on the diverse bacterial community. Here, we characterized the diversity and potential function of bacteria inhabiting in the fruiting body of the most widely applied S. vaninii using a combination method of high-throughput sequencing with pure culturing for the first time, and tested the biological activities of bacterial isolates, of which Illumina NovaSeq provided a more comprehensive results on the bacterial community structure. Total 33 phyla, 82 classes, 195 orders, 355 families, 601 genera and 679 species were identified in the fruiting body, and our results revealed that the community was predominated by the common Proteobacteria, Gammaproteobacteria, Burkholderiales, Methylophilaceae (partly consistent with pure-culturing findings), and was dominated by the genera of distinctive Methylotenera and Methylomonas (yet-uncultured taxa). Simultaneously, the functional analysis showed that companion bacteria were involved in the pathways of carbohydrate transport and metabolism, metabolism of terpenoids and polyketides, cell wall/membrane/envelope biogenesis, etc. Hence, it was inferred that bacteria associated with fruiting body may have the potential to adjust the growth, development and active metabolite production of host S. vaninii combined with the tested results of indole-3-acetic acid and total antioxidant capacity. Altogether, this report first provided new findings which can be inspiring for further in-depth studies to exploit bioactive microbial resources for increased production of Sanghuangporus, as well as to explore the relationship between medicinal macrofungi and their associated endophytes.

Low Input Genome-Wide DNA Methylation Analysis with Minimal Library Amplification.

Whole genome bisulfite sequencing (WGBS) is a high-throughput DNA sequencing-based technique that is used to determine genome-wide DNA methylation patterns at base resolution. Library construction by post-bisulfite adaptor tagging (PBAT ) extends the application of WGBS to several hundred cells and minimizes the required number of library amplification cycles. We herein describe a PBAT protocol to prepare WGBS libraries from 200 cells and introduce the outline of a downstream bioinformatic analysis. The prepared library can typically generate 800 million sequencing reads, which is sufficient to cover the human and mouse genomes approximately 15 times, using the Illumina NovaSeq 6000 sequencing system.