A field guide to whole-genome sequencing, assembly and annotation

Systems biology of Haemonchus contortus - Advancing biotechnology for parasitic nematode control

Parasitic nematodes represent a substantial global burden, impacting animal health, agriculture and economies worldwide. Of these worms, Haemonchus contortus - a blood-feeding nematode of ruminants - is a major pathogen and a model for molecular and applied parasitology research. This review synthesises some key advances in understanding the molecular biology, genetic diversity and host-parasite interactions of H. contortus, highlighting its value for comparative studies with the free-living nematode Caenorhabditis elegans. Key themes include recent developments in genomic, transcriptomic and proteomic technologies and resources, which are illuminating critical molecular pathways, including the ubiquitination pathway, protease/protease inhibitor systems and the secretome of H. contortus. Some of these insights are providing a foundation for identifying essential genes and exploring their potential as targets for novel anthelmintics or vaccines, particularly in the face of widespread anthelmintic resistance. Advanced bioinformatic tools, such as machine learning (ML) algorithms and artificial intelligence (AI)-driven protein structure prediction, are enhancing annotation capabilities, facilitating and accelerating analyses of gene functions, and biological pathways and processes. This review also discusses the integration of these tools with cutting-edge single-cell sequencing and spatial transcriptomics to dissect host-parasite interactions at the cellular level. The discussion emphasises the importance of curated databases, improved culture systems and functional genomics platforms to translate molecular discoveries into practical outcomes, such as novel interventions. New research findings and resources not only advance research on H. contortus and related nematodes but may also pave the way for innovative solutions to the global challenges with anthelmintic resistance.

Establishing genome sequencing and assembly for non-model and emerging model organisms: a brief guide

Reference genome assemblies are the basis for comprehensive genomic analyses and comparisons. Due to declining sequencing costs and growing computational power, genome projects are now feasible in smaller labs. De novo genome sequencing for non-model or emerging model organisms requires knowledge about genome size and techniques for extracting high molecular weight DNA. Next to quality, the amount of DNA obtained from single individuals is crucial, especially, when dealing with small organisms. While long-read sequencing technologies are the methods of choice for creating high quality genome assemblies, pure short-read assemblies might bear most of the coding parts of a genome but are usually much more fragmented and do not well resolve repeat elements or structural variants. Several genome initiatives produce more and more non-model organism genomes and provide rules for standards in genome sequencing and assembly. However, sometimes the organism of choice is not part of such an initiative or does not meet its standards. Therefore, if the scientific question can be answered with a genome of low contiguity in intergenic parts, missing the high standards of chromosome scale assembly should not prevent publication. This review describes how to set up an animal genome sequencing project in the lab, how to estimate costs and resources, and how to deal with suboptimal conditions. Thus, we aim to suggest optimal strategies for genome sequencing that fulfil the needs according to specific research questions, e.g. "How are species related to each other based on whole genomes?" (phylogenomics), "How do genomes of populations within a species differ?" (population genomics), "Are differences between populations relevant for conservation?" (conservation genomics), "Which selection pressure is acting on certain genes?" (identification of genes under selection), "Did repeats expand or contract recently?" (repeat dynamics).

SNP-RFLP Markers for the Study of Arabidopsis lyrata

Arabidopsis lyrata has become a useful system for the study of comparative genomics, hybridization, polyploidization, and evolutionary transitions from outcrossing to selfing. Previous studies of its mating system have used microsatellite loci, but low allelic diversity, particularly in self-compatible populations characterized by low levels of outcrossing, reduces the utility of these markers for more detailed studies. Here, we aimed to develop population-level SNP markers for A. lyrata ssp. lyrata sampled from a self-compatible population at Rondeau Provincial Park, Ontario, Canada. We performed de novo SNP discovery and identified 6808 putative SNPs from genome-wide sequences of 22 individuals originating from a highly selfing population. Further filtering and marker validation enabled the development of 17 SNP marker loci that can be visualized using standard PCR-RFLP protocols. These markers had average minor-allele frequencies of 0.40 in the target population, and four of seven markers were variable in a small sample from nine other A. lyrata populations. These PCR-RFLP markers have the potential to be useful for the analysis of mating patterns within and beyond the inbred self-compatible populations of A. lyrata studied here and enable the continued development of A. lyrata as a model for studying evolutionary transitions from outcrossing to selfing.

Insights into the Genomic Background of Nine Common Chinese Medicinal Plants by Flow Cytometry and Genome Survey

Medicinal plants have long played a crucial role in healthcare systems, but limited genomic information on these species has impeded the integration of modern biological technologies into medicinal plant research. In this study, we selected nine common medicinal plants, each belonging to a different plant family, including Sarcandra glabra (Chloranthaceae), Nekemias grossedentata (Vitaceae), Uraria crinita (Fabaceae), Gynostemma pentaphyllum (Cucurbitaceae), Reynoutria japonica (Polygonaceae), Pseudostellaria heterophylla (Caryophyllaceae), Morinda officinalis (Rubiaceae), Vitex rotundifolia (Lamiaceae), and Gynura formosana (Asteraceae), to estimate their genome sizes and conduct preliminary genomic surveys. The estimated genome sizes by flow cytometry were 3.66 Gb, 0.65 Gb, 0.58 Gb, 1.02 Gb, 3.96 Gb, 2.99 Gb, 0.43 Gb, 0.78 Gb, and 7.27 Gb, respectively. The genome sizes of M. officinalis, R. japonica, and G. pentaphyllum have been previously reported. Comparative analyses suggest that variations in genome size may arise due to differences in measurement methods and sample sources. Therefore, employing multiple approaches to assess genome size is necessary to provide more reliable information for further genomic research. Based on the genome survey, species with considerable genome size variation or polyploidy, such as G. pentaphyllum, should undergo a ploidy analysis in conjunction with population genomics studies to elucidate the development of the diversified genome size. Additionally, a genome survey of U. crinita, a medicinal plant with a relatively small genome size (509.08 Mb) and of considerable interest in southern China, revealed a low heterozygosity rate (0.382%) and moderate repeat content (51.24%). Given the limited research costs, this species represents a suitable candidate for further genomic studies on Leguminous medicinal plants characteristic of southern China. This foundational genomic information will serve as a critical reference for the sustainable development and utilization of these medicinal plants.

Identification of Laccase Family of Auricularia auricula-judae and Structural Prediction Using Alphafold

Laccase is an enzyme that plays an important role in fungi, including lignin degradation, stress defense, and formation of fruiting bodies. Auricularia auricula-judae is a white-rot fungus in the Basidiomycota phylum, capable of delignifying wood. In this study, seven genes belonging to the laccase family were identified through de novo sequencing, containing Cu-Oxidase, Cu-Oxidase_2, and Cu-Oxidase_3 domains. Subsequently, the physical characteristics, phylogenetic relationships, protein secondary structure, and tertiary structure of the laccase family (AaLac1-AaLac7) were analyzed. Prediction of N-glycosylation sites identified 2 to 10 sites in the laccase family, with AaLac7 having the highest number of sites at 10. Sequence alignment and analysis of the laccase family showed high consistency in signature sequences. Phylogenetic analysis confirmed the relationship among laccases within the family, with AaLac3-AaLac4 and AaLac5-AaLac6 being closely positioned on the tree, exhibiting high similarity in tertiary structure predictions. This study identified and analyzed laccase family genes in Auricularia auricula-judae using de novo sequencing, offering a simple method for identifying and analyzing the laccase family in organisms with unknown genetic information.

Museomics help resolving the phylogeny of snowfinches (Aves, Passeridae, Montifringilla and allies)

Historical specimens from museum collections provide a valuable source of material also from remote areas or regions of conflict that are not easily accessible to scientists today. With this study, we are providing a taxon-complete phylogeny of snowfinches using historical DNA from whole skins of an endemic species from Afghanistan, the Afghan snowfinch, Pyrgilauda theresae. To resolve the strong conflict between previous phylogenetic hypotheses, we generated novel mitogenome sequences for selected taxa and genome-wide SNP data using ddRAD sequencing for all extant snowfinch species endemic to the Qinghai-Tibet Plateau (QTP) and for an extended intraspecific sampling of the sole Central and Western Palearctic snowfinch species (Montifringilla nivalis). Our phylogenetic reconstructions unanimously refuted the previously suggested paraphyly of genus Pyrgilauda. Misplacement of one species-level taxon (Onychostruthus tazcanowskii) in previous snowfinch phylogenies was undoubtedly inferred from chimeric mitogenomes that included heterospecific sequence information. Furthermore, comparison of novel and previously generated sequence data showed that the presumed sister-group relationship between M. nivalis and the QTP endemic M. henrici was suggested based on flawed taxonomy. Our phylogenetic reconstructions based on genome-wide SNP data and on mitogenomes were largely congruent and supported reciprocal monophyly of genera Montifringilla and Pyrgilauda with monotypic Onychostruthus being sister to the latter. The Afghan endemic P. theresae likely originated from a rather ancient Pliocene out-of-Tibet dispersal probably from a common ancestor with P. ruficollis. Our extended trans-Palearctic sampling for the white-winged snowfinch, M. nivalis, confirmed strong lineage divergence between an Asian and a European clade dated to 1.5 - 2.7 million years ago (mya). Genome-wide SNP data suggested subtle divergence among European samples from the Alps and from the Cantabrian mountains.

Unraveling genomic features and phylogenomics through the analysis of three Mexican endemic Myotis genomes

Background

Genomic resource development for non-model organisms is rapidly progressing, seeking to uncover molecular mechanisms and evolutionary adaptations enabling thriving in diverse environments. Limited genomic data for bat species hinder insights into their evolutionary processes, particularly within the diverse Myotis genus of the Vespertilionidae family. In Mexico, 15 Myotis species exist, with three-M. vivesi, M. findleyi, and M. planiceps-being endemic and of conservation concern.

Methods

We obtained samples of Myotis vivesi, M. findleyi, and M. planiceps for genomic analysis. Each of three genomic DNA was extracted, sequenced, and assembled. The scaffolding was carried out utilizing the M. yumanensis genome via a genome-referenced approach within the ntJoin program. GapCloser was employed to fill gaps. Repeat elements were characterized, and gene prediction was done via ab initio and homology methods with MAKER pipeline. Functional annotation involved InterproScan, BLASTp, and KEGG. Non-coding RNAs were annotated with INFERNAL, and tRNAscan-SE. Orthologous genes were clustered using Orthofinder, and a phylogenomic tree was reconstructed using IQ-TREE.

Results

We present genome assemblies of these endemic species using Illumina NovaSeq 6000, each exceeding 2.0 Gb, with over 90% representing single-copy genes according to BUSCO analyses. Transposable elements, including LINEs and SINEs, constitute over 30% of each genome. Helitrons, consistent with Vespertilionids, were identified. Values around 20,000 genes from each of the three assemblies were derived from gene annotation and their correlation with specific functions. Comparative analysis of orthologs among eight Myotis species revealed 20,820 groups, with 4,789 being single copy orthogroups. Non-coding RNA elements were annotated. Phylogenomic tree analysis supported evolutionary chiropterans' relationships. These resources contribute significantly to understanding gene evolution, diversification patterns, and aiding conservation efforts for these endangered bat species.

In vitro, in planta, and comparative genomic analyses of Pseudomonas syringae pv. syringae strains of pepper (Capsicum annuum var. annuum)

Pseudomonas syringae pv. syringae (Pss) is an emerging phytopathogen that causes Pseudomonas leaf spot (PLS) disease in pepper plants. Pss can cause serious economic damage to pepper production, yet very little is known about the virulence factors carried by Pss that cause disease in pepper seedlings. In this study, Pss strains isolated from pepper plants showing PLS symptoms in Ohio between 2013 and 2021 (n = 16) showed varying degrees of virulence (Pss populations and disease symptoms on leaves) on 6-week-old pepper seedlings. In vitro studies assessing growth in nutrient-limited conditions, biofilm production, and motility also showed varying degrees of virulence, but in vitro and in planta variation in virulence between Pss strains did not correlate. Comparative whole-genome sequencing studies identified notable virulence genes including 30 biofilm genes, 87 motility genes, and 106 secretion system genes. Additionally, a total of 27 antimicrobial resistance genes were found. A multivariate correlation analysis and Scoary analysis based on variation in gene content (n = 812 variable genes) and single nucleotide polymorphisms within virulence genes identified no significant correlations with disease severity, likely due to our limited sample size. In summary, our study explored the virulence and antimicrobial gene content of Pss in pepper seedlings as a first step toward understanding the virulence and pathogenicity of Pss in pepper seedlings. Further studies with additional pepper Pss strains will facilitate defining genes in Pss that correlate with its virulence in pepper seedlings, which can facilitate the development of effective measures to control Pss in pepper and other related P. syringae pathovars.

IMPORTANCE

Pseudomonas leaf spot (PLS) caused by Pseudomonas syringae pv. syringae (Pss) causes significant losses to the pepper industry. Highly virulent Pss strains under optimal environmental conditions (cool-moderate temperatures, high moisture) can cause severe necrotic lesions on pepper leaves that consequently can decrease pepper yield if the disease persists. Hence, it is important to understand the virulence mechanisms of Pss to be able to effectively control PLS in peppers. In our study, in vitro, in planta, and whole-genome sequence analyses were conducted to better understand the virulence and pathogenicity characteristics of Pss strains in peppers. Our findings fill a knowledge gap regarding potential virulence and pathogenicity characteristics of Pss in peppers, including virulence and antimicrobial gene content. Our study helps pave a path to further identify the role of specific virulence genes in causing disease in peppers, which can have implications in developing strategies to effectively control PLS in peppers.

European Food Safety Authority (EFSA), César‐Razquin A, Casacuberta J, Dalmay T, Federici S, Jacchia S, Kagkli DM, Moxon S, Papadopoulou N. Technical Note on the quality of DNA sequencing for the molecular characterisation of genetically modified plants

As part of the risk assessment (RA) requirements for genetically modified (GM) plants, according to Regulation (EU) No 503/2013 and the EFSA guidance on the RA of food and feed from GM plants (EFSA GMO Panel 2011), applicants need to perform a molecular characterisation of the DNA sequences inserted in the GM plant genome. This Technical Note to the applicants puts together requirements and recommendations for the quality assessment of the methodology, analysis and reporting when DNA sequencing is used for the molecular characterisation of GM plants. In particular, it applies to the use of Sanger sequencing and next-generation sequencing for the characterisation of the inserted genetic material and its flanking regions at each insertion site, the determination of the copy number of all detectable inserts and the analysis of the genetic stability of the inserts. This updated document replaces the EFSA 2018 Technical Note and reflects the current knowledge in scientific-technical methods for generating and verifying, in a standardised manner, DNA sequencing data in the context of RA of GM plants. It does not take into consideration the verification and validation of the detection method which remains under the remit of the Joint Research Centre (JRC).

Molecular characterization of Campylobacter spp. isolates obtained from commercial broilers and native chickens in Southern Thailand using whole genome sequencing

Chickens are the primary reservoirs of Campylobacter spp., mainly C. jejuni and C. coli, that cause human bacterial gastrointestinal infections. However, genomic characteristics and antimicrobial resistance of Campylobacter spp. in low- to middle-income countries need more comprehensive exploration. This study aimed to characterize 21 C. jejuni and 5 C. coli isolates from commercial broilers and native chickens using whole genome sequencing and compare them to 28 reference Campylobacter sequences. Among the 26 isolates, 13 sequence types (ST) were identified in C. jejuni and 5 ST in C. coli. The prominent ST was ST 2274 (5 isolates, 19.2%), followed by ST 51, 460, 2409, and 6455 (2 isolates in each ST, 7.7%), while all remaining ST (464, 536, 595, 2083, 6736, 6964, 8096, 10437, 828, 872, 900, 8237, and 13540) had 1 isolate per ST (3.8%). Six types of antimicrobial resistance genes (ant(6)-Ia, aph(3')-III, blaOXA, cat, erm(B), and tet(O)) and one point mutations in the gyrA gene (Threonine-86-Isoleucine) and another in the rpsL gene (Lysine-43-Arginine) were detected. The blaOXA resistance gene was present in all isolates, the gyrA mutations was in 95.2% of C. jejuni and 80.0% of C. coli, and the tet(O) resistance gene in 76.2% of C. jejuni and 80.0% of C. coli. Additionally, 203 virulence-associated genes linked to 16 virulence factors were identified. In terms of phenotypic resistance, the C. jejuni isolates were all resistant to ciprofloxacin, enrofloxacin, and nalidixic acid, with lower levels of resistance to tetracycline (76.2%), tylosin (52.3%), erythromycin (23.8%), azithromycin (22.2%), and gentamicin (11.1%). Most C. coli isolates were resistant to all tested antimicrobials, while 1 C. coli was pan-susceptible except for tylosin. Single-nucleotide polymorphisms concordance varied widely, with differences of up to 13,375 single-nucleotide polymorphisms compared to the reference Campylobacter isolates, highlighting genetic divergence among comparative genomes. This study contributes to a deeper understanding of the molecular epidemiology of Campylobacter spp. in Thai chicken production systems.

The updated genome of the Hungarian population of Aedes koreicus

Vector-borne diseases pose a potential risk to human and animal welfare, and understanding their spread requires genomic resources. The mosquito Aedes koreicus is an emerging vector that has been introduced into Europe more than 15 years ago but only a low quality, fragmented genome was available. In this study, we carried out additional sequencing and assembled and characterized the genome of the species to provide a background for understanding its evolution and biology. The updated genome was 1.1 Gbp long and consisted of 6099 contigs with an N50 value of 329,610 bp and a BUSCO score of 84%. We identified 22,580 genes that could be functionally annotated and paid particular attention to the identification of potential insecticide resistance genes. The assessment of the orthology of the genes indicates a high turnover at the terminal branches of the species tree of mosquitoes with complete genomes, which could contribute to the adaptation and evolutionary success of the species. These results could form the basis for numerous downstream analyzes to develop targets for the control of mosquito populations.

The draft genome of Spiraea crenata L. (Rosaceae) - the first complete genome in tribe Spiraeeae

Spiraea crenata L. is a deciduous shrub distributed across the Eurasian steppe zone. The species is of cultural and horticultural importance and occurs in scattered populations throughout its westernmost range. Currently, there is no genomic information on the tribe of Spiraeeae. Therefore we sequenced and assembled the whole genome of S. crenata using second- and third-generation sequencing and a hybrid assembly approach to expand genomic resources for conservation and support research on this horticulturally important lineage. In addition to the organellar genomes (the plastome and the mitochondrion), we present the first draft genome of the species with an estimated size of 220 Mbp, an N50 value of 7.7 Mbp, and a BUSCO score of 96.0%. Being the first complete genome in tribe Spiraeeae, this may not only be the first step in the genomic study of a rare plant but also a contribution to genomic resources supporting the study of biodiversity and evolutionary history of Rosaceae.

Acipenserid herpesvirus 2 genome and partial validation of a qPCR for its detection in white sturgeon Acipenser transmontanus

White sturgeon Acipenser transmontanus is the primary species used for caviar and sturgeon meat production in the USA. An important pathogen of white sturgeon is acipenserid herpesvirus 2 (AciHV-2). In this study, 4 archived isolates from temporally discrete natural outbreaks spanning the past 30 yr were sequenced via Illumina and Oxford Nanopore Technologies platforms. Assemblies of approximately 134 kb were obtained for each isolate, and the putative ATPase subunit of the terminase gene was selected as a potential quantitative PCR (qPCR) target based on sequence conservation among AciHV-2 isolates and low sequence homology with other important viral pathogens. The qPCR was repeatable and reproducible, with a linear dynamic range covering 5 orders of magnitude, an efficiency of approximately 96%, an R2 of 0.9872, and an analytical sensitivity of 103 copies per reaction after 35 cycles. There was no cross-reaction with other known viruses or closely related sturgeon species, and no inhibition by sturgeon DNA. Clinical accuracy was assessed from white sturgeon juveniles exposed to AciHV-2 by immersion. Viral culture (gold standard) and qPCR were in complete agreement for both cell culture negative and cell culture positive samples, indicating that this assay has 100% relative accuracy compared to cell culture during an active outbreak. The availability of a whole-genome sequence for AciHV-2 and a highly specific and sensitive qPCR assay for detection of AciHV-2 in white sturgeon lays a foundation for further studies on host-pathogen interactions while providing a specific and rapid test for AciHV-2 in captive and wild populations.

Genotype and phenotype data standardization, utilization and integration in the big data era for agricultural sciences

Large-scale genotype and phenotype data have been increasingly generated to identify genetic markers, understand gene function and evolution and facilitate genomic selection. These datasets hold immense value for both current and future studies, as they are vital for crop breeding, yield improvement and overall agricultural sustainability. However, integrating these datasets from heterogeneous sources presents significant challenges and hinders their effective utilization. We established the Genotype-Phenotype Working Group in November 2021 as a part of the AgBioData Consortium (https://www.agbiodata.org) to review current data types and resources that support archiving, analysis and visualization of genotype and phenotype data to understand the needs and challenges of the plant genomic research community. For 2021-22, we identified different types of datasets and examined metadata annotations related to experimental design/methods/sample collection, etc. Furthermore, we thoroughly reviewed publicly funded repositories for raw and processed data as well as secondary databases and knowledgebases that enable the integration of heterogeneous data in the context of the genome browser, pathway networks and tissue-specific gene expression. Based on our survey, we recommend a need for (i) additional infrastructural support for archiving many new data types, (ii) development of community standards for data annotation and formatting, (iii) resources for biocuration and (iv) analysis and visualization tools to connect genotype data with phenotype data to enhance knowledge synthesis and to foster translational research. Although this paper only covers the data and resources relevant to the plant research community, we expect that similar issues and needs are shared by researchers working on animals. Database URL: https://www.agbiodata.org.

Comparative De Novo and Pan-Genome Analysis of MDR Nosocomial Bacteria Isolated from Hospitals in Jeddah, Saudi Arabia

Multidrug-resistant (MDR) bacteria are one of the most serious threats to public health, and one of the most important types of MDR bacteria are those that are acquired in a hospital, known as nosocomial. This study aimed to isolate and identify MDR bacteria from selected hospitals in Jeddah and analyze their antibiotic-resistant genes. Bacteria were collected from different sources and wards of hospitals in Jeddah City. Phoenix BD was used to identify the strains and perform susceptibility testing. Identification of selected isolates showing MDR to more than three classes on antibiotics was based on 16S rRNA gene and whole genome sequencing. Genes conferring resistance were characterized using de novo and pan-genome analyses. In total, we isolated 108 bacterial strains, of which 75 (69.44%) were found to be MDR. Taxonomic identification revealed that 24 (32%) isolates were identified as Escherichia coli, 19 (25.3%) corresponded to Klebsiella pneumoniae, and 17 (22.67%) were methicillin-resistant Staphylococcus aureus (MRSA). Among the Gram-negative bacteria, K. pneumoniae isolates showed the highest resistance levels to most antibiotics. Of the Gram-positive bacteria, S. aureus (MRSA) strains were noticed to exhibit the uppermost degree of resistance to the tested antibiotics, which is higher than that observed for K. pneumoniae isolates. Taken together, our results illustrated that MDR Gram-negative bacteria are the most common cause of nosocomial infections, while MDR Gram-positive bacteria are characterized by a wider antibiotic resistance spectrum. Whole genome sequencing found the appearance of antibiotic resistance genes, including SHV, OXA, CTX-M, TEM-1, NDM-1, VIM-1, ere(A), ermA, ermB, ermC, msrA, qacA, qacB, and qacC.

Genome sequencing and application of Taiwanese macaque Macaca cyclopis

Formosan macaque (Macaca cyclopis) is the only non-human primate in Taiwan Island. We performed de novo hybrid assembly for M. cyclopis using Illumina paired-end short reads, mate-pair reads and Nanopore long reads and obtained 5065 contigs with a N50 of 2.66 megabases. M. cyclopis contigs >  = 10 kb were assigned to chromosomes using Indian rhesus macaque (Macaca mulatta mulatta) genome assembly Mmul_10 as reference, resulting in a draft of M. cyclopis genome of 2,846,042,475 bases, distributed in 21 chromosomes. The draft genome contains 23,462 transcriptional origins (genes), capable of expressing 716,231 exons in 59,484 transcripts. Genome-based phylogenetic study using the assembled M. cyclopis genome together with genomes of four other macaque species, human, orangutan and chimpanzee showed similar result as previously reported. However, the M. cyclopis species was found to diverge from Chinese M. mulatta lasiota about 1.8 million years ago. Fossil gene analysis detected the presence of gap and pol endogenous viral elements of simian retrovirus in all macaques tested, including M. fascicularis, M. m. mulatta and M. cyclopis. However, M. cyclopis showed ~ 2 times less in number and more uniform in chromosomal locations. The constrain in foreign genome disturbance, presumably due to geographical isolation, should be able to simplify genomics-related investigations, making M. cyclopis an ideal primate species for medical research.

Khairuddin RF, Mahmood Y. Genome assembly composition of the String "ACGT" array: a review of data structure accuracy and performance challenges

Background

The development of sequencing technology increases the number of genomes being sequenced. However, obtaining a quality genome sequence remains a challenge in genome assembly by assembling a massive number of short strings (reads) with the presence of repetitive sequences (repeats). Computer algorithms for genome assembly construct the entire genome from reads in two approaches. The de novo approach concatenates the reads based on the exact match between their suffix-prefix (overlapping). Reference-guided approach orders the reads based on their offsets in a well-known reference genome (reads alignment). The presence of repeats extends the technical ambiguity, making the algorithm unable to distinguish the reads resulting in misassembly and affecting the assembly approach accuracy. On the other hand, the massive number of reads causes a big assembly performance challenge.

Method

The repeat identification method was introduced for misassembly by prior identification of repetitive sequences, creating a repeat knowledge base to reduce ambiguity during the assembly process, thus enhancing the accuracy of the assembled genome. Also, hybridization between assembly approaches resulted in a lower misassembly degree with the aid of the reference genome. The assembly performance is optimized through data structure indexing and parallelization. This article's primary aim and contribution are to support the researchers through an extensive review to ease other researchers' search for genome assembly studies. The study also, highlighted the most recent developments and limitations in genome assembly accuracy and performance optimization.

Results

Our findings show the limitations of the repeat identification methods available, which only allow to detect of specific lengths of the repeat, and may not perform well when various types of repeats are present in a genome. We also found that most of the hybrid assembly approaches, either starting with de novo or reference-guided, have some limitations in handling repetitive sequences as it is more computationally costly and time intensive. Although the hybrid approach was found to outperform individual assembly approaches, optimizing its performance remains a challenge. Also, the usage of parallelization in overlapping and reads alignment for genome assembly is yet to be fully implemented in the hybrid assembly approach.

Conclusion

We suggest combining multiple repeat identification methods to enhance the accuracy of identifying the repeats as an initial step to the hybrid assembly approach and combining genome indexing with parallelization for better optimization of its performance.

Bioinformatics approaches and big data analytics opportunities in improving fisheries and aquaculture

Aquaculture has witnessed an excellent growth rate during the last two decades and offers huge potential to provide nutritional as well as livelihood security. Genomic research has contributed significantly toward the development of beneficial technologies for aquaculture. The existing high throughput technologies like next-generation technologies generate oceanic data which requires extensive analysis using appropriate tools. Bioinformatics is a rapidly evolving science that involves integrating gene based information and computational technology to produce new knowledge for the benefit of aquaculture. Bioinformatics provides new opportunities as well as challenges for information and data processing in new generation aquaculture. Rapid technical advancements have opened up a world of possibilities for using current genomics to improve aquaculture performance. Understanding the genes that govern economically relevant characteristics, necessitates a significant amount of additional research. The various dimensions of data sources includes next-generation DNA sequencing, protein sequencing, RNA sequencing gene expression profiles, metabolic pathways, molecular markers, and so on. Appropriate bioinformatics tools are developed to mine the biologically relevant and commercially useful results. The purpose of this scoping review is to present various arms of diverse bioinformatics tools with special emphasis on practical translation to the aquaculture industry.

Next-generation sequencing approach to investigate genome variability of Parapoxvirus in Canadian muskoxen (Ovibos moschatus)

In 2016, the first orf virus, a double-stranded DNA (dsDNA) virus of the genus parapoxvirus, from a muskox was isolated on Victoria Island, Nunavut (NU), Canada. We used deep sequencing on DNA extracted from orf virus-positive tissues from wild muskoxen from locations on Victoria Island and the adjacent mainland. Orf virus sequence reads derived from four samples were nearly identical. The consensus sequences generated from pooled reads of MxOV comprises of a large contiguous sequence (contig) of 131,759 bp and a smaller right terminal contig of 3552 bp, containing all coding sequences identified as Parapoxvirus. Individual gene comparisons reveal that MxOV shares genetic characteristics with reference strains from both sheep and goat origin. Recombination analysis using Bootscan, MAXCHI, GENECONV, CHIMAERA, SISCAN, and RDP algorithms within the RDP4 software predicted recombination events in two virulence factors, and a large 3000 bp segment of the MxOV genome. Partial B2L nucleotide sequences from strains around the world and other North American isolates were compared to MxOV using MUSCLE alignments and RAxML phylogenetic trees. MxOV was identical to our previously characterized isolate, and shared similarity with orf virus isolated from sheep and goats. The phylogenetic grouping of partial B2L nucleotide sequences did not follow the sample geographic distribution. More full genomes of orf virus, or at least full B2L gene squences, in wildlife are needed especially in North America to better understand the epidemiology of the disease in muskoxen.

Multicenter Preparedness Exercise Enables Rapid Development of Cluster-Specific PCR-Based Screening Assays from Bacterial Genomic Data

PCR-based screening assays targeting strain-specific genetic markers allow the timely detection and specific differentiation of bacterial strains. Especially in situations where an infection cluster occurs, fast assay development is crucial for supporting targeted control measures. However, the turnaround times (TATs) for assay setup may be high due to insufficient knowledge about screening assay methods, workflows, and software tools. Here, two blind-coded and quality-controlled ring trials were performed in which five German laboratories established PCR-based screening assays from genomic data that specifically target selected bacterial clusters within two bacterial monospecies sample panels. While the first ring trial was conducted without a time limit to train the participants and assess assay feasibility, in the second ring trial, a challenging time limit of 2 weeks was set to force fast assay development as soon as genomic data were available. During both ring trials, we detected high interlaboratory variability regarding the screening assay methods and targets, the TATs for assay setup, and the number of screening assays. The participants designed between one and four assays per cluster that targeted cluster-specific unique genetic sequences, genes, or single nucleotide variants using conventional PCRs, high-resolution melting assays, or TaqMan PCRs. Assays were established within the 2-week time limit, with TATs ranging from 4 to 13 days. TaqMan probe delivery times strongly influenced TATs. In summary, we demonstrate that a specific exercise improved the preparedness to develop functional cluster-specific PCR-based screening assays from bacterial genomic data. Furthermore, the parallel development of several assays enhances assay availability.

Highly accurate long reads are crucial for realizing the potential of biodiversity genomics

BACKGROUND

Generating the most contiguous, accurate genome assemblies given available sequencing technologies is a long-standing challenge in genome science. With the rise of long-read sequencing, assembly challenges have shifted from merely increasing contiguity to correctly assembling complex, repetitive regions of interest, ideally in a phased manner. At present, researchers largely choose between two types of long read data: longer, but less accurate sequences, or highly accurate, but shorter reads (i.e., >Q20 or 99% accurate). To better understand how these types of long-read data as well as scale of data (i.e., mean length and sequencing depth) influence genome assembly outcomes, we compared genome assemblies for a caddisfly, Hesperophylax magnus, generated with longer, but less accurate, Oxford Nanopore (ONT) R9.4.1 and highly accurate PacBio HiFi (HiFi) data. Next, we expanded this comparison to consider the influence of highly accurate long-read sequence data on genome assemblies across 6750 plant and animal genomes. For this broader comparison, we used HiFi data as a surrogate for highly accurate long-reads broadly as we could identify when they were used from GenBank metadata.

RESULTS

HiFi reads outperformed ONT reads in all assembly metrics tested for the caddisfly data set and allowed for accurate assembly of the repetitive ~ 20 Kb H-fibroin gene. Across plants and animals, genome assemblies that incorporated HiFi reads were also more contiguous. For plants, the average HiFi assembly was 501% more contiguous (mean contig N50 = 20.5 Mb) than those generated with any other long-read data (mean contig N50 = 4.1 Mb). For animals, HiFi assemblies were 226% more contiguous (mean contig N50 = 20.9 Mb) versus other long-read assemblies (mean contig N50 = 9.3 Mb). In plants, we also found limited evidence that HiFi may offer a unique solution for overcoming genomic complexity that scales with assembly size.

CONCLUSIONS

Highly accurate long-reads generated with HiFi or analogous technologies represent a key tool for maximizing genome assembly quality for a wide swath of plants and animals. This finding is particularly important when resources only allow for one type of sequencing data to be generated. Ultimately, to realize the promise of biodiversity genomics, we call for greater uptake of highly accurate long-reads in future studies.

Improvements to the Gulf pipefish Syngnathus scovelli genome

The Gulf pipefish Syngnathus scovelli has emerged as an important species for studying sexual selection, development, and physiology. Comparative evolutionary genomics research involving fishes from Syngnathidae depends on having a high-quality genome assembly and annotation. However, the first S. scovelli genome assembled using short-read sequences and a smaller RNA-sequence dataset has limited contiguity and a relatively poor annotation. Here, using PacBio long-read high-fidelity sequences and a proximity ligation library, we generate an improved assembly to obtain 22 chromosome-level scaffolds. Compared to the first assembly, the gaps in the improved assembly are smaller, the N75 is larger, and our genome is ~95% BUSCO complete. Using a large body of RNA-Seq reads from different tissue types and NCBI's Eukaryotic Annotation Pipeline, we discovered 28,162 genes, of which 8,061 are non-coding genes. Our new genome assembly and annotation are tagged as a RefSeq genome by NCBI and provide enhanced resources for research work involving S. scovelli..

Transfer of knowledge from model organisms to evolutionarily distant non-model organisms: The coral Pocillopora damicornis membrane signaling receptome

With the ease of gene sequencing and the technology available to study and manipulate non-model organisms, the extension of the methodological toolbox required to translate our understanding of model organisms to non-model organisms has become an urgent problem. For example, mining of large coral and their symbiont sequence data is a challenge, but also provides an opportunity for understanding functionality and evolution of these and other non-model organisms. Much more information than for any other eukaryotic species is available for humans, especially related to signal transduction and diseases. However, the coral cnidarian host and human have diverged over 700 million years ago and homologies between proteins in the two species are therefore often in the gray zone, or at least often undetectable with traditional BLAST searches. We introduce a two-stage approach to identifying putative coral homologues of human proteins. First, through remote homology detection using Hidden Markov Models, we identify candidate human homologues in the cnidarian genome. However, for many proteins, the human genome alone contains multiple family members with similar or even more divergence in sequence. In the second stage, therefore, we filter the remote homology results based on the functional and structural plausibility of each coral candidate, shortlisting the coral proteins likely to have conserved some of the functions of the human proteins. We demonstrate our approach with a pipeline for mapping membrane receptors in humans to membrane receptors in corals, with specific focus on the stony coral, P. damicornis. More than 1000 human membrane receptors mapped to 335 coral receptors, including 151 G protein coupled receptors (GPCRs). To validate specific sub-families, we chose opsin proteins, representative GPCRs that confer light sensitivity, and Toll-like receptors, representative non-GPCRs, which function in the immune response, and their ability to communicate with microorganisms. Through detailed structure-function analysis of their ligand-binding pockets and downstream signaling cascades, we selected those candidate remote homologues likely to carry out related functions in the corals. This pipeline may prove generally useful for other non-model organisms, such as to support the growing field of synthetic biology.

GALA: a computational framework for de novo chromosome-by-chromosome assembly with long reads

High-quality genome assembly has wide applications in genetics and medical studies. However, it is still very challenging to achieve gap-free chromosome-scale assemblies using current workflows for long-read platforms. Here we report on GALA (Gap-free long-read Assembly tool), a computational framework for chromosome-based sequencing data separation and de novo assembly implemented through a multi-layer graph that identifies discordances within preliminary assemblies and partitions the data into chromosome-scale scaffolding groups. The subsequent independent assembly of each scaffolding group generates a gap-free assembly likely free from the mis-assembly errors which usually hamper existing workflows. This flexible framework also allows us to integrate data from various technologies, such as Hi-C, genetic maps, and even motif analyses to generate gap-free chromosome-scale assemblies. As a proof of principle we de novo assemble the C. elegans genome using combined PacBio and Nanopore sequencing data and a rice cultivar genome using Nanopore sequencing data from publicly available datasets. We also demonstrate the proposed method's applicability with a gap-free assembly of the human genome using PacBio high-fidelity (HiFi) long reads. Thus, our method enables straightforward assembly of genomes with multiple data sources and overcomes barriers that at present restrict the application of de novo genome assembly technology.

Elucidating the genomic history of commercially used Bacillus thuringiensis subsp. tenebrionis strain NB176

Bacillus thuringiensis subsp. tenebrionis (Btt) produces a coleopteran-specific crystal protoxin protein (Cry3Aa δ-endotoxin). After its discovery in 1982, the strain NB125 (DSM 5526) was eventually registered in 1990 to control the Colorado potato beetle (Leptinotarsa decemlineata). Gamma-irradiation of NB125 resulted in strain NB176-1 (DSM 5480) that exhibited higher cry3Aa production and became the active ingredient of the plant protection product Novodor^® FC. Here, we report a comparative genome analysis of the parental strain NB125, its derivative NB176-1 and the current commercial production strain NB176. The entire genome sequences of the parental and derivative strains were deciphered by a hybrid de novo approach using short (Illumina) and long (Nanopore) read sequencing techniques. Genome assembly revealed a chromosome of 5.4 to 5.6 Mbp and six plasmids with a size range from 14.9 to 250.5 kbp for each strain. The major differences among the original NB125 and the derivative strains NB176-1 and NB176 were an additional copy of the cry3Aa gene, which translocated to another plasmid as well as a chromosomal deletion (~ 178 kbp) in NB176. The assembled genome sequences were further analyzed in silico for the presence of virulence and antimicrobial resistance (AMR) genes.

Immunological Studies to Understand Hybrid/Recombinant Variants of SARS-CoV-2

The zoonotic SARS-CoV-2 virus was present before the onset of the pandemic. It undergoes evolution, adaptation, and selection to develop variants that gain high transmission rates and virulence, resulting in the pandemic. Structurally, the spike protein of the virus is required for binding to ACE2 receptors of the host cells. The gene coding for the spike is known to have a high propensity of mutations, as a result generating numerous variants. The variants can be generated by random point mutations or recombination during replication. However, SARS-CoV-2 can also produce hybrid variants on co-infection of the host by two distinct lineages of the virus. The genomic sequences of the two variants undergo recombination to produce the hybrid variants. Additionally, these sub-variants also contain numerous mutations from both the parent variants, as well as some novel mutations unique to the hybrids. The hybrid variants (XD, XE, and XF) can be identified through numerous techniques, such as peak PCR, NAAT, and hybrid capture SARS-CoV-2 NGS (next generation sequencing) assay, etc., but the most accurate approach is genome sequencing. There are numerous immunological diagnostic assays, such as ELISA, chemiluminescence immunoassay, flow-cytometry-based approaches, electrochemiluminescence immunoassays, neutralization assays, etc., that are also designed and developed to provide an understanding of the hybrid variants, their pathogenesis, and other reactions. The objective of our study is to comprehensively analyze the variants of SARS-CoV-2, especially the hybrid variants. We have also discussed the techniques available for the identification of hybrids, as well as the immunological assays and studies for analyzing the hybrid variants.

Alkaloid production and response to natural adverse conditions in Peganum harmala: in silico transcriptome analyses

Peganum harmala is a valuable wild plant that grows and survives under adverse conditions and produces pharmaceutical alkaloid metabolites. Using different assemblers to develop a transcriptome improves the quality of assembled transcriptome. In this study, a concrete and accurate method for detecting stress-responsive transcripts by comparing stress-related gene ontology (GO) terms and public domains was designed. An integrated transcriptome for P. harmala including 42 656 coding sequences was created by merging de novo assembled transcriptomes. Around 35 000 transcripts were annotated with more than 90% resemblance to three closely related species of Citrus, which confirmed the robustness of the assembled transcriptome; 4853 stress-responsive transcripts were identified. CYP82 involved in alkaloid biosynthesis showed a higher number of transcripts in P. harmala than in other plants, indicating its diverse alkaloid biosynthesis attributes. Transcription factors (TFs) and regulatory elements with 3887 transcripts comprised 9% of the transcriptome. Among the TFs of the integrated transcriptome, cystein2/histidine2 (C2H2) and WD40 repeat families were the most abundant. The Kyoto Encyclopedia of Genes and Genomes (KEGG) MAPK (mitogen-activated protein kinase) signaling map and the plant hormone signal transduction map showed the highest assigned genes to these pathways, suggesting their potential stress resistance. The P. harmala whole-transcriptome survey provides important resources and paves the way for functional and comparative genomic studies on this plant to discover stress-tolerance-related markers and response mechanisms in stress physiology, phytochemistry, ecology, biodiversity, and evolution. P. harmala can be a potential model for studying adverse environmental cues and metabolite biosynthesis and a major source for the production of various alkaloids.

Comparative Genomics for Evolutionary Cell Biology Using AMOEBAE: Understanding the Golgi and Beyond

Taking an evolutionary approach to cell biology can yield important new information about how the cell works and how it evolved to do so. This is true of the Golgi apparatus, as it is of all systems within the cell. Comparative genomics is one of the crucial first steps to this line of research, but comes with technical challenges that must be overcome for rigor and robustness. We here introduce AMOEBAE, a workflow for mid-range scale comparative genomic analyses. It allows for customization of parameters, queries, and taxonomic sampling of genomic and transcriptomics data. This protocol article covers the rationale for an evolutionary approach to cell biological study (i.e., when would AMOEBAE be useful), how to use AMOEBAE, and discussion of limitations. It also provides an example dataset, which demonstrates that the Golgi protein AP4 Epsilon is present as the sole retained subunit of the AP4 complex in basidiomycete fungi. AMOEBAE can facilitate comparative genomic studies by balancing reproducibility and speed with user-input and interpretation. It is hoped that AMOEBAE or similar tools will encourage cell biologists to incorporate an evolutionary context into their research.

Domestication-associated PHYTOCHROME C is a flowering time repressor and a key factor determining Setaria as a short-day plant

Phytochromes play vital roles in the regulation of flowering time, but little is known in Panicoideae species, especially the C₄ model Setaria. Here, genomic variations of PHYTOCHROME C (PHYC) between wild and cultivated Setaria gene pools were analysed and three SiphyC mutants were identified. The function of SiPHYC was verified by CRISPR-Cas9 approach and transcriptome sequencing. Furthermore, efficiency of indoor cultivation of SiphyC mutants were systematically evaluated. An extreme purified selection of PHYC was detected in wild to cultivated domestication process of Setaria. SiphyC mutants and knockout transgenic plants showed an early heading date and a loss of response to short-day photoperiod. Furthermore, variable expression of SiFTa, SiMADS14 and SiMADS15 might be responsible for promoting flowering of SiphyC mutants. Moreover, SiphyC mutant was four times that of the indoor plot ratio of wild-type and produced over 200 seeds within 45 d per individual. Our results suggest that domestication-associated SiPHYC repressed flowering and determined Setaria as a short-day plant, and SiphyC mutants possess the potential for creating efficient indoor cultivation system suitable for research on Setaria as a model, and either for maize or sorghum as well.

Characterization of driver mutations in Chinese non-small cell lung cancer patients using a novel targeted sequencing panel

Background

The identification of driver mutations has greatly promoted the precise diagnosis and treatment of non-small cell lung cancer (NSCLC), but there is lack of targeted sequencing panels specifically designed and applied to Chinese NSCLC patients. This study aimed to design and validate of a novel sequencing panel for comprehensive characterization of driver mutations in Chinese NSCLC patients, facilitating further exploration of downstream pathway alterations and therapeutic utility.

Methods

A novel target sequencing panel including 21 driver genes was designed and examined in a cohort of 260 Chinese NSCLC patients who underwent surgery in Peking Union Medical College Hospital (PUMCH). Genetic alterations were identified and further analyzed for driver mutations, downstream pathways and therapeutic utilities.

Results

The most frequently identified driver mutations in PUMCH NSCLC cohort were on genes TP53 (28%), EGFR (27%) and PIK3CA (19%) for lung adenocarcinoma (LUAD), and TP53 (41%), PIK3CA (14%) and CDKN2A (13%) for lung squamous cell carcinoma (LUSC), respectively. Downstream pathway analysis revealed common pathways like G1_AND_S1_PHASES pathway were shared not only between LUAD and LUSC patients, but also among three different NSCLC cohorts, while other pathways were subtype-specific, like the unique enrichment of SHC1_EVENT_IN_EGFR_SIGNALING pathway in LUAD patients, and P38_ALPHA_BETA_DOWNSTREAM pathway in LUSC patients, respectively. About 60% of both LUAD and LUSC patients harbored driver mutations as sensitive biomarkers for different targeted therapies, covering not only frequent mutations like EGFR L858R mutation, but also rare mutations like BRAF D594N mutation.

Conclusions

Our study provides a novel target sequencing panel suitable for Chinese NSCLC patients, which can effectively identify driver mutations, analyze downstream pathway alterations and predict therapeutic utility. Overall it is promising to further optimize and apply this panel in clinic with convenience and effectiveness.

Evaluating Illumina-, Nanopore-, and PacBio-based genome assembly strategies with the bald notothen, Trematomus borchgrevinki

For any genome-based research, a robust genome assembly is required. De novo assembly strategies have evolved with changes in DNA sequencing technologies and have been through at least 3 phases: (1) short-read only, (2) short- and long-read hybrid, and (3) long-read only assemblies. Each of the phases has its own error model. We hypothesized that hidden short-read scaffolding errors and erroneous long-read contigs degrade the quality of short- and long-read hybrid assemblies. We assembled the genome of Trematomus borchgrevinki from data generated during each of the 3 phases and assessed the quality problems we encountered. We developed strategies such as k-mer-assembled region replacement, parameter optimization, and long-read sampling to address the error models. We demonstrated that a k-mer-based strategy improved short-read assemblies as measured by Benchmarking Universal Single-Copy Ortholog while mate-pair libraries introduced hidden scaffolding errors and perturbed Benchmarking Universal Single-Copy Ortholog scores. Furthermore, we found that although hybrid assemblies can generate higher contiguity they tend to suffer from lower quality. In addition, we found long-read-only assemblies can be optimized for contiguity by subsampling length-restricted raw reads. Our results indicate that long-read contig assembly is the current best choice and that assemblies from phase I and phase II were of lower quality.

Leveraging an existing whole-genome resequencing population data set to characterize toll-like receptor gene diversity in a threatened bird

Species recovery programs are increasingly using genomic data to measure neutral genetic diversity and calculate metrics like relatedness. While these measures can inform conservation management, determining the mechanisms underlying inbreeding depression requires information about functional genes associated with adaptive or maladaptive traits. Toll-like receptors (TLRs) are one family of functional genes, which play a crucial role in recognition of pathogens and activation of the immune system. Previously, these genes have been analysed using species-specific primers and PCR. Here, we leverage an existing short-read reference genome, whole-genome resequencing population data set, and bioinformatic tools to characterize TLR gene diversity in captive and wild tchūriwat'/tūturuatu/shore plover (Thinornis novaeseelandiae), a threatened bird endemic to Aotearoa New Zealand. Our results show that TLR gene diversity in tchūriwat'/tūturuatu is low, and forms two distinct captive and wild genetic clusters. The bioinformatic approach presented here has broad applicability to other threatened species with existing genomic resources in Aotearoa New Zealand and beyond.

Development of Genomic Resources in Mexican Bursera (Section: Bullockia: Burseraceae): Genome Assembly, Annotation, and Marker Discovery for Three Copal Species

Bursera comprises ~100 tropical shrub and tree species, with the center of the species diversification in Mexico. The genomic resources developed for the genus are scarce, and this has limited the study of the gene flow, local adaptation, and hybridization dynamics. In this study, based on ~155 million Illumina paired-end reads per species, we performed a de novo genome assembly and annotation of three Bursera species of the Bullockia section: Bursera bipinnata, Bursera cuneata, and Bursera palmeri. The total lengths of the genome assemblies were 253, 237, and 229 Mb for B. cuneata, B. palmeri, and B. bipinnata, respectively. The assembly of B. palmeri retrieved the most complete and single-copy BUSCOs (87.3%) relative to B. cuneata (86.5%) and B. bipinnata (76.6%). The ab initio gene prediction recognized between 21,000 and 32,000 protein-coding genes. Other genomic features, such as simple sequence repeats (SSRs), were also detected. Using the de novo genome assemblies as a reference, we identified single-nucleotide polymorphisms (SNPs) for a set of 43 Bursera individuals. Moreover, we mapped the filtered reads of each Bursera species against the chloroplast genomes of five Burseraceae species, obtaining consensus sequences ranging from 156 to 160 kb in length. Our work contributes to the generation of genomic resources for an important but understudied genus of tropical-dry-forest species.

Widespread false gene gains caused by duplication errors in genome assemblies

Background

False duplications in genome assemblies lead to false biological conclusions. We quantified false duplications in popularly used previous genome assemblies for platypus, zebra finch, and Anna’s Hummingbird, and their new counterparts of the same species generated by the Vertebrate Genomes Project, of which the Vertebrate Genomes Project pipeline attempted to eliminate false duplications through haplotype phasing and purging. These assemblies are among the first generated by the Vertebrate Genomes Project where there was a prior chromosomal level reference assembly to compare with.

Results

Whole genome alignments revealed that 4 to 16% of the sequences are falsely duplicated in the previous assemblies, impacting hundreds to thousands of genes. These lead to overestimated gene family expansions. The main source of the false duplications is heterotype duplications, where the haplotype sequences were relatively more divergent than other parts of the genome leading the assembly algorithms to classify them as separate genes or genomic regions. A minor source is sequencing errors. Ancient ATP nucleotide binding gene families have a higher prevalence of false duplications compared to other gene families. Although present in a smaller proportion, we observe false duplications remaining in the Vertebrate Genomes Project assemblies that can be identified and purged.

Conclusions

This study highlights the need for more advanced assembly methods that better separate haplotypes and sequence errors, and the need for cautious analyses on gene gains.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02764-1.

Low carbon-to-nitrogen ratio digestate from high-rate anaerobic baffled reactor facilitates heterotrophic/autotrophic nitrifiers involved in nitrogen removal

In this study, baffled anaerobic-aerobic reactors (AOBRs) with modified basalt fiber (MBF) carriers and felt were used to treat domestic wastewater (DWW). The influent was first treated in anaerobic compartments, with the NH₄⁺-N containing digestate refluxed into aerobic compartment for nitrification. The nitrified liquid was channeled to the anaerobic compartments for further denitrification. Under optimal conditions, AOBR with MBF carriers could remove 91% chemical oxygen demand (COD) and 81% total nitrogen (TN), with biomass production increased by 7.6%, 4.5% and 8.7% in three successive anaerobic compartments compared to the control. Biological viability analysis showed that live cells outnumbered dead cells in bio-nests. Metagenomics analysis showed that multiple metabolic pathways accounted for nitrogen conversion in anaerobic and aerobic compartments. More importantly, low COD/TN ratio digestate facilitated heterotrophic nitrification-aerobic denitrification (HN-AD) species growth in aerobic compartment. This study provides a promising strategy to source treatment of DWW from urban communities.

Validation of reference-assisted assembly using existing and novel Heliothine genomes

Chloridea subflexa and Chloridea virescens are a pair of closely related noctuid species exhibiting pheromone-based sexual isolation and divergent host plant preferences. We produced a novel Illumina short read C. subflexa genome assembly and an improved C. virescens genome assembly, which offer opportunities to study the genomic basis for evolutionarily important traits in this lepidopteran family with few genomic resources. We then examined the feasibility of reference-assisted assembly, an approach that leverages existing high quality genomic resources for genome improvement in closely related taxa and applied it to our Heliothine genomes. Our work demonstrates that reference-assisted assembly has the potential to enhance contiguity and completeness of existing insect genomic resources with minimal additional laboratory costs. We conclude by discussing both the potential and pitfalls of reference-assisted assembly according to the intended downstream assembly application.

Insights from cyanobacterial genomic and transcriptomic analyses into adaptation strategies in terrestrial environments

Phylogenomic analysis of Nostoc sp. MG11, a terrestrial cyanobacterium, and some terrestrial and freshwater Nostoc strains showed that the terrestrial strains grouped together in a distinctive clade, which reveals the effect of habitat on shaping Nostoc genomes. Terrestrial strains showed larger genomes and had higher predicted CDS contents than freshwater strains. Comparative genomic analysis demonstrated that genome expansion in the terrestrial Nostoc is supported by an increase in copy number of the core genes and acquisition of shared genes. Transcriptomic profiling analysis under desiccation stress revealed that Nostoc sp. MG11 protected its cell by induction of catalase, proteases, sucrose synthase, trehalose biosynthesis and maltodextrin utilization genes and maintained its normal metabolism during this condition by up-regulation of genes related to phycobilisomes and light reactions of photosynthesis, CO₂ fixation and protein metabolism. These results provide insights into the strategies related to survival and adaptation of Nostoc strains to terrestrial environments.

CRISPR-CasB technology in forensic DNA analysis: challenges and solutions

CRISPR-Cas technology has revolutionized the field of biotechnology with its precise therapeutic use from genetic as well as infectious diseases point of view. This technology is rapidly evolving to single tool enabling site-directed cut in the genome and highly specific activation or inhibition of gene expression or the exchange of single bases. Besides clinical applications, CRISPR-Cas technology has also shown promising use in the field of forensic DNA analysis. Enrichment of targeted genetic marker for identification followed by sequencing and non-PCR-dependent technique ensures the use of CRISPR-Cas technology in challenging forensic biological samples. The use of this advanced technology is also deemed helpful in mixed profile attribution, mostly in LCN contributors and the generation of a useful DNA profile in degraded samples. Besides its useful applications in forensic DNA analysis, CRISPR-Cas technology poses a huge threat from the generation of ghost DNA profiles by modification/alteration of target genetic markers. Forensic DNA analysts should carry out analysis of additional markers such as non-CODIS markers, Y-, X-chromosome markers, and mitochondrial DNA sequencing in a suspected ghost DNA profile case. KEY POINTS: • CRISPR-Cas9 technique is useful in analyzing LCN, mixed and degraded samples • Alteration of DNA using this technique can lead to generation of ghost DNA profiles • Alternative genetic markers and methylation pattern may detect a ghost DNA profile.

Intra-specific copy number variation of MHC class II genes in the Siamese fighting fish

Duplicates of genes for major histocompatibility complex (MHC) molecules can be subjected to selection independently and vary markedly in their evolutionary rates, sequence polymorphism, and functional roles. Therefore, without a thorough understanding of their copy number variation (CNV) in the genome, the MHC-dependent fitness consequences within a species could be misinterpreted. Studying the intra-specific CNV of this highly polymorphic gene, however, has long been hindered by the difficulties in assigning alleles to loci and the lack of high-quality genomic data. Here, using the high-quality genome of the Siamese fighting fish (Betta splendens), a model for mate choice studies, and the whole-genome sequencing (WGS) data of 17 Betta species, we achieved locus-specific amplification of their three classical MHC class II genes - DAB1, DAB2, and DAB3. By performing quantitative PCR and depth-of-coverage analysis using the WGS data, we revealed intra-specific CNV at the DAB3 locus. We identified individuals that had two allelic copies (i.e., heterozygous or homozygous) or one allele (i.e., hemizygous) and individuals without this gene. The CNV was due to the deletion of a 20-kb-long genomic region harboring both the DAA3 and DAB3 genes. We further showed that the three DAB genes were under different modes of selection, which also applies to their corresponding DAA genes that share similar pattern of polymorphism. Our study demonstrates a combined approach to study CNV within a species, which is crucial for the understanding of multigene family evolution and the fitness consequences of CNV.

Genome Assembly and Evolutionary Analysis of the Mandarin Duck Aix galericulata Reveal Strong Genome Conservation among Ducks

The mandarin duck, Aix galericulata, is popular in East Asian cultures and displays exaggerated sexual dimorphism, especially in feather traits during breeding seasons. We generated and annotated the first mandarin duck de novo assembly, which was 1.08 Gb in size and encoded 16,615 proteins. Using a phylogenomic approach calibrated with fossils and molecular divergences, we inferred that the last common ancestor of ducks occurred 13.3-26.7 Ma. The majority of the mandarin duck genome repetitive sequences belonged to the chicken repeat 1 (CR1) retroposon CR1-J2_Pass, which underwent a duck lineage-specific burst. Synteny analyses among ducks revealed infrequent chromosomal rearrangements in which breaks were enriched in LINE retrotransposons and DNA transposons. The calculation of the dN/dS ratio revealed that the majority of duck genes were under strong purifying selection. The expanded gene families in the mandarin duck are primarily involved in olfactory perception as well as the development and morphogenesis of feather and branching structures. This new reference genome will improve our understanding of the morphological and physiological characteristics of ducks and provide a valuable resource for functional genomics studies to investigate the feather traits of the mandarin duck.

Analysis and Performance Assessment of the Whole Genome Bisulfite Sequencing Data Workflow: Currently Available Tools and a Practical Guide to Advance DNA Methylation Studies

DNA methylation is associated with transcriptional repression, genomic imprinting, stem cell differentiation, embryonic development, and inflammation. Aberrant DNA methylation can indicate disease states, including cancer and neurological disorders. Therefore, the prevalence and location of 5-methylcytosine in the human genome is a topic of interest. Whole-genome bisulfite sequencing (WGBS) is a high-throughput method for analyzing DNA methylation. This technique involves library preparation, alignment, and quality control. Advancements in epigenetic technology have led to an increase in DNA methylation studies. This review compares the detailed experimental methodology of WGBS using accessible and up-to-date analysis tools. Practical codes for WGBS data processing are included as a general guide to assist progress in DNA methylation studies through a comprehensive case study.

Whole-genome sequence analysis through online web interfaces: a review

The recent development of whole-genome sequencing technologies paved the way for understanding the genomes of microorganisms. Every whole-genome sequencing (WGS) project requires a considerable cost and a massive effort to address the questions at hand. The final step of WGS is data analysis. The analysis of whole-genome sequence is dependent on highly sophisticated bioinformatics tools that the research personal have to buy. However, many laboratories and research institutions do not have the bioinformatics capabilities to analyze the genomic data and therefore, are unable to take maximum advantage of whole-genome sequencing. In this aspect, this study provides a guide for research personals on a set of bioinformatics tools available online that can be used to analyze whole-genome sequence data of bacterial genomes. The web interfaces described here have many advantages and, in most cases exempting the need for costly analysis tools and intensive computing resources.

Sequencing and Reconstructing Helminth Mitochondrial Genomes Directly from Genomic Next-Generation Sequencing Data

We present a detailed method for extraction of high-molecular weight genomic DNA suitable for numerous DNA sequencing applications, and a straightforward in silico approach for reconstructing novel mitochondrial (mt) genomes directly from total genomic DNA extracts derived from next-generation sequencing (NGS) data sets. The in silico post-sequencing pipeline described is fast, accurate, and highly efficient, with modest memory requirements that can be performed using a standard desktop computer. The approach is particularly effective for obtaining mitochondrial genomes for species with little or no mitochondrial sequence information currently available and overcomes many of the limitations of traditional strategies. The described methodologies are also applicable for metagenomics sequencing from mixed or pooled samples containing multiple species and subsequent specific assembly of specific mitochondrial genomes.

Actinobacteria From Desert: Diversity and Biotechnological Applications

Deserts, as an unexplored extreme ecosystem, are known to harbor diverse actinobacteria with biotechnological potential. Both multidrug-resistant (MDR) pathogens and environmental issues have sharply raised the emerging demand for functional actinobacteria. From 2000 to 2021, 129 new species have been continuously reported from 35 deserts worldwide. The two largest numbers are of the members of the genera Streptomyces and Geodermatophilus, followed by other functional extremophilic strains such as alkaliphiles, halotolerant species, thermophiles, and psychrotolerant species. Improved isolation strategies for the recovery of culturable and unculturable desert actinobacteria are crucial for the exploration of their diversity and offer a better understanding of their survival mechanisms under extreme environmental stresses. The main bioprospecting processes involve isolation of target actinobacteria on selective media and incubation and selection of representatives from isolation plates for further investigations. Bioactive compounds obtained from desert actinobacteria are being continuously explored for their biotechnological potential, especially in medicine. To date, there are more than 50 novel compounds discovered from these gifted actinobacteria with potential antimicrobial activities, including anti-MDR pathogens and anti-inflammatory, antivirus, antifungal, antiallergic, antibacterial, antitumor, and cytotoxic activities. A range of plant growth-promoting abilities of the desert actinobacteria inspired great interest in their agricultural potential. In addition, several degradative, oxidative, and other functional enzymes from desert strains can be applied in the industry and the environment. This review aims to provide a comprehensive overview of desert environments as a remarkable source of diverse actinobacteria while such rich diversity offers an underexplored resource for biotechnological exploitations.

Fourth generation biofuel from genetically modified algal biomass: Challenges and future directions

Genetic engineering applications in the field of biofuel are rapidly expanding due to their potential to boost biomass productivity while lowering its cost and enhancing its quality. Recently, fourth-generation biofuel (FGB), which is biofuel obtained from genetically modified (GM) algae biomass, has gained considerable attention from academic and industrial communities. However, replacing fossil resources with FGB is still beset with many challenges. Most notably, technical aspects of genetic modification operations need to be more fully articulated and elaborated. However, relatively little attention has been paid to GM algal biomass. There is a limited number of reviews on the progress and challenges faced in the algal genetics of FGB. Therefore, the present review aims to fill this gap in the literature by recapitulating the findings of recent studies and achievements on safe and efficient genetic manipulation in the production of FGB. Then, the essential issues and parameters related to genome editing in algal strains are highlighted. Finally, the main challenges to FGB pertaining to the diffusion risk and regulatory frameworks are addressed. This review concluded that the technical and biosafety aspects of FGB, as well as the complexity and diversity of the related regulations, legitimacy concerns, and health and environmental risks, are among the most important challenges that require a strong commitment at the national/international levels to reach a global consensus.

Ancient hybridization patterns between bighorn and thinhorn sheep

Whole-genome sequencing has advanced the study of species evolution, including the detection of genealogical discordant events such as ancient hybridization and incomplete lineage sorting (ILS). The evolutionary history of bighorn (Ovis canadensis) and thinhorn (Ovis dalli) sheep present an ideal system to investigate evolutionary discordance due to their recent and rapid radiation and putative secondary contact between bighorn and thinhorn sheep subspecies, specifically the dark pelage Stone sheep (O. dalli stonei) and predominately white Dall sheep (O. dalli dalli), during the last ice age. Here, we used multiple genomes of bighorn and thinhorn sheep, together with snow (O. nivicola) and the domestic sheep (O. aries) as outgroups, to assess their phylogenomic history, potential introgression patterns and their adaptive consequences. Among the Pachyceriforms (snow, bighorn and thinhorn sheep) a consistent monophyletic species tree was retrieved; however, many genealogical discordance patterns were observed. Alternative phylogenies frequently placed Stone and bighorn as sister clades. This relationship occurred more often and was less divergent than that between Dall and bighorn. We also observed many blocks containing introgression signal between Stone and bighorn genomes in which coat colour genes were present. Introgression signals observed between Dall and bighorn were more random and less frequent, and therefore probably due to ILS or intermediary secondary contact. These results strongly suggest that Stone sheep originated from a complex series of events, characterized by multiple, ancient periods of secondary contact with bighorn sheep.

Comparative Analysis of Genotyping by Sequencing and Whole-Genome Sequencing Methods in Diversity Studies of Olea europaea L

Olive, Olea europaea L., is a tree of great economic and cultural importance in the Mediterranean basin. Thousands of cultivars have been described, of which around 1200 are conserved in the different olive germplasm banks. The genetic characterisation of these cultivars can be performed in different ways. Whole-genome sequencing (WGS) provides more information than the reduced representation methods such as genotype by sequencing (GBS), but at a much higher cost. This may change as the cost of sequencing continues to drop, but, currently, genotyping hundreds of cultivars using WGS is not a realistic goal for most research groups. Our aim is to systematically compare both methodologies applied to olive genotyping and summarise any possible recommendations for the geneticists and molecular breeders of the olive scientific community. In this work, we used a selection of 24 cultivars from an olive core collection from the World Olive Germplasm Collection of the Andalusian Institute of Agricultural and Fisheries Research and Training (WOGBC), which represent the most of the cultivars present in cultivated fields over the world. Our results show that both methodologies deliver similar results in the context of phylogenetic analysis and popular population genetic analysis methods such as clustering. Furthermore, WGS and GBS datasets from different experiments can be merged in a single dataset to perform these analytical methodologies with proper filtering. We also tested the influence of the different olive reference genomes in this type of analysis, finding that they have almost no effect when estimating genetic relationships. This work represents the first comparative study between both sequencing techniques in olive. Our results demonstrate that the use of GBS is a perfectly viable option for replacing WGS and reducing research costs when the goal of the experiment is to characterise the genetic relationship between different accessions. Besides this, we show that it is possible to combine variants from GBS and WGS datasets, allowing the reuse of publicly available data.