De novo genome assembly and functional insights of the first commercial pink Auricularia cornea.

A pioneering pink cultivar of Auricularia cornea, first commercially cultivated in 2022, lacks genomic data, hindering research in genetic breeding, gene discovery, and product development. Here, we report the de novo assembly of the pink A. cornea Fen-A1 genome and provide a detailed functional annotation. The genome is 73.17 Mb in size, contains 86 scaffolds (N50 âˆ¼ 5.49 Mb), 59.09% GC content and encodes 19,120 predicted genes with a BUSCO completeness of 92.60%. Comparative genomic analysis reveals the phylogenetic relatedness of Fen-A1 and remarkable gene family dynamics. Putative genes were found mapped to 3 antibiotic-related, 36 light-dependent and 25 terpene metabolites. In addition, 789 CAZymes genes were classified, revealing the dynamics of quality loss due to postharvest refrigeration. Overall, our work is the first report on a pink A. cornea genome and provides a comprehensive insight into its complex functions.

Transcriptomic Analysis Under Drought and Salt Stress Provides Insight into Genes Putatively Involved in Ginsenoside Biosynthesis in Panax japonicus Meyer.

Panax japonicus Meyer, a perennial herb of the dicotyledonaceae family Araliaceae, is a rare folk traditional Chinese medicine, known as "the king of herbal medicine" in China. To understand the genes involved in secondary pathways under drought and salt stress, the transcriptomic analysis of P. japonicus is of vital importance. The transcriptome of underground rhizomes, stems, and leaves under drought and salt stress in P. japonicus were performed using the Illumina HiSeq platform. After de novo assembly of transcripts, expression profiling and identified differentially expressed genes (DEGs) were performed. Furthermore, putative functions of identified DEGs correlated with ginsenoside in P. japonicus were explored using Gene Ontology terms and Kyoto Encyclopedia of Genes and Genome (KEGG) pathway enrichment analysis. A total of 221,804 unigenes were obtained from the transcriptome of P. japonicus. The further analysis revealed that 10,839 unigenes were mapped to 91 KEGG pathways. Furthermore, a total of two metabolic pathways of P. japonicus in response to drought and salt stress related to triterpene saponin synthesis were screened. The sesquiterpene and triterpene metabolic pathways were annotated and finally putatively involved in ginsenoside content and correlation analysis of the expression of these genes were analyzed to identify four genes, ß-amyrin synthase, isoprene synthase, squalene epoxidase, and 1-deoxy-D-ketose-5-phosphate synthase, respectively. Our results paves the way for screening highly expressed genes and mining genes related to triterpenoid saponin synthesis. It also provides valuable references for the study of genes involved in ginsenoside biosynthesis and signal pathway of P. japonicus.

MHConstructor: A high-throughput, haplotype-informed solution to the MHC assembly challenge.

The extremely high levels of genetic polymorphism within the human major histocompatibility complex (MHC) limit the usefulness of reference-based alignment methods for sequence assembly. We incorporate a short read de novo assembly algorithm into a workflow for novel application to the MHC. MHConstructor is a containerized pipeline designed for high-throughput, haplotype-informed, reproducible assembly of both whole genome sequencing and target-capture short read data in large, population cohorts. To-date, no other self-contained tool exists for the generation of de novo MHC assemblies from short read data. MHConstructor facilitates wide-spread access to high quality, alignment-free MHC sequence analysis.

De novo assembly of transcriptomes and differential gene expression analysis using short-read data from emerging model organisms - a brief guide.

Many questions in biology benefit greatly from the use of a variety of model systems. High-throughput sequencing methods have been a triumph in the democratization of diverse model systems. They allow for the economical sequencing of an entire genome or transcriptome of interest, and with technical variations can even provide insight into genome organization and the expression and regulation of genes. The analysis and biological interpretation of such large datasets can present significant challenges that depend on the 'scientific status' of the model system. While high-quality genome and transcriptome references are readily available for well-established model systems, the establishment of such references for an emerging model system often requires extensive resources such as finances, expertise and computation capabilities. The de novo assembly of a transcriptome represents an excellent entry point for genetic and molecular studies in emerging model systems as it can efficiently assess gene content while also serving as a reference for differential gene expression studies. However, the process of de novo transcriptome assembly is non-trivial, and as a rule must be empirically optimized for every dataset. For the researcher working with an emerging model system, and with little to no experience with assembling and quantifying short-read data from the Illumina platform, these processes can be daunting. In this guide we outline the major challenges faced when establishing a reference transcriptome de novo and we provide advice on how to approach such an endeavor. We describe the major experimental and bioinformatic steps, provide some broad recommendations and cautions for the newcomer to de novo transcriptome assembly and differential gene expression analyses. Moreover, we provide an initial selection of tools that can assist in the journey from raw short-read data to assembled transcriptome and lists of differentially expressed genes.

Chromosome-level genome assembly provides insights into the genetic diversity, evolution, and flower development of Prunus conradinae.

Prunus conradinae, a valuable flowering cherry belonging to the Rosaceae family subgenus Cerasus and endemic to China, has high economic and ornamental value. However, a high-quality P. conradinae genome is unavailable, which hinders our understanding of its genetic relationships and phylogenesis, and ultimately, the possibility of mining of key genes for important traits. Herein, we have successfully assembled a chromosome-scale P. conradinae genome, identifying 31,134 protein-coding genes, with 98.22% of them functionally annotated. Furthermore, we determined that repetitive sequences constitute 46.23% of the genome. Structural variation detection revealed some syntenic regions, inversions, translocations, and duplications, highlighting the genetic diversity and complexity of Cerasus. Phylogenetic analysis demonstrated that P. conradinae is most closely related to P. campanulata, from which it diverged ~ 19.1 million years ago (Mya). P. avium diverged earlier than P. cerasus and P. conradinae. Similar to the other Prunus species, P. conradinae underwent a common whole-genome duplication event at ~ 138.60 Mya. Furthermore, 79 MADS-box members were identified in P. conradinae, accompanied by the expansion of the SHORT VEGETATIVE PHASE subfamily. Our findings shed light on the complex genetic relationships, and genome evolution of P. conradinae and will facilitate research on the molecular breeding and functions of key genes related to important horticultural and economic characteristics of subgenus Cerasus.

First Chromosome-Level Genome Assembly of a Ribbon Worm from the Hoplonemertea Clade, Emplectonema gracile, and Its Structural Annotation.

Genome-wide information has so far been unavailable for ribbon worms of the clade Hoplonemertea, the most species-rich class within the phylum Nemertea. While species within Pilidiophora, the sister clade of Hoplonemertea, possess a pilidium larval stage and lack stylets on their proboscis, Hoplonemertea species have a planuliform larva and are armed with stylets employed for the injection of toxins into their prey. To further compare these developmental, physiological, and behavioral differences from a genomic perspective, the availability of a reference genome for a Hoplonemertea species is crucial. Such data will be highly useful for future investigations toward a better understanding of molecular ecology, venom evolution, and regeneration not only in Nemertea but also in other marine invertebrate phyla. To this end, we herein present the annotated chromosome-level genome assembly for Emplectonema gracile (Nemertea; Hoplonemertea; Monostilifera; Emplectonematidae), an easily collected nemertean well suited for laboratory experimentation. The genome has an assembly size of 157.9 Mb. Hi-C scaffolding yielded chromosome-level scaffolds, with a scaffold N50 of 10.0 Mb and a score of 95.1% for complete BUSCO genes found as a single copy. Annotation predicted 20,684 protein-coding genes. The high-quality reference genome reaches an Earth BioGenome standard level of 7.C.Q50.

Transcriptome profiling of barnyard millet (Echinochloa frumentacea L.) during grain development to reveal the genomic insights into iron accumulation.

In the realm of food nutritional security, the development of mineral-rich grains assumes a pivotal role in combating malnutrition. Within the scope of the current investigation, we endeavoured to discern the transcripts accountable for the improved accumulation of grain-Fe within Indian barnyard millet. This pursuit entailed transcriptome sequencing of genotypes BAR-1433 (with high Fe content) and BAR-1423 (with low Fe content) during two distinct stages of spike development-spike emergence and milking stage. In the context of spike emergence, we identified a cohort of 895 up-regulated transcripts and 126 down-regulated transcripts that delineated the difference between the high and low grain-Fe genotypes. In contrast, during the milking stage, the tally of up-regulated transcripts reached 436, while down-regulated transcripts numbered 285. The transcripts that consistently ascended in both developmental stages underwent functional annotation, aligning their roles with nucleolar proteins, metal-nicotianamine transporters, ribonucleoprotein complexes, vinorine synthases, cellulose synthases, auxin response factors, embryogenesis abundant proteins, cytochrome c oxidases, and zinc finger BED domain-containing proteins. Meanwhile, a heterogeneous spectrum of transcripts exhibited differential expression and upregulation throughout the distinct stages. These transcripts encompassed various facets, such as ABC Transporter family proteins, Calcium-dependent kinase family, Ferritin, Metal ion binding, Iron-sulfur cluster binding, Cytochrome family, Zinc finger transcription factor family, Ferredoxin-NADP reductase type 1 family, Putative laccase, Multicopper oxidase family, and Terpene synthase family. To authenticate the reliability of these transcripts, six contigs representing probable functions, including metal transporters, iron sulfur coordination, metal ion binding, auxin-responsive GH3-like protein 2, and cytochrome P450 71B16, were harnessed for primer design. Subsequently, these primers were utilized in the validation process through qRT-PCR, with the outcomes aligning harmoniously with the transcriptome results. This study chronicles a constellation of genes linked to elevated iron content within barnyard millet, showcasing a proof of concept for leveraging transcriptome insights in marker-assisted selection to fortify barnyard millet with iron. This marks the inaugural comprehensive transcriptome analysis delineating transcripts associated with varying levels of grain-iron content during the panicle developmental stages within the barnyard millet paradigm.

Genome assembly of Ottelia alismoides, a multiple-carbon utilisation aquatic plant.

OBJECTIVES: Ottelia Pers. is in the Hydrocharitaceae family. Species in the genus are aquatic, and China is their centre of origin in Asia. Ottelia alismoides (L.) Pers., which is distributed worldwide, is a distinguishing element in China, while other species of this genus are endemic to China. However, O. alismoides is also considered endangered due to habitat loss and pollution in some Asian countries. Ottelia alismoides is the only submerged macrophyte that contains three carbon dioxide-concentrating mechanisms, i.e. bicarbonate (HCO3-) use, crassulacean acid metabolism and the C4 pathway. In this study, we present its first genome assembly to help illustrate the various carbon metabolism mechanisms and to enable genetic conservation in the future. DATA DESCRIPTION: Using DNA and RNA extracted from one O. alismoides leaf, this work produced ∼ 73.4 Gb HiFi reads, ∼ 126.4 Gb whole genome sequencing short reads and ∼ 21.9 Gb RNA-seq reads. The de novo genome assembly was 6,455,939,835 bp in length, with 11,923 scaffolds/contigs and an N50 of 790,733 bp. Genome assembly completeness assessment with Benchmarking Universal Single-Copy Orthologs revealed a score of 94.4%. The repetitive sequence in the assembly was 4,875,817,144 bp (75.5%). A total of 116,176 genes were predicted. The protein sequences were functionally annotated against multiple databases, facilitating comparative genomic analysis.

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.

A metagenome-wide association study of gut microbiome in patients with AMD, ASD, RA, T2D & VKH diseases.

Clinical studies have already illustrated the associations between gut microbes and diseases, yet fundamental questions remain unclear that how we can universalize this knowledge. Considering the important role of human gut microbial composition in maintaining overall health, it is important to understand the microbial diversity and altered disease conditions of the human gut. Metagenomics provides a way to analyze and understand the microbes and their role in a community manner. It provides qualitative as well as quantitative measurements, in terms of relative abundance. Various studies are already going on to find out the association between microbes and diseases; still, the mined knowledge is limited. Considering the current scenario, using the targeted metagenomics approach, we analyzed the gut microbiome of 99 samples from healthy and diseased individuals. Our metagenomic analysis mainly targeted five diseased microbiomes (i.e., Age-related macular degeneration, Autism spectrum disorder, Rheumatoid arthritis, Type 2 diabetes and Vogt-Koyanagi harada), with compare to healthy microbiome, and reported disease-associated microbiome shift in different conditions.

Phasing Diploid Genome Assembly Graphs with Single-Cell Strand Sequencing.

Haplotype information is crucial for biomedical and population genetics research. However, current strategies to produce de-novo haplotype-resolved assemblies often require either difficult-to-acquire parental data or an intermediate haplotype-collapsed assembly. Here, we present Graphasing, a workflow which synthesizes the global phase signal of Strand-seq with assembly graph topology to produce chromosome-scale de-novo haplotypes for diploid genomes. Graphasing readily integrates with any assembly workflow that both outputs an assembly graph and has a haplotype assembly mode. Graphasing performs comparably to trio-phasing in contiguity, phasing accuracy, and assembly quality, outperforms Hi-C in phasing accuracy, and generates human assemblies with over 18 chromosome-spanning haplotypes.

The first complete mitochondrial genome data of the Afghan pika Ochotona rufescens (Lagomorpha, Ochotonidae), near the type locality.

The Afghan pika Ochotona rufescens (Gray, 1842) is widely distributed across the mountains of Afghanistan, Iran, Pakistan, and southwestern Turkmenistan, most often at elevations between 2,000 and 3,000 m. Here we present, for the first time, the complete mitochondrial genomes of two specimens of the nominotypical subspecies Ochotona rufescens rufescens, de novo assembled from Illumina short reads of fragmented probe-enriched DNA. The lengths of the circular mitogenomes are 16,408 bp and 16,407 bp, respectively. Both mitogenomes contain 13 protein-coding genes (PCGs), two ribosomal RNAs (16S rRNA and 12S rRNA), 22 transfer RNA genes, and a control region. The gene NAD6 and the tRNA (Gln), tRNA (Ala), tRNA (Asn), tRNA (Cys), tRNA (Tyr), tRNA (Ser), tRNA (Glu), and tRNA (Pro) are encoded on the light strand while the rest are encoded on the heavy strand. The overall nucleotide composition was ∼30% for A, 25% for T, 15% for G, and 29% for C. The mitogenome data are available in the GenBank under the accession numbers ON859136 and ON859137.

De novo genome assembly of a high-protein soybean variety HJ117.

OBJECTIVES: Soybean is an important feed and oil crop in the world due to its high protein and oil content. China has a collection of more than 43,000 soybean germplasm resources, which provides a rich genetic diversity for soybean breeding. However, the rich genetic diversity poses great challenges to the genetic improvement of soybean. This study reports on the de novo genome assembly of HJ117, a soybean variety with high protein content of 52.99%. These data will prove to be valuable resources for further soybean quality improvement research, and will aid in the elucidation of regulatory mechanisms underlying soybean protein content. DATA DESCRIPTION: We generated a contiguous reference genome of 1041.94 Mb for HJ117 using a combination of Illumina short reads (23.38 Gb) and PacBio long reads (25.58 Gb), with high-quality sequence coverage of approximately 22.44× and 24.55×, respectively. HJ117 was developed through backcross breeding, using Jidou 12 as the recurrent parent and Chamoshidou as the donor parent. The assembly was further assisted by 114.5 Gb Hi-C data (109.9×), resulting in a contig N50 of 19.32 Mb and scaffold N50 of 51.43 Mb. Notably, Core Eukaryotic Genes Mapping Approach (CEGMA) assessment and Benchmarking Universal Single-Copy Orthologs (BUSCO) assessment results indicated that most core eukaryotic genes (97.18%) and genes in the BUSCO dataset (99.4%) were identified, and 96.44% of the genomic sequences were anchored onto twenty pseudochromosomes.

The Complete Mitochondrial Genome of Paeonia lactiflora Pall. (Saxifragales: Paeoniaceae): Evidence of Gene Transfer from Chloroplast to Mitochondrial Genome.

Paeonia lactiflora (P. lactiflora), a perennial plant renowned for its medicinal roots, provides a unique case for studying the phylogenetic relationships of species based on organelle genomes, as well as the transference of DNA across organelle genomes. In order to investigate this matter, we sequenced and characterized the mitochondrial genome (mitogenome) of P. lactiflora. Similar to the chloroplast genome (cpgenome), the mitogenome of P. lactiflora extends across 181,688 base pairs (bp). Its unique quadripartite structure results from a pair of extensive inverted repeats, each measuring 25,680 bp in length. The annotated mitogenome includes 27 protein-coding genes, 37 tRNAs, 8 rRNAs, and two pseudogenes (rpl5, rpl16). Phylogenetic analysis was performed to identify phylogenetic trees consistent with Paeonia species phylogeny in the APG Ⅳ system. Moreover, a total of 12 MTPT events were identified and 32 RNA editing sites were detected during mitogenome analysis of P. lactiflora. Our research successfully compiled and annotated the mitogenome of P. lactiflora. The study provides valuable insights regarding the taxonomic classification and molecular evolution within the Paeoniaceae family.

Assembly and characterization of the complete mitochondrial genome of Ventilago leiocarpa.

KEY MESSAGE: We reported the mitochondrial genome of Ventilago leiocarpa for the first time. Two and one sites lead to the generation of stop and stat codon through editing were verified. Ventilago leiocarpa, a member of the Rhamnaceae family, is frequently utilized in traditional medicine due to the medicinal properties of its roots. In this study, we successfully assembled the mitogenome of V. leiocarpa using both BGI short reads and Nanopore long reads. This mitogenome has a total length of 331,839 bp. The annotated results showed 36 unique protein-coding, 16 tRNA and 3 rRNA genes in this mitogenome. Furthermore, we confirmed the presence of a branched structure through the utilization of long reads mapping, PCR amplification, and Sanger sequencing. Specifically, the ctg1 can form a single circular molecule or combine with ctg4 to form a linear molecule. Likewise, ctg2 can form a single circular molecule or can be connected to ctg4 to form a linear molecule. Subsequently, through a comparative analysis of the mitogenome and cpgenome sequences, we identified ten mitochondrial plastid sequences (MTPTs), including two complete protein-coding genes and five complete tRNA genes. The existence of MTPTs was verified by long reads. Colinear analysis showed that the mitogenomes of Rosales were highly divergent in structure. Finally, we identified 545 RNA editing sites involving 36 protein-coding genes by Deepred-mt. To validate our findings, we conducted PCR amplification and Sanger sequencing, which confirmed the generation of stop codons in atp9-223 and rps10-391, as well as the generation of a start codon in nad4L-2. This project reported the complex structure and RNA editing event of the V. Leiocarpa mitogenome, which will provide valuable information for the study of mitochondrial gene expression.

A review on complete silk gene sequencing and de novo assembly of artificial silk.

Silk, especially spider and insect silk, is a highly versatile biomaterial with potential applications in biomedicine, materials science, and biomimetic engineering. The primary structure of silk proteins is the basis for the mechanical properties of silk fibers. Biotechnologies such as single-molecule sequencing have facilitated an increasing number of reports on new silk genes and assembled silk proteins. Therefore, this review aims to provide a comprehensive overview of the recent advances in representative spider and insect silk proteins, focusing on identification methods, sequence characteristics, and de novo design and assembly. The review discusses three identification methods for silk genes: polymerase chain reaction (PCR)-based sequencing, PCR-free cloning and sequencing, and whole-genome sequencing. Moreover, it reveals the main spider and insect silk proteins and their sequences. Subsequent de novo assembly of artificial silk is covered and future research directions in the field of silk proteins, including new silk genes, customizable artificial silk, and the expansion of silk production and applications are discussed. This review provides a basis for the genetic aspects of silk production and the potential applications of artificial silk in material science and biomedical engineering.

De novo genome assembly of a probiotic Lacticaseibacillus rhamnosus ISO20, isolated from raw milk in South Africa.

Lactic acid bacteria are known to exhibit probiotic properties through various mechanisms including production of antimicrobial substances and bile salts tolerance. Here, we report a draft genome sequence of Lacticaseibacillus rhamnosus ISO20, a lactic acid bacterium isolated from raw goat's milk to provide genomic insight into its strategies as probiotic strain.

Mining of disease-resistance genes in Crocus sativus based on transcriptome sequencing.

Introduction: Crocus sativus L. has an important medicinal and economic value in traditional perennial Chinese medicine. However, due to its unique growth characteristics, during cultivation it is highly susceptible to disease. The absence of effective resistance genes restricts us to breed new resistant varieties of C. sativus. Methods: In present study, comprehensive transcriptome sequencing was introduced to explore the disease resistance of the candidate gene in healthy and corm rot-infected C. sativus. Results and discussion: Totally, 43.72 Gb of clean data was obtained from the assembly to generate 65,337 unigenes. By comparing the gene expression levels, 7,575 differentially expressed genes (DEGs) were primarily screened. A majority of the DEGs were completely in charge of defense and metabolism, and 152 of them were annotated as pathogen recognition genes (PRGs) based on the PGRdb dataset. The expression of some transcription factors including NAC, MYB, and WRKY members, changed significantly based on the dataset of transcriptome sequencing. Therefore, this study provides us some valuable information for exploring candidate genes involved in the disease resistance in C. sativus.

Reference transcriptome assembly of a protogynous sex change fish, harlequin sandsmelt (Parapercis pulchella).

The harlequin sandsmelt (Parapercis pulchella) is a female-to-male sex change fish in which functional females possess ovotestes that consist of both ovarian and testicular tissues. These features indicate that this species could be an excellent model for studying the flexibility of sex differentiation in vertebrates. However, genetic resources in this species have so far been limited. Therefore, in this study, the reference transcriptome of this fish was constructed through RNA-sequencing, de novo transcriptome assembly, superTranscripts construction, and functional annotations. To obtain as many genes as possible, RNA was extracted from various tissues (brains, gills, hearts, livers, guts, and gonads) and various sexual stages (females, individuals during sex change, and males) and then subjected to sequencing and downstream analyses. As a result, 91,884 representative transcripts with 32,627 protein-coding sequences were generated. 72.2% of protein-coding sequences (23,566 sequences) were functionally annotated. Also, our analysis shows that the superTranscripts method effectively removes redundant sequences from raw-assembled data compared with other strategies. The resultant dataset is a valuable resource for future molecular developmental studies on sex change in P. pulchella.

Genome assembly provides insights into the genome evolution of Baccaurea ramiflora Lour.

Baccaurea ramiflora Lour., an evergreen tree of the Baccaurea genus of the Phyllanthaceae family, is primarily distributed in South Asia, Southeast Asia, and southern China, including southern Yunnan Province. It is a wild or semi-cultivated tree species with ornamental, edible, and medicinal value, exhibiting significant development potential. In this study, we present the whole-genome sequencing of B. ramiflora, employing a combination of PacBio SMRT and Illumina HiSeq 2500 sequencing techniques. The assembled genome size was 975.8 Mb, with a contig N50 of 509.33 kb and the longest contig measuring 7.74 Mb. The genome comprises approximately 73.47% highly repetitive sequences, of which 52.1% are long terminal repeat-retrotransposon sequences. A total of 29,172 protein-coding genes were predicted, of which 25,980 (89.06%) have been annotated, Additionally, 3452 non-coding RNAs were identified. Comparative genomic analysis revealed a close relationship between B. ramiflora and the Euphorbiaceae family, with both being sister groups that diverged approximately 59.9 million years ago. During the evolutionary process, B. ramiflora exhibited positive selection in 278 candidate genes. Synonymous substitution rate and collinearity analysis demonstrated that B. ramiflora underwent a single ancient genome-wide triploidization event, without recent genome-wide duplication events. This high-quality B. ramiflora genome provides a valuable resource for basic research and tree improvement programs focusing on the Phyllanthaceae family.