Long walk to genomics: History and current approaches to genome sequencing and assembly

Next-generation data filtering in the genomics era

Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering - removing sequencing bases, reads, genetic variants and/or individuals from a dataset - to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy-Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima's D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne).

Genetic structure, selective characterization and specific molecular identity cards of high-yielding Houdan chickens based on genome-wide SNP

The high-yielding Houdan chicken (GGF) is characterized by high egg production and disease resistance. This study conducted whole genome resequencing of the GGF population and compared it to data from other breeds. Genetic diversity analysis revealed higher observed heterozygosity (Ho), Polymorphism information content (PIC), number of runs of homozygosity (ROH), and inbreeding coefficient (FROH) in GGF. Linkage disequilibrium (LD) decay was slowest in GGF, indicating intensive inbreeding and strong selection. These findings suggest a need for appropriate strategies to enhance genetic diversity conservation in this breed. Population structure analysis demonstrated that GGF was genetically distinct from both the red jungle fowl (RJF) and Chinese indigenous chicken (CIC) populations, highlighting GGF as a unique genetic resource warranting intensive protection and utilization. Selective sweep analysis identified genes under selection in GGF, primarily enriched in signaling pathways related to oocyte meiosis and progesterone-mediated oocyte maturation. Key candidate genes included: CCNE1, SKP1, CDC20, CDK2, ADCY8, RPS6KA6, PPP3CB, PDE3B, HSP90AB1, and AKT3. These findings provide a theoretical foundation for their potential application in poultry breeding. Additionally, this study combined bioinformatics analysis with PCR amplification and Sanger sequencing to identify 4 SNPs that can serve as a molecular identity card (ID) for GGF: SNP1 (Chr2: 136130976), SNP3 (Chr4:11705164), SNP4 (Chr4: 63255588), and SNP5 (Chr24: 3271008). This study provides a scientific basis for effective management and conservation of GGF genetic resources, and establishes a simple, economical, and accurate set of molecular IDs to combat the proliferation of inferior breeds and protect genetic resources.

Exploring the molecular aspect and updating evolutionary approaches to the DNA polymerase enzymes for biotechnological needs: A comprehensive review

DNA polymerases are essential enzymes that play a key role in living organisms, as they participate in the synthesis and maintenance of the DNA molecule. The intrinsic properties of these enzymes have been widely observed and studied to understand their functions, activities, and behavior, which has allowed their natural power in DNA synthesis to be exploited in modern biotechnology, to the point of making them true pillars of the field. In this context, the laboratory evolution of these enzymes, either by directed evolution or rational design, has led to the generation of a wide range of new DNA polymerases with novel properties, suitable for a variety of biotechnological needs. In this review, we examine DNA polymerases at the molecular level, their biotechnological use, and their evolutionary methods in relation to the novel properties sought, providing a chronological selection of evolved DNA polymerases cited in the literature that we consider to be of great interest. To our knowledge, this work is the first to bring together the molecular, functional and evolutionary aspects of the DNA polymerase enzyme. We believe it will be of great interest to researchers whose aim is to produce new lines of evolved DNA polymerases.

Advancement of Next-Generation DNA Sequencing through Ionic Blockade and Transverse Tunneling Current Methods

DNA sequencing is transforming the field of medical diagnostics and personalized medicine development by providing a pool of genetic information. Recent advancements have propelled solid-state material-based sequencing into the forefront as a promising next-generation sequencing (NGS) technology, offering amplification-free, cost-effective, and high-throughput DNA analysis. Consequently, a comprehensive framework for diverse sequencing methodologies and a cross-sectional understanding with meticulous documentation of the latest advancements is of timely need. This review explores a broad spectrum of progress and accomplishments in the field of DNA sequencing, focusing mainly on electrical detection methods. The review delves deep into both the theoretical and experimental demonstrations of the ionic blockade and transverse tunneling current methods across a broad range of device architectures, nanopore, nanogap, nanochannel, and hybrid/heterostructures. Additionally, various aspects of each architecture are explored along with their strengths and weaknesses, scrutinizing their potential applications for ultrafast DNA sequencing. Finally, an overview of existing challenges and future directions is provided to expedite the emergence of high-precision and ultrafast DNA sequencing with ionic and transverse current approaches.

Genome size and chromosome number are critical metrics for accurate genome assembly assessment in Eukaryota

The number of genome assemblies has rapidly increased in recent history, with NCBI databases reaching over 41,000 eukaryotic genome assemblies across about 2,300 species. Increases in read length and improvements in assembly algorithms have led to increased contiguity and larger genome assemblies. While this number of assemblies is impressive, only about a third of these assemblies have corresponding genome size estimations for their respective species on publicly available databases. In this paper, genome assemblies are assessed regarding their total size compared to their respective publicly available genome size estimations. These deviations in size are assessed related to genome size, kingdom, sequencing platform, and standard assembly metrics, such as N50 and BUSCO values. A large proportion of assemblies deviate from their estimated genome size by more than 10%, with increasing deviations in size with increased genome size, suggesting nonprotein coding and structural DNA may be to blame. Modest differences in performance of sequencing platforms are noted as well. While standard metrics of genome assessment are more likely to indicate an assembly approaching the estimated genome size, much of the variation in this deviation in size is not explained with these raw metrics. A new, proportional N50 metric is proposed, in which N50 values are made relative to the average chromosome size of each species. This new metric has a stronger relationship with complete genome assemblies and, due to its proportional nature, allows for a more direct comparison across assemblies for genomes with variation in sizes and architectures.

Sequence confirmation of synthetic DNA exceeding 100 nucleotides using restriction enzyme mediated digestion combined with high-resolution tandem mass spectrometry

Oligonucleotides have emerged as important therapeutic options for inherited diseases. In recent years, RNA therapeutics, especially mRNA, have been pushed to the market. Analytical methods for these molecules have been published extensively in the last few years. Notably, mass spectrometry has proven as a state-of-the-art quality control method. For RNA based therapeutics, numerous methods are available, while DNA therapeutics lack of suitable MS-based methods when it comes to molecules exceeding approximately 60 nucleotides. We present a method which combines the use of common restriction enzymes and short enzyme-directing oligonucleotides to generate DNA digestion products with the advantages of high-resolution tandem mass spectrometry. The instrumentation includes ion pair reverse phase chromatography coupled to a time-of-flight mass spectrometer with a collision induced dissociation (CID) for sequence analysis. Utilizing this approach, we increased the sequence coverage from 23.3% for a direct CID-MS/MS experiment of a 100 nucleotide DNA molecule to 100% sequence coverage using the restriction enzyme mediated approach presented in this work. This approach is suitable for research and development and quality control purposes in a regulated environment, which makes it a versatile tool for drug development.

Introducing 'identification probability' for automated and transferable assessment of metabolite identification confidence in metabolomics and related studies

Methods for assessing compound identification confidence in metabolomics and related studies have been debated and actively researched for the past two decades. The earliest effort in 2007 focused primarily on mass spectrometry and nuclear magnetic resonance spectroscopy and resulted in four recommended levels of metabolite identification confidence - the Metabolite Standards Initiative (MSI) Levels. In 2014, the original MSI Levels were expanded to five levels (including two sublevels) to facilitate communication of compound identification confidence in high resolution mass spectrometry studies. Further refinement in identification levels have occurred, for example to accommodate use of ion mobility spectrometry in metabolomics workflows, and alternate approaches to communicate compound identification confidence also have been developed based on identification points schema. However, neither qualitative levels of identification confidence nor quantitative scoring systems address the degree of ambiguity in compound identifications in context of the chemical space being considered, are easily automated, or are transferable between analytical platforms. In this perspective, we propose that the metabolomics and related communities consider identification probability as an approach for automated and transferable assessment of compound identification and ambiguity in metabolomics and related studies. Identification probability is defined simply as 1/N, where N is the number of compounds in a reference library or chemical space that match to an experimentally measured molecule within user-defined measurement precision(s), for example mass measurement or retention time accuracy, etc. We demonstrate the utility of identification probability in an in silico analysis of multi-property reference libraries constructed from the Human Metabolome Database and computational property predictions, provide guidance to the community in transparent implementation of the concept, and invite the community to further evaluate this concept in parallel with their current preferred methods for assessing metabolite identification confidence.

A practical introduction to holo-omics

Holo-omics refers to the joint study of non-targeted molecular data layers from host-microbiota systems or holobionts, which is increasingly employed to disentangle the complex interactions between the elements that compose them. We navigate through the generation, analysis, and integration of omics data, focusing on the commonalities and main differences to generate and analyze the various types of omics, with a special focus on optimizing data generation and integration. We advocate for careful generation and distillation of data, followed by independent exploration and analyses of the single omic layers to obtain a better understanding of the study system, before the integration of multiple omic layers in a final model is attempted. We highlight critical decision points to achieve this aim and flag the main challenges to address complex biological questions regarding the integrative study of host-microbiota relationships.

LSTrAP-denovo: Automated Generation of Transcriptome Atlases for Eukaryotic Species Without Genomes

Despite the abundance of species with transcriptomic data, a significant number of species still lack sequenced genomes, making it difficult to study gene function and expression in these organisms. While de novo transcriptome assembly can be used to assemble protein-coding transcripts from RNA-sequencing (RNA-seq) data, the datasets used often only feature samples of arbitrarily selected or similar experimental conditions, which might fail to capture condition-specific transcripts. We developed the Large-Scale Transcriptome Assembly Pipeline for de novo assembled transcripts (LSTrAP-denovo) to automatically generate transcriptome atlases of eukaryotic species. Specifically, given an NCBI TaxID, LSTrAP-denovo can (1) filter undesirable RNA-seq accessions based on read data, (2) select RNA-seq accessions via unsupervised machine learning to construct a sample-balanced dataset for download, (3) assemble transcripts via over-assembly, (4) functionally annotate coding sequences (CDS) from assembled transcripts and (5) generate transcriptome atlases in the form of expression matrices for downstream transcriptomic analyses. LSTrAP-denovo is easy to implement, written in Python, and is freely available at https://github.com/pengkenlim/LSTrAP-denovo/.

Escherichia coli non-coding regulatory regions are highly conserved

Microbial genome sequences are rapidly accumulating, enabling large-scale studies of sequence variation. Existing studies primarily focus on coding regions to study amino acid substitution patterns in proteins. However, non-coding regulatory regions also play a distinct role in determining physiologic responses. To investigate intergenic sequence variation on a large-scale, we identified non-coding regulatory region alleles across 2350 Escherichia coli strains. This 'alleleome' consists of 117 781 unique alleles for 1169 reference regulatory regions (transcribing 1975 genes) at single base-pair resolution. We find that 64% of nucleotide positions are invariant, and variant positions vary in a median of just 0.6% of strains. Additionally, non-coding alleles are sufficient to recover E. coli phylogroups. We find that core promoter elements and transcription factor binding sites are significantly conserved, especially those located upstream of essential or highly-expressed genes. However, variability in conservation of transcription factor binding sites is significant both within and across regulons. Finally, we contrast mutations acquired during adaptive laboratory evolution with wild-type variation, finding that the former preferentially alter positions that the latter conserves. Overall, this analysis elucidates the wealth of information found in E. coli non-coding sequence variation and expands pangenomic studies to non-coding regulatory regions at single-nucleotide resolution.

Application of Next-Generation Sequencing (NGS) Techniques for Selected Companion Animals

Next-Generation Sequencing (NGS) techniques have revolutionized veterinary medicine for cats and dogs, offering insights across various domains. In veterinary parasitology, NGS enables comprehensive profiling of parasite populations, aiding in understanding transmission dynamics and drug resistance mechanisms. In infectious diseases, NGS facilitates rapid pathogen identification, characterization of virulence factors, and tracking of outbreaks. Moreover, NGS sheds light on metabolic processes by elucidating gene expression patterns and metabolic pathways, essential for diagnosing metabolic disorders and designing tailored treatments. In autoimmune diseases, NGS helps identify genetic predispositions and molecular mechanisms underlying immune dysregulation. Veterinary oncology benefits from NGS through personalized tumor profiling, mutation analysis, and identification of therapeutic targets, fostering precision medicine approaches. Additionally, NGS plays a pivotal role in veterinary genetics, unraveling the genetic basis of inherited diseases and facilitating breeding programs for healthier animals. Physiological investigations leverage NGS to explore complex biological systems, unraveling gene-environment interactions and molecular pathways governing health and disease. Application of NGS in treatment planning enhances precision and efficacy by enabling personalized therapeutic strategies tailored to individual animals and their diseases, ultimately advancing veterinary care for companion animals.

Comparative Genomics Unveils Functional Diversity, Pangenome Openness, and Underlying Biological Drivers among Bacillus subtilis Group

The Bacillus subtilis group (Bs group), with Bacillus subtilis as its core species, holds significant research and economic value in various fields, including science, industrial production, food, and pharmaceuticals. However, most studies have been confined to comparative genomics analyses and exploration within individual genomes at the level of species, with few conducted within groups across different species. This study focused on Bacillus subtilis, the model of Gram-positive bacteria, and 14 other species with significant research value, employing comparative pangenomics as well as population enrichment analysis to ascertain the functional enrichment and diversity. Through the quantification of pangenome openness, this work revealed the underlying biological drivers and significant correlation between pangenome openness and various factors, including the distribution of toxin-antitoxin- and integrase-related genes, as well as the number of endonucleases, recombinases, repair system-related genes, prophages, integrases, and transfer mobile elements. Furthermore, the functional enrichment results indicated the potential for secondary metabolite, probiotic, and antibiotic exploration in Bacillus licheniformis, Bacillus paralicheniformis, and Bacillus spizizenii, respectively. In general, this work systematically exposed the quantification of pangenome openness, biological drivers, the pivotal role of genomic instability factors, and mobile elements, providing targeted exploration guidance for the Bs group.

Genetic Characterization of Lumpy Skin Disease Viruses Circulating in Lesotho Cattle

Lumpy skin disease is one of the fast-spreading viral diseases of cattle and buffalo that can potentially cause severe economic impact. Lesotho experienced LSD for the first time in 1947 and episodes of outbreaks occurred throughout the decades. In this study, eighteen specimens were collected from LSD-clinically diseased cattle between 2020 and 2022 from Mafeteng, Leribe, Maseru, Berea, and Mohales' Hoek districts of Lesotho. A total of 11 DNA samples were analyzed by PCR and sequencing of the extracellular enveloped virus (EEV) glycoprotein, G-protein-coupled chemokine receptor (GPCR), 30 kDa RNA polymerase subunit (RPO30), and B22R genes. All nucleotide sequences of the above-mentioned genes confirmed that the PCR amplicons of clinical samples are truly LSDV, as they were identical to respective LSDV isolates on the NCBI GenBank. Two of the elevem samples were further characterized by whole-genome sequencing. The analysis, based on both CaPV marker genes and complete genome sequences, revealed that the LSDV isolates from Lesotho cluster with the NW-like LSDVs, which includes the commonly circulating LSDV field isolates from Africa, the Middle East, the Balkans, Turkey, and Eastern Europe.

Modification of coronary artery disease clinical risk factors by coronary artery disease polygenic risk score

BACKGROUND

The extent to which the relationships between clinical risk factors and coronary artery disease (CAD) are altered by CAD polygenic risk score (PRS) is not well understood. Here, we determine whether the interactions between clinical risk factors and CAD PRS further explain risk for incident CAD.

METHODS

Participants were of European ancestry from the UK Biobank without prevalent CAD. An externally trained genome-wide CAD PRS was generated and then applied. Clinical risk factors were ascertained at baseline. Cox proportional hazards models were fitted to examine the incident CAD effects of CAD PRS, risk factors, and their interactions. Next, the PRS and risk factors were stratified to investigate the attributable risk of clinical risk factors.

FINDINGS

A total of 357,144 individuals of European ancestry without prevalent CAD were included. During a median of 11.1 years of follow-up (interquartile range 10.4-14.1 years), CAD PRS was associated with 1.35-fold (95% confidence interval [CI] 1.332-1.368) risk per SD for incident CAD. The prognostic relevance of the following risk factors was relatively diminished for those with high CAD PRS on a continuous scale: type 2 diabetes (hazard ratio [HR]interaction 0.91, 95% CIinteraction 0.88-0.94), increased body mass index (HRinteraction 0.97, 95% CIinteraction 0.96-0.98), and increased C-reactive protein (HRinteraction 0.98, 95% CIinteraction 0.96-0.99). However, a high CAD PRS yielded joint risk increases with low-density lipoprotein cholesterol (HRinteraction 1.05, 95% CIinteraction 1.04-1.06) and total cholesterol (HRinteraction 1.05, 95% CIinteraction 1.03-1.06).

CONCLUSION

The CAD PRS is associated with incident CAD, and its application improves the prognostic relevance of several clinical risk factors.

FUNDING

P.N. (R01HL127564, R01HL151152, and U01HG011719) is supported by the National Institutes of Health.

Translating genomic advances into biodiversity conservation

A key action of the new Global Biodiversity Framework is the maintenance of genetic diversity in all species to safeguard their adaptive potential. To achieve this goal, a translational mindset, which aims to convert results of basic research into direct practical benefits, needs to be applied to biodiversity conservation. Despite much discussion on the value of genomics to conservation, a disconnect between those generating genomic resources and those applying it to biodiversity management remains. As global efforts to generate reference genomes for non-model species increase, investment into practical biodiversity applications is critically important. Applications such as understanding population and multispecies diversity and longitudinal monitoring need support alongside education for policymakers on integrating the data into evidence-based decisions. Without such investment, the opportunity to revolutionize global biodiversity conservation using genomics will not be fully realized.

The promise and pitfalls of synteny in phylogenomics

Reconstructing the tree of life remains a central goal in biology. Early methods, which relied on small numbers of morphological or genetic characters, often yielded conflicting evolutionary histories, undermining confidence in the results. Investigations based on phylogenomics, which use hundreds to thousands of loci for phylogenetic inquiry, have provided a clearer picture of life's history, but certain branches remain problematic. To resolve difficult nodes on the tree of life, 2 recent studies tested the utility of synteny, the conserved collinearity of orthologous genetic loci in 2 or more organisms, for phylogenetics. Synteny exhibits compelling phylogenomic potential while also raising new challenges. This Essay identifies and discusses specific opportunities and challenges that bear on the value of synteny data and other rare genomic changes for phylogenomic studies. Synteny-based analyses of highly contiguous genome assemblies mark a new chapter in the phylogenomic era and the quest to reconstruct the tree of life.

Detection accuracy and clinical applications of DP-TOF mass spectrometry

OBJECTIVE

Matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS) is currently used in clinical microbiology laboratories. This study aimed to determine whether dual-polarity time-of-flight mass spectrometry (DP-TOF MS) could be applied to clinical nucleotide detection.

METHODS

This prospective study included 40 healthy individuals and 110 patients diagnosed with cardiovascular diseases. We used DP-TOF MS and Sanger sequencing to evaluate 17 loci across 11 genes associated with cardiovascular drug responses. In addition, we used DP-TOF MS to test 998 retrospectively collected clinical DNA samples with known results.

RESULTS

A, T, and G nucleotide detection by DP-TOF MS and Sanger sequencing revealed 100% concordance, whereas the C nucleotide concordance was 99.86%. Genotyping based on the results of the two methods showed 99.96% concordance. Regarding clinical applications, DP-TOF MS yielded a 99.91% concordance rate for known loci. The minimum detection limit for DNA was 0.4 ng; the inter-assay and intra-assay precision rates were both 100%. Anti-interference analysis showed that aerosol contamination greater than 1013 copies/µL in the laboratory environment could influence the results of DP-TOF MS.

CONCLUSIONS

The DP-TOF MS platform displayed good detection performance, as demonstrated by its 99.96% concordance rate with Sanger sequencing. Thus, it may be applied to clinical nucleotide detection.

Genomic and phenotypic analyses reveal Paenibacillus polymyxa PJH16 is a potential biocontrol agent against cucumber fusarium wilt

In recent years, bacterial-based biocontrol agents (BCA) have become a new trend for the control of fungal diseases such as fusarium wilt that seriously threaten the yield and quality of cucumber, which are transmitted through infested soil and water. This study was set out with the aim of figuring the mechanism of the isolated rhizobacterial strain Paenibacillus polymyxa PJH16 in preventing Fusarium oxysporum f. sp. cucumerinum (Foc). Biocontrol and growth-promoting experiments revealed that bacterial strain causes effective inhibition of the fungal disease through a significant growth-promoting ability of plants, and had activities of β-1,3-glucanase, cellulase, amylase and protease. It could produce siderophore and indole-3-acetic acid, too. Using the high-throughput sequencing tool PacBio Sequel II system and the database annotation, the bacterial strain was identified as P. polymyxa PJH16 and contained genes encoding for presence of biofilm formation, antimicrobial peptides, siderophores and hydrolyases. From comparing data between the whole genome of P. polymyxa PJH16 with four closely related P. polymyxa strains, findings revealed markedly the subtle differences in their genome sequences and proposed new antifungal substances present in P. polymyxa PJH16. Therefore, P. polymyxa PJH16 could be utilized in bioengineering a microbial formulation for application as biocontrol agent and bio-stimulant, in the future.

Implementing Whole Genome Sequencing (WGS) in Clinical Practice: Advantages, Challenges, and Future Perspectives

The integration of whole genome sequencing (WGS) into all aspects of modern medicine represents the next step in the evolution of healthcare. Using this technology, scientists and physicians can observe the entire human genome comprehensively, generating a plethora of new sequencing data. Modern computational analysis entails advanced algorithms for variant detection, as well as complex models for classification. Data science and machine learning play a crucial role in the processing and interpretation of results, using enormous databases and statistics to discover new and support current genotype-phenotype correlations. In clinical practice, this technology has greatly enabled the development of personalized medicine, approaching each patient individually and in accordance with their genetic and biochemical profile. The most propulsive areas include rare disease genomics, oncogenomics, pharmacogenomics, neonatal screening, and infectious disease genomics. Another crucial application of WGS lies in the field of multi-omics, working towards the complete integration of human biomolecular data. Further technological development of sequencing technologies has led to the birth of third and fourth-generation sequencing, which include long-read sequencing, single-cell genomics, and nanopore sequencing. These technologies, alongside their continued implementation into medical research and practice, show great promise for the future of the field of medicine.

Recent advances in plant translational genomics for crop improvement

The growing population, climate change, and limited agricultural resources put enormous pressure on agricultural systems. A plateau in crop yields is occurring and extreme weather events and urbanization threaten the livelihood of farmers. It is imperative that immediate attention is paid to addressing the increasing food demand, ensuring resilience against emerging threats, and meeting the demand for more nutritious, safer food. Under uncertain conditions, it is essential to expand genetic diversity and discover novel crop varieties or variations to develop higher and more stable yields. Genomics plays a significant role in developing abundant and nutrient-dense food crops. An alternative to traditional breeding approach, translational genomics is able to improve breeding programs in a more efficient and precise manner by translating genomic concepts into practical tools. Crop breeding based on genomics offers potential solutions to overcome the limitations of conventional breeding methods, including improved crop varieties that provide more nutritional value and are protected from biotic and abiotic stresses. Genetic markers, such as SNPs and ESTs, contribute to the discovery of QTLs controlling agronomic traits and stress tolerance. In order to meet the growing demand for food, there is a need to incorporate QTLs into breeding programs using marker-assisted selection/breeding and transgenic technologies. This chapter primarily focuses on the recent advances that are made in translational genomics for crop improvement and various omics techniques including transcriptomics, metagenomics, pangenomics, single cell omics etc. Numerous genome editing techniques including CRISPR Cas technology and their applications in crop improvement had been discussed.

Data Storage Using DNA

The exponential growth of global data has outpaced the storage capacities of current technologies, necessitating innovative storage strategies. DNA, as a natural medium for preserving genetic information, has emerged as a highly promising candidate for next-generation storage medium. Storing data in DNA offers several advantages, including ultrahigh physical density and exceptional durability. Facilitated by significant advancements in various technologies, such as DNA synthesis, DNA sequencing, and DNA nanotechnology, remarkable progress has been made in the field of DNA data storage over the past decade. However, several challenges still need to be addressed to realize practical applications of DNA data storage. In this review, the processes and strategies of in vitro DNA data storage are first introduced, highlighting recent advancements. Next, a brief overview of in vivo DNA data storage is provided, with a focus on the various writing strategies developed to date. At last, the challenges encountered in each step of DNA data storage are summarized and promising techniques are discussed that hold great promise in overcoming these obstacles.

Wildlife endogenous retroviruses: colonization, consequences, and cooption

Endogenous retroviruses (ERVs) are inherited genomic remains of past germline retroviral infections. Research on human ERVs has focused on medical implications of their dysregulation on various diseases. However, recent studies incorporating wildlife are yielding remarkable perspectives on long-term retrovirus-host interactions. These initial forays into broader taxonomic analysis, including sequencing of multiple individuals per species, show the incredible plasticity and variation of ERVs within and among wildlife species. This demonstrates that stochastic processes govern much of the vertebrate genome. In this review, we elaborate on discoveries pertaining to wildlife ERV origins and evolution, genome colonization, and consequences for host biology.

De Novo Assembly of the Genome of the Sea Urchin Paracentrotus lividus (Lamarck 1816)

The Mediterranean purple sea urchin Paracentrotus lividus (Lamarck 1816) is a remarkable model system for molecular, evolutionary and cell biology studies, particularly in the field of developmental biology. We sequenced the genome, performed a de novo assembly, and analysed the assembly content. The genome of P. lividus was sequenced using Illumina NextSeq 500 System (Illumina) in a 2 × 150 paired-end format. More than 30,000 open reading frames (ORFs), (more than 8000 are unique), were identified and analysed to provide molecular tools accessible for the scientific community. In particular, several genes involved in complex innate immune responses, oxidative metabolism, signal transduction, and kinome, as well as genes regulating the membrane receptors, were identified in the P. lividus genome. In this way, the employment of the Mediterranean sea urchin for investigations and comparative analyses was empowered, leading to the explanation of cis-regulatory networks and their evolution in a key developmental model occupying an important evolutionary position with respect to vertebrates and humans.

Unraveling metagenomics through long-read sequencing: a comprehensive review

The study of microbial communities has undergone significant advancements, starting from the initial use of 16S rRNA sequencing to the adoption of shotgun metagenomics. However, a new era has emerged with the advent of long-read sequencing (LRS), which offers substantial improvements over its predecessor, short-read sequencing (SRS). LRS produces reads that are several kilobases long, enabling researchers to obtain more complete and contiguous genomic information, characterize structural variations, and study epigenetic modifications. The current leaders in LRS technologies are Pacific Biotechnologies (PacBio) and Oxford Nanopore Technologies (ONT), each offering a distinct set of advantages. This review covers the workflow of long-read metagenomics sequencing, including sample preparation (sample collection, sample extraction, and library preparation), sequencing, processing (quality control, assembly, and binning), and analysis (taxonomic annotation and functional annotation). Each section provides a concise outline of the key concept of the methodology, presenting the original concept as well as how it is challenged or modified in the context of LRS. Additionally, the section introduces a range of tools that are compatible with LRS and can be utilized to execute the LRS process. This review aims to present the workflow of metagenomics, highlight the transformative impact of LRS, and provide researchers with a selection of tools suitable for this task.

Recent insights into crosstalk between genetic parasites and their host genome

The bulk of higher order organismal genomes is comprised of transposable element (TE) copies, i.e. genetic parasites. The host-parasite relation is multi-faceted, varying across genomic region (genic versus intergenic), life-cycle stages, tissue-type and of course in health versus pathological state. The reach of functional genomics though, in investigating genotype-to-phenotype relations, has been limited when TEs are involved. The aim of this review is to highlight recent progress made in understanding how TE origin biochemical activity interacts with the central dogma stages of the host genome. Such interaction can also bring about modulation of the immune context and this could have important repercussions in disease state where immunity has a role to play. Thus, the review is to instigate ideas and action points around identifying evolutionary adaptations that the host genome and the genetic parasite have evolved and why they could be relevant.

ToxCodAn-Genome: an automated pipeline for toxin-gene annotation in genome assembly of venomous lineages

BACKGROUND

The rapid development of sequencing technologies resulted in a wide expansion of genomics studies using venomous lineages. This facilitated research focusing on understanding the evolution of adaptive traits and the search for novel compounds that can be applied in agriculture and medicine. However, the toxin annotation of genomes is a laborious and time-consuming task, and no consensus pipeline is currently available. No computational tool currently exists to address the challenges specific to toxin annotation and to ensure the reproducibility of the process.

RESULTS

Here, we present ToxCodAn-Genome, the first software designed to perform automated toxin annotation in genomes of venomous lineages. This pipeline was designed to retrieve the full-length coding sequences of toxins and to allow the detection of novel truncated paralogs and pseudogenes. We tested ToxCodAn-Genome using 12 genomes of venomous lineages and achieved high performance on recovering their current toxin annotations. This tool can be easily customized to allow improvements in the final toxin annotation set and can be expanded to virtually any venomous lineage. ToxCodAn-Genome is fast, allowing it to run on any personal computer, but it can also be executed in multicore mode, taking advantage of large high-performance servers. In addition, we provide a guide to direct future research in the venomics field to ensure a confident toxin annotation in the genome being studied. As a case study, we sequenced and annotated the toxin repertoire of Bothrops alternatus, which may facilitate future evolutionary and biomedical studies using vipers as models.

CONCLUSIONS

ToxCodAn-Genome is suitable to perform toxin annotation in the genome of venomous species and may help to improve the reproducibility of further studies. ToxCodAn-Genome and the guide are freely available at https://github.com/pedronachtigall/ToxCodAn-Genome.

How to make DNA data storage more applicable

The storage of digital data is becoming a worldwide problem. DNA has been recognized as a biological solution due to its ability to store genetic information without alteration over long periods. The first demonstrations of high-capacity long-lasting DNA digital data storage have been shown. However, high storage costs and slow retrieval of the data must be overcome to make DNA data storage more applicable and marketable. Herein, we discuss the issues and recent advances in DNA data storage methods and highlight pathways to make this technology more applicable to real-world digital data storage. We envision that a combination of molecular biology, nanotechnology, novel polymers, electronics, and automation with systematic development will allow DNA data storage sufficient for everyday use.

Whole-genome sequencing provides novel insights into the evolutionary history and genetic adaptation of reindeer populations in northern Eurasia

Domestic reindeer (Rangifer tarandus) play a vital role in the culture and livelihoods of indigenous people across northern Eurasia. These animals are well adapted to harsh environmental conditions, such as extreme cold, limited feed availability and long migration distances. Therefore, understanding the genomics of reindeer is crucial for improving their management, conservation and utilisation. In this study, we have generated a new genome assembly for the Fennoscandian domestic reindeer with high contiguity, making it the most complete reference genome for reindeer to date. The new genome assembly was utilised to explore genetic diversity, population structure and selective sweeps in Eurasian Rangifer tarandus populations which was based on the largest population genomic dataset for reindeer, encompassing 58 individuals from diverse populations. Phylogenetic analyses revealed distinct genetic clusters, with the Finnish wild forest reindeer (Rangifer tarandus fennicus) standing out as a unique subspecies. Divergence time estimates suggested a separation of ~ 52 thousand years ago (Kya) between the northern European Rangifer tarandus fennicus and Rangifer tarandus tarandus. Our study identified four main genetic clusters: Fennoscandian, the eastern/northern Russian and Alaskan group, the Finnish forest reindeer, and the Svalbard reindeer. Furthermore, two independent reindeer domestication processes were inferred, suggesting separate origins for the domestic Fennoscandian and eastern/northern Russian reindeer. Notably, shared genes under selection, including retroviral genes, point towards molecular domestication processes that aided adaptation of this species to diverse environments.

Advancing pathogen genomics in resource-limited settings

Genomic sequencing has emerged as a powerful tool to enhance early pathogen detection and characterization with implications for public health and clinical decision making. Although widely available in developed countries, the application of pathogen genomics among low-resource, high-disease burden settings remains at an early stage. In these contexts, tailored approaches for integrating pathogen genomics within infectious disease control programs will be essential to optimize cost efficiency and public health impact. We propose a framework for embedding pathogen genomics within national surveillance plans across a spectrum of surveillance and laboratory capacities. We adopt a public health approach to genomics and examine its application to high-priority diseases relevant in resource-limited settings. For each grouping, we assess the value proposition for genomics to inform public health and clinical decision-making, alongside its contribution toward research and development of novel diagnostics, therapeutics, and vaccines.

Current technical advancements in plant epitranscriptomic studies

The growth and development of plants are the result of the interplay between the internal developmental programming and plant-environment interactions. Gene expression regulations in plants are made up of multi-level networks. In the past few years, many studies were carried out on co- and post-transcriptional RNA modifications, which, together with the RNA community, are collectively known as the "epitranscriptome." The epitranscriptomic machineries were identified and their functional impacts characterized in a broad range of physiological processes in diverse plant species. There is mounting evidence to suggest that the epitranscriptome provides an additional layer in the gene regulatory network for plant development and stress responses. In the present review, we summarized the epitranscriptomic modifications found so far in plants, including chemical modifications, RNA editing, and transcript isoforms. The various approaches to RNA modification detection were described, with special emphasis on the recent development and application potential of third-generation sequencing. The roles of epitranscriptomic changes in gene regulation during plant-environment interactions were discussed in case studies. This review aims to highlight the importance of epitranscriptomics in the study of gene regulatory networks in plants and to encourage multi-omics investigations using the recent technical advancements.

Harnessing synthetic biology for advancing RNA therapeutics and vaccine design

Recent global events have drawn into focus the diversity of options for combatting disease across a spectrum of prophylactic and therapeutic approaches. The recent success of the mRNA-based COVID-19 vaccines has paved the way for RNA-based treatments to revolutionize the pharmaceutical industry. However, historical treatment options are continuously updated and reimagined in the context of novel technical developments, such as those facilitated through the application of synthetic biology. When it comes to the development of genetic forms of therapies and vaccines, synthetic biology offers diverse tools and approaches to influence the content, dosage, and breadth of treatment with the prospect of economic advantage provided in time and cost benefits. This can be achieved by utilizing the broad tools within this discipline to enhance the functionality and efficacy of pharmaceutical agent sequences. This review will describe how synthetic biology principles can augment RNA-based treatments through optimizing not only the vaccine antigen, therapeutic construct, therapeutic activity, and delivery vector. The enhancement of RNA vaccine technology through implementing synthetic biology has the potential to shape the next generation of vaccines and therapeutics.

Unveiling DNA Translocation in Pristine Graphene Nanopores: Understanding Pore Clogging via Polarizable Simulations

Graphene has garnered remarkable attention in recent years as an attractive nanopore membrane for rapid and accurate sequencing of DNA. The inherent characteristics of graphene offer exquisite experimental control over pore dimensions, encompassing both the width (pore diameter) and height. Despite these promising prospects, the practical deployment of pristine graphene nanopores for DNA sequencing has encountered a formidable challenge in the form of pore clogging, which is primarily attributed to hydrophobic interactions. However, a comprehensive understanding of the atomistic origins underpinning this clogging phenomenon and the nuanced impact of individual nucleobase identities on clogging dynamics remain an underexplored domain. Elucidating the atomistic intricacies governing pore clogging is pivotal to devising strategies for its mitigation and advancing our understanding of graphene nanopore behavior. We harness Drude polarizable simulations to systematically dissect the nucleobase-dependent mechanisms that play a pivotal role in nanopore clogging. We unveil nucleobase-specific interactions that illuminate the multifaceted roles played by both hydrophobic and electrostatic forces in driving nanopore clogging events. Notably, the Drude simulations also unveil the bias-dependent translocation dynamics and its pivotal role in alleviating pore clogging─a facet that remains significantly underestimated in conventional additive (nonpolarizable) simulations. Our findings underscore the indispensability of incorporating polarizability to faithfully capture the intricate dynamics governing graphene nanopore translocation phenomena, thus deepening our insights into this crucial field.

Genome Assembly of Cordia subcordata, a Coastal Protection Species in Tropical Coral Islands

Cordia subcordata trees or shrubs, belonging to the Boraginaceae family, have strong resistance and have adapted to their habitat on a tropical coral island in China, but the lack of genome information regarding its genetic background is unclear. In this study, the genome was assembled using both short/long whole genome sequencing reads and Hi-C reads. The assembled genome was 475.3 Mb, with 468.7 Mb (99.22%) of the sequences assembled into 16 chromosomes. Repeat sequences accounted for 54.41% of the assembled genome. A total of 26,615 genes were predicted, and 25,730 genes were functionally annotated using different annotation databases. Based on its genome and the other 17 species, phylogenetic analysis using 336 single-copy genes obtained from ortholog analysis showed that C. subcordata was a sister to Coffea eugenioides, and the divergence time was estimated to be 77 MYA between the two species. Gene family evolution analysis indicated that the significantly expanded gene families were functionally related to chemical defenses against diseases. These results can provide a reference to a deeper understanding of the genetic background of C. subcordata and can be helpful in exploring its adaptation mechanism on tropical coral islands in the future.

A chromosome-level genome assembly of the Korean crossbred pig Nanchukmacdon (Sus scrofa)

As plentiful high-quality genome assemblies have been accumulated, reference-guided genome assembly can be a good approach to reconstruct a high-quality assembly. Here, we present a chromosome-level genome assembly of the Korean crossbred pig called Nanchukmacdon (the NCMD assembly) using the reference-guided assembly approach with short and long reads. The NCMD assembly contains 20 chromosome-level scaffolds with a total size of 2.38 Gbp (N50: 138.77 Mbp). Its BUSCO score is 93.1%, which is comparable to the pig reference assembly, and a total of 20,588 protein-coding genes, 8,651 non-coding genes, and 996.14 Mbp of repetitive elements are annotated. The NCMD assembly was also used to close many gaps in the pig reference assembly. This NCMD assembly and annotation provide foundational resources for the genomic analyses of pig and related species.

Construction of heat stress regulation networks based on Illumina and SMRT sequencing data in potato

Potato (Solanum tuberosum L.) is one of the most important tuber food crops in the world; however, the cultivated potatoes are susceptible to high temperature, by which potato production is adversely affected. Understanding the coping mechanism of potato to heat stress is essential to secure yield and expand adaptability under environmental conditions with rising temperature. However, the lack of heat-related information has significantly limited the identification and application of core genes. To gain deeper insights into heat tolerance genes, next-generation sequencing and single-molecule real-time sequencing were used to learn the transcriptional response of potato to heat stress and 13,159 differentially expressed genes (DEGs) were identified in this study. All DEGs were grouped into 12 clusters using the K-means clustering algorithm. Gene Ontology enrichment analysis revealed that they were involved in temperature signaling, phytohormone, and protein modification. Among them, there were 950 differentially expressed transcription factors (DETFs). According to the network analysis of DETFs at the sixth hour under heat stress, we found some genes that were previously reported to be associated with photoperiodic tuberization, StCO (CONSTANS), tuber formation, StBEL11 (BEL1-LIKE 11), and earliness in potato, StCDF1 (CYCLING DOF FACTOR 1) responding to temperature. Furthermore, we verified the relative expression levels using quantitative real-time polymerase chain reaction, and the results were consistent with the inferences from transcriptomes. In addition, there were 22,125 alternative splicing events and 2,048 long non-coding RNAs. The database and network established in this study will extend our understanding of potato response to heat stress. It ultimately provided valuable resources for molecular analysis of heat stress response in potato and cultivation of potato varieties with heat tolerance.

The haplotype-resolved T2T genome of teinturier cultivar Yan73 reveals the genetic basis of anthocyanin biosynthesis in grapes

Teinturier grapes are characterized by the typical accumulation of anthocyanins in grape skin, flesh, and vegetative tissues, endowing them with high utility value in red wine blending and nutrient-enriched foods developing. However, due to the lack of genome information, the mechanism involved in regulating teinturier grape coloring has not yet been elucidated and their genetic utilization research is still insufficient. Here, the cultivar 'Yan73' was used for assembling the telomere-to-telomere (T2T) genome of teinturier grapes by combining the High Fidelity (HiFi)， Hi-C and ultralong Oxford Nanopore Technologies (ONT) reads. Two haplotype genomes were assembled, at the sizes of 501.68 Mb and 493.38 Mb, respectively. In the haplotype 1 genome, the transposable elements (TEs) contained 32.77% of long terminal repeats (LTRs), while in the haplotype 2 genome, 31.53% of LTRs were detected in TEs. Furthermore, obvious inversions were identified in chromosome 18 between the two haplotypes. Transcriptome profiling suggested that the gene expression patterns in 'Cabernet Sauvignon' and 'Yan73' were diverse depending on tissues, developmental stages, and varieties. The transcription program of genes in the anthocyanins biosynthesis pathway between the two cultivars exhibited high similarity in different tissues and developmental stages, whereas the expression levels of numerous genes showed significant differences. Compared with other genes, the expression levels of VvMYBA1 and VvUFGT4 in all samples, VvCHS2 except in young shoots and VvPAL9 except in the E-L23 stage of 'Yan73' were higher than those of 'Cabernet Sauvignon'. Further sequence alignments revealed potential variant gene loci and structure variations of anthocyanins biosynthesis related genes and a 816 bp sequence insertion was found in the promoter of VvMYBA1 of 'Yan73' haplotype 2 genome. The 'Yan73' T2T genome assembly and comparative analysis provided valuable foundations for further revealing the coloring mechanism of teinturier grapes and the genetic improvement of grape coloring traits.

Microfluidic delivery of cutting enzymes for fragmentation of surface-adsorbed DNA molecules

We describe a method for fragmenting, in-situ, surface-adsorbed and immobilized DNAs on polymethylmethacrylate(PMMA)-coated silicon substrates using microfluidic delivery of the cutting enzyme DNase I. Soft lithography is used to produce silicone elastomer (Sylgard 184) gratings which form microfluidic channels for delivery of the enzyme. Bovine serum albumin (BSA) is used to reduce DNase I adsorption to the walls of the microchannels and enable diffusion of the cutting enzyme to a distance of 10mm. Due to the DNAs being immobilized, the fragment order is maintained on the surface. Possible methods of preserving the order for application to sequencing are discussed.

A practical approach to genome assembly and annotation of Basidiomycota using the example of Armillaria

Technological advancements in genome sequencing, assembly and annotation platforms and algorithms that resulted in several genomic studies have created an opportunity to further our understanding of the biology of phytopathogens, including Armillaria species. Most Armillaria species are facultative necrotrophs that cause root- and stem-rot, usually on woody plants, significantly impacting agriculture and forestry worldwide. Genome sequencing, assembly and annotation in terms of samples used and methods applied in Armillaria genome projects are evaluated in this review. Infographic guidelines and a database of resources to facilitate future Armillaria genome projects were developed. Knowledge gained from genomic studies of Armillaria species is summarized and prospects for further research are provided. This guide can be applied to other diploid and dikaryotic fungal genomics.

Complete assembly, annotation of virulence genes and CRISPR editing of the genome of Leishmania amazonensis PH8 strain

We report the sequencing and assembly of the PH8 strain of Leishmania amazonensis one of the etiological agents of leishmaniasis. After combining data from long Pacbio reads, short Illumina reads and synteny with the Leishmania mexicana genome, the sequence of 34 chromosomes with 8317 annotated genes was generated. Multigene families encoding three virulence factors, A2, amastins and the GP63 metalloproteases, were identified and compared to their annotation in other Leishmania species. As they have been recently recognized as virulence factors essential for disease establishment and progression of the infection, we also identified 14 genes encoding proteins involved in parasite iron and heme metabolism and compared to genes from other Trypanosomatids. To follow these studies with a genetic approach to address the role of virulence factors, we tested two CRISPR-Cas9 protocols to generate L. amazonensis knockout cell lines, using the Miltefosine transporter gene as a proof of concept.

OGRE: calculate, visualize, and analyze overlap between genomic input regions and public annotations

BACKGROUND

Modern genome sequencing leads to an ever-growing collection of genomic annotations. Combining these elements with a set of input regions (e.g. genes) would yield new insights in genomic associations, such as those involved in gene regulation. The required data are scattered across different databases making a manual approach tiresome, unpractical, and prone to error. Semi-automatic approaches require programming skills in data parsing, processing, overlap calculation, and visualization, which most biomedical researchers lack. Our aim was to develop an automated tool providing all necessary algorithms, benefiting both bioinformaticians and researchers without bioinformatic training.

RESULTS

We developed overlapping annotated genomic regions (OGRE) as a comprehensive tool to associate and visualize input regions with genomic annotations. It does so by parsing regions of interest, mining publicly available annotations, and calculating possible overlaps between them. The user can thus identify location, type, and number of associated regulatory elements. Results are presented as easy to understand visualizations and result tables. We applied OGRE to recent studies and could show high reproducibility and potential new insights. To demonstrate OGRE's performance in terms of running time and output, we have conducted a benchmark and compared its features with similar tools.

CONCLUSIONS

OGRE's functions and built-in annotations can be applied as a downstream overlap association step, which is compatible with most genomic sequencing outputs, and can thus enrich pre-existing analyses pipelines. Compared to similar tools, OGRE shows competitive performance, offers additional features, and has been successfully applied to two recent studies. Overall, OGRE addresses the lack of tools for automatic analysis, local genomic overlap calculation, and visualization by providing an easy to use, end-to-end solution for both biologists and computational scientists.

The State of Squamate Genomics: Past, Present, and Future of Genome Research in the Most Speciose Terrestrial Vertebrate Order

Squamates include more than 11,000 extant species of lizards, snakes, and amphisbaenians, and display a dazzling diversity of phenotypes across their over 200-million-year evolutionary history on Earth. Here, we introduce and define squamates (Order Squamata) and review the history and promise of genomic investigations into the patterns and processes governing squamate evolution, given recent technological advances in DNA sequencing, genome assembly, and evolutionary analysis. We survey the most recently available whole genome assemblies for squamates, including the taxonomic distribution of available squamate genomes, and assess their quality metrics and usefulness for research. We then focus on disagreements in squamate phylogenetic inference, how methods of high-throughput phylogenomics affect these inferences, and demonstrate the promise of whole genomes to settle or sustain persistent phylogenetic arguments for squamates. We review the role transposable elements play in vertebrate evolution, methods of transposable element annotation and analysis, and further demonstrate that through the understanding of the diversity, abundance, and activity of transposable elements in squamate genomes, squamates can be an ideal model for the evolution of genome size and structure in vertebrates. We discuss how squamate genomes can contribute to other areas of biological research such as venom systems, studies of phenotypic evolution, and sex determination. Because they represent more than 30% of the living species of amniote, squamates deserve a genome consortium on par with recent efforts for other amniotes (i.e., mammals and birds) that aim to sequence most of the extant families in a clade.

Continuous synthesis of E. coli genome sections and Mb-scale human DNA assembly

Whole-genome synthesis provides a powerful approach for understanding and expanding organism function1-3. To build large genomes rapidly, scalably and in parallel, we need (1) methods for assembling megabases of DNA from shorter precursors and (2) strategies for rapidly and scalably replacing the genomic DNA of organisms with synthetic DNA. Here we develop bacterial artificial chromosome (BAC) stepwise insertion synthesis (BASIS)-a method for megabase-scale assembly of DNA in Escherichia coli episomes. We used BASIS to assemble 1.1 Mb of human DNA containing numerous exons, introns, repetitive sequences, G-quadruplexes, and long and short interspersed nuclear elements (LINEs and SINEs). BASIS provides a powerful platform for building synthetic genomes for diverse organisms. We also developed continuous genome synthesis (CGS)-a method for continuously replacing sequential 100 kb stretches of the E. coli genome with synthetic DNA; CGS minimizes crossovers1,4 between the synthetic DNA and the genome such that the output for each 100 kb replacement provides, without sequencing, the input for the next 100 kb replacement. Using CGS, we synthesized a 0.5 Mb section of the E. coli genome-a key intermediate in its total synthesis1-from five episomes in 10  days. By parallelizing CGS and combining it with rapid oligonucleotide synthesis and episome assembly5,6, along with rapid methods for compiling a single genome from strains bearing distinct synthetic genome sections1,7,8, we anticipate that it will be possible to synthesize entire E. coli genomes from functional designs in less than 2 months.

Scalable, accessible, and reproducible reference genome assembly and evaluation in Galaxy

Improvements in genome sequencing and assembly are enabling high-quality reference genomes for all species. However, the assembly process is still laborious, computationally and technically demanding, lacks standards for reproducibility, and is not readily scalable. Here we present the latest Vertebrate Genomes Project assembly pipeline and demonstrate that it delivers high-quality reference genomes at scale across a set of vertebrate species arising over the last ~500 million years. The pipeline is versatile and combines PacBio HiFi long-reads and Hi-C-based haplotype phasing in a new graph-based paradigm. Standardized quality control is performed automatically to troubleshoot assembly issues and assess biological complexities. We make the pipeline freely accessible through Galaxy, accommodating researchers even without local computational resources and enhanced reproducibility by democratizing the training and assembly process. We demonstrate the flexibility and reliability of the pipeline by assembling reference genomes for 51 vertebrate species from major taxonomic groups (fish, amphibians, reptiles, birds, and mammals).

The complete reference genome for grapevine (Vitis vinifera L.) genetics and breeding

Grapevine is one of the most economically important crops worldwide. However, the previous versions of the grapevine reference genome tipically consist of thousands of fragments with missing centromeres and telomeres, limiting the accessibility of the repetitive sequences, the centromeric and telomeric regions, and the study of inheritance of important agronomic traits in these regions. Here, we assembled a telomere-to-telomere (T2T) gap-free reference genome for the cultivar PN40024 using PacBio HiFi long reads. The T2T reference genome (PN_T2T) is 69 Mb longer with 9018 more genes identified than the 12X.v0 version. We annotated 67% repetitive sequences, 19 centromeres and 36 telomeres, and incorporated gene annotations of previous versions into the PN_T2T assembly. We detected a total of 377 gene clusters, which showed associations with complex traits, such as aroma and disease resistance. Even though PN40024 derives from nine generations of selfing, we still found nine genomic hotspots of heterozygous sites associated with biological processes, such as the oxidation-reduction process and protein phosphorylation. The fully annotated complete reference genome therefore constitutes an important resource for grapevine genetic studies and breeding programs.

Rapid survey of de novo mutations in naturally growing tree species following the March 2011 disaster in Fukushima: The effect of low-dose-rate radiation

The impact of low-dose-rate radiation on genetics is largely unknown, particularly in natural environments. The Fukushima Dai-ich Nuclear Power Plant disaster resulted in the creation of contaminated natural lands. In this study, de novo mutations (DNMs) in germ line cells were surveyed from double-digest RADseq fragments in Japanese cedar and flowering cherry trees exposed to ambient dose rates ranging from 0.08 to 6.86 μGy h^-1. These two species are among the most widely cultivated Japanese gymnosperm and angiosperm trees for forestry and horticultural purpose, respectively. For Japanese flowering cherry, open crossings were performed to produce seedlings, and only two candidate DNMs were detected from uncontaminated area. For Japanese cedar, the haploid megagametophytes were used as next generation samples. The use of megagametophytes from open crossing for next generation mutation screening had many advantages such as reducing exposure to radiation in contaminated areas because artificial crossings are not needed and the ease of data analysis owing to the haploid nature of megagametophytes. A direct comparison of the nucleotide sequences of parents and megagametophytes revealed an average of 1.4 candidate DNMs per megagametophyte sample (range: 0-40) after filtering procedures were optimized based on the validation of DNMs via Sanger sequencing. There was no relationship between the observed mutations and the ambient dose rate in the growing area or the concentration of ¹³⁷Cs in cedar branches. The present results also suggest that mutation rates differ among lineages and that the growing environment has a relatively large influence on these mutation rates. These results suggested there was no significant increase in the mutation rate of the germplasm of Japanese cedar and flowering cherry trees growing in the contaminated areas.

A practical comparison of the next-generation sequencing platform and assemblers using yeast genome

Assembling fragmented whole-genomic information from the sequencing data is an inevitable process for further genome-wide research. However, it is intricate to select the appropriate assembly pipeline for unknown species because of the species-specific genomic properties. Therefore, our study focused on relatively more static proclivities of sequencing platforms and assembly algorithms than the fickle genome sequences. A total of 212 draft and polished de novo assemblies were constructed under the different sequencing platforms and assembly algorithms with the repetitive yeast genome. Our comprehensive data indicated that sequencing reads from Oxford Nanopore with R7.3 flow cells generated more continuous assemblies than those derived from the PacBio Sequel, although the homopolymer-based assembly errors and chimeric contigs exist. In addition, the comparison between two second-generation sequencing platforms showed that Illumina NovaSeq 6000 provides more accurate and continuous assembly in the second-generation-sequencing-first pipeline, but MGI DNBSEQ-T7 provides a cheap and accurate read in the polishing process. Furthermore, our insight into the relationship among the computational time, read length, and coverage depth provided clues to the optimal pipelines of yeast assembly.

Identification of Genetic Risk Factors for Monogenic and Complex Canine Diseases

Advances in DNA sequencing and other technologies have greatly facilitated the identification of genetic risk factors for inherited diseases in dogs. We review recent technological developments based on selected examples from canine disease genetics. The identification of disease-causing variants in dogs with monogenic diseases may become a widely employed diagnostic approach in clinical veterinary medicine in the not-too-distant future. Diseases with complex modes of inheritance continue to pose challenges to researchers but have also become much more tangible than in the past. In addition to strategies for identifying genetic risk factors, we provide some thoughts on the interpretation of sequence variants that are largely inspired by developments in human clinical genetics.

Genetic and epigenetic control of the plant metabolome

Plant metabolites are mainly produced through chemical reactions catalysed by enzymes encoded in the genome. Mutations in enzyme-encoding or transcription factor-encoding genes can alter the metabolome by changing the enzyme's catalytic activity or abundance, respectively. Insertion of transposable elements into non-coding regions has also been reported to affect transcription and ultimately metabolite content. In addition to genetic mutations, transgenerational epigenetic variations have also been found to affect metabolic content by controlling the transcription of metabolism-related genes. However, the majority of cases reported so far, in which epigenetic mechanisms are associated with metabolism, are non-transgenerational, and are triggered by developmental signals or environmental stress. Although, accumulating research has provided evidence of strong genetic control of the metabolome, epigenetic control has been largely untouched. Here, we provide a review of the genetic and epigenetic control of metabolism with a focus on epigenetics. We discuss both transgenerational and non-transgenerational epigenetic marks regulating metabolism as well as prospects of the field of metabolic control where intricate interactions between genetics and epigenetics are involved.

The gap between development and manufacturing in gene therapy: Strategic options for overcoming traps

Gene therapy has been one of the most promising therapeutic approaches in recent years. This study analyzed a research and development (R&D) system for adeno-associated virus (AAV)-based gene therapies, and confirmed that there was a gap between the development and manufacturing capabilities. Although a start-up company that has no academic or manufacturing facilities can begin the clinical development process, it cannot successfully continue development activities without forming alliances and capital investment or, at a certain stage, without appropriate manufacturing and marketing strategies. We reviewed a series of case studies to categorize the acquisition patterns of pharmaceutical companies that are engaged in AAV gene therapy. These results provide insights into the R&D structures for AAV gene therapies from a technological management perspective.

In it for the long run: perspectives on exploiting long-read sequencing in livestock for population scale studies of structural variants

Studies have demonstrated that structural variants (SV) play a substantial role in the evolution of species and have an impact on Mendelian traits in the genome. However, unlike small variants (< 50 bp), it has been challenging to accurately identify and genotype SV at the population scale using short-read sequencing. Long-read sequencing technologies are becoming competitively priced and can address several of the disadvantages of short-read sequencing for the discovery and genotyping of SV. In livestock species, analysis of SV at the population scale still faces challenges due to the lack of resources, high costs, technological barriers, and computational limitations. In this review, we summarize recent progress in the characterization of SV in the major livestock species, the obstacles that still need to be overcome, as well as the future directions in this growing field. It seems timely that research communities pool resources to build global population-scale long-read sequencing consortiums for the major livestock species for which the application of genomic tools has become cost-effective.