BiGM-lncLoc: Bi-level Multi-Graph Meta-Learning for Predicting Cell-Specific Long Noncoding RNAs Subcellular Localization

The precise spatiotemporal expression of long noncoding RNAs (lncRNAs) plays a pivotal role in biological regulation, and aberrant expression of lncRNAs in different subcellular localizations has been intricately linked to the onset and progression of a variety of cancers. Computational methods provide effective means for predicting lncRNA subcellular localization, but current studies either ignore cell line and tissue specificity or the correlation and shared information among cell lines. In this study, we propose a novel approach, BiGM-lncLoc, treating the prediction of lncRNA subcellular localization across cell lines as a multi-graph meta-learning task. Our investigation involves two categories of data: the localization data of nucleotide sequences in different cell lines and cell line expression data. BiGM-lncLoc comprises a cell line-specific optimization network learning specific knowledge from cell line expression data and a graph neural network optimized across cell lines. Subsequently, the specific and shared knowledge acquired through bi-level optimization is applied to a new cell-line prediction task without the need for re-training or fine-tuning. Additionally, through key feature analysis of the impact of different nucleotide combinations on the model, we confirm the necessity of cell line-specific studies based on correlation analysis. Finally, experiments conducted on various cell lines with different data sizes indicate that BiGM-lncLoc outperforms other methods in terms of prediction accuracy, with an average accuracy of 97.7%. After removing overlapping samples to ensure data independence for each cell line, the accuracy ranged from 82.4% to 94.7%, still surpassing existing models. Our code can be found at https://github.com/BioCL1/BiGM-lncLoc .

Mapping-based genome size estimation

While the size of chromosomes can be measured under a microscope, obtaining the exact size of a genome remains a challenge. Biochemical methods and k-mer distribution-based approaches allow only estimations. An alternative approach to estimate the genome size based on high contiguity assemblies and read mappings is presented here. Analyses of Arabidopsis thaliana and Beta vulgaris data sets are presented to show the impact of different parameters. Oryza sativa, Brachypodium distachyon, Solanum lycopersicum, Vitis vinifera, and Zea mays were also analyzed to demonstrate the broad applicability of this approach. Further, MGSE was also used to analyze Escherichia coli, Saccharomyces cerevisiae, and Caenorhabditis elegans datasets to show its utility beyond plants. Mapping-based Genome Size Estimation (MGSE) and additional scripts are available on GitHub: https://github.com/bpucker/MGSE . MGSE predicts genome sizes based on short reads or long reads requiring a minimal coverage of 5-fold.

Integrated genomic and transcriptomic analyses reveal the genetic and molecular mechanisms underlying hawthorn peel color and seed hardness diversity

Hawthorn (Crataegus pinnatifida) fruit peel color and seed hardness are key traits that significantly impact economic value. We present here the high-quality chromosome-scale genomes of two cultivars, including the hard-seed, yellow-peel C. pinnatifida "Jinruyi" (JRY) and the soft-seed, red-peel C. pinnatifida "Ruanzi" (RZ). The assembled genomes comprising 17 chromosomes are 809.1 Mb and 760.5 Mb in size, achieving scaffold N50 values of 48.5 Mb and 46.8 Mb for JRY and RZ, respectively. Comparative genomic analysis identifies 3.6-3.8 million single nucleotide polymorphisms, 8.5-9.3 million insertions/deletions, and approximately 30 Mb of presence/absence variations across different hawthorn genomes. Through integrating differentially expressed genes and accumulated metabolites, we filter candidate genes CpMYB114 and CpMYB44 associated with differences in hawthorn fruit peel color and seed hardness, respectively. Functional validation confirms that the CpMYB114-CpANS regulates anthocyanin biosynthesis in hawthorn peels, contributing to the observed variation in peel color. CpMYB44-CpCOMT is significantly upregulated in JRY and is verified to promote lignin biosynthesis, resulting in the distinction in seed hardness. Overall, this study reveals the new insights into understanding of distinct peel pigmentation and seed hardness in hawthorn and provides an abundant resource for molecular breeding.

Draft genome sequence of benzo[a]pyrene degrading Bacillus altitudinis strain AR19 isolated from Digboi oil refinery (India)

Bacillus altitudinis AR19, isolated from the Digboi oil refinery (India), has a genome size of 3,630,000 bp with a G+C content of 42.45%. The genome encodes 3,755 protein-coding genes, including those for ring-cleaving dioxygenases and biofilm formation. These genes likely play crucial roles in the bacterium's survival in hydrocarbon-enriched environments.

Genomic garden: From societal and scientific impacts to biodiversity conservation

Because of urbanization, deforestation, pollution, climate change, and natural disasters, the loss of biodiversity is a pressing concern globally. As part of our efforts toward biodiversity conservation, we propose the establishment of a genomic garden, where the genome of each plant in the garden is elucidated. Combining science, horticulture, and a digital content hub accessible with any handheld device, the genomic garden serves multiple purposes, from enhancing urban landscapes, facilitating biomedical research, and improving population health to providing entertainment and education for visitors.

Chromosome-scale genome assembly of Zoysia japonica uncovers cold tolerance candidate genes

Zoysiagrass stands out as a crucial native turfgrass due to its exceptional abiotic stress tolerance, extensive adaptability, and high ornamental value. In this study, we generated a high-quality chromosome-level genome assembly of Compadre (COM) zoysiagrass, leveraging PacBio SMRT sequencing and Hi-C scaffolding technologies. The resulting genome assembly (312.42 Mb) is anchored on 20 chromosomes, with a Scaffold N50 of 18.72 Mb. In total, 49,074 genes and 306,768 repeat sequences were annotated in the assembled genome. The first chromosome-scale genome of Zoysia japonica 'Compadre' provides a critical genetic resource for cold-tolerant turfgrass breeding through identifying stress-responsive candidate genes. Additionally, we have successfully established a cell nucleus extraction and library construction protocol tailored for zoysiagrass ATAC-seq technology, and a total of 80 low temperature tolerance candidate genes were preliminarily identified via ATAC-seq and RNA-seq profiling, thereby initiating the exploration of turfgrass epigenomics.

A taxonomic revision of the genus Raorchestes in Myanmar and Thailand with the description of two new species from Myanmar (Amphibia, Anura, Rhacophoridae)

We revise the frogs of the genus Raorchestes from Myanmar and Thailand based on data of external morphology, bioacoustics, and molecular genetics. The results of this integrative study provide evidence for the recognition of seven species, two of which we describe as new: Raorchestes mindat sp. nov. from Mindat District, Chin State, western Myanmar, and Raorchestes leiktho sp. nov. from Hpa-an District, Kayin State, southeastern Myanmar. The other species that we recognize in Myanmar and Thailand are R. cangyuanensis, R. huanglianshan, R. longchuanensis, R. menglaensis, and R. parvulus. We have compared the external morphology of the lectotype and four paralectotypes of Ixalus parvulus Boulenger, 1893 with the species of the Raorchestes parvulus group currently recognized from South-east Asia. Although the type series of Ixalus parvulus is morphologically most similar to specimens of R. cangyuanensis from Thailand, we refrain from formally synonymizing these two taxa until genetic data for I. parvulus are available that would allow this hypothesis to be tested. Thus, R. parvulus remains an enigmatic taxon still only known from the original type series. As now defined, R. cangyuanensis is distributed across most of Myanmar except for the Malayan Peninsula, and also in adjacent Yunnan Province, China, and adjacent northeastern Bangladesh. Raorchestes longchuanensis occurs in northwestern Thailand as well as in eastern Myanmar and western Yunnan, China. Raorchestes menglaensis ranges from southern Yunnan, China, across Thailand, Laos, and Cambodia to northern Western Malaysia. Raorchestes huanglianshan is distributed in southern Yunnan, China, and northwestern Thailand. Often two, at some places even three species of this genus occur sympatrically (e.g., R. leiktho sp. nov., R. longchuanensis and R. parvulus near Leiktho, Kayin State, Myanmar; R. longchuanensis and R. huanglianshan at Doi Inthanon). We provide new bioacoustic data for R. longchuanensis, R. menglaensis, and R. leiktho sp. nov., and compare these with data of R. cangyuanensis and R. rezakhani.

Chromosome-level de novo genome assembly of wild, anoxia-tolerant crucian carp, Carassius carassius

Crucian carp (Carassius carassius), a member of the carp family (Cyprinidae), is known for its remarkable anoxia tolerance. The physiological responses and adaptations to anoxia are well documented, but there is a need for better understanding of the molecular regulation and evolutionary mechanisms behind these adaptations. Here we present a high-quality, functionally annotated, chromosome-level genome assembly that can facilitate such further studies. Genomic DNA was obtained from a wild-caught crucian carp specimen and used for PacBio long-read, Illumina short-read and Hi-C sequencing. Short-read mRNA data were used for structural annotation using the BRAKER3 pipeline, while PacBio long-read RNA sequencing data were used for annotation of untranslated regions and refinement of gene-isoform relationships, using the PASA pipeline. The full assembly had a contig-level N50 of 15Mbp in 290 scaffolds and 98.6% of the total length (1.65Gbp) placed in 50 chromosomes. Structural annotation resulted in 82,557 protein-coding transcripts (in 45,667 genes), with a BUSCO completeness of 99.6% and of which 77,370 matched a protein in the UniProtKB/Swiss-Prot database.

K-mer-based Approaches to Bridging Pangenomics and Population Genetics

Many commonly studied species now have more than one chromosome-scale genome assembly, revealing a large amount of genetic diversity previously missed by approaches that map short reads to a single reference. However, many species still lack multiple reference genomes and correctly aligning references to build pangenomes can be challenging for many species, limiting our ability to study this missing genomic variation in population genetics. Here, we argue that k-mers are a very useful but underutilized tool for bridging the reference-focused paradigms of population genetics with the reference-free paradigms of pangenomics. We review current literature on the uses of k-mers for performing three core components of most population genetics analyses: identifying, measuring, and explaining patterns of genetic variation. We also demonstrate how different k-mer-based measures of genetic variation behave in population genetic simulations according to the choice of k, depth of sequencing coverage, and degree of data compression. Overall, we find that k-mer-based measures of genetic diversity scale consistently with pairwise nucleotide diversity (π) up to values of about π=0.025 (R2=0.97) for neutrally evolving populations. For populations with even more variation, using shorter k-mers will maintain the scalability up to at least π=0.1. Furthermore, in our simulated populations, k-mer dissimilarity values can be reliably approximated from counting bloom filters, highlighting a potential avenue to decreasing the memory burden of k-mer-based genomic dissimilarity analyses. For future studies, there is a great opportunity to further develop methods to identifying selected loci using k-mers.

Evidence for evolution of a new sex chromosome within the haploid-dominant Marchantiales plant lineage

Sex chromosomes have evolved independently in numerous lineages across the Tree of Life, in both diploid-dominant species, including many animals and plants, and the less studied haploid-dominant plants and algae. Strict genetic sex determination ensures that individuals reproduce by outcrossing. However, species with separate sexes (termed dioecy in diploid plants, and dioicy in haploid plants) may sometimes evolve different sex systems, and become monoicous, with the ability to self-fertilize. Here, we studied dioicy-monoicy transitions in the ancient liverwort haploid-dominant plant lineage, using three telomere-to-telomere gapless chromosome-scale reference genome assemblies from the Ricciaceae group of Marchantiales. Ancestral liverworts are believed to have been dioicous, with U and V chromosomes (chromosome 9) determining femaleness and maleness, respectively. We confirm the finding that monoicy in Ricciocarpos natans evolved from a dioicous ancestor, and most ancestrally U chromosomal genes have been retained on autosomes in this species. We also describe evidence suggesting the possible re-evolution of dioicy in the genus Riccia, with probable de novo establishment of a sex chromosome from an autosome (chromosome 5), and further translocations of genes from the new sex chromosome to autosomes. Our results also indicated that micro-chromosomes are consistent genomic features, and may have evolved independently from sex chromosomes in Ricciocarpos and Riccia lineages.

IMA GENOME - F20 A draft genome assembly of Agroatheliarolfsii, Ceratobasidiumpapillatum, Pyrenopezizabrassicae, Neopestalotiopsismacadamiae, Sphaerellopsisfilum and genomic resources for Colletotrichumspaethianum and Colletotrichumfructicola

This is a genome announacment there is no abstract

A comparative study of antibiotic resistance patterns in Mycobacterium tuberculosis

This study leverages the Bacterial and Viral Bioinformatics Resource Center (BV-BRC) to analyze over 27,000 Mycobacterium tuberculosis (MTB) genomic strains, providing a comprehensive and large-scale overview of antibiotic resistance (AMR) prevalence and resistance patterns. We used MTB++, which is the newest and most comprehensive AI-based MTB drug resistance profiler tool, to predict the resistance profile of each of the 27,000 MTB isolates and then used feature analysis to identify key genes that were associated with the resistance. There are three main contributions to this study. Firstly, it provides a detailed picture of the prevalence of specific AMR genes in the BV-BRC dataset as well as their biological implications, providing critical insight into MTB's resistance mechanisms that can help identify genes of high priority for further investigation. The second aspect of this study is to compare the prevalence of antibiotic resistance across previous studies that have addressed both the temporal and geographical evolution of MTB drug resistance. Lastly, this study emphasizes the need for targeted diagnostics and personalized treatment plans. In addition to these contributions, the study acknowledges the limitations of computational prediction and recommends future experimental validation.

X-Mapper: fast and accurate sequence alignment via gapped x-mers

Sequence alignment is foundational to many bioinformatic analyses. Many aligners start by splitting sequences into contiguous, fixed-length seeds, called k-mers. Alignment is faster with longer, unique seeds, but more accurate with shorter seeds avoiding mutations. Here, we introduce X-Mapper, aiming to offer high speed and accuracy via dynamic-length seeds containing gaps, called gapped x-mers. We observe 11-24-fold fewer suboptimal alignments analyzing a human reference and 3-579-fold lower inconsistency across bacterial references than other aligners, improving on 53% and 30% of reads aligned to non-target strains and species, respectively. Other seed-based analysis algorithms might benefit from gapped x-mers too.

Adaptive evolution of stress response genes in parasites aligns with host niche diversity

BACKGROUND

Stress responses are key the survival of parasites and, consequently, also the evolutionary success of these organisms. Despite this importance, our understanding of the evolution of molecular pathways dealing with environmental stressors in parasitic animals remains limited. Here, we tested the link between adaptive evolution of parasite stress response genes and their ecological diversity and species richness. We comparatively investigated antioxidant, heat shock, osmoregulatory, and behaviour-related genes (foraging) in two model parasitic flatworm lineages with contrasting ecological diversity, Cichlidogyrus and Kapentagyrus (Platyhelminthes: Monopisthocotyla), through whole-genome sequencing of 11 species followed by in silico exon bait capture as well as phylogenetic and codon analyses.

RESULTS

We assembled the sequences of 48 stress-related genes and report the first foraging (For) gene orthologs in flatworms. We found duplications of heat shock (Hsp) and oxidative stress genes in Cichlidogyrus compared to Kapentagyrus. We also observed positive selection patterns in genes related to mitochondrial protein import (Hsp) and behaviour (For) in species of Cichlidogyrus infecting East African cichlids-a host lineage under adaptive radiation. These patterns are consistent with a potential adaptation linked to a co-radiation of these parasites and their hosts. Additionally, the absence of cytochrome P450 and kappa and sigma-class glutathione S-transferases in monogenean flatworms is reported, genes considered essential for metazoan life.

CONCLUSIONS

This study potentially identifies the first molecular function linked to a flatworm radiation. Furthermore, the observed gene duplications and positive selection indicate the potentially important role of stress responses for the ecological adaptation of parasite species.

Population sequencing of cherry accessions unravels the evolution of Cerasus species and the selection of genetic characteristics in edible cherries

Cerasus is a subgenus of Prunus in the family Rosaceae that is popular owing to its ornamental, edible, and medicinal properties. Understanding the evolution of the Cerasus subgenus and identifying selective trait loci in edible cherries are crucial for the improvement of cherry cultivars to meet producer and consumer demands. In this study, we performed a de novo assembly of a chromosome-scale genome for the sweet cherry (Prunus avium L.) cultivar 'Burlat', covering 297.55 Mb and consisting of eight chromosomes with 33,756 protein-coding genes. The resequencing and population structural analysis of 384 Cerasus representative accessions revealed that they could be divided into four groups (Group 1, Group 2, Group 3, and Group 4). We inferred that Group 1 was the oldest population and Groups 2, 3, and 4 were clades derived from it. In addition, we found selective sweeps for fruit flavor and improved stress resistance in different varieties of edible cherries (P. avium, P. cerasus, and P. pseudocerasus). Transcriptome analysis revealed significant differential expression of genes associated with key pathways, such as sucrose starch and sucrose metabolism, fructose and mannose metabolism, and the pentose phosphate pathway, between the leaves and fruits of P. avium. This study enhances the understanding of the evolutionary processes of the Cerasus subgenus and provides resources for functional genomics research and the improvement of edible cherries.

The haplotype-resolved T2T genome for Bauhinia × blakeana sheds light on the genetic basis of flower heterosis

BACKGROUND

The Hong Kong orchid tree Bauhinia × blakeana Dunn has long been proposed to be a sterile interspecific hybrid exhibiting flower heterosis when compared to its likely parental species, Bauhinia purpurea L. and Bauhinia variegata L. Here, we report comparative genomic and transcriptomic analyses of the 3 Bauhinia species.

FINDINGS

We generated chromosome-level assemblies for the parental species and applied a trio-binning approach to construct a haplotype-resolved telomere-to-telomere (T2T) genome for B. blakeana. Comparative chloroplast genome analysis confirmed B. purpurea as the maternal parent. Transcriptome profiling of flower tissues highlighted a closer resemblance of B. blakeana to its maternal parent. Differential gene expression analyses revealed distinct expression patterns among the 3 species, particularly in biosynthetic and metabolic processes. To investigate the genetic basis of flower heterosis observed in B. blakeana, we focused on gene expression patterns within pigment biosynthesis-related pathways. High-parent dominance and overdominance expression patterns were observed, particularly in genes associated with carotenoid biosynthesis. Additionally, allele-specific expression analysis revealed a balanced contribution of maternal and paternal alleles in shaping the gene expression patterns in B. blakeana.

CONCLUSIONS

Our study offers valuable insights into the genome architecture of hybrid B. blakeana, establishing a comprehensive genomic and transcriptomic resource for future functional genetics research within the Bauhinia genus. It also serves as a model for exploring the characteristics of hybrid species using T2T haplotype-resolved genomes, providing a novel approach to understanding genetic interactions and evolutionary mechanisms in complex genomes with high heterozygosity.

Complete Mitochondrial Genome of Niphon spinosus (Perciformes: Niphonidae): Genome Characterization and Phylogenetic Analysis

The species Niphon spinosus (Cuvier, 1829) is the only representative of the family Niphonidae and the genus Niphon, and its taxonomic history is complicated; it is still unclear in a phylogenetic sense. In this study, we report the complete mitochondrial genome of N. spinosus (OP391482), which was determined to be 16,503 bp long with biased A + T contents (53.8%) using next-generation technology. The typical set of 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and one control region (D-loop) are included in the mitochondrial genome. The H-strand encoded 28 genes (14 tRNA, 2 rRNA, and 12 PCGs), and D-loop, whereas the L-strand encoded the remaining 9 genes (8 tRNA and ND6). Its nucleotide composition, gene arrangement, codon usage patterns, and tRNA secondary structures are identical with other members of the Percoidei suborder. Furthermore, we reconstructed phylogenetic trees based on the 13 PCGs. The resulting phylogenetic trees showed N. spinosus placing as a separate lineage within the family Niphonidae, its close relationship to Trachinus draco (Trachinidae), and the clustering of major subfamilies like Luciopercinae and Percinae of the Percoidei suborder. These findings will contribute to future studies on the evolutionary history, population genetics, molecular taxonomy, and phylogeny of N. spinosus and related species.

A chromosome-level genome of Lobelia seguinii provides insights into the evolution of Campanulaceae and the lobeline biosynthesis

Lobelia seguinii is a plant with great ecological and medicinal value and belongs to Campanulaceae. Lobelia contains lobeline, a well-known compound used to treat respiratory diseases. Nevertheless, lobeline biosynthesis needs further exploration. Moreover, whole-genome duplication (WGD) and karyotype evolution within Campanulaceae still need to be better understood. In this study, we obtained a chromosome-level genome of L. seguinii with a size of 1.4 Gb and 38253 protein-coding genes. Analyses revealed two WGDs within Campanulaceae, one at the most recent common ancestor (MRCA) of Campanula and Adenophora, and another at the MRCA of Lobelioideae. Analyses further revealed that the karyotype of Platycodon grandiflorus represents the ancient type within Asterales. We proposed eight enzymes involved in the lobeline biosynthesis pathway of L. seguinii. Molecular cloning and heterologous expression of phenylalanine ammonia-lyase (PAL), a candidate enzyme involved in the first step of lobeline biosynthesis, verified its function to catalyze the deamination of phenylalanine to cinnamic acid. This study sheds light on the evolution of Campanulaceae and lobeline biosynthesis.

Prokrustean Graph: A substring index for rapid k-mer size analysis

Despite the widespread adoption of k-mer-based methods in bioinformatics, understanding the influence of k-mer sizes remains a persistent challenge. Selecting an optimal k-mer size or employing multiple k-mer sizes is often arbitrary, application-specific, and fraught with computational complexities. Typically, the influence of k-mer size is obscured by the outputs of complex bioinformatics tasks, such as genome analysis, comparison, assembly, alignment, and error correction. However, it is frequently overlooked that every method is built above a well-defined k-mer-based object like Jaccard Similarity, de Bruijn graphs, k-mer spectra, and Bray-Curtis Dissimilarity. Despite these objects offering a clearer perspective on the role of k-mer sizes, the dynamics of k-mer-based objects with respect to k-mer sizes remain surprisingly elusive. This paper introduces a computational framework that generalizes the transition of k-mer-based objects across k-mer sizes, utilizing a novel substring index, the Prokrustean graph. The primary contribution of this framework is to compute quantities associated with k-mer-based objects for all k-mer sizes, where the computational complexity depends solely on the number of maximal repeats and is independent of the range of k-mer sizes. For example, counting vertices of compacted de Bruijn graphs for k = 1 , … , 100 can be accomplished in mere seconds with our substring index constructed on a gigabase-sized read set. Additionally, we derive a space-efficient algorithm to extract the Prokrustean graph from the Burrows-Wheeler Transform. It becomes evident that modern substring indices, mostly based on longest common prefixes of suffix arrays, inherently face difficulties at exploring varying k-mer sizes due to their limitations at grouping co-occurring substrings. We have implemented four applications that utilize quantities critical in modern pangenomics and metagenomics. The code for these applications and the construction algorithm is available at https://github.com/KoslickiLab/prokrustean.

Neo-Sex Chromosome Evolution in Treehoppers Despite Long-Term X Chromosome Conservation

Sex chromosomes follow distinct evolutionary trajectories compared to the rest of the genome. In many cases, sex chromosomes (X and Y or Z and W) significantly differentiate from one another resulting in heteromorphic sex chromosome systems. Such heteromorphic systems are thought to act as an evolutionary trap that prevents subsequent turnover of the sex chromosome system. For old, degenerated sex chromosome systems, chromosomal fusion with an autosome may be one way that sex chromosomes can "refresh" their sequence content. We investigated these dynamics using treehoppers (hemipteran insects of the family Membracidae), which ancestrally have XX/X0 sex chromosomes. We assembled the most complete reference assembly for treehoppers to date for Umbonia crassicornis and employed comparative genomic analyses of 12 additional treehopper species to analyze X chromosome variation across different evolutionary timescales. We find that the X chromosome is largely conserved, with one exception being an X-autosome fusion in Calloconophora caliginosa. We also compare the ancestral treehopper X with other X chromosomes in Auchenorrhyncha (the clade containing treehoppers, leafhoppers, spittlebugs, cicadas, and planthoppers), revealing X conservation across more than 300 million years. These findings shed light on chromosomal evolution dynamics in treehoppers and the role of chromosomal rearrangements in sex chromosome evolution.

De novo sequencing allows genome-wide identification of genes involved in galactomannan synthesis in locust bean (Ceratonia siliqua)

Locust bean (Ceratonia siliqua) accumulates the galactomannan (GM) locust bean gum (LBG) in its seeds. LBG is a major industrial raw material used as a food thickener and gelling agent, whose unique properties mean that it cannot be readily replaced by other GMs. Whereas much is known about GM accumulation and the genes associated with GM biosynthesis in legumes, the genes involved in GM biosynthesis in C. siliqua are largely unknown. Here, we present a genome-wide list of genes predicted to be associated with the GM biosynthesis pathway in C. siliqua. We confirmed high GM accumulation in endosperm using a newly established GM quantification method involving LC-MS/MS. Through de novo draft genome assembly, we comprehensively identified genes predicted to be related to the GM biosynthesis pathway in C. siliqua by identifying orthologous groups. In particular, we identified all genes predicted to encode mannan synthase (ManS) and galactomannan galactosyltransferase (GMGT), enzymes functioning in the final step of GM biosynthesis, from the C. siliqua draft genome. ManS and the GMGT paralogs were predominantly expressed in endosperm. The genome and transcriptome produced in this study should facilitate research examining why C. siliqua produces LBG, unlike other legumes.

A survey of k-mer methods and applications in bioinformatics

The rapid progression of genomics and proteomics has been driven by the advent of advanced sequencing technologies, large, diverse, and readily available omics datasets, and the evolution of computational data processing capabilities. The vast amount of data generated by these advancements necessitates efficient algorithms to extract meaningful information. K-mers serve as a valuable tool when working with large sequencing datasets, offering several advantages in computational speed and memory efficiency and carrying the potential for intrinsic biological functionality. This review provides an overview of the methods, applications, and significance of k-mers in genomic and proteomic data analyses, as well as the utility of absent sequences, including nullomers and nullpeptides, in disease detection, vaccine development, therapeutics, and forensic science. Therefore, the review highlights the pivotal role of k-mers in addressing current genomic and proteomic problems and underscores their potential for future breakthroughs in research.

Long-term methylome changes after experimental seed demethylation and their interaction with recurrent water stress in Erodium cicutarium (Geraniaceae)

The frequencies and lengths of drought periods are increasing in subtropical and temperate regions worldwide. Epigenetic responses to water stress could be key for plant resilience to these largely unpredictable challenges. Experimental DNA demethylation, together with application of a stress factor is an appropriate strategy to reveal the contribution of epigenetics to plant responses to stress. We analysed leaf cytosine methylation changes in adult plants of the annual Mediterranean herb, Erodium cicutarium, in a greenhouse, after seed demethylation with 5-Azacytidine and/or recurrent water stress. We used bisulfite RADseq (BsRADseq) and a newly reported reference genome for E. cicutarium to characterize methylation changes in a 2 × 2 factorial design, controlling for plant relatedness. In the long term, 5-Azacytidine treatment alone caused both hypo- and hyper-methylation at individual cytosines, with substantial hypomethylation in CG contexts. In control conditions, drought resulted in a decrease in methylation in all but CHH contexts. In contrast, the genome of plants that experienced recurrent water stress and had been treated with 5-Azacytidine increased DNA methylation level by ca. 5%. Seed demethylation and recurrent drought produced a highly significant interaction in terms of global and context-specific cytosine methylation. Most methylation changes occurred around genic regions and within Transposable Elements. The annotation of these Differentially Methylated Regions associated with genes included several with a potential role in stress responses (e.g., PAL, CDKC, and ABCF), confirming an epigenetic contribution in response to stress at the molecular level.

Databases and computational methods for the identification of piRNA-related molecules: A survey

Piwi-interacting RNAs (piRNAs) are a class of small non-coding RNAs (ncRNAs) that plays important roles in many biological processes and major cancer diagnosis and treatment, thus becoming a hot research topic. This study aims to provide an in-depth review of computational piRNA-related research, including databases and computational models. Herein, we perform literature analysis and use comparative evaluation methods to summarize and analyze three aspects of computational piRNA-related research: (i) computational models for piRNA-related molecular identification tasks, (ii) computational models for piRNA-disease association prediction tasks, and (iii) computational resources and evaluation metrics for these tasks. This study shows that computational piRNA-related research has significantly progressed, exhibiting promising performance in recent years, whereas they also suffer from the emerging challenges of inconsistent naming systems and the lack of data. Different from other reviews on piRNA-related identification tasks that focus on the organization of datasets and computational methods, we pay more attention to the analysis of computational models, algorithms, and performances that aim to provide valuable references for computational piRNA-related identification tasks. This study will benefit the theoretical development and practical application of piRNAs by better understanding computational models and resources to investigate the biological functions and clinical implications of piRNA.

Chromosome-level genome assemblies and genetic maps reveal heterochiasmy and macrosynteny in endangered Atlantic Acropora

BACKGROUND

Over their evolutionary history, corals have adapted to sea level rise and increasing ocean temperatures, however, it is unclear how quickly they may respond to rapid change. Genome structure and genetic diversity contained within may highlight their adaptive potential.

RESULTS

We present chromosome-scale genome assemblies and linkage maps of the critically endangered Atlantic acroporids, Acropora palmata and A. cervicornis. Both assemblies and linkage maps were resolved into 14 chromosomes with their gene content and colinearity. Repeats and chromosome arrangements were largely preserved between the species. The family Acroporidae and the genus Acropora exhibited many phylogenetically significant gene family expansions. Macrosynteny decreased with phylogenetic distance. Nevertheless, scleractinians shared six of the 21 cnidarian ancestral linkage groups as well as numerous fission and fusion events compared to other distantly related cnidarians. Genetic linkage maps were constructed from one A. palmata family and 16 A. cervicornis families using a genotyping array. The consensus maps span 1,013.42 cM and 927.36 cM for A. palmata and A. cervicornis, respectively. Both species exhibited high genome-wide recombination rates (3.04 to 3.53 cM/Mb) and pronounced sex-based differences, known as heterochiasmy, with 2 to 2.5X higher recombination rates estimated in the female maps.

CONCLUSIONS

Together, the chromosome-scale assemblies and genetic maps we present here are the first detailed look at the genomic landscapes of the critically endangered Atlantic acroporids. These data sets revealed that adaptive capacity of Atlantic acroporids is not limited by their recombination rates. The sister species maintain macrosynteny with few genes with high sequence divergence that may act as reproductive barriers between them. In the Atlantic Acropora, hybridization between the two sister species yields an F1 hybrid with limited fertility despite the high levels of macrosynteny and gene colinearity of their genomes. Together, these resources now enable genome-wide association studies and discovery of quantitative trait loci, two tools that can aid in the conservation of these species.

A new species of giant gecko of the genus Rhacodactylus from New Caledonia (Squamata, Gekkota, Diplodactylidae)

We describe a new species of giant gecko, Rhacodactylus willihenkeli sp. nov. from New Caledonia. The new species is most similar in external appearance and molecular data (16S and ND2 sequences) with R. leachianus from which it differs in coloration as well as having a genetic distance of 5.0% in the 16S gene fragment and 8.5% in the ND2 fragment, respectively.

LKLPDA: A Low-Rank Fast Kernel Learning Approach for Predicting piRNA-Disease Associations

Piwi-interacting RNAs (piRNAs) are increasingly recognized as potential biomarkers for various diseases. Investig-ating the complex relationship between piRNAs and diseases through computational methods can reduce the costs and risks associated with biological experiments. Fast kernel learning (FKL) is a classical method for multi-source data fusion that is widely employed in association prediction research. However, biological networks are noisy due to the limitations of measurement technology and inherent natural variation, which can hamper the effectiveness of the network-based ideal kernel. The conventional FKL method does not address this issue. In this study, we propose a low-rank fast kernel learning (LRFKL) algorithm, which consists of low-rank representation (LRR) and the FKL algorithm. The LRFKL algorithm is designed to mitigate the effects of noise on the network-based ideal kernel. Using LRFKL, we propose a novel approach for predicting piRNA-disease associations called LKLPDA. Specifically, we first compute the similarity matrices for piRNAs and diseases. Then we use the LRFKL to fuse the similarity matrices for piRNAs and diseases separately. Finally, the LKLPDA employs AutoGluon-Tabular for predictive analysis. Computational results show that LKLPDA effectively predicts piRNA-disease associations with higher accuracy compared to previous methods. In addition, case studies confirm the reliability of the model in predicting piRNA-disease associations.

findGSEP: estimating genome size of polyploid species using k-mer frequencies

SUMMARY

Estimating genome size using k-mer frequencies, which plays a fundamental role in designing genome sequencing and analysis projects, has remained challenging for polyploid species, i.e., ploidy p > 2. To address this, we introduce "findGSEP," which is designed based on iterative curve fitting of k-mer frequencies. Precisely, it first disentangles up to p normal distributions by analyzing k-mer frequencies in whole genome sequencing of the focal species. Second, it computes the sizes of genomic regions related to 1∼p (homologous) chromosome(s) using each respective curve fitting, from which it infers the full polyploid and average haploid genome size. "findGSEP" can handle any level of ploidy p, and infer more accurate genome size than other well-known tools, as shown by tests using simulated and real genomic sequencing data of various species including octoploids.

AVAILABILITY AND IMPLEMENTATION

"findGSEP" was implemented as a web server, which is freely available at http://146.56.237.198:3838/findGSEP/. Also, "findGSEP" was implemented as an R package for parallel processing of multiple samples. Source code and tutorial on its installation and usage is available at https://github.com/sperfu/findGSEP.

Comparative Genomics of the World's Smallest Mammals Reveals Links to Echolocation, Metabolism, and Body Size Plasticity

Originating 30 million years ago, shrews (Soricidae) have diversified into around 400 species worldwide. Shrews display a wide array of adaptations, with some species having developed distinctive traits such as echolocation, underwater diving, and venomous saliva. Accordingly, these tiny insectivores are ideal to study the genomic mechanisms of evolution and adaptation. We conducted a comparative genomic analysis of four shrew species and 16 other mammals to identify genomic variations unique to shrews. Using two existing shrew genomes and two de novo assemblies for the maritime (Sorex maritimensis) and smoky (Sorex fumeus) shrews, we identified mutations in conserved regions of the genomes, also known as accelerated regions, gene families that underwent significant expansion, and positively selected genes. Our analyses unveiled shrew-specific genomic variants in genes associated with the nervous, metabolic, and auditory systems, which can be linked to unique traits in shrews. Notably, genes suggested to be under convergent evolution in echolocating mammals exhibited accelerated regions in shrews, and pathways linked to putative body size plasticity were detected. These findings provide insight into the evolutionary mechanisms shaping shrew species, shedding light on their adaptation and divergence over time.

Genome sequencing of Caridina pseudogracilirostris and its comparative analysis with malacostracan crustaceans

The Caridina pseudogracilirostris is commonly found in the brackish waters of the southwestern coastal regions of India. This study provides a comprehensive genomic investigation of the shrimp species C. pseudogracilirostris, offering insights into its genetic makeup, evolutionary dynamics, and functional annotations. The genomic DNA was isolated from tissue samples, sequenced using next-generation sequencing (NGS), and stored in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) database (Accession No: PRJNA847710). De novo sequencing indicated a genome size of 1.31 Gbp with a low heterozygosity of about 0.81%. Repeat masking and annotation revealed that repeated elements constitute 24.60% of the genome, with simple sequence repeats (SSRs) accounting for 7.26%. Gene prediction identified 14,101 genes, with functional annotations indicating involvement in critical biological processes such as development, cellular function, immunological responses, and reproduction. Furthermore, phylogenetic analysis revealed genomic links among Malacostraca species, indicating gene duplication as a strategy for genetic diversity and adaptation. C. pseudogracilirostris has 1,856 duplicated genes, reflecting a distinct genomic architecture and evolutionary strategy within the Malacostraca branch. These findings enhance our understanding of the genetic characteristics and evolutionary relationships of C. pseudogracilirostris, providing significant insights into the overall evolutionary dynamics of the Malacostraca group.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04121-4.

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.

Whole genome phylogenomics helps to resolve the phylogenetic position of the Zygothrica genus group (Diptera, Drosophilidae) and the causes of previous incongruences

Incomplete Lineage Sorting (ILS) and introgression are among the two main factors causing incongruence between gene and species trees. Advances in phylogenomic studies have allowed us to overcome most of these issues, providing reliable phylogenetic hypotheses while revealing the underlying evolutionary scenario. Across the last century, many incongruent phylogenetic reconstructions were recovered for Drosophilidae, employing a limited sampling of genetic markers or species. In these studies, the monophyly and the phylogenetic positioning of the Zygothrica genus group stood out as one of the most controversial questions. Thus, here, we addressed these issues using a phylogenomic approach, while accessing the influence of ILS and introgressions on the diversification of these species and addressing the spatio-temporal scenario associated with their evolution. For this task, the genomes of nine specimens from six Neotropical species belonging to the Zygothrica genus group were sequenced and evaluated in a phylogenetic framework encompassing other 39 species of Drosophilidae. Nucleotide and amino acid sequences recovered for a set of 2,534 single-copy genes by BUSCO were employed to reconstruct maximum likelihood (ML) concatenated and multi-species coalescent (MSC) trees. Likelihood mapping, quartet sampling, and reticulation tests were employed to infer the level and causes of incongruence. Lastly, a penalized-likelihood molecular clock strategy with fossil calibrations was performed to infer divergence times. Taken together, our results recovered the subdivision of Drosophila into six different lineages, one of which clusters species of the Zygothrica genus group (except for H. duncani). The divergence of this lineage was dated to Oligocene ∼ 31 Mya and seems to have occurred in the same timeframe as other key diversification within Drosophila. According to the concatenated and MSC strategies, this lineage is sister to the clade joining Drosophila (Siphlodora) with the Hawaiian Drosophila and Scaptomyza. Likelihood mapping, quartet sampling, reticulation reconstructions as well as introgression tests revealed that this lineage was the target of several hybridization events involving the ancestors of different Drosophila lineages. Thus, our results generally show introgression as a major source of previous incongruence. Nevertheless, the similar diversification times recovered for several of the Neotropical Drosophila lineages also support the scenario of multiple and simultaneous diversifications taking place at the base of Drosophilidae phylogeny, at least in the Neotropics.

Range-limited Heaps' law for functional DNA words in the human genome

Heaps' or Herdan-Heaps' law is a linguistic law describing the relationship between the vocabulary/dictionary size (type) and word counts (token) to be a power-law function. Its existence in genomes with certain definition of DNA words is unclear partly because the dictionary size in genome could be much smaller than that in a human language. We define a DNA word as a coding region in a genome that codes for a protein domain. Using human chromosomes and chromosome arms as individual samples, we establish the existence of Heaps' law in the human genome within limited range. Our definition of words in a genomic or proteomic context is different from other definitions such as over-represented k-mers which are much shorter in length. Although an approximate power-law distribution of protein domain sizes due to gene duplication and the related Zipf's law is well known, their translation to the Heaps' law in DNA words is not automatic. Several other animal genomes are shown herein also to exhibit range-limited Heaps' law with our definition of DNA words, though with various exponents. When tokens were randomly sampled and sample sizes reach to the maximum level, a deviation from the Heaps' law was observed, but a quadratic regression in log-log type-token plot fits the data perfectly. Investigation of type-token plot and its regression coefficients could provide an alternative narrative of reusage and redundancy of protein domains as well as creation of new protein domains from a linguistic perspective.

A chromosome-level genome assembly of Drosophila madeirensis, a fruit fly species endemic to the island of Madeira

Drosophila subobscura is distributed across Europe, the Near East, and the Americas, while its sister species, Drosophila madeirensis, is endemic to the island of Madeira in the Atlantic Ocean. D. subobscura is known for its strict light-dependence in mating and its unique courtship displays, including nuptial gift-giving. D. subobscura has also attracted the interest of researchers because of its abundant variations in chromosomal polymorphisms correlated to the latitude and season, which have been used as a tool to track global climate warming. Although D. madeirensis can be an important resource for understanding the evolutionary underpinning of these genetic characteristics of D. subobscura, little work has been done on the biology of this species. Here, we used a HiFi long-read sequencing data set to produce a de novo genome assembly for D. madeirensis. This assembly comprises a total of 111 contigs spanning 135.5 Mb and has an N50 of 24.2 Mb and a BUSCO completeness score of 98.6%. Each of the 6 chromosomes of D. madeirensis consisted of a single contig except for some centromeric regions. Breakpoints of the chromosomal inversions between D. subobscura and D. madeirensis were characterized using this genome assembly, updating some of the previously identified locations.

Comparative genomics of dusky kob (Argyrosomus japonicus, Sciaenidae) conspecifics: Evidence for speciation and the genetic mechanisms underlying traits

Dusky kob (Argyrosomus japonicus) is a commercially important finfish, indigenous to South Africa, Australia, and China. Previous studies highlighted differences in genetic composition, life history, and morphology of the species across geographic regions. A draft genome sequence of 0.742 Gb (N50 = 5.49 Mb; BUSCO completeness = 97.8%) and 22,438 predicted protein-coding genes was generated for the South African (SA) conspecific. A comparison with the Chinese (CN) conspecific revealed a core set of 32,068 orthologous protein clusters across both genomes. The SA genome exhibited 440 unique clusters compared to 1928 unique clusters in the CN genome. Transportation and immune response processes were overrepresented among the SA accessory genome, whereas the CN accessory genome was enriched for immune response, DNA transposition, and sensory detection (FDR-adjusted p < 0.01). These unique clusters may represent an adaptive component of the species' pangenome that could explain population divergence due to differential environmental specialisation. Furthermore, 700 single-copy orthologues (SCOs) displayed evidence of positive selection between the SA and CN genomes, and globally these genomes shared only 92% similarity, suggesting they might be distinct species. These genes primarily play roles in metabolism and digestion, illustrating the evolutionary pathways that differentiate the species. Understanding these genomic mechanisms underlying adaptation and evolution within and between species provides valuable insights into growth and maturation of kob, traits that are particularly relevant to commercial aquaculture.

Graphite: painting genomes using a colored de Bruijn graph

The recent growth of microbial sequence data allows comparisons at unprecedented scales, enabling the tracking of strains, mobile genetic elements, or genes. Querying a genome against a large reference database can easily yield thousands of matches that are tedious to interpret and pose computational challenges. We developed Graphite that uses a colored de Bruijn graph (cDBG) to paint query genomes, selecting the local best matches along the full query length. By focusing on the best genomic match of each query region, Graphite reduces the number of matches while providing the most promising leads for sequence tracking or genomic forensics. When applied to hundreds of Campylobacter genomes we found extensive gene sharing, including a previously undetected C. coli plasmid that matched a C. jejuni chromosome. Together, genome painting using cDBGs as enabled by Graphite, can reveal new biological phenomena by mitigating computational hurdles.

Chromosome-level genome assemblies reveal genome evolution of an invasive plant Phragmites australis

Biological invasions pose a significant threat to ecosystems, disrupting local biodiversity and ecosystem functions. The genomic underpinnings of invasiveness, however, are still largely unknown, making it difficult to predict and manage invasive species effectively. The common reed (Phragmites australis) is a dominant grass species in wetland ecosystems and has become particularly invasive when transferred from Europe to North America. Here, we present a high-quality gap-free, telomere-to-telomere genome assembly of Phragmites australis consisting of 24 pseudochromosomes and a B chromosome. Fully phased subgenomes demonstrated considerable subgenome dominance and revealed the divergence of diploid progenitors approximately 30.9 million years ago. Comparative genomics using chromosome-level scaffolds for three other lineages and a previously published draft genome assembly of an invasive lineage revealed that gene family expansions in the form of tandem duplications may have contributed to the invasiveness of the lineage. This study sheds light on the genome evolution of Arundinoideae grasses and suggests that genetic drivers, such as gene family expansions and tandem duplications, may underly the processes of biological invasion in plants. These findings provide a crucial step toward understanding and managing the genetic basis of invasiveness in plant species.

Taking advantage of reference-guided assembly in a slowly-evolving lineage: Application to Testudo graeca

BACKGROUND

Obtaining de novo chromosome-level genome assemblies greatly enhances conservation and evolutionary biology studies. For many research teams, long-read sequencing technologies (that produce highly contiguous assemblies) remain unaffordable or unpractical. For the groups that display high synteny conservation, these limitations can be overcome by a reference-guided assembly using a close relative genome. Among chelonians, tortoises (Testudinidae) are considered one of the most endangered taxa, which calls for more genomic resources. Here we make the most of high synteny conservation in chelonians to produce the first chromosome-level genome assembly of the genus Testudo with one of the most iconic tortoise species in the Mediterranean basin: Testudo graeca.

RESULTS

We used high-quality, paired-end Illumina sequences to build a reference-guided assembly with the chromosome-level reference of Gopherus evgoodei. We reconstructed a 2.29 Gb haploid genome with a scaffold N50 of 107.598 Mb and 5.37% gaps. We sequenced 25,998 protein-coding genes, and identified 41.2% of the assembly as repeats. Demographic history reconstruction based on the genome revealed two events (population decline and recovery) that were consistent with previously suggested phylogeographic patterns for the species. This outlines the value of such reference-guided assemblies for phylogeographic studies.

CONCLUSIONS

Our results highlight the value of using close relatives to produce de novo draft assemblies in species where such resources are unavailable. Our annotated genome of T. graeca paves the way to delve deeper into the species' evolutionary history and provides a valuable resource to enhance direct conservation efforts on their threatened populations.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

Genomics of ecological adaptation in Canary Island Descurainia (Brassicaceae) and comparisons with other Brassicaceae

Oceanic archipelagos provide striking examples of lineages that have radiated over pronounced ecological gradients. Accompanying this diversification, lineages have evolved adaptations allowing survival in extreme environments. Here, we investigate the genomic basis of ecological adaptation in Canary Island Descurainia (Brassicaceae), an island relative of Arabidopsis. The seven endemic species have diversified in situ along an elevational and ecological gradient, from low-elevation scrub to high-elevation sub-alpine desert. We first generated a reference genome for Descurainia millefolia, phylogenetic analysis of which placed it as sister to D. sophioides. Ninety-six gene families were found to be specific to D. millefolia and a further 1087 and 1469 gene families have expanded or contracted in size, respectively, along the D. millefolia branch. We then employed genome re-sequencing to sample 14 genomes across the seven species of Canary Island Descurainia and an outgroup. Phylogenomic analyses were consistent with previous reconstructions of Canary Island Descurainia in resolving low- and high-elevation clades. Using the branch-site dN/dS method, we detected positive selection for 275 genes on the branch separating the low- and high-elevation species and these positively selected genes (PSGs) were significantly enriched for functions related to reproduction and stress tolerance. Comparing PSGs to those in analyses of adaptation to elevation and/or latitude in other Brassicaceae, we found little evidence of widespread convergence and gene reuse, except for two examples, one of which was a significant overlap between Descurainia and Draba nivalis, a species restricted to high latitudes. The study of Canary Island Descurainia suggests that the transition to high-elevation environments such as that found in the high mountains of the Canary Islands involves selection on genes related to reproduction and stress tolerance but that repeated evolution across different lineages that have evolved into similar habitats is limited, indicating substantially different molecular trajectories to adaptation.

Development of a Reference Transcriptome and Identification of Differentially Expressed Genes Linked to Salt Stress in Salt Marsh Grass (Sporobolus alterniflorus) along Delaware Coastal Regions

Salt marsh grass (Sporobolus alterniflorus) plays a crucial role in Delaware coastal regions by serving as a physical barrier between land and water along the inland bays and beaches. This vegetation helps to stabilize the shoreline and prevent erosion, protecting the land from the powerful forces of the waves and tides. In addition to providing a physical barrier, salt marsh grass is responsible for filtering nutrients in the water, offering an environment for aquatic species and presenting a focal point of study for high salt tolerance in plants. As seawater concentrations vary along the Delaware coast from low to medium to high salinity, our study seeks to identify the impact of salt tolerance in marsh grass and to identify genes associated with salt tolerance levels. We developed more than 211,000 next-generation-sequencing (Illumina) transcriptomic reads to create a reference transcriptome from low-, medium-, and high-salinity marsh grass leaf samples collected from the Delaware coastline. Contiguous sequences were annotated based on a homology search using BLASTX against rice (Oryza sativa), foxtail millet (Setaria italica), and non-redundant species within the Viridiplantae database. Additionally, we identified differentially expressed genes related to salinity stress as candidates for salt stress qPCR analysis. The data generated from this study may help to elucidate the genetic signatures and physiological responses of plants to salinity stress, thereby offering valuable insight into the use of innovative approaches for gene expression studies in crops that are less salt tolerant.

A high-quality chromosome-level genome assembly of the endangered tree Kmeria septentrionalis

Kmeria septentrionalis is a critically endangered tree endemic to Guangxi, China, and is listed on the International Union for Conservation of Nature's Red List. The lack of genetic information and high-quality genome data has hindered conservation efforts and studies on this species. In this study, we present a chromosome-level genome assembly of K. septentrionalis. The genome was initially assembled to be 2.57 Gb, with a contig N50 of 11.93 Mb. Hi-C guided genome assembly allowed us to anchor 98.83% of the total length of the initial contigs onto 19 pseudochromosomes, resulting in a scaffold N50 of 135.08 Mb. The final chromosome-level genome, spaning 2.54 Gb, achieved a BUSCO completeness of 98.9% and contained 1.67 Gb repetitive elements and 35,927 coding genes. This high-quality genome assembly provides a valuable resource for understanding the genetic basis of conservation-related traits and biological properties of this endangered tree species. Furthermore, it lays a critical foundation for evolutionary studies within the Magnoliaceae family.

Enhancing SNV identification in whole-genome sequencing data through the incorporation of known genetic variants into the minimap2 index

MOTIVATION

Alignment of reads to a reference genome sequence is one of the key steps in the analysis of human whole-genome sequencing data obtained through Next-generation sequencing (NGS) technologies. The quality of the subsequent steps of the analysis, such as the results of clinical interpretation of genetic variants or the results of a genome-wide association study, depends on the correct identification of the position of the read as a result of its alignment. The amount of human NGS whole-genome sequencing data is constantly growing. There are a number of human genome sequencing projects worldwide that have resulted in the creation of large-scale databases of genetic variants of sequenced human genomes. Such information about known genetic variants can be used to improve the quality of alignment at the read alignment stage when analysing sequencing data obtained for a new individual, for example, by creating a genomic graph. While existing methods for aligning reads to a linear reference genome have high alignment speed, methods for aligning reads to a genomic graph have greater accuracy in variable regions of the genome. The development of a read alignment method that takes into account known genetic variants in the linear reference sequence index allows combining the advantages of both sets of methods.

RESULTS

In this paper, we present the minimap2_index_modifier tool, which enables the construction of a modified index of a reference genome using known single nucleotide variants and insertions/deletions (indels) specific to a given human population. The use of the modified minimap2 index improves variant calling quality without modifying the bioinformatics pipeline and without significant additional computational overhead. Using the PrecisionFDA Truth Challenge V2 benchmark data (for HG002 short-read data aligned to the GRCh38 linear reference (GCA_000001405.15) with parameters k = 27 and w = 14) it was demonstrated that the number of false negative genetic variants decreased by more than 9500, and the number of false positives decreased by more than 7000 when modifying the index with genetic variants from the Human Pangenome Reference Consortium.

The genomic and cellular basis of biosynthetic innovation in rove beetles

How evolution at the cellular level potentiates macroevolutionary change is central to understanding biological diversification. The >66,000 rove beetle species (Staphylinidae) form the largest metazoan family. Combining genomic and cell type transcriptomic insights spanning the largest clade, Aleocharinae, we retrace evolution of two cell types comprising a defensive gland-a putative catalyst behind staphylinid megadiversity. We identify molecular evolutionary steps leading to benzoquinone production by one cell type via a mechanism convergent with plant toxin release systems, and synthesis by the second cell type of a solvent that weaponizes the total secretion. This cooperative system has been conserved since the Early Cretaceous as Aleocharinae radiated into tens of thousands of lineages. Reprogramming each cell type yielded biochemical novelties enabling ecological specialization-most dramatically in symbionts that infiltrate social insect colonies via host-manipulating secretions. Our findings uncover cell type evolutionary processes underlying the origin and evolvability of a beetle chemical innovation.

De-novo assembly of four rail (Aves: Rallidae) genomes: A resource for comparative genomics

Rails are a phenotypically diverse family of birds that includes 130 species and displays a wide distribution around the world. Here we present annotated genome assemblies for four rails from Aotearoa New Zealand: two native volant species, pūkeko Porphyrio melanotus and mioweka Gallirallus philippensis, and two endemic flightless species takahē Porphyrio hochstetteri and weka Gallirallus australis. Using the sequence read data, heterozygosity was found to be lowest in the endemic flightless species and this probably reflects their relatively small populations. The quality checks and comparison with other rallid genomes showed that the new assemblies were of good quality. This study significantly increases the number of available rallid genomes and will enable future genomic studies on the evolution of this family.

Scalable de novo classification of antibiotic resistance of Mycobacterium tuberculosis

MOTIVATION

World Health Organization estimates that there were over 10 million cases of tuberculosis (TB) worldwide in 2019, resulting in over 1.4 million deaths, with a worrisome increasing trend yearly. The disease is caused by Mycobacterium tuberculosis (MTB) through airborne transmission. Treatment of TB is estimated to be 85% successful, however, this drops to 57% if MTB exhibits multiple antimicrobial resistance (AMR), for which fewer treatment options are available.

RESULTS

We develop a robust machine-learning classifier using both linear and nonlinear models (i.e. LASSO logistic regression (LR) and random forests (RF)) to predict the phenotypic resistance of Mycobacterium tuberculosis (MTB) for a broad range of antibiotic drugs. We use data from the CRyPTIC consortium to train our classifier, which consists of whole genome sequencing and antibiotic susceptibility testing (AST) phenotypic data for 13 different antibiotics. To train our model, we assemble the sequence data into genomic contigs, identify all unique 31-mers in the set of contigs, and build a feature matrix M, where M[i, j] is equal to the number of times the ith 31-mer occurs in the jth genome. Due to the size of this feature matrix (over 350 million unique 31-mers), we build and use a sparse matrix representation. Our method, which we refer to as MTB++, leverages compact data structures and iterative methods to allow for the screening of all the 31-mers in the development of both LASSO LR and RF. MTB++ is able to achieve high discrimination (F-1 >80%) for the first-line antibiotics. Moreover, MTB++ had the highest F-1 score in all but three classes and was the most comprehensive since it had an F-1 score >75% in all but four (rare) antibiotic drugs. We use our feature selection to contextualize the 31-mers that are used for the prediction of phenotypic resistance, leading to some insights about sequence similarity to genes in MEGARes. Lastly, we give an estimate of the amount of data that is needed in order to provide accurate predictions.

AVAILABILITY

The models and source code are publicly available on Github at https://github.com/M-Serajian/MTB-Pipeline.

First whole genome sequencing data of a Mycobacterium tuberculosis STB-T1A strain isolated from a spinal tuberculosis patient in Sabah, Malaysia

Spinal tuberculosis, also referred to as Pott's disease, presents a significant risk of severe paralysis if not promptly detected and treated, owing to complications such as spinal cord compression and deformity. This article presents the genetic analysis of a Mycobacterium tuberculosis STB-T1A strain, isolated from the spine of a 29-year-old female diagnosed with spinal tuberculosis. Genomic DNA was extracted from pure culture and subjected to sequencing using the Illumina NovaSeq 6000 sequencing system. The genome of the M. tuberculosis STB-T1A strain spans 4,367,616 base pairs with a G+C content of 65.56 % and 4174 protein-coding genes. Comparative genomic analysis, conducted via single nucleotide polymorphism (SNP)-based phylogenetic analysis using the Maximum Likelihood method, revealed that the strain falls within the Indo-Oceanic lineage (Lineage 1). It clusters with the M. tuberculosis 43-16836 strain, which was isolated from the cerebrospinal fluid of a patient with tuberculous meningitis in Thailand. The complete genome sequence has been deposited at the National Center for Biotechnology Information (NCBI) GenBank database with the accession number JBBMVZ000000000.

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.

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.

Genomic factors shape carbon and nitrogen metabolic niche breadth across Saccharomycotina yeasts

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Two general paradigms have been proposed to explain this variation: (i) trade-offs between performance efficiency and breadth and (ii) the joint influence of extrinsic (environmental) and intrinsic (genomic) factors. We assembled genomic, metabolic, and ecological data from nearly all known species of the ancient fungal subphylum Saccharomycotina (1154 yeast strains from 1051 species), grown in 24 different environmental conditions, to examine niche breadth evolution. We found that large differences in the breadth of carbon utilization traits between yeasts stem from intrinsic differences in genes encoding specific metabolic pathways, but we found limited evidence for trade-offs. These comprehensive data argue that intrinsic factors shape niche breadth variation in microbes.