De Novo Assembly Research Articles

Strain-resolved de-novo metagenomic assembly of viral genomes and microbial 16S rRNAs

BackgroundMetagenomics is a powerful approach to study environmental and human-associated microbial communities and, in particular, the role of viruses in shaping them. Viral genomes are challenging to assemble from metagenomic samples due to their genomic diversity caused by high mutation rates. In the standard de Bruijn graph assemblers, this genomic diversity leads to complex k-mer assembly graphs with a plethora of loops and bulges that are challenging to resolve into strains or haplotypes because variants more than the k-mer size apart cannot be phased. In contrast, overlap assemblers can phase variants as long as they are covered by a single read.ResultsHere, we present PenguiN, a software for strain resolved assembly of viral DNA and RNA genomes and bacterial 16S rRNA from shotgun metagenomics. Its exhaustive detection of all read overlaps in linear time combined with a Bayesian model to select strain-resolved extensions allow it to assemble longer viral contigs, less fragmented genomes, and more strains than existing assembly tools, on both real and simulated datasets. We show a 3–40-fold increase in complete viral genomes and a 6-fold increase in bacterial 16S rRNA genes.ConclusionPenguiN is the first overlap-based assembler for viral genome and 16S rRNA assembly from large and complex metagenomic datasets, which we hope will facilitate studying the key roles of viruses in microbial communities.C22dm7AjjRZTtiuc33W1snVideo

De-novo long-read genome assembly and annotation of the luna moth (Actias luna) fully resolves repeat-rich silk genes.

We present the first long-read de-novo -assembly and annotation of the luna moth (Actias luna) and provide the full characterization of heavy chain fibroin (h-fibroin)--, a long and highly repetitive gene (>20 Kbp) essential in silk fiber production. There are more than 160,000 described species of moths and butterflies (Lepidoptera), but only within the last five years have we begun to recover high-quality annotated whole genomes across the order which capture h-fibroin. Using PacBio HiFi reads, we produce the first high-quality long-read reference genome for this species. The assembled genome has a length of 532 Mbp, a contig N50 of 16.8 Mbp, an L50 of 14 contigs, and 99.4% completeness (BUSCO). Our annotation using Bombyx mori protein and A.luna RNAseq evidence captured a total of 20,866 genes at 98.9% completeness with 10,267 functionally annotated proteins and a full-length h-fibroin annotation of 2,679 amino acid residues.

Optical Genome Mapping for Applications in Repeat Expansion Disorders.

Short tandem repeat (STR) expansions are associated with more than 60 genetic disorders. The size and stability of these expansions correlate with the severity and age of onset of the disease. Therefore, being able to accurately detect the absolute length of STRs is important. Current diagnostic assays include laborious lab experiments, including repeat-primed PCR and Southern blotting, that still cannot precisely determine the exact length of very long repeat expansions. Optical genome mapping (OGM) is a cost-effective and easy-to-use alternative to traditional cytogenetic techniques and allows the comprehensive detection of chromosomal aberrations and structural variants >500 bp in length, including insertions, deletions, duplications, inversions, translocations, and copy number variants. Here, we provide methodological guidance for preparing samples and performing OGM as well as running the analysis pipelines and using the specific repeat expansion workflows to determine the exact repeat length of repeat expansions expanded beyond 500 bp. Together these protocols provide all details needed to analyze the length and stability of any repeat expansion with an expected repeat size difference from the expected wild-type allele of >500 bp. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Genomic ultra-high-molecular-weight DNA isolation, labeling, and staining Basic Protocol 2: Data generation and genome mapping using the Bionano Saphyr® System Basic Protocol 3: Manual De Novo Assembly workflow Basic Protocol 4: Local guided assembly workflow Basic Protocol 5: EnFocus Fragile X workflow Basic Protocol 6: Molecule distance script workflow.

Complete mitogenome and intra-family comparative mitogenomics showed distinct position of Pama Croaker Otolithoides pama

The Pama Croaker, Otolithoides pama, is an economically important fish species in Bangladesh. Intra-family similarities in morphology and typical barcode sequences of cox1 create ambiguities in its identification. Therefore, morphology and the complete mitochondrial genome of O. pama, and comparative mitogenomics within the family Sciaenidae have been studied. Extracted genomic DNA was subjected to Illumina-based short read sequencing for De-Novo mitogenome assembly. The complete mitogenome of O. pama (Accession: OQ784575.1) was 16,513 bp, with strong AC biasness and strand asymmetry. Relative synonymous codon usage (RSCU) among 13 protein-coding genes (PCGs) of O. pama was also analyzed. The studied mitogenomes including O. pama exhibited consistent sizes and gene orders, except for the genus Johnius which possessed notably longer mitogenomes with unique gene rearrangements. Different genetic distance metrics across 30 species of Sciaenidae family demonstrated 12S rRNA and the control region (CR) as the most conserved and variable regions, respectively, while most of the PCGs undergone a purifying selection. Different phylogenetic trees were congruent with one another, where O. pama was distinctly placed. This study would contribute to distinguishing closely related fish species of Sciaenidae family and can be instrumental in conserving the genetic diversity of O. pama.

De Novo Assembly and Annotation of the Siganus fuscescens (Houttuyn, 1782) Genome: Marking a Pioneering Advance for the Siganidae Family.

This study presents the first draft genome of Siganus fuscescens, and thereby establishes the first whole-genome sequence for a species in the Siganidae family. Leveraging both long and short read sequencing technologies, i.e., Oxford Nanopore and Illumina sequencing, we successfully assembled a mitogenome spanning 16.494 Kb and a first haploid genome encompassing 498Mb. The assembled genome accounted for a 99.6% of the estimated genome size and was organized into 164 contigs with an N50 of 7.2Mb. This genome assembly showed a GC content of 42.9% and a high Benchmarking Universal Single-Copy Orthologue (BUSCO) completeness score of 99.5% using actinopterygii_odb10 lineage, thereby meeting stringent quality standards. In addition to its structural aspects, our study also examined the functional genomics of this species, including the intricate capacity to biosynthesize long-chain polyunsaturated fatty acids (LC-PUFAs) and secrete venom. Notably, our analyses revealed various repeats elements, which collectively constituted 17.43% of the genome. Moreover, annotation of 28,351 genes uncovered both shared genetic signatures and those that are unique to S. fuscescens. Our assembled genome also displayed a moderate prevalence of gene duplication compared to other fish species, which suggests that this species has a distinctive evolutionary trajectory and potentially unique functional constraints. Taken altogether, this genomic resource establishes a robust foundation for future research on the biology, evolution, and the aquaculture potential of S. fuscescens.

Mitochondrial genome complexity in Stemona sessilifolia: nanopore sequencing reveals chloroplast gene transfer and DNA rearrangements.

Mitochondria are semi-autonomous organelles in eukaryotic cells with their own genome. Plant mitogenomes differ from animal mitogenomes in size, structure, and repetitive DNA sequences. Despite larger sizes, plant mitogenomes do not have significantly more genes. They exhibit diverse structures due to variations in size, repetitive DNA, recombination frequencies, low gene densities, and reduced nucleotide substitution rates. In this study, we analyzed the mitochondrial genome of Stemona sessilifolia using Nanopore and Illumina sequencing. De-novo assembly and annotation were conducted using Unicycler, Geseq, tRNAscan-SE and BLASTN, followed by codon usage, repeat sequence, RNA-editing, synteny, and phylogenetic analyses. S. sessilifolia's mitogenome consisted of one linear contig and six circular contigs totaling 724,751bp. It had 39 protein-coding genes, 27 tRNA genes, and 3 rRNA genes. Transfer of chloroplast sequences accounted for 13.14% of the mitogenome. Various analyses provided insights into genetic characteristics, evolutionary dynamics, and phylogenetic placement. Further investigations can explore transferred genes' functions and RNA-editing's role in mitochondrial gene expression in S. sessilifolia.

De Novo Assembly and Comparative Analysis of Mitochondrial Genomes of Two Pueraria montana Varieties.

Pueraria montana is a species with important medicinal value and a complex genetic background. In this study, we sequenced and assembled the mitochondrial (mt) genomes of two varieties of P. montana. The mt genome lengths of P. montana var. thomsonii and P. montana var. montana were 457,390 bp and 456,731 bp, respectively. Both P. montana mitogenomes showed a multi-branched structure consisting of two circular molecules, with 56 genes annotated, comprising 33 protein-coding genes, 18 tRNA genes (trnC-GCA and trnM-CAU are multi-copy genes), and 3 rRNA genes. Then, 207 pairs of long repeats and 96 simple sequence repeats (SSRs) were detected in the mt genomes of P. montana, and 484 potential RNA-editing sites were found across the 33 mitochondrial protein-coding genes of each variety. Additionally, a syntenic sequence analysis showed a high collinearity between the two mt genomes. This work is the first to analyze the mt genomes of P. montana. It can provide information that can be used to analyze the structure of mt genomes of higher plants and provide a foundation for future comparative genomic studies and evolutionary biology research in related species.

De Novo Assembly, Characterization, and Comparative Transcriptome Analysis of Mature Male and Female Gonads of Rabbitfish (Siganus oramin) (Bloch & Schneider, 1801).

The rabbitfish, Siganus oramin, is a commercially important table fish in southeastern China. However, there have been few studies on its gonad development and reproduction regulation. Comparative transcriptome analysis was first performed on adult male and female gonads of S. oramin. In total, 47,070 unigenes were successfully assembled and 22,737 unigenes were successfully annotated. Through comparative transcriptome analysis of male and female gonads, a total of 6722 differentially expressed genes were successfully identified, with 3528 upregulated genes and 3154 downregulated genes in the testes. In addition, 39 differentially expressed reproduction-related genes were identified. Finally, quantitative real-time PCR was used to validate the expression levels of several differentially expressed genes. These results provide important data for further studying the function of reproduction-related genes and the molecular mechanism regulating gonad development and reproduction in S. oramin.

Mitogenomic analysis of Rüppell’s fox (Vulpes rueppellii) confirms phylogenetic placement within the Palaearctic clade shared with its sister species, the red fox (Vulpes vulpes)

The Rüppell’s fox (Vulpes rueppellii) inhabits desert regions across North Africa, the Arabian Peninsula and southwestern Asia. Its phylogenetic relationship with other fox species, especially within the phylogeographic context of its sister species, V. vulpes, remain unclear. We here report the sequencing and de-novo assembly of the first annotated mitogenome of V. rueppellii, analysed with data from other foxes (tribe Vulpini, subfamily Caninae). We used four bioinformatic approaches to reconstruct the V. rueppellii mitogenome, obtaining identical sequences except for the incompletely assembled tandem-repeat region within the D-loop. The mitogenome displayed an identical organization, number and length of genes as V. vulpes. We found high support for clustering of both known subclades of V. rueppellii within the Palearctic clade of V. vulpes, rendering the latter species paraphyletic, consistent with previous analyses of shorter mtDNA fragments. More work is needed for a full understanding of the evolutionary drivers and consequences of hybridization in foxes.

Transcriptome profiles reveal response mechanisms and key role of PsNAC1 in Pinus sylvestris var. mongolica to drought stress

BackgroundDrought stress severely impedes plant growth, and only a limited number of species exhibit long-term resistance to such conditions. Pinus sylvestris var. mongolica, a dominant tree species in arid and semi-arid regions of China, exhibits strong drought resistance and plays a crucial role in the local ecosystem. However, the molecular mechanisms underlying this resistance remain poorly understood.ResultsHere, we conducted transcriptome sequence and physiological indicators analysis of needle samples during drought treatment and rehydration stages. De-novo assembly yielded approximately 114,152 unigenes with an N50 length of 1,363 bp. We identified 6,506 differentially expressed genes (DEGs), with the majority being concentrated in the heavy drought stage (4,529 DEGs). Functional annotation revealed enrichment of drought-related GO terms such as response to water (GO:0009415: enriched 108 genes) and response to water deprivation (GO:0009414: enriched 106 genes), as well as KEGG categories including MAPK signaling pathway (K04733: enriched 35 genes) and monoterpenoid biosynthesis (K21374: enriched 27 genes). Multiple transcription factor families and functional protein families were differentially expressed during drought treatment. Co-expression network analysis identified a potential drought regulatory network between cytochrome P450 genes (Unigene4122_c1_g1) and a core regulatory transcription factor Unigene9098_c3_g1 (PsNAC1) with highly significant expression differences. We validated PsNAC1 overexpression in Arabidopsis and demonstrated enhanced drought resistance.ConclusionsThese findings provide insight into the molecular basis of drought resistance in P. sylvestris var. mongolica and lay the foundation for further exploration of its regulatory network.

Annotated genome and transcriptome of the endangered Caribbean mountainous star coral (Orbicella faveolata) using PacBio long-read sequencing

Long-read sequencing is revolutionizing de-novo genome assemblies, with continued advancements making it more readily available for previously understudied, non-model organisms. Stony corals are one such example, with long-read de-novo genome assemblies now starting to be publicly available, opening the door for a wide array of ‘omics-based research. Here we present a new de-novo genome assembly for the endangered Caribbean star coral, Orbicella faveolata, using PacBio circular consensus reads. Our genome assembly improved the contiguity (51 versus 1,933 contigs) and complete and single copy BUSCO orthologs (93.6% versus 85.3%, database metazoa_odb10), compared to the currently available reference genome generated using short-read methodologies. Our new de-novo assembled genome also showed comparable quality metrics to other coral long-read genomes. Telomeric repeat analysis identified putative chromosomes in our scaffolded assembly, with these repeats at either one, or both ends, of scaffolded contigs. We identified 32,172 protein coding genes in our assembly through use of long-read RNA sequencing (ISO-seq) of additional O. faveolata fragments exposed to a range of abiotic and biotic treatments, and publicly available short-read RNA-seq data. With anthropogenic influences heavily affecting O. faveolata, as well as its increasing incorporation into reef restoration activities, this updated genome resource can be used for population genomics and other ‘omics analyses to aid in the conservation of this species.

First Report of Melon Chlorotic Spot Virus in Cultivated Sorrel (Rumex acetosa) in Belgium

In 2020, symptoms of putative viral origin were observed on 7% of tomatoes in an organic vegetable farm in Belgium (deformed uneven ripened fruits, vein clearing, mosaic and purple leaves, stunted plants). The leaves of twenty symptomatic plants were collected, pooled and screened for viruses using high throughput sequencing technologies (HTS) on Illumina NextSeq500 following a virion-associated nucleic acid (VANA) protocol (Temple et al., 2021, Be_SL1). In total, 3,665,498 reads (PE150) were generated. Bioinformatic analyses (denovo assembly, tblastx search on NCBI and mapping) using Geneious Prime® 2020.1.2 revealed the presence of three viruses known to infect tomatoes: Physostegia chlorotic mottle virus (PhCMoV), 547,142 reads map on NC_055466, potato virus Y (PVY), 4056 reads map on MW595184, and melon chlorotic spot virus (MeCSV), 55 reads mapped to six out of the eight different MeCSV segments (NC_040448-55). Tomato plants have already been artificially inoculated by MeCSV (Lecoq et al., 2019) but this detection (confirmed by independent RT-PCR on the pooled sample) is the first one in natural condition on farm. The high prevalence of symptoms triggered the research of alternative perennial hosts that can serve as a reservoir during inter-cropping season. One plant of Rumex acetosa showing vein clearing (CT-122) was collected in the same greenhouse the year after. Total RNA was extracted, followed by ribodepletion, and Illumina HTS using the protocol described in Temple et al., (2021) for Be_GP1. In total, 4,549,721 PE150 reads were obtained and bioinformatic analyses confirmed the presence of MeCSV (8,816 reads mapped on eight RNA segments NC_040448-55 with an average 96,52% coverage of the reference sequences, supplementary table 1) and suggested the presence of an unclassified partitivirus. Consensus sequences were extracted for each segment of MeCSV (OQ818038-45) and showed between 83% and 87% of nucleotide identity with the reference sequences NC_040448-55. RNA1 segment was used to design MeCSV-specific RT-PCR primers for detection (MeCSV-125F 5’-TTTAAGGCCAGATCCAGAGGTTC-3’/ MeCSV-498R 5’-TGGATGTGACAACCTGGTAGTAC-3’). Thereafter, in July 2022, 42 R. acetosa plants were collected in the same greenhouse. Among them, seven plants showed vein clearing, two showed yellowing with necrosis, two exhibited yellowing and vein clearing (Supplementary figure 1), and one showed mosaic. The 42 plants were subjected to RNA extraction and RT-PCR for MeCSV (Supplementary figure 2) and PhCMoV detection. MeCSV was detected in 13 plants (two asymptomatic plants and all the symptomatic plants except the one exhibiting mosaic where PhCMoV was detected). PhCMoV was also detected in three plants with vein clearing, one with yellowing and one of the two asymptomatic plants infected by MeCSV. Our results report the first detection of MeCSV in R. acetosa and the first detection of MeCSV in Belgium. In addition, according to the hierarchical approach for assessing causal relationships in plant virology (Fox et al., 2020), a preliminary association was observed between symptoms and MeCSV detection [6% prevalence on asymptomatic plants and 92% prevalence on diseased plants (from which seven symptomatic samples were not co-infected by PhCMoV)]. Symptom causality should be further investigated but this results are important for disease management because they suggested that cultivated perennial R. acetosa may serve as a reservoir for two emergent plant viruses (PhCMoV and MeCSV) (Lecoq et al., 2019, Temple et al., 2021).

De novo transcriptome analysis of Argemone mexicana plant (SRR287649) and annotation for secondary metabolite producing genes

Argemone mexicana (family: Papaveraceae) is an important medicinal plant that produces a wide range of secondary metabolites. In this study, the RNA-seq data of Argemone mexicana plant stem was retrieved from the NCBI SRA database (SRR287649). The quality checking and trimming of the raw data were performed by using FastQC, Cutadapt and Trimmomatic 0.39 software. Further, de-novo sequence assembly by Trinity software package generated a total 29591 number of contig sequences. As subsequent work, BLASTX and annotation of all the contig by using Blast2GO and WEGO2 tool produced 19780 annotated genes with 44440 gene ontology (GO) terms. Further, the genes involved in secondary metabolite production of the Argemone mexicana plant were predicted by functional assignment of the contig sequences.

Whole Genome Sequences, De Novo Assembly, and Annotation of Antibiotic Resistant Campylobacter jejuni Strains S27, S33, and S36 Newly Isolated from Chicken Meat.

Campylobacter is a leading bacterial cause of gastrointestinal infections in humans and has imposed substantial medical and public health burdens worldwide. Among a total of 39 species in the Campylobacter genus, C. jejuni is the most important species responsible for approx. 90% of human Campylobacter illness. Most cases of the infection were acquired by ingesting undercooked poultry meat due to the high prevalence of Campylobacter in the products. Here, we reported the dataset of raw sequences, de novo assembled and annotated genomes of C. jejuni strains S27, S33, and S36 recently isolated from retail chicken by using PacBio highly accurate long-read sequencing technology combined with bioinformatics tools. Our data revealed several virulence and antibiotic resistance genes in each of the chromosomes, a type IV secretion system in the plasmid (pCjS33) of C. jejuni S33, and a type VI secretion system and a phage in the plasmid (pCjS36) of C. jejuni S36. This study not only provides new sequence data but also extends the knowledge pertaining to the genomic and functional aspects of this important foodborne pathogen, including the genetic determinants of virulence and antibiotic resistance.

De Novo Assembly of Venom Gland Transcriptome of Tropidolaemus wagleri (Temple Pit Viper, Malaysia) and Insights into the Origin of Its Major Toxin, Waglerin.

The venom proteome of Temple Pit Viper (Tropidolaemus wagleri) is unique among pit vipers, characterized by a high abundance of a neurotoxic peptide, waglerin. To further explore the genetic diversity of its toxins, the present study de novo assembled the venom gland transcriptome of T. wagleri from west Malaysia. Among the 15 toxin gene families discovered, gene annotation and expression analysis reveal the dominating trend of bradykinin-potentiating peptide/angiotensin-converting enzyme inhibitor-C-type natriuretic peptide (BPP/ACEI-CNP, 76.19% of all-toxin transcription) in the transcriptome, followed by P-III snake venom metalloproteases (13.91%) and other toxins. The transcript TwBNP01 of BPP/ACEI-CNP represents a large precursor gene (209 amino acid residues) containing the coding region for waglerin (24 residues). TwBNP01 shows substantial sequence variations from the corresponding genes of its sister species, Tropidolaemus subannulatus of northern Philippines, and other viperid species which diversely code for proline-rich small peptides such as bradykinin-potentiating peptides (BPPs). The waglerin/waglerin-like peptides, BPPs and azemiopsin are proline-rich, evolving de novo from multiple highly diverged propeptide regions within the orthologous BPP/ACEI-CNP genes. Neofunctionalization of the peptides results in phylogenetic constraints consistent with a phenotypic dichotomy, where Tropidolaemus spp. and Azemiops feae convergently evolve a neurotoxic trait while vasoactive BPPs evolve only in other species.

Adaptation to the High-Arctic island environment despite long-term reduced genetic variation in Svalbard reindeer

Typically much smaller in number than their mainland counterparts, island populations are ideal systems to investigate genetic threats to small populations. The Svalbard reindeer (Rangifer tarandus platyrhynchus) is an endemic subspecies that colonized the Svalbard archipelago ca. 6,000-8,000 years ago and now shows numerous physiological and morphological adaptations to its arctic habitat. Here, we report a de-novo chromosome-level assembly for Svalbard reindeer and analyze 133 reindeer genomes spanning Svalbard and most of the species' Holarctic range, to examine the genomic consequences of long-term isolation and small population size in this insular subspecies. Empirical data, demographic reconstructions, and forward simulations show that long-term isolation and high inbreeding levels may have facilitated the reduction of highly deleterious-and to a lesser extent, moderately deleterious-variation. Our study indicates that long-term reduced genetic diversity did not preclude local adaptation to the High Arctic, suggesting that even severely bottlenecked populations can retain evolutionary potential.

De Novo Assembly, Characterization and Comparative Transcriptome Analysis of the Gonads of Jade Perch (Scortum barcoo).

Due to the high meat yield and rich nutritional content, jade perch (Scortum barcoo) has become an important commercial aquaculture species in China. Jade perch has a slow growth rate, taking 3-4 years to reach sexual maturity, and has almost no difference in body size between males and females. However, the study of its gonad development and reproduction regulation is still blank, which limited the yield increase. Herein, the gonad transcriptomes of juvenile males and females of S. barcoo were identified for the first time. A total of 107,060 unigenes were successfully annotated. By comparing male and female gonad transcriptomes, a total of 23,849 differentially expressed genes (DEGs) were identified, of which 9517 were downregulated, and 14,332 were upregulated in the testis. In addition, a large number of DEGs involved in sex differentiation, gonadal development and differentiation and gametogenesis were identified, and the differential expression patterns of some genes were further verified using real-time fluorescence quantitative PCR. The results of this study will provide a valuable resource for further studies on sex determination and gonadal development of S. barcoo.

Complete genome sequence of Treponema pedis GNW45 isolated from dairy cattle with active bovine digital dermatitis in Korea.

Treponema pedis, a fastidious anaerobic spirochete, is one of the main pathogens involved in the development and progression of bovine digital dermatitis (BDD), a lameness-causing hoof infection in cattle. Here, the complete genome sequencing of T. pedis GNW45 isolated from a dairy cow infected with BDD, was presented. Libraries for long and short reads were sequenced using PacBioRSII and Illimuna HiSeqXTen platforms, respectively. De-novo assembly was done using the long reads, producing a circular contig, by which the short reads were aligned to generate a more accurate genome sequence. The genome has a total size of 3,077,465 base pairs, with 36.84% guanine-cytosine content. A total of 2,749 protein-coding sequences, seven ribosomal RNA's, and 45 transfer RNA's were annotated. Functional analysis revealed genes associated with pathogenicity and survivability in the complex pathobiome of BDD. This study provided novel insights into the survival and pathogenic mechanisms of T. pedis GNW45.

Genome Sequencing and Organization of Three Geographically Different Isolates of Nucleopolyhedrovirus from the Gypsy Moth Reveal Significant Genomic Differences.

The gypsy moth (Lymantria dispar L., Lepidoptera: Erebidae) is a worldwide pest of trees and forests. Lymantria dispar nucleopolyhedrovirus (LdMNPV) belongs to the Baculoviridae family and is an insect virus specific to gypsy moth larvae. In this study, we describe the complete genome sequences of three geographically diverse isolates, H2 (China), J2 (Japan), and T3 (Turkey), of Lymantria dispar multiple nucleopolyhedrovirus (LdMNPV). The genomes of isolates H2, J2, and T3 were subjected to shotgun pyrosequencing using Roche 454 FLX and assembled using Roche GS De Novo Assembler. Comparative analysis of all isolates was performed using bioinformatics methods. The genomes of LdMNPV-H2, J2, and T3 were 164,746, 162,249, and 162,614 bp in size, had GC content of 57.25%, 57.30%, and 57.46%, and contained 162, 165, and 164 putative open reading frames (ORFs ≥ 150 nt), respectively. Comparison between the reference genome LdMNPV-5/6 (AF081810) and the genomes of LdMNPV-H2, J2, and T3 revealed differences in gene content. Compared with LdMNPV-5/6, ORF5, 6, 8, 10, 31, and 67 were absent in LdMNPV-H2, ORF5, 13, and 66 were absent in LdMNPV-J2, and ORF10, 13, 31, and 67 were absent in LdMNPV-T3. In addition, the gene encoding the mucin-like protein (ORF4) was split into two parts in isolates H2 and T3 and designated ORF4a and ORF4b. Phylogenetic analysis grouped isolates H2 and J2 in a different cluster than isolate T3, which is more closely related to the Turkish and Polish isolates. In addition, H2 was found to be closely related to a South Korean LdMNPV isolate. This study provided a more detailed overview of the relationships between different geographic LdMNPV isolates. The results showed remarkable differences between groups at the genome level.

Integrative analysis of metabolite and transcriptome reveals biosynthetic pathway and candidate genes for eupatilin and jaceosidin biosynthesis in Artemisia argyi.

Artemisia argyi (A. argyi) is a medicinal plant belonging to the Asteraceae family and Artemisia genus. Flavonoids abundant in A. argyi are associated with anti-inflammatory, anticancer, and antioxidative effects. Eupatilin and jaceosidin are representative polymethoxy flavonoids with medicinal properties significant enough to warrant the development of drugs using their components. However, the biosynthetic pathways and related genes of these compounds have not been fully explored in A. argyi. This study comprehensively analyzed the transcriptome data and flavonoids contents from four different tissues of A. argyi (young leaves, old leaves, trichomes collected from stems, and stems without trichomes) for the first time. We obtained 41,398 unigenes through the de-novo assembly of transcriptome data and mined promising candidate genes involved in the biosynthesis of eupatilin and jaceosidin using differentially expressed genes, hierarchical clustering, phylogenetic tree, and weighted gene co-expression analysis. Our analysis led to the identification of a total of 7,265 DEGs, among which 153 genes were annotated as flavonoid-related genes. In particular, we were able to identify eight putative flavone-6-hydroxylase (F6H) genes, which were responsible for providing a methyl group acceptor into flavone basic skeleton. Furthermore, five O-methyltransferases (OMTs) gene were identified, which were required for the site-specific O-methylation during the biosynthesis of eupatilin and jaceosidin. Although further validation would be necessary, our findings pave the way for the modification and mass-production of pharmacologically important polymethoxy flavonoids through genetic engineering and synthetic biological approaches.