Reference Genome Research Articles

Characterization of the genetic variation and evolutionary divergence of the CLEC18 family.

The C-type lectin family 18 (CLEC18) with lipid and glycan binding capabilities is important to metabolic regulation and innate immune responses against viral infection. However, human CLEC18 comprises three paralogous genes with highly similar sequences, making it challenging to distinguish genetic variations, expression patterns, and biological functions of individual CLEC18 paralogs. Additionally, the evolutionary relationship between human CLEC18 and its counterparts in other species remains unclear. To identify the sequence variation and evolutionary divergence of human CLEC18 paralogs, we conducted a comprehensive analysis using various resources, including human and non-human primate reference genome assemblies, human pangenome assemblies, and long-read-based whole-genome and -transcriptomesequencing datasets. We uncovered paralogous sequence variants (PSVs) and polymorphic variants (PVs) of human CLEC18 proteins, and identified distinct signatures specific to each CLEC18 paralog. Furthermore, we unveiled a novel segmental duplication for human CLEC18A gene. By comparing CLEC18 across human and non-human primates, our research showed that the CLEC18 paralogy probably occurred in the common ancestor of human and closely related non-human primates, and the lipid-binding CAP/SCP/TAPS domain of CLEC18 is more diverse than its glycan-binding CTLD. Moreover, we found that certain amino acids alterations at variant positions are exclusive to human CLEC18 paralogs. Our findings offer a comprehensive profiling of the intricate variations and evolutionary characteristics of human CLEC18.

Patterns of genomic variation reveal a single evolutionary origin of the wild allotetraploid Mimulus sookensis.

Polyploidy occurs across the tree of life and is especially common in plants. Because newly formed cytotypes are often incompatible with their progenitors, polyploidy is also said to trigger "instantaneous" speciation. If a polyploid can self-fertilize or reproduce asexually, it is even possible for one individual to produce an entirely new lineage, but how often this scenario occurs is unclear. Here, we investigate the evolutionary history of the wild allotetraploid Mimulus sookensis, which was formed through hybridization between self-compatible, diploid species in the Mimulus guttatus complex. We generate a chromosome-scale reference assembly for M. sookensis and define its distinct subgenomes. Despite previous reports suggesting multiple origins of this highly selfing polyploid, we discover patterns of population genomic variation that provide unambiguous support for a single origin. One M. sookensis subgenome is clearly derived from the selfer M. nasutus, which organellar variation suggests is the maternal progenitor. The ancestor of the other subgenome is less certain, but it shares variation with both M. decorus and M. guttatus, two outcrossing diploids with geographic ranges that overlap broadly with M. sookensis. This study establishes M. sookensis as an example of instantaneous speciation, likely facilitated by the polyploid's predisposition to self-fertilize.

Graph-Based Pangenome of Actinidia chinensis Reveals Structural Variations Mediating Fruit Degreening.

Fruit ripening is associated with the degreening process (loss of chlorophyll) that occurs in most fruit species. Kiwifruit is one of the special species whose fruits may maintain green flesh by accumulating a large amount of chlorophyll even after ripening. However, little is known about the genetic variations related to the fruit degreening process. Here, a graph-based kiwifruit pangenome by analyzing 14 chromosome-scale haplotype-resolved genome assemblies from seven representative cultivars or lines in Actinidia chinensis is built. A total of 49,770 non-redundant gene families are identified, with core genes constituting 46.6%, and dispensable genes constituting 53.4%. A total of 84,591 non-redundant structural variations (SVs) are identified. The pangenome graph integrating both reference genome sequences and variant information facilitates the identification of SVs related to fruit color. The SV in the promoter of the AcBCM gene determines its high expression in the late developmental stage of fruits, which causes chlorophyll accumulation in the green-flesh fruits by post-translationally regulating AcSGR2, a key enzyme of chlorophyll catabolism. Taken together, a high-quality pangenome is constructed, unraveled numerous genetic variations, and identified a novel SV mediating fruit coloration and fruit quality, providing valuable information for further investigating genome evolution and domestication, QTL genes function, and genomics-assisted breeding.

Application and Comparison of Machine Learning and Database-Based Methods in Taxonomic Classification of High-Throughput Sequencing Data.

The advent of high-throughput sequencing technologies has not only revolutionized the field of bioinformatics but has also heightened the demand for efficient taxonomic classification. Despite technological advancements, efficiently processing and analysing the deluge of sequencing data for precise taxonomic classification remains a formidable challenge. Existing classification approaches primarily fall into two categories, database-based methods and machine learning methods, each presenting its own set of challenges and advantages. On this basis, the aim of our study was to conduct a comparative analysis between these two methods while also investigating the merits of integrating multiple database-based methods. Through an in-depth comparative study, we evaluated the performance of both methodological categories in taxonomic classification by utilizing simulated datasets. Our analysis revealed that database-based methods excel in classification accuracy when backed by a rich and comprehensive reference database. Conversely, while machine learning methods show superior performance in scenarios where reference sequences are sparse or lacking, they generally show inferior performance compared to database methods under most conditions. Moreover, our study confirms that integrating multiple database-based methods does, in fact, enhance classification accuracy. These findings shed new light on the taxonomic classification of high-throughput sequencing data and bear substantial implications for the future development of computational biology. For those interested in further exploring our methods, the source code of this study is publicly available on https://github.com/LoadStar822/Genome-Classifier-Performance-Evaluator. Additionally, a dedicated webpage showcasing our collected database, datasets, and various classification software can be found at: http://lab.malab.cn/∼tqz/project/taxonomic/.

Phylum-level studies of bacterial cutinases for unravelling enzymatic specificity toward PET degradation: an in silico approach.

The overwhelming use of PET plastic in various day-to-day activities led to the voluminous increase in PET waste and growing environmental hazards. A plethora of methods have been used that are associated with secondary pollutants. Therefore, microbial degradation of PET provides a sustainable approach due to its versatile metabolic diversity and capacity. The present work highlights the cutinase enzyme's role in PET degradation. This study focuses on the bacterial cutinases homologs screened from 43 reported phylum of bacteria. The reported bacterial cutinases for plastic degradation have been chosen as reference sequences, and 917 sequences have shown homology across the bacterial phyla. The dienelactone hydrolase (DLH) domain was identified for attaining specificity towards PET binding in 196 of 917 sequences. Various computational tools have been used for the physicochemical characterization of 196 sequences. The analysis revealed that most selected sequences are hydrophilic, extracellular, and thermally stable. Based on this analysis, 17 sequences have been further pursued for three-dimensional structure prediction and validation. The molecular docking studies of 17 selected sequences revealed efficient PET binding with the three sequences derived from the phylum Bacteroidota, the lowest binding energy of -5.9kcal/mol, Armatimonadota, and Nitrososphaerota with -5.8kcal/mol. The two enzyme sequences retrieved from the phylum Bacteroidota and Armatimonadota are metagenomically derived. Therefore, the present studies concluded that there is a high probability of finding cutinase homologs in various environmental resources that can be further explored for PET degradation.

Determination of hepatitis C virus subtype prevalent in Sindh, Pakistan: a phylogenetic analysis

Hepatitis is a major public health issue, affecting 10–17 million people worldwide, with its prevalence continuously increasing. The Hepatitis C virus (HCV) is responsible for liver related diseases, which include liver cirrhosis, hepatocellular carcinoma, and chronic hepatitis. Pakistan is experiencing a serious rise in HCV cases. This study aimed to assess the prevalence and distribution of HCV genotypes in Sindh, Pakistan. Serum samples from HCV-positive patients were collected from various local hospitals in Sindh. These samples were first screened for HCV antibodies using ELISA. Samples that tested positive for HCV RNA underwent further genotyping through sequencing using the standard Sanger method. The genotypes were identified by comparing the sequences with those available in the National Center for Biotechnology Information (NCBI) database, and a phylogenetic tree was constructed. The phylogenetic analysis showed that all isolates in this study were clustered with genotypes 3a and 3b, except for one sequence that was clustered with genotype 1a. No isolates were found to be clustered with reference genomes of genotypes 2, 4, 5, 6, and 7 suggesting that genotype 3a is endemic in this region. The analyzed sequences demonstrated a 98% similarity with reference and isolated sequences. In summary, sequencing of the HCV 5′ UTR essential for identifying the predominant genotype of HCV RNA in the Sindh region Further research on the distribution of HCV genotypes in other regions of Pakistan could aid in improving screening processes, identifying more effective treatment options, and developing suitable prevention strategies.

Genome-wide identification analysis in wild-type Solanum pinnatisectum reveals some genes defending against Phytophthora infestans

Solanum pinnatisectum exhibits strong resistance to late blight caused by Phytophthora infestans but only an incomplete genome assembly based on short Illumina reads has been published. In this study, we generated the first chromosome-level draft genome for the wild-type potato species S. pinnatisectum in China using Oxford Nanopore technology sequencing and Hi-C technology. The high-quality assembled genome size is 664 Mb with a scaffold N50 value of 49.17 Mb, of which 65.87% was occupied by repetitive sequences, and predominant long terminal repeats (42.51% of the entire genome). The genome of S. pinnatisectum was predicted to contain 34,245 genes, of which 99.34% were functionally annotated. Moreover, 303 NBS-coding disease resistance (R) genes were predicted in the S. pinnatisectum genome to investigate the potential mechanisms of resistance to late blight disease. The high-quality chromosome-level reference genome of S. pinnatisectum is expected to provide potential valuable resources for intensively and effectively investigating molecular breeding and genetic research in the future.

Phased Telomere-to-Telomere Reference Genome and Pangenome Reveal an Expansion of Resistance Genes during Apple Domestication.

The cultivated apple (Malus domestica Borkh.) is a cross-pollinated perennial fruit tree of great economic importance. Previous versions of apple reference genomes were unphased, fragmented, and lacked comprehensive insights into the highly heterozygous genome, which impeded genetic studies and breeding programs in apple. In this study, we assembled a haplotype-resolved telomere-to-telomere reference genome for the diploid apple cultivar Golden Delicious. Subsequently, we constructed a pangenome based on twelve assemblies from wild and cultivated apples to investigate different types of resistance gene analogs (RGAs). Our results revealed the dynamics of the gene gain and loss events during apple domestication. Compared with cultivated species, more gene families in wild species were significantly enriched in oxidative phosphorylation, pentose metabolic process, responses to salt, and abscisic acid biosynthesis process. Interestingly, our analyses demonstrated a higher prevalence of RGAs in cultivated apples than their wild relatives, partially attributed to segmental and tandem duplication events in certain RGAs classes. Other types of structural variations, mainly deletions and insertions, have affected the presence and absence of TIR-NB-ARC-LRR (TNL), NB-ARC-LRR (NL), and CC-NB-ARC-LRR (CNL) genes. Additionally, hybridization/introgression from wild species has also contributed to the expansion of resistance genes in domesticated apples. Our haplotype-resolved T2T genome and pangenome provide important resources for genetic studies of apples, emphasizing the need to study the evolutionary mechanisms of resistance genes in apple breeding programs.

Identification of novel single nucleotide variants in the drug resistance mechanism of Mycobacterium tuberculosis isolates by whole-genome analysis

BackgroundTuberculosis (TB) represents a major global health challenge. Drug resistance in Mycobacterium tuberculosis (MTB) poses a substantial obstacle to effective TB treatment. Identifying genomic mutations in MTB isolates holds promise for unraveling the underlying mechanisms of drug resistance in this bacterium.MethodsIn this study, we investigated the roles of single nucleotide variants (SNVs) in MTB isolates resistant to four antibiotics (moxifloxacin, ofloxacin, amikacin, and capreomycin) through whole-genome analysis. We identified the drug-resistance-associated SNVs by comparing the genomes of MTB isolates with reference genomes using the MuMmer4 tool.ResultsWe observed a strikingly high proportion (94.2%) of MTB isolates resistant to ofloxacin, underscoring the current prevalence of drug resistance in MTB. An average of 3529 SNVs were detected in a single ofloxacin-resistant isolate, indicating a mutation rate of approximately 0.08% under the selective pressure of ofloxacin exposure. We identified a set of 60 SNVs associated with extensively drug-resistant tuberculosis (XDR-TB), among which 42 SNVs were non-synonymous mutations located in the coding regions of nine key genes (ctpI, desA3, mce1R, moeB1, ndhA, PE_PGRS4, PPE18, rpsA, secF). Protein structure modeling revealed that SNVs of three genes (PE_PGRS4, desA3, secF) are close to the critical catalytic active sites in the three-dimensional structure of the coding proteins.ConclusionThis comprehensive study elucidates novel resistance mechanisms in MTB against antibiotics, paving the way for future design and development of anti-tuberculosis drugs.

Comprehensive analysis of silk proteins and gland compartments in Limnephilus lunatus, a case-making trichopteran

Caddisfly larvae produce silk containing heavy and light fibroins, similar to the silk of Lepidoptera, for the construction of underwater structures. We analyzed the silk of Limnephilus lunatus belonging to the case-forming suborder Integripalpia. We analyzed the transcriptome, mapped the transcripts to a reference genome and identified over 80 proteins using proteomic methods, and checked the specificity of their expression. For comparison, we also analyzed the transcriptome and silk proteome of Limnephilus flavicornis. Our results show that fibroins and adhesives are produced together in the middle and posterior parts of the silk glands, while the anterior part produces enzymes and an unknown protein AT24. The number of silk proteins of L. lunatus far exceeds that of the web-spinning Plectrocnemia conspersa, a previously described species from the suborder Annulipalpia. Our results support the idea of increasing the structural complexity of silk in rigid case builders compared to trap web builders.

Epidemic characteristics of local HIV-2 transmission across Hunan province, China

ObjectiveTo elucidate the epidemiological features of HIV-2 in Hunan Province, China, utilizing sequence analysis. MethodsThirteen individuals diagnosed with HIV-2 infection in Hunan Province, China, from 2017 to 2023 were included in this study. Amplification of HIV-2 env and pol regions was conducted, followed by Sanger sequencing. Phylogenetic and molecular transmission network analyses were performed to delineate molecular features and transmission dynamics. ResultsAll 14 individuals contracted HIV-2 through heterosexual intercourse, comprising 7 males and 7 females, with a median age of 58 years. Among them, three couples (HN001 and HN013, HN010 and HN011, HN008 and HN009) were identified, along with commercial sexual activity engagement reported for subject HN004. Notably, subjects HN001, HN003, HN008, and HN010 engaged in commercial sexual activities at the same location as subject HN004. Phylogenetic analysis of the pol gene revealed close proximity of sequences from all subjects to reference sequences from Gambia (Sub-type A). Employing a genetic distance threshold of 1.5 %, eight out of the 14 subjects formed a molecular transmission network, with HN002 and HN004 identified as central nodes. ConclusionFrom 2017 to 2023, all HIV-2-infected individuals in Hunan Province, China, acquired the virus through identifiable routes, indicating transmission of similar HIV-2 strains among them.

The chromosome-scale reference genome for the pinfish (Lagodon rhomboides) provides insights into their evolutionary and demographic history.

The pinfish (Lagodon rhomboides) is an ecologically, economically, and culturally relevant member of family Sparidae, playing crucial roles in marine food webs of the western Atlantic Ocean and Gulf of Mexico. Despite their high abundance and ecological importance, there is a scarcity of genomic resources for this species. We assembled and annotated a chromosome-scale genome for the pinfish, resulting in a highly contiguous 785 Mb assembly of 24 scaffolded chromosomes. The high-quality assembly contains 98.9% complete BUSCOs and shows strong synteny to other chromosome-scale genomes of fish in the family Sparidae, with a limited number of large-scale genomic rearrangements. Leveraging this new genomic resource, we found evidence of significant expansions of dietary gene families over the evolutionary history of the pinfish, which may be associated with an ontogenetic shift from carnivory to herbivory seen in this species. Estimates of historical patterns of population demography using this new reference genome identified several periods of population growth and contraction which were associated with ancient climatic shifts and sea level changes. This genome serves as a valuable reference for future studies of population genomics and differentiation and provides a much-needed genomic resource for this western Atlantic sparid.

Identification and validation of stable reference genes for RT-qPCR analyses of Kobresia littledalei seedlings

BackgroundKobreisa littledalei, belonging to the Cyperaceae family is the first Kobresia species with a reference genome and the most dominant species in Qinghai-Tibet Plateau alpine meadows. It has several resistance genes which could be used to breed improved crop varieties. Reverse Transcription Quantitative Real-Time Polymerase Chain Reaction (RT-qPCR) is a popular and accurate gene expression analysis method. Its reliability depends on the expression levels of reference genes, which vary by species, tissues and environments. However, K.littledalei lacks a stable and normalized reference gene for RT-qPCR analysis.ResultsThe stability of 13 potential reference genes was tested and the stable reference genes were selected for RT-qPCR normalization for the expression analysis in the different tissues of K. littledalei under two abiotic stresses (salt and drought) and two hormonal treatments (abscisic acid (ABA) and gibberellin (GA)). Five algorithms were used to assess the stability of putative reference genes. The results showed a variation amongst the methods, and the same reference genes showed tissue expression differences under the same conditions. The stability of combining two reference genes was better than a single one. The expression levels of ACTIN were stable in leaves and stems under normal conditions, in leaves under drought stress and in roots under ABA treatment. The expression of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) expression was stable in the roots under the control conditions and salt stress and in stems exposed to drought stress. Expression levels of superoxide dismutase (SOD) were stable in stems of ABA-treated plants and in the roots under drought stress. Moreover, RPL6 expression was stable in the leaves and stems under salt stress and in the stems of the GA-treated plants. EF1-alpha expression was stable in leaves under ABA and GA treatments. The expression levels of 28 S were stable in the roots under GA treatment. In general, ACTIN and GAPDH could be employed as housekeeping genes for K. littledalei under different treatments.ConclusionThis study identified the best RT-qPCR reference genes for different K. littledalei tissues under five experimental conditions. ACTIN and GAPDH genes can be employed as the ideal housekeeping genes for expression analysis under different conditions. This is the first study to investigate the stable reference genes for normalized gene expression analysis of K. littledalei under different conditions. The results could aid molecular biology and gene function research on Kobresia and other related species.

A set of nematode rRNA cistron databases and a primer assessment tool to enable more flexible and comprehensive metabarcoding.

The ITS-2-rRNA has been particularly useful for nematode metabarcoding but does not resolve all phylogenetic relationships, and reference sequences are not available for many nematode species. This is a particular issue when metabarcoding complex communities such as wildlife parasites or terrestrial and aquatic free-living nematode communities. We have used markerDB to produce four databases of distinct regions of the rRNA cistron: the 18S rRNA gene, the 28S rRNA gene, the ITS-1 intergenic spacer and the region spanning ITS-1_5.8S_ITS-2. These databases comprise 2645, 254, 13,461 and 10,107 unique full-length sequences representing 1391, 204, 1837 and 1322 nematode species, respectively. The comparative analysis illustrates the complementary value but also reveals a better representation of Clade III, IV and V than Clade I and Clade II nematodes in each case. Although the ITS-1 database includes the largest number of unique full-length sequences, the 18S rRNA database provides the widest taxonomic coverage. We also developed PrimerTC, a tool to assess primer sequence conservation across any reference sequence database, and have applied it to evaluate a large number of previously published rRNA cistron primers. We identified sets of primers that currently provide the broadest taxonomic coverage for each rRNA marker across the nematode phylum. These new resources will facilitate more comprehensive metabarcoding of nematode communities using either short-read or long-read sequencing platforms. Further, PrimerTC is available as a simple WebApp to guide or assess PCR primer design for any genetic marker and/or taxonomic group beyond the nematode phylum.

A chromosome-level genome assembly of an avivorous bat species (Nyctalus aviator)

Currently, three carnivorous bat species, namely Ia io, Nyctalus lasiopterus, and Nyctalus aviator, are known to actively prey on seasonal migratory birds (hereinafter referred to as “avivorous bats”). However, the absence of reference genomes impedes a thorough comprehension of the molecular adaptations of avivorous bat species. Herein, we present the high-quality chromosome-scale reference genome of N. aviator based on PacBio subreads, DNBSEQ short-reads and Hi-C sequencing data. The genome assembly size of N. aviator is 1.77 Gb, with a scaffold N50 of 102 Mb, of which 99.8% assembly was anchored into 21 pseudo-chromosomes. After masking 635.1 Mb repetitive sequences, a total of 19,412 protein-coding genes were identified, of which 99.3% were functionally annotated. The genome assembly and gene prediction reached 96.1% and 96.1% completeness of Benchmarking Universal Single-Copy Orthologs (BUSCO), respectively. This chromosome-level reference genome of N. aviator fills a gap in the existing information on the genomes of carnivorous bats, especially avivorous ones, and will be valuable for mechanism of adaptations to dietary niche expansion in bat species.

Machine learning and metagenomics identifies uncharacterized taxa inferred to drive biogeochemical cycles in a subtropical hypereutrophic estuary

Abstract Anthropogenic influences have drastically increased nutrient concentrations in many estuaries globally, and microbial communities have adapted to the resulting hypereutrophic ecosystems. However, our knowledge of the dominant microbial taxa and their potential functions in these ecosystems has remained sparse. Here, we study prokaryotic community dynamics in a temporal–spatial dataset, from a subtropical hypereutrophic estuary. Screening 54 water samples across brackish to marine sites revealed that nutrient concentrations and salinity best explained spatial community variations, whereas temperature and dissolved oxygen likely drive seasonal shifts. By combining short and long read sequencing data, we recovered 2459 metagenome-assembled genomes (MAGs), proposed new taxon names for previously uncharacterised lineages, and created an extensive, habitat specific genome reference database. Community profiling based on this genome reference database revealed a diverse prokaryotic community comprising 61 bacterial and 18 archaeal phyla, and resulted in an improved taxonomic resolution at lower ranks down to genus level. We found that the vast majority (61 out of 73) of abundant genus level taxa (>1% average) represented unnamed and novel lineages, and that all genera could be clearly separated into brackish and marine ecotypes with inferred habitat specific functions. Applying supervised machine learning and metabolic reconstruction, we identified several microbial indicator taxa responding directly or indirectly to elevated nitrate and total phosphorus concentrations. In conclusion, our analysis highlights the importance of improved taxonomic resolution, sheds light on the role of previously uncharacterised lineages in estuarine nutrient cycling, and identifies microbial indicators for nutrient levels crucial in estuary health assessments.

Detection of mpox virus in wastewater provides forewarning of clinical cases in Canadian cities

Wastewater-based surveillance (WBS) has shown to be an effective tool in monitoring the spread of SARS-CoV-2 and has helped guide public health actions. Consequently, WBS has expanded to now include the monitoring of mpox virus (MPXV) to contribute to its mitigation efforts. In this study, we demonstrate a unique sample processing and a molecular diagnostic strategy for MPXV detection that can inform on the epidemiological situation of mpox outbreaks through WBS.We conducted WBS for MPXV in 22 Canadian wastewater treatment plants (WWTPs) for 14 weeks. Three MPXV qPCR assays were assessed in this study for the detection of MPXV which include the G2R assays (G2R_WA and G2R_G) developed by the Centers for Disease Control and Prevention (CDC) in 2010, and an in-house-developed assay that we have termed G2R_NML. The G2R_NML assay was designed using reference genomes from the 2022 MPXV outbreak and provides a larger qPCR amplicon size to facilitate Sanger sequencing.Results show that all three assays have similar limits of detection and are able to detect the presence of MPXV in wastewater. The G2R_NML assay produced a significantly greater number of Sanger sequence-confirmed MPXV results compared to the CDC G2R assays. Detection of MPXV was possible where provincial surveillance indicated overall low caseloads, and in some sites forewarning of up to several weeks was observed.Overall, this study proposes that WBS of MPXV provides additional information to help fill knowledge gaps in clinical case-surveillance and is potentially an essential component to the management of mpox.

Parsnp 2.0: Scalable Core-Genome alignment for massive microbial datasets.

Since 2016, the number of microbial species with available reference genomes in NCBI has more than tripled. Multiple genome alignment, the process of identifying nucleotides across multiple genomes which share a common ancestor, is used as the input to numerous downstream comparative analysis methods. Parsnp is one of the few multiple genome alignment methods able to scale to the current era of genomic data; however, there has been no major release since its initial release in 2014. To address this gap, we developed Parsnp v2, which significantly improves on its original release. Parsnp v2 provides users with more control over executions of the program, allowing Parsnp to be better tailored for different use-cases. We introduce a partitioning option to Parsnp, which allows the input to be broken up into multiple parallel alignment processes which are then combined into a final alignment. The partitioning option can reduce memory usage by over 4x and reduce runtime by over 2x, all while maintaining a precise core-genome alignment. The partitioning workflow is also less susceptible to complications caused by assembly artifacts and minor variation, as alignment anchors only need to be conserved within their partition and not across the entire input set. We highlight the performance on datasets involving thousands of bacterial and viral genomes. Parsnp v2 is available at https://github.com/marbl/parsnp.

Transcriptome dynamics of Gossypium purpurascens in response to abiotic stresses by Iso-seq and RNA-seq data

Gossypium purpurascens is a member of the Malvaceae family, holds immense economic significance as a fiber crop worldwide. Abiotic stresses harm cotton crops, reduce yields, and cause economic losses. Generating high-quality reference genomes and large-scale transcriptomic datasets across diverse conditions can offer valuable insights into identifying preferred agronomic traits for crop breeding. The present research used leaf tissues to conduct PacBio Iso-seq and RNA-seq analysis. We carried out an in-depth analysis of DEGs using both correlations with cluster analysis and principal component analysis. Additionally, the study also involved the identification of both lncRNAs and CDS. We have prepared RNA-seq libraries from 75 RNA samples to study the effects of drought, salinity, alkali, and saline-alkali stress, as well as control conditions. A total of 454.06 Gigabytes of transcriptome data were effectively validated through the identification of differentially expressed genes and KEGG and GO analysis. Overwhelmingly, gene expression profiles and full-length transcripts from cotton tissues will aid in understanding the genetic mechanism of abiotic stress tolerance in G. purpurascens.

Chromosomal-level genome assembly and annotation of the tropical sea cucumber Holothuria scabra

Holothuria scabra, a commercially valuable yet ecologically vulnerable tropical holothuroid, has experienced a severe decline in its wild populations, especially in China. Genomic resources are crucial for the development of effective genomic breeding projects and stock conservation strategies to restore these natural populations. Until now, a high-quality, chromosome-level reference genome for H. scabra has not been available. Here, we employed Oxford Nanopore and Hi-C sequencing technologies to assemble and annotate a high-quality, chromosome-level reference genome of H. scabra. The final genome comprised 31 scaffolds with a total length of 1.19 Gb and a scaffold N50 length of 53.52 Mb. Remarkably, 1,191.67 Mb (99.95%) of the sequences were anchored to 23 pseudo-chromosomes, with the longest one spanning 79.75 Mb. A total of 34,418 protein-coding genes were annotated in the final genome, with BUSCO analysis revealing 98.01% coverage of metazoa_odb10 genes, marking a significant improvement compared to the previous report. These chromosome-level sequences and annotations will provide an essential genomic basis for further investigation into molecular breeding and conservation management of H. scabra.