Single-molecule Real-time Research Articles

SARS-CoV-2 Co-Infections and Recombinations Identified by Long-Read Single-Molecule Real-Time Sequencing.

Co-infection with at least 2 strains of virus is the prerequisite for recombination, one of the means of genetic diversification. Little is known about the prevalence of these events in SARS-CoV-2, partly because it is difficult to detect them. We used long-read PacBio single-molecule real-time (SMRT) sequencing technology to sequence whole genomes and targeted regions for haplotyping. We identified 17 co-infections with SARS-CoV-2 strains belonging to different clades in 6829 samples sequenced between January and October, 2022 (prevalence 0.25%). There were 3 Delta/Omicron co-infections and 14 Omicron/Omicron co-infections (4 cases of 21K/21L, 1 case of 21L/22A, 2 cases of 21L/22B, 4 cases of 22A/22B, 2 cases of 22B/22C and 1 case of 22B/22E). Four of these patients (24%) also harbored recombinant minor haplotypes, including one with a recombinant virus that was selected in the viral quasispecies over the course of his chronic infection. While co-infections remain rare among SARS-CoV-2-infected individuals, long-read SMRT sequencing is a useful tool for detecting them as well as recombinant events, providing the basis for assessing their clinical impact, and a precise indicator of epidemic evolution. IMPORTANCE SARS-CoV-2 variants have been responsible for the successive waves of infection over the 3 years of pandemic. While co-infection followed by recombination is one driver of virus evolution, there have been few reports of co-infections, mainly between Delta and Omicron variants or between the first 2 Omicron variants 21K_BA.1 and 21L_BA.2. The 17 co-infections we detected during 2022 included cases with the recent clades of Omicron 22A, 22B, 22C, and 22E; 24% harbored recombinant variants. This study shows that long-read SMRT sequencing is well suited to SARS-CoV-2 genomic surveillance.

The frequency of HKαα allele in silent deletional α-thalassemia carriers in the Yulin region of southern China using the third-generation sequencing

Objectivesα-thalassemia is relatively prevalent in Yulin Region in southern China. In order to accurately detect α-globin gene aberrations for genetic counseling, the prevalence of HKαα (Hong Kong αα) allele in this subpopulation of silent deletional α-thalassemia were examined. Materials and methodsA total of 1845 subjects were selected in Yulin Region from January 2021 to March 2021. Peripheral blood was collected from each participant for routine genetic analysis of thalassemia. The HKαα allele was determined using the Single-molecule real-time (SMRT) technology for samples with -α3.7/αα, βN/βN genotype. ResultsTwo samples were identified with HKαα allele from 100 samples with -α3.7/αα, βN/βN genotype. The frequency of HKαα allele was 2.0 % (2/100) in -α3.7/αα, βN/βN carriers in Yulin Region. One sample was identified with a novel variant of the α-globin gene cluster named αHKαα by SMRT technology. One rare HBA2 variant and six HBB variants were found by SMRT technology, including -α3.7/HBA2:c.300 + 34G > A, HBB:c.316-45G > C/βN, HBB:c.315 + 180 T > C/βN, HBB:c.316-179A > C/βN. ConclusionA certain proportion of HKαα allele had been detected in Yulin Region. SMRT technology plays a crucial role for improving the diagnostic accuracy and positive detection rate of thalassemia. The completion of this study has great meaning for strengthening the prevention and control of thalassemia in Yulin Region.

Chromosome-level assembly of Dictyophora rubrovolvata genome using third-generation DNA sequencing and Hi-C analysis.

Dictyophora rubrovolvata, a rare edible mushroom with both nutritional and medicinal values, was regarded as the "queen of the mushroom" for its attractive appearance. Dictyophora rubrovolvata has been widely cultivated in China in recent years, and many researchers were focusing on its nutrition, culture condition, and artificial cultivation. Due to a lack of genomic information, research on bioactive substances, cross breeding, lignocellulose degradation, and molecular biology is limited. In this study, we report a chromosome-level reference genome of D. rubrovolvata using the PacBio single-molecule real-time-sequencing technique and high-throughput chromosome conformation capture (Hi-C) technologies. A total of 1.83 Gb circular consensus sequencing reads representing ∼983.34 coverage of the D. rubrovolvata genome were generated. The final genome was assembled into 136 contigs with a total length of 32.89 Mb. The scaffold and contig N50 length were 2.71 and 2.48 Mb, respectively. After chromosome-level scaffolding, 11 chromosomes with a total length of 28.24 Mb were constructed. Genome annotation further revealed that 9.86% of the genome was composed of repetitive sequences, and a total of 508 noncoding RNA (rRNA: 329, tRNA: 150, ncRNA: 29) were annotated. In addition, 9,725 protein-coding genes were predicted, among which 8,830 (90.79%) genes were predicted using homology or RNA-seq. Benchmarking Universal Single-Copy Orthologs results further revealed that there were 80.34% complete single-copy fungal orthologs. In this study, a total of 360 genes were annotated as belonging to the carbohydrate-active enzymes family. Further analysis also predicted 425 cytochromes P450 genes, which can be classified into 41 families. This highly accurate, chromosome-level reference genome of D. rubrovolvata will provide essential genomic information for understanding the molecular mechanism in its fruiting body formation during morphological development and facilitate the exploitation of medicinal compounds produced by this mushroom.

Full-length transcriptome analysis of Zanthoxylum nitidum (Roxb.) DC.

Zanthoxylum nitidum (Roxb.) DC. (Z. nitidum) is a type of Chinese Dao-di herb, also called Liangmianzhen, which is widely used to treat arthralgia, rheumatic arthralgia, and stomach pain. However, genomic resources for Z. nitidum are still scarce. This study provides transcriptomic resources for Z. nitidum by applying single-molecule real-time (SMRT) sequencing technology. In total, 456,109 circular consensus sequencing (CCS) reads were generated with a mean length of 2,216 bp from Z. nitidum roots, old stems, young branches, leaves, flowers, and fruits. Of these total reads, 353,932 were full-length nonchimeric (FLNC) reads with an average length of 1,996 bp. A total of 16,163 transcripts with a mean length of 1,171 bp were acquired. Of these transcripts, 14,231 (88%) were successfully annotated using public databases. Across all the 16,163 transcripts, we identified 6,255 long non-coding RNAs (lncRNAs) and 22,780 simple sequence repeats (SSRs). Furthermore, 3,482 transcription factors were identified. Among the SSR loci, 1-3 nucleotide repeats were dominant, occupying 99.36% of the total SSR loci, with mono-, di-, and tri-nucleotide repeats accounting for 61.80%, 19.89%, and 5.02% of the total SSR loci, respectively. A total of 36 out of 100 randomly selected primer pairs were verified to be positive, 20 of which showed polymorphism. These findings enrich the genetic resources available for facilitating future studies and research on relevant topics such as population genetics in Z. nitidum.

MaxDEL: Accurate and Efficient Calling of Genomic Deletions from Single Molecular Real-time Sequencing Using Integrated Method

Background: Single-molecule real-time (SMRT) sequencing data are characterized by long read length and high read depth. Compared to next-generation sequencing (NGS), SMRT sequencing data can present more structural variations (SVs) and have greater advantages in calling variation. However, there are high sequencing errors and noises in SMRT sequencing data, which causes inaccuracy in calling SVs from sequencing data. Most existing tools cannot overcome sequencing errors and detect genomic deletions. Objective: In this investigation, we propose a new method for calling deletions from SMRT sequencing data called MaxDEL. Methods: Firstly, MaxDEL uses a machine learning method to calibrate the deletion regions from the variant call format (VCF) file. Secondly, it develops a novel feature visualization method to convert the variant features to images and uses these images to accurately call the deletions based on a convolutional neural network (CNN). Results: The result shows that MaxDEL performs better in terms of accuracy and recall for calling variants when compared to existing methods in both real data and simulative data. Conclusion: MaxDEL can effectively overcome SMRT sequencing data's noise and integrate new machine learning and deep learning technologies. The method can capture the variant features of the deletions and establish the learning model between images and gene data. In our experiment, the MaxDEL method is superior to NextSV, SVIM, Sniffles, Picky and SMRT-SV, especially in recall and F1-score.

Identification of four novel large deletions and complex variants in the α-globin locus in Chinese population

BackgroundAt present, the methods generally used to detect α-thalassemia mutations are confined to detecting common mutations, which may lead to misdiagnosis or missed diagnosis. The single-molecule real-time (SMRT) sequencing enables long-read single-molecule sequencing with high detection accuracy, and long-length DNA chain reads in high-fidelity read mode. This study aimed to identify novel large deletions and complex variants in the α-globin locus in Chinese population.MethodsWe used SMRT sequencing to detect rare and complex variants in the α-globin locus in four individuals whose hematological data indicated microcytic hypochromic anemia. However, the conventional thalassemia detection result was negative. Multiplex ligation-dependent probe amplification and droplet digital polymerase chain reaction were used to confirm SMRT sequencing results.ResultsFour novel large deletions were observed ranging from 23 to 81 kb in the α-globin locus. One patient also had a duplication of upstream of HBZ in the deletional region, while another, with a 27.31-kb deletion on chromosome 16 (hg 38), had abnormal hemoglobin Siriraj (Hb Siriraj).ConclusionWe first identified the four novel deletions in the α-globin locus using SMRT sequencing. Considering that the conventional methods might lead to misdiagnosis or missed diagnosis, SMRT sequencing proved to be an excellent method to discover rare and complex variants in thalassemia, especially in prenatal diagnosis.

Bacterial microbiota and metabolic characteristics of traditional Zhangjiajie Jiuqu aromatic ester vinegar in China

Aged vinegars are popular acidic condiments worldwide, used for their flavour and nutritional value; they also have a rich microbial diversity. This study combined the ultra-performance liquid chromatography-quadrupole-time of flight mass spectrometry (UPLC-Q-TOF-MS) technology with Pacific Biosciences (PacBio) single-molecule real-time (SMRT) sequencing to analyze the microbial population and associated metabolism in six traditional Zhangjiajie Jiuqu aromatic ester vinegars (JAEVs) in China. A total of 265 species and subspecies from 170 bacterial genera and 14 phyla were identified; among these Lactococcus lactis (18.23%) and Acinetobacter johnsonii (17.18%) were the predominant species. A total of 1352 metabolites were detected in the metabolomes of all samples. In addition, occurrence of 40 metabolites detected were significantly different among six samples. Based on correlation analyses between the prevalent bacterial species and main differentially recorded metabolites, we conclude that there were significant correlations (P < 0.05) between 20 metabolites and five bacterial species. Limosilactobacillus fermentum and Hyphomicrobium zavarzinii were the main regulatory strains caused differential metabolites. Some potential functional metabolites were also identified, mainly chlorogenic acid and trans-piceid. This study revealed the relationship between the unique flavour of JAEVs and the microbial community, and provided a basis for industrial production of traditional aged vinegar.

Division of developmental phases of freshwater leech Whitmania pigra and key genes related to neurogenesis revealed by whole genome and transcriptome analysis

The freshwater leech Whitmania pigra (W. pigra) Whitman (Annelida phylum) is a model organism for neurodevelopmental studies. However, molecular biology research on its embryonic development is still scarce. Here, we described a series of developmental stages of the W. pigra embryos and defined five broad stages of embryogenesis: cleavage stages, blastocyst stage, gastrula stage, organogenesis and refinement, juvenile. We obtained a total of 239.64 Gb transcriptome data of eight representative developmental phases of embryos (from blastocyst stage to maturity), which was then assembled into 21,482 unigenes according to our reference genome sequenced by single-molecule real-time (SMRT) long-read sequencing. We found 3114 genes differentially expressed during the eight phases with phase-specific expression pattern. Using a comprehensive transcriptome dataset, we demonstrated that 57, 49 and 77 DEGs were respectively related to morphogenesis, signal pathways and neurogenesis. 49 DEGs related to signal pathways included 30 wnt genes, 14 notch genes, and 5 hedgehog genes. In particular, we found a cluster consisting of 7 genes related to signal pathways as well as synapses, which were essential for regulating embryonic development. Eight genes cooperatively participated in regulating neurogenesis. Our results reveal the whole picture of W. pigra development mechanism from the perspective of transcriptome and provide new clues for organogenesis and neurodevelopmental studies of Annelida species.

Advances in Nanopore and Photoelectron-Based High-Throughput Sequencing Technology for Single-Molecule Sequencing

Next-Generation sequencing is a widespread technology of sequencing which is used in genomics research and biotechnology, although it comes with many shortcomings such as short read length or polymerase chain reactions (PCR) bias. In this review, a brief introduction is provided to the Helicos Biosciences true single-molecule sequencing (tSMS), Pacific Biosciences single-molecule real-time (SMRT) sequencing and Oxford Nanopore single-molecule sequencing technologies. In comparison with both first- and second-generation sequencing technologies, third-generation sequencing such as nanopore and photoelectron high-throughput sequencing technologies have the advantages of ultra-long read length and no PCR bias, which help to study repetitive regions of the genome and complex structural variants, thus correctly identifying potential gene mutations in various diseases, and are important for diagnosis and treatment of various diseases. This paper focuses on the sequencing principles of nanopore and photoelectron high-throughput sequencing technologies as well as their applications in cancer, pathogen detection and genetic diseases.

Distribution characteristics and clinical phenotype analyses of hemoglobin variants in the Z region of Central Guangxi, Southern China

ABSTRACTObjectiveTo investigate the molecular diagnosis of hemoglobin variants in Z region by Capillary electrophoresis in Central Guangxi, Southern China, and analyze their distribution and phenotypic characteristics, to provide a reference for clinical consultation and prenatal diagnosis for couples.MethodsA total of 23,709 subjects were collected for blood routine analysis, hemoglobin analysis, and common α- and β-globin gene loci in Chinese population. The hemoglobin electrophoresis components were divided into Zone 1-Zone 15 (Z1-Z15) by Capillary zone electrophoresis (CE). For samples not clearly detected by the conventional technology, Sanger sequencing, and multiplex ligation-dependent probe amplification (MLPA) were used. Single-molecule real-time (SMRT) sequencing technology was used to analyze rare-type genes in a sample with a structural variation.ResultsTen rare hemoglobin variants distributed in Z region were detected in 23,709 samples, including Hb Cibeles, which was reported for the first time in Asia; Hb J-Broussais, Hb G-Honolulu and J-Wenchang-Wuming, they were first reported in Guangxi; 1 case of Hb Anti-Lepore Liuzhou, which was a newly discovered hemoglobin variant; hemoglobin variants Hb G-Siriraj, Hb Handsworth, Hb Q-Thailand, Hb Ube-2, Hb NewYork were also detected.ConclusionThere are a few studies on rare hemoglobin variants in Z region in Southern China. Ten rare hemoglobin variants were found in this study. The hematological phenotype and component content of hemoglobin variants are related to the occurrence of thalassemia. This study enriched the data of rare hemoglobin variants in Southern China and provided a comprehensive data basis for prenatal diagnosis of hemoglobin variants in this area.

A Chinese SCA36 pedigree analysis of NOP56 expansion region based on long-read sequencing.

Introduction: Spinocerebellar ataxias 36 (SCA36) is the neurodegenerative disease caused by the GGCCTG Hexanucleotide repeat expansions in NOP56, which is too long to sequence using short-read sequencing. Single molecule real time (SMRT) sequencing can sequence across disease-causing repeat expansion. We report the first long-read sequencing data across the expansion region in SCA36. Methods: We collected and described the clinical manifestations and imaging features of Han Chinese pedigree with three generations of SCA36. Also, we focused on structural variation analysis for intron 1 of the NOP56 gene by SMRT sequencing in the assembled genome. Results: The main clinical features of this pedigree are late-onset ataxia symptoms, with a presymptomatic presence of affective and sleep disorders. In addition, the results of SMRT sequencing showed the specific repeat expansion region and demonstrated that the region was not composed of single GGCCTG hexanucleotides and there were random interruptions. Discussion: We extended the phenotypic spectrum of SCA36. We applied SMRT sequencing to reveal the correlation between genotype and phenotype of SCA36. Our findings indicated that long-read sequencing is well suited to characterize known repeat expansion.

Complete Genome Sequence of Staphylococcus epidermidis CCSM0287, Isolated from Healthy Facial Skin.

Staphylococcus epidermidis strain CCSM0287 was isolated from healthy facial skin. The complete genome of CCSM0287 was sequenced using a combination of Pacific Biosciences (PacBio) RS II single-molecule real-time (SMRT) and Illumina sequencing. The assembled 2.5-Mbp genome consisted of one chromosome and three plasmids.

Complete Genome Sequence of Staphylococcus capitis CCSM0123, Isolated from Healthy Facial Skin

Staphylococcus capitis strain CCSM0123 was isolated from healthy facial skin. The complete genome of CCSM0123 was sequenced using a combination of Pacific Biosciences (PacBio) RS II single-molecule real-time (SMRT) and Illumina sequencing. The assembled 2.5-Mbp genome consisted of one chromosome and four plasmids.

Hb Q-Thailand heterozygosity unlinked with the (–α4.2/) α+-thalassemia deletion allele identified by long-read SMRT sequencing: hematological and molecular analyses

ABSTRACTObjectiveIn the present study, two unrelated cases of Hb Q-Thailand heterozygosity unlinked with the (–α4.2/) α+-thalassemia deletion allele were identified by long-read single molecule real-time (SMRT) sequencing in southern China. The aim of this study was to report the hematological and molecular features as well as diagnostic aspects of the rare manifestation.MethodsHematological parameters and hemoglobin analysis results were recorded. A suspension array system for routine thalassemia genetic analysis and long-read SMRT sequencing were applied in parallel for thalassemia genotyping. Traditional methods, including Sanger sequencing, multiplex gap-polymerase chain reaction (gap-PCR) and multiplex ligation-dependent probe amplification (MLPA), were used together to confirm the thalassemia variants.ResultsLong-read SMRT sequencing was used to diagnose two Hb Q-Thailand heterozygous patients for whom the hemoglobin variant was unlinked to the (–α4.2/) allele for the first time. The hitherto undescribed genotypes were verified by traditional methods. Hematological parameters were compared with those of Hb Q-Thailand heterozygosity linked with the (–α4.2/) deletion allele in our study. For the positive control samples, long-read SMRT sequencing revealed a linkage relationship between the Hb Q-Thailand allele and the (–α4.2/) deletion allele.ConclusionsIdentification of the two patients confirms that the linkage relationship between the Hb Q-Thailand allele and the (−α4.2/) deletion allele is a common possibility but not a certainty. Remarkably, as it is superior to traditional methods, SMRT technology may eventually serve as a more comprehensive and precise method that holds promising prospects in clinical practice, especially for rare variants.

Profiling DNA Ligase Substrate Specificity with a Pacific Biosciences Single-Molecule Real-Time Sequencing Assay.

DNA ligases catalyze the joining of breaks in nucleic acid backbones and are essential enzymes for in vivo genome replication and repair across all domains of life. These enzymes are also critically important to in vitro manipulation of DNA in applications such as cloning, sequencing, and molecular diagnostics. DNA ligases generally catalyze the formation of a phosphodiester bond between an adjacent 5'-phosphate and 3'-hydroxyl in DNA, but they exhibit different substrate structure preferences, sequence-dependent biases in reaction kinetics, and variable tolerance for mismatched base pairs. Information on substrate structure and sequence specificity can inform both biological roles and molecular biology applications of these enzymes. Given the high complexity of DNA sequence space, testing DNA ligase substrate specificity on individual nucleic acid sequences in parallel rapidly becomes impractical when a large sequence space is investigated. Here, we describe methods for investigating DNA ligase sequence bias and mismatch discrimination using Pacific Biosciences Single-Molecule Real-Time (PacBio SMRT) sequencing technology. Through its rolling-circle amplification methodology, SMRT sequencing can give multiple reads of the same insert. This feature permits high-quality top- and bottom-strand consensus sequences to be determined while preserving information on top-bottom strand mismatches that can be obfuscated or lost when using other sequencing methods. Thus, PacBio SMRT sequencing is uniquely suited to measuring substrate bias and enzyme fidelity through multiplexing a diverse set of sequences in a single reaction. The protocols describe substrate synthesis, library preparation, and data analysis methods suitable for measuring fidelity and bias of DNA ligases. The methods are easily adapted to different nucleic acid substrate structures and can be used to characterize many enzymes under a variety of reaction conditions and sequence contexts in a rapid and high-throughput manner. © 2023 New England Biolabs and The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Preparation of overhang DNA substrates for ligation Basic Protocol 2: Preparation of ligation fidelity libraries Support Protocol 1: Preparation of ligation libraries for PacBio Sequel II sequencing Support Protocol 2: Loading and sequencing of a prepared library on the Sequel II instrument Basic Protocol 3: Computational processing of ligase fidelity sequencing data.

Characterization of Extensively Drug-Resistant Salmonella enterica Serovar Kentucky Sequence Type 198 Isolates from Chicken Meat Products in Xuancheng, China.

The purpose of this study was to characterize extensively drug-resistant Salmonella enterica serovar Kentucky sequence type 198 (ST198) isolates from chicken meat products. Ten S. Kentucky strains obtained from chicken meat products in Xuancheng, China, carried 12 to 17 resistance genes, such as blaCTX-M-55, rmtB, tet(A), floR, and fosA3, combined with mutations within gyrA (S83F and D87N) and parC (S80I), resulting in resistance to numerous antimicrobial agents, including the clinically important antibiotics cephalosporin, ciprofloxacin, tigecycline, and fosfomycin. These S. Kentucky isolates shared a close phylogenetic relationship (21 to 36 single-nucleotide polymorphisms [SNPs]) and showed close genetic relatedness to two human clinical isolates from China. Three S. Kentucky strains were subjected to whole-genome sequencing using Pacific Biosciences (PacBio) single-molecule real-time (SMRT) technology. All antimicrobial resistance genes were located on their chromosomes and clustered in one multiresistance region (MRR) and Salmonella genomic island (SGI) SGI1-K. The MRRs in three S. Kentucky strains were bounded by IS26 at both ends and were inserted downstream of the bcfABCDEFG cluster with 8-bp direct repeats. The MRRs were related to those of IncHI2 plasmids but differed by insertions, deletions, and rearrangements of multiple segments involving resistance genes and plasmid backbones. This finding suggests that the MRR fragment possibly originates from IncHI2 plasmids. Four SGI1-K variants with slight differences were identified in 10 S. Kentucky strains. Mobile elements, particularly IS26, play an essential role in forming distinct MRRs and SGI1-K structures. In conclusion, the emergence of extensively drug-resistant S. Kentucky ST198 strains containing numerous chromosomally located resistance genes is alarming and needs continued surveillance. IMPORTANCE Salmonella spp. are important foodborne pathogens, and multidrug-resistant (MDR) Salmonella strains have become a serious threat to clinical therapy. MDR S. Kentucky ST198 strains have been increasingly reported from various sources and have become a global risk. In this study, we described extensively drug-resistant S. Kentucky ST198 strains from chicken meat products from a city in China. Numerous resistance genes are clustered in the chromosomes of S. Kentucky ST198 strains, possibly acquired with the help of mobile elements. This would facilitate the spread of numerous resistance genes as intrinsic chromosomal genes within this global epidemic clone, with the potential to capture more resistance genes. The emergence and dissemination of extensively drug-resistant S. Kentucky ST198 pose a severe clinical and public health threat; therefore, continuous surveillance is warranted.

Actinobacillus pleuropneumoniae encodes multiple phase-variable DNA methyltransferases that control distinct phasevarions.

Actinobacillus pleuropneumoniae is the cause of porcine pleuropneumonia, a severe respiratory tract infection that is responsible for major economic losses to the swine industry. Many host-adapted bacterial pathogens encode systems known as phasevarions (phase-variable regulons). Phasevarions result from variable expression of cytoplasmic DNA methyltransferases. Variable expression results in genome-wide methylation differences within a bacterial population, leading to altered expression of multiple genes via epigenetic mechanisms. Our examination of a diverse population of A. pleuropneumoniae strains determined that Type I and Type III DNA methyltransferases with the hallmarks of phase variation were present in this species. We demonstrate that phase variation is occurring in these methyltransferases, and show associations between particular Type III methyltransferase alleles and serovar. Using Pacific BioSciences Single-Molecule, Real-Time (SMRT) sequencing and Oxford Nanopore sequencing, we demonstrate the presence of the first ever characterised phase-variable, cytosine-specific Type III DNA methyltransferase. Phase variation of distinct Type III DNA methyltransferase in A. pleuropneumoniae results in the regulation of distinct phasevarions, and in multiple phenotypic differences relevant to pathobiology. Our characterisation of these newly described phasevarions in A. pleuropneumoniae will aid in the selection of stably expressed antigens, and direct and inform development of a rationally designed subunit vaccine against this major veterinary pathogen.

Illumina RNA and SMRT Sequencing Reveals the Mechanism of Uptake and Transformation of Selenium Nanoparticles in Soybean Seedlings

Selenium (Se) is an essential element for mammals, and its deficiency in the diet is a global problem. Agronomic biofortification through exogenous Se provides a valuable strategy to enhance human Se intake. Selenium nanoparticles (SeNPs) have been regarded to be higher bioavailability and less toxicity in comparison with selenite and selenate. Still, little has been known about the mechanism of their metabolism in plants. Soybean (Glycine max L.) can enrich Se, providing an ideal carrier for Se biofortification. In this study, soybean sprouts were treated with SeNPs, and a combination of next-generation sequencing (NGS) and single-molecule real-time (SMRT) sequencing was applied to clarify the underlying molecular mechanism of SeNPs metabolism. A total of 74,662 nonredundant transcripts were obtained, and 2109 transcription factors, 9687 alternative splice events, and 3309 long non-coding RNAs (lncRNAs) were predicted, respectively. KEGG enrichment analysis of the DEGs revealed that metabolic pathways, biosynthesis of secondary metabolites, and peroxisome were most enriched both in roots and leaves after exposure to SeNPs. A total of 117 transcripts were identified to be putatively involved in SeNPs transport and biotransformation in soybean. The top six hub genes and their closely coexpressed Se metabolism-related genes, such as adenylylsulfate reductase (APR3), methionine-tRNA ligase (SYM), and chloroplastic Nifs-like cysteine desulfurases (CNIF1), were screened by WGCNA and identified to play crucial roles in SeNPs accumulation and tolerance in soybean. Finally, a putative metabolism pathway of SeNPs in soybean was proposed. These findings have provided a theoretical foundation for future elucidation of the mechanism of SeNPs metabolism in plants.

Identification of thyroid hormone response genes in the remodeling of dorsal muscle during Microhyla fissipes metamorphosis.

Extensive morphological, biochemical, and cellular changes occur during anuran metamorphosis, which is triggered by a single hormone, thyroid hormone (TH). The function of TH is mainly mediated through thyroid receptor (TR) by binding to the specific thyroid response elements (TREs) of direct response genes, in turn regulating the downstream genes in the cascade. The remodeling of dorsal skeletal muscle during anuran metamorphosis provides the perfect model to identify the immediate early and direct response genes that are important during apoptosis, proliferation, and differentiation of the muscle. In our current study, we performed Illumina sequencing combined with single-molecule real-time (SMRT) sequencing in the dorsal muscle of Microhyla fissipes after TH, cycloheximide (CHX), and TH_CHX treatment. We first identified 1,245 differentially expressed transcripts (DETs) after TH exposure, many of which were involved in DNA replication, protein processing in the endoplasmic reticulum, cell cycle, apoptosis, p53 signaling pathway, and protein digestion and absorption. In the comparison of the TH group vs. control group and TH_CHX group vs. CHX group overlapping gene, 39 upregulated and 6 downregulated genes were identified as the TH directly induced genes. Further analysis indicated that AGGTCAnnTnAGGTCA is the optimal target sequence of target genes for TR/RXR heterodimers in M. fissipes. Future investigations on the function and regulation of these genes and pathways should help to reveal the mechanisms governing amphibian dorsal muscle remodeling. These full-length and high-quality transcriptomes in this study also provide an important foundation for future studies in M. fissipes metamorphosis.

WHOLE-GENOME SINGLE MOLECULE REAL-TIME SEQUENCING OF SARS-CoV-2 OMICRON.

New variants and genetic mutations of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome can only be identified using accurate sequencing methods. Single molecule real-time (SMRT) sequencing has been used to characterize Alpha and Delta variants, but not Omicron variants harbouring numerous mutations in the SARS-CoV-2 genome. This study assesses the performance of a target capture SMRT sequencing protocol for whole genome sequencing (WGS) of SARS-CoV-2 Omicron variants and compared it to that of an amplicon SMRT sequencing protocol optimized for Omicron variants. The failure rate of the target capture protocol (6%) was lower than that of the amplicon protocol (34%, p<0.001) on our dataset, and the median genome coverage with the target capture protocol (98.6% [IQR: 86-99.4]) was greater than that with the amplicon protocol (76.6% [IQR: 66-89.6], (p<0.001)). The percentages of samples with >95% whole genome coverage were 64% with the target capture protocol and 19% with the amplicon protocol (p<0.05). The clades of 96 samples determined with both protocols were 93% concordant and the lineages of 59 samples were 100% concordant. Thus, target capture SMRT sequencing appears to be an efficient method for WGS, genotyping and detecting mutations of SARS-CoV-2 Omicron variants. This article is protected by copyright. All rights reserved.