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

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.

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.

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.

Unraveling Plant-Pathogen Interactions in Cereals Using RNA-seq

Over the past two decades, there have been significant advancements in the realm of transcriptomics, or the study of genes and their expression. Modern RNA sequencing technologies and high-performance computing are creating a "big data" revolution that provides new opportunities to explore the interactions between cereals and pathogens that affect grain yield and food safety. These data are being used to annotate genes and gene variants, as well as identify differentially expressed genes and create global gene co-expression networks. Moreover, these data can unravel the complex interactions between pathogen and host and identify genes and pathways involved in these interactions. This information can then be used for disease mitigation and the development of crops with superior resistance.

Chromosome-level genome and the identification of sex chromosomes in Uloborus diversus

The orb web is a remarkable example of animal architecture that is observed in families of spiders that diverged over 200 million years ago. While several genomes exist for araneid orb-weavers, none exist for other orb-weaving families, hampering efforts to investigate the genetic basis of this complex behavior. Here we present a chromosome-level genome assembly for the cribellate orb-weaving spider Uloborus diversus. The assembly reinforces evidence of an ancient arachnid genome duplication and identifies complete open reading frames for every class of spidroin gene, which encode the proteins that are the key structural components of spider silks. We identified the 2 X chromosomes for U. diversus and identify candidate sex-determining loci. This chromosome-level assembly will be a valuable resource for evolutionary research into the origins of orb-weaving, spidroin evolution, chromosomal rearrangement, and chromosomal sex determination in spiders.

Near-chromosomal de novo assembly of Bengal tiger genome reveals genetic hallmarks of apex predation

The tiger, a poster child for conservation, remains an endangered apex predator. Continued survival and recovery will require a comprehensive understanding of genetic diversity and the use of such information for population management. A high-quality tiger genome assembly will be an important tool for conservation genetics, especially for the Indian tiger, the most abundant subspecies in the wild. Here, we present high-quality near-chromosomal genome assemblies of a female and a male wild Indian tiger (Panthera tigris tigris). Our assemblies had a scaffold N50 of >140 Mb, with 19  scaffolds corresponding to the 19 numbered chromosomes, containing 95% of the genome. Our assemblies also enabled detection of longer stretches of runs of homozygosity compared to previous assemblies, which will help improve estimates of genomic inbreeding. Comprehensive genome annotation identified 26,068 protein-coding genes, including several gene families involved in key morphological features such as the teeth, claws, vision, olfaction, taste, and body stripes. We also identified 301 microRNAs, 365 small nucleolar RNAs, 632 transfer RNAs, and other noncoding RNA elements, several of which are predicted to regulate key biological pathways that likely contribute to the tiger's apex predatory traits. We identify signatures of positive selection in the tiger genome that are consistent with the Panthera lineage. Our high-quality genome will enable use of noninvasive samples for comprehensive assessment of genetic diversity, thus supporting effective conservation and management of wild tiger populations.

Advances in experimental and computational methodologies for the study of microbial-surface interactions at different omics levels

The study of the biological response of microbial cells interacting with natural and synthetic interfaces has acquired a new dimension with the development and constant progress of advanced omics technologies. New methods allow the isolation and analysis of nucleic acids, proteins and metabolites from complex samples, of interest in diverse research areas, such as materials sciences, biomedical sciences, forensic sciences, biotechnology and archeology, among others. The study of the bacterial recognition and response to surface contact or the diagnosis and evolution of ancient pathogens contained in archeological tissues require, in many cases, the availability of specialized methods and tools. The current review describes advances in in vitro and in silico approaches to tackle existing challenges (e.g., low-quality sample, low amount, presence of inhibitors, chelators, etc.) in the isolation of high-quality samples and in the analysis of microbial cells at genomic, transcriptomic, proteomic and metabolomic levels, when present in complex interfaces. From the experimental point of view, tailored manual and automatized methodologies, commercial and in-house developed protocols, are described. The computational level focuses on the discussion of novel tools and approaches designed to solve associated issues, such as sample contamination, low quality reads, low coverage, etc. Finally, approaches to obtain a systems level understanding of these complex interactions by integrating multi omics datasets are presented.

Hybrid assembly using long reads resolves repeats and completes the genome sequence of a laboratory strain of Staphylococcus aureus subsp. Aureus RN4220

Staphylococcus aureus RN4220 has been extensively used by staphylococcal researchers as an intermediate strain for genetic manipulation due to its ability to accept foreign DNA. Despite its wide use in laboratories, its complete genome is not available. In this study, we used a hybrid genome assembly approach using minION long reads and Illumina short reads to sequence the complete genome of S. aureus RN4220. The comparative analysis of the annotated complete genome showed the presence of 39 genes fragmented in the previous assembly, many of which were located near the repeat regions. Using RNA-Seq reads, we showed that a higher number of reads could be mapped to the complete genome than the draft genome and the gene expression profile obtained using the complete genome also differs from that obtained from the draft genome. Furthermore, by comparative transcriptomic analysis, we showed the correlation between expression levels of staphyloxanthin biosynthetic genes and the production of yellow pigment. This study highlighted the importance of long reads in completing microbial genomes, especially those possessing repetitive elements.

•

S. aureus RN4220 is used as an intermediate strain for genetic manipulation.

•

We completed its genome and found 39 fragmented genes in previous genome assembly.

•

RNA-Seq analysis improved mapping of the reads with the use of complete genome.

•

Expression of staphyloxanthin biosynthetic genes was correlated with its production.

Semi-automated assembly of high-quality diploid human reference genomes

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

Comparative genomics reveals putative evidence for high-elevation adaptation in the American pika (Ochotona princeps)

High-elevation environments have lower atmospheric oxygen content, reduced temperatures, and higher levels of UV radiation than found at lower elevations. As such, species living at high elevations must overcome these challenges to survive, grow, and reproduce. American pikas (Ochotona princeps) are alpine lagomorphs that are habitat specialists typically found at elevations >2,000 m. Previous research has shown putative evidence for high-elevation adaptation; however, investigations to date have been limited to a fraction of the genome. Here, we took a comparative genomics approach to identify putative regions under selection using a chromosomal reference genome assembly for the American pika relative to 8 other mammalian species targeted based on phylogenetic relatedness and (dis)similarity in ecology. We first identified orthologous gene groups across species and then extracted groups containing only American pika genes as well as unclustered pika genes to inform functional enrichment analyses; among these, we found 141 enriched terms with many related to hypoxia, metabolism, mitochondrial function/development, and DNA repair. We identified 15 significantly expanded gene families within the American pika across all orthologous gene groups that displayed functionally enriched terms associated with hypoxia adaptation. We further detected 196 positively selected genes, 41 of which have been associated with putative adaptation to hypoxia, cold tolerance, and response to UV following a literature review. In particular, OXNAD1, NRDC, and those genes critical in DNA repair represent important targets for future research to examine their functional implications in the American pika, especially as they may relate to adaptation to rapidly changing environments.

Identification of interaction partners using protein aggregation and NMR spectroscopy

The interaction among proteins is one of the most fundamental methods of information transfer in the living system. Many methods have been developed in order to identify the interaction pairs or groups either in vivo or in vitro. The in vitro pulldown/coprecipitation assay directly observes the protein that binds to the target. This method involves electrophoresis, which is a technique of a low resolution as well as a low throughput. As a better alternative, we wish to propose a new method that is based on the NMR spectroscopy. This method utilizes the aggregation of the target protein and the concomitant signal disappearance of the interacting partner. The aggregation is accomplished by the elastin-like polypeptide, which is fused to the target. If a protein binds to this supramolecular complex, its NMR signal then becomes too broadened in order to be observed, which is the basic phenomenon of the NMR spectroscopy. Thus, the protein that loses its signal is the one that binds to the target. A compound that interferes with these types of bindings among the proteins can be identified by observing the reappearance of the protein signals with the simultaneous disappearance of the signals of the compound. This technique will be applied in order to find an interaction pair in the information transfer pathway as well as a compound that disrupts it. This proposed method should be able to work with a mixture of proteins and provide a higher resolution in order to find the binding partner in a higher throughput fashion.