Twelve quick steps for genome assembly and annotation in the classroom

Evolutionary insights into adaptation of hemocyanins from deep-sea hydrothermal vent shrimps

Deep-sea hydrothermal vent shrimps inhabit environments with low oxygen levels and may even be exposed to hypoxic conditions. In response, their respiratory pigment, hemocyanin (Hc) may undergo molecular adaptations to enable them to survive in such extreme ecosystems. Therefore, we sampled four Alvinocarididae species from hydrothermal vents in the northern Central Indian Ridge and two types of Hc genes (α and γ) were observed. Employing the branch model, we detected positive selection for the deep-sea hydrothermal vent lineage, including 11 Decapoda species. Furthermore, using the branch-site model, we identified a putative mutant residue (Leu226, Ser377, and Ile390) close to the active site of Hc. Moreover, our results suggested potential molecular docking between two α-type Hc proteins. Thus, this study provides valuable and novel perspectives on the functional significance of the Hc gene in deep-sea hydrothermal vent shrimps, laying the foundation for future investigations in this intriguing area of research.

Collaborative metabolic curation of an emerging model marine bacterium, Alteromonas macleodii ATCC 27126

Inferring the metabolic capabilities of an organism from its genome is a challenging process, relying on computationally-derived or manually curated metabolic networks. Manual curation can correct mistakes in the draft network and add missing reactions based on the literature, but requires significant expertise and is often the bottleneck for high-quality metabolic reconstructions. Here, we present a synopsis of a community curation workshop for the model marine bacterium Alteromonas macleodii ATCC 27126 and its genome database in BioCyc, focusing on pathways for utilizing organic carbon and nitrogen sources. Due to the scarcity of biochemical information or gene knock-outs, the curation process relied primarily on published growth phenotypes and bioinformatic analyses, including comparisons with related Alteromonas strains. We report full pathways for the utilization of the algal polysaccharides alginate and pectin in contrast to inconclusive evidence for one-carbon metabolism and mixed acid fermentation, in accordance with the lack of growth on methanol and formate. Pathways for amino acid degradation are ubiquitous across Alteromonas macleodii strains, yet enzymes in the pathways for the degradation of threonine, tryptophan and tyrosine were not identified. Nucleotide degradation pathways are also partial in ATCC 27126. We postulate that demonstrated growth on nitrate as sole nitrogen source proceeds via a nitrate reductase pathway that is a hybrid of known pathways. Our evidence highlights the value of joint and interactive curation efforts, but also shows major knowledge gaps regarding Alteromonas metabolism. The manually-curated metabolic reconstruction is available as a "Tier-2" database on BioCyc.

Evaluation of strategies for evidence-driven genome annotation using long-read RNA-seq

While the production of a draft genome has become more accessible due to long-read sequencing, the annotation of these new genomes has not been developed at the same pace. Long-read RNA sequencing offers a promising solution for enhancing gene annotation. In this study, we explore how sequencing platforms, Oxford Nanopore R9.4.1 chemistry or Pacific Biosciences (PacBio) Sequel II CCS, and data processing methods influence evidence-driven genome annotation using long reads. Incorporating PacBio transcripts into our annotation pipeline significantly outperformed traditional methods, such as ab initio predictions and short-read-based annotations. We applied this strategy to a nonmodel species, the Florida manatee, and compared our results to existing short-read-based annotation. At the loci level, both annotations were highly concordant, with 90% agreement. However, at the transcript level, the agreement was only 35%. We identified 4906 novel loci, represented by 5707 isoforms, with 64% of these isoforms matching known sequences in other mammalian species. Overall, our findings underscore the importance of using high-quality curated transcript models in combination with ab initio methods for effective genome annotation.

Genome assembly, gene content, and plastic gene expression responses to salinity changes in the Brackishwater Clam (Corbicula japonica) from a dynamic estuarine environment

Estuaries are dynamic transition zones between marine and freshwater environments, where salinity varies greatly on spatial and temporal scales. The temporal salinity fluctuations of these habitats require organisms to rapidly regulate ionic concentrations and osmotic pressure to survive in these dynamic conditions. Understanding the extent of plasticity of euryhaline animals is vital for predicting their responses and resilience to salinity change. We generated the first high-resolution genome and transcriptome sequences of C. japonica. In comparison with 11 other molluscan genomes, the C. japonica genome displayed striking expansions of putative neuron-related genes and gene families. The involvement of these genes in the glutamate/GABA-glutamine and glycine cycle suggests a possible contribution to the excitation of neuronal networks, particularly under high salinity conditions. This study contributes to our understanding of mechanisms underlying the rapid responses of estuarine species to changing conditions and raises many intriguing hypotheses and questions for future investigation.

Genome-Wide Identification and Biochemical Characterization of Alcohol Acyltransferases for Aroma Generation in Wickerhamomyces subpelliculosus Isolates from Fermented Food

The importance of nonconventional yeasts has increasingly been highlighted, particularly for aroma formation in fermented foods. Here, we performed de novo whole-genome sequencing of Wickerhamomyces subpelliculosus, which produces a variety of volatile flavor compounds, leading to the identification of the alcohol acyltransferase (AATase) family of genes. The genome of W. subpelliculosus contains seven AATase genes, encoding alcohol-O-acetyltransferases (ATFs) and ethanol acetyltransferase 1 (EAT1) for acetate ester formation, along with ethanol hexanoyl transferase 1 (EHT1) for ethyl ester formation. Among five ATF homologues, only WsATF5 showed acetyltransferase activity toward myriocin, a structural analogue of sphingosine. In contrast, heterologous expression of WsEHT1 and WsEAT1 in Saccharomyces cerevisiae promoted the production of ethyl decanoate and ethyl acetate, respectively, supporting their AATase activity. The enzymatic activity analyses revealed the additional alcoholysis activity of WsEAT1 and the thioesterase activity of WsEHT1. Subcellular localization analysis indicated that WsEAT1 was localized in the mitochondria, WsEHT1 in the endoplasmic reticulum and lipid droplets (LDs), and WsATF5 in the LDs. The novel W. subpelliculosus AATases could be usefully applied to produce flavor components in various food industries.

Evaluation of Currently Available Laboratory Methods to Detect Terbinafine Resistant Dermatophytes Including a Gradient Strip for Terbinafine, EUCAST Microdilution E.Def 11.0, a Commercial Real-Time PCR Assay, Squalene Epoxidase Sequencing and Whole Genome Sequencing

BACKGROUND

Terbinafine resistance in dermatophytes is an increasing problem worldwide. Several outbreaks of terbinafine-resistant dermatophytosis are currently occurring in India and surrounding countries, and these recent years, European countries have also been affected by this issue. Currently, antifungal susceptibility testing of dermatophytes is not routinely performed in clinical laboratories.

OBJECTIVES

Given the current situation and associated public health concerns, there is an urgent need for accurate and rapid detection of terbinafine resistance in laboratories. Therefore, we evaluated different methods currently available for the detection of terbinafine resistance in dermatophytes.

METHODS

Twenty-eight strains previously identified as T. indotineae/mentagrophytes/interdigitale were concurrently characterised using terbinafine gradient strips (HiMedia), EUCAST E.Def 11.0 microdilution, the DermaGenius resistance PCR assay (PathoNostics), and SQLE sequencing. These four methods were compared to terbinafine resistance characterisation obtained by whole genome sequencing (WGS).

RESULTS

All four evaluated methods were able to detect terbinafine resistant strains either by showing high MICs (> 0.125 μg/mL) or by detecting SQLE substitutions.

CONCLUSIONS

The gradient strips, despite questionable essential agreement with EUCAST E.Def 11.0, can be an easy, fast and cheap method to screen terbinafine resistance among dermatophytes in clinical laboratories. The DermaGenius resistance PCR assay enables rapid detection of the most common substitutions in SQLE associated with terbinafine resistance. However, its inability to precisely determine specific substitutions on SQLE or identify new ones may pose a problem in the future. These limitations can be addressed by using SQLE sequencing or whole genome sequencing (WGS).

Comparative genomic analysis of copepod humoral immunity genes with sex-biased expression in Labidocera rotunda

Studies of innate immune system function in invertebrates have contributed significantly to our understanding of the mammalian innate immune system. However, in-depth research on innate immunity in marine invertebrates remains sparse. We generated the first de novo genome and transcriptome sequences of copepod Labidocera rotunda using Illumina paired-end data and conducted a comparative genome analysis including five crustaceans (four copepods and one branchiopod species). We cataloged the presence of Toll, Imd, JAK/STAT, and JNK pathway components among them and compared them with 17 previously reported diverse arthropod species representative of insects, myriapods, chelicerates, and malacostracans. Our results indicated that copepod Gram-negative binding proteins may function in direct digestion or pathogen killing. The phylogenetic analysis of arthropod TEP and copepod-specific GCGEQ motif patterns suggested that the evolutionary history of copepod TEPs may have diverged from that of other arthropods. We classified the copepod Toll-like receptors identified in our analysis as either vertebrate or protostome types based on their cysteine motifs and the tree built with their Toll/interleukin-1 receptor domains. LrotCrustin, the first copepod AMP, was identified based on the structure of its WAP domain and deep-learning AMP predictors. Gene expression level analysis of L. rotunda innate immunity-related transcripts in each sex showed higher Toll pathway-related expression in male L. rotunda than in females, which may reflect an inverse correlation between allocation of reproductive investment and elevated immune response in males. Taken together, the results of our study provide insight into copepod innate immunity-related gene families and illuminate the evolutionary potential of copepods relative to other crustaceans.

Limited Migration From Physiological Refugia Constrains the Rescue of Native Gastropods Facing an Invasive Predator

Biological invasions have caused the loss of freshwater biodiversity worldwide. The interplay between adaptive responses and demographic characteristics of populations impacted by invasions is expected to be important for their resilience, but the interaction between these factors is poorly understood. The freshwater gastropod Amnicola limosus is native to the Upper St. Lawrence River and distributed along a water calcium concentration gradient within which high-calcium habitats are impacted by an invasive predator fish (Neogobius melanostomus, round goby), whereas low-calcium habitats provide refuges for the gastropods from the invasive predator. Our objectives were to (1) test for adaptation of A. limosus to the invasive predator and the low-calcium habitats, and (2) investigate if migrant gastropods could move from refuge populations to declining invaded populations (i.e., demographic rescue), which could also help maintain genetic diversity through gene flow (i.e., genetic rescue). We conducted a laboratory reciprocal transplant of wild F0 A. limosus sourced from the two habitat types (high calcium/invaded and low calcium/refuge) to measure adult survival and fecundity in home and transplant treatments of water calcium concentration (low/high) and round goby cue (present/absent). We then applied pooled whole-genome sequencing of 12 gastropod populations from across the calcium/invasion gradient. We identified patterns of life-history traits and genetic differentiation across the habitats that are consistent with local adaptation to low-calcium concentrations in refuge populations and to round goby predation in invaded populations. We also detected restricted gene flow from the low-calcium refugia towards high-calcium invaded populations, implying that the potential for demographic and genetic rescue is limited by natural dispersal. Our study highlights the importance of considering the potentially conflicting effects of local adaptation and gene flow for the resilience of populations coping with invasive predators.

The occurrence of positive selection on BicA transporter of Microcystis aeruginosa

The rapid growth of cyanobacteria, particularly Microcystis aeruginosa, poses a significant threat to global water security. The proliferation of toxic Microcystis aeruginosa raises concerns due to its potential harm to human health and socioeconomic impacts. Dense blooms contribute to spatiotemporal inorganic carbon depletion, promoting interest in the roles of carbon-concentrating mechanisms (CCMs) for competitive carbon uptake. Despite the importance of HCO3- transporters, genetic evaluations and functional predictions in M. aeruginosa remain insufficient. In this study, we explored the diversity of HCO3- transporters in the genomes of 46 strains of M. aeruginosa, assessing positive selection for each. Intriguingly, although the Microcystis BicA transporter became a partial gene in 23 out of 46 genomic strains, we observed significant positive sites. Structural analyses, including predicted 2D and 3D models, confirmed the structural conservation of the Microcystis BicA transporter. Our findings suggest that the Microcystis BicA transport likely plays a crucial role in competitive carbon uptake, emphasizing its ecological significance. The ecological function of the Microcystis BicA transport in competitive growth during cyanobacterial blooms raises important questions. Future studies require experimental confirmation to better understand the role of the Microcysits BicA transporter in cyanobacterial blooms dynamics.

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.

CD59 gene: 143 haplotypes of 22,718 nucleotides length by computational phasing in 113 individuals from different ethnicities

BACKGROUND

CD59 deficiency due to rare germline variants in the CD59 gene causes disabilities, ischemic strokes, neuropathy, and hemolysis. CD59 deficiency due to common somatic variants in the PIG-A gene in hematopoietic stem cells causes paroxysmal nocturnal hemoglobinuria. The ISBT database lists one nonsense and three missense germline variants that are associated with the CD59-null phenotype. To analyze the genetic diversity of the CD59 gene, we determined long-range CD59 haplotypes among individuals from different ethnicities.

METHODS

We determined a 22.7 kb genomic fragment of the CD59 gene in 113 individuals using next-generation sequencing (NGS), which covered the whole NM_203330.2 mRNA transcript of 7796 base pairs. Samples came from an FDA reference repository and our Ethiopia study cohorts. The raw genotype data were computationally phased into individual haplotype sequences.

RESULTS

Nucleotide sequencing of the CD59 gene of 226 chromosomes identified 216 positions with single nucleotide variants. Only three haplotypes were observed in homozygous form, which allowed us to assign them unambiguously as experimentally verified CD59 haplotypes. They were also the most frequent haplotypes among both cohorts. An additional 140 haplotypes were imputed computationally.

DISCUSSION

We provided a large set of haplotypes and proposed three verified long-range CD59 reference sequences, based on a population approach, using a generalizable rationale for our choice. Correct long-range haplotypes are useful as template sequences for allele calling in high-throughput NGS and precision medicine approaches, thus enhancing the reliability of clinical diagnostics. Long-range haplotypes can also be used to evaluate the influence of genetic variation on the risk of transfusion reactions or diseases.

Transcriptome and functional analyses of phenotypic plasticity in sea grape Caulerpa okamurae

Caulerpa is a marine green macroalga distinguished by a large single cell with multiple nuclei. It also exhibits remarkable morphological intraspecies variations, in response to diverse environmental types. However, the molecular mechanisms underlying this phenotypic plasticity remain poorly understood. In this work, we compare the transcriptomes of Caulerpa okamurae Weber Bosse, 1897 displaying altered phenotypes of cultivation and natural phenotypes and investigate significantly regulated genes and their biological functions using differential expression analyses. We observe light-harvesting complex upregulation and cellular framework stability downregulation in altered phenotypes compared to the natural phenotypes. Intertidal macrophytes reduce light capture to avoid photodamage and regulate their morphology to protect against wave damage. In contrast, the lower light conditions and the cultivation environment augment light capture and increase a morphology prioritizing light trapping. Moreover, the addition of simulated wave-sweeping stimuli induces a return to the natural morphology under high-light conditions, showing how mechanical stress affects morphological organization in C. okamurae. We provide detailed gene expression patterns in C. okamurae under varying light intensities and water conditions, suggesting a distinct influence on its morphological traits.

Dietary exposure of the water flea Daphnia galeata to microcystin-LR

Harmful substances like the cyanotoxin microcystin-leucine-arginine (MC-LR) are commonly found in eutrophic freshwater environments, posing risks to aquatic organisms. The water flea, Daphnia, is a well-established model organism for environmental toxicology research. Nevertheless, there is currently insufficient research on the genes that respond to MC-LR in Daphnia galeata. This study aimed to gain insights into the notable genes that react significantly to MC-LR. In this study, we generated an extensive RNA-Seq sequences isolated from the D. galeata HK strain, Han River in Korea. This strain was nourished with a diet of the green microalga Chlorella vulgaris and treated with pure MC-LR at a concentration of 36 ug/L. The transcriptome profile in response to the MC-LR treatment was obtained and 336 differentially expressed genes were subjected to Gene Ontology (GO) and euKaryotic Orthologous Groups of proteins analyses. GO enrichment analysis showed that chemical stimulus, amino sugar metabolic and catabolic process, oxidative stress, and detoxification were highly enriched, in reverse, proteolysis and fucosylation were underpresented. Detoxification process related genes such as peroxidase-like, chorion, and thyroid peroxidase-like were enriched for eliminating or neutralizing MC_LR from an organism's body. Furthermore, functional protein classification revealed an upregulation of lipid and inorganic ion transport processes, while amino acid and carbohydrate transport processes were found to be downregulated. These findings offer insights into how organisms respond to ecotoxic stimuli, providing valuable information for understanding adaptation or defense pathways.

Recent advances in genome annotation and synthetic biology for the development of microbial chassis

This article provides an overview of microbial host selection, synthetic biology, genome annotation, metabolic modeling, and computational methods for predicting gene essentiality for developing a microbial chassis. This article focuses on lactic acid bacteria (LAB) as a microbial chassis and strategies for genome annotation of the LAB genome. As a case study, Lactococcus lactis is chosen based on its well-established therapeutic applications such as probiotics and oral vaccine development. In this article, we have delineated the strategies for genome annotations of lactic acid bacteria. These strategies also provide insights into streamlining genome reduction without compromising the functionality of the chassis and the potential for minimal genome chassis development. These insights underscore the potential for the development of efficient and sustainable synthetic biology systems using streamlined microbial chassis with minimal genomes.

Contrasting new and available reference genomes to highlight uncertainties in assemblies and areas for future improvement: an example with monodontid species

BACKGROUND

Reference genomes provide a foundational framework for evolutionary investigations, ecological analysis, and conservation science, yet uncertainties in the assembly of reference genomes are difficult to assess, and by extension rarely quantified. Reference genomes for monodontid cetaceans span a wide spectrum of data types and analytical approaches, providing the context to derive broader insights related to discrepancies and regions of uncertainty in reference genome assembly. We generated three beluga (Delphinapterus leucas) and one narwhal (Monodon monoceros) reference genomes and contrasted these with published chromosomal scale assemblies for each species to quantify discrepancies associated with genome assemblies.

RESULTS

The new reference genomes achieved chromosomal scale assembly using a combination of PacBio long reads, Illumina short reads, and Hi-C scaffolding data. For beluga, we identified discrepancies in the order and orientation of contigs in 2.2-3.7% of the total genome depending on the pairwise comparison of references. In addition, unsupported higher order scaffolding was identified in published reference genomes. In contrast, we estimated 8.2% of the compared narwhal genomes featured discrepancies, with inversions being notably abundant (5.3%). Discrepancies were linked to repetitive elements in both species.

CONCLUSIONS

We provide several new reference genomes for beluga (Delphinapterus leucas), while highlighting potential avenues for improvements. In particular, additional layers of data providing information on ultra-long genomic distances are needed to resolve persistent errors in reference genome construction. The comparative analyses of monodontid reference genomes suggested that the three new reference genomes for beluga are more accurate compared to the currently published reference genome, but that the new narwhal genome is less accurate than one published. We also present a conceptual summary for improving the accuracy of reference genomes with relevance to end-user needs and how they relate to levels of assembly quality and uncertainty.

Integrated genomics and phenotype microarray analysis of Saccharomyces cerevisiae industrial strains for rice wine fermentation and recombinant protein production

The industrial potential of Saccharomyces cerevisiae has extended beyond its traditional use in fermentation to various applications, including recombinant protein production. Herein, comparative genomics was performed with three industrial S. cerevisiae strains and revealed a heterozygous diploid genome for the 98-5 and KSD-YC strains (exploited for rice wine fermentation) and a haploid genome for strain Y2805 (used for recombinant protein production). Phylogenomic analysis indicated that Y2805 was closely associated with the reference strain S288C, whereas KSD-YC and 98-5 were grouped with Asian and European wine strains, respectively. Particularly, a single nucleotide polymorphism (SNP) in FDC1, involved in the biosynthesis of 4-vinylguaiacol (4-VG, a phenolic compound with a clove-like aroma), was found in KSD-YC, consistent with its lack of 4-VG production. Phenotype microarray (PM) analysis showed that KSD-YC and 98-5 displayed broader substrate utilization than S288C and Y2805. The SNPs detected by genome comparison were mapped to the genes responsible for the observed phenotypic differences. In addition, detailed information on the structural organization of Y2805 selection markers was validated by Sanger sequencing. Integrated genomics and PM analysis elucidated the evolutionary history and genetic diversity of industrial S. cerevisiae strains, providing a platform to improve fermentation processes and genetic manipulation.

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.

Welcome to the big leaves: Best practices for improving genome annotation in non-model plant genomes

Premise

Robust standards to evaluate quality and completeness are lacking in eukaryotic structural genome annotation, as genome annotation software is developed using model organisms and typically lacks benchmarking to comprehensively evaluate the quality and accuracy of the final predictions. The annotation of plant genomes is particularly challenging due to their large sizes, abundant transposable elements, and variable ploidies. This study investigates the impact of genome quality, complexity, sequence read input, and method on protein-coding gene predictions.

Methods

The impact of repeat masking, long-read and short-read inputs, and de novo and genome-guided protein evidence was examined in the context of the popular BRAKER and MAKER workflows for five plant genomes. The annotations were benchmarked for structural traits and sequence similarity.

Results

Benchmarks that reflect gene structures, reciprocal similarity search alignments, and mono-exonic/multi-exonic gene counts provide a more complete view of annotation accuracy. Transcripts derived from RNA-read alignments alone are not sufficient for genome annotation. Gene prediction workflows that combine evidence-based and ab initio approaches are recommended, and a combination of short and long reads can improve genome annotation. Adding protein evidence from de novo assemblies, genome-guided transcriptome assemblies, or full-length proteins from OrthoDB generates more putative false positives as implemented in the current workflows. Post-processing with functional and structural filters is highly recommended.

Discussion

While the annotation of non-model plant genomes remains complex, this study provides recommendations for inputs and methodological approaches. We discuss a set of best practices to generate an optimal plant genome annotation and present a more robust set of metrics to evaluate the resulting predictions.

South Africa's indigenous microbial diversity for industrial applications: A review of the current status and opportunities

The unique metagenomic, metaviromic libraries and indigenous micro diversity within Southern Africa have the potential for global beneficiation in academia and industry. Microorganisms that flourish at high temperatures, adverse pH conditions, and high salinity are likely to have enzyme systems that function efficiently under those conditions. These attributes afford researchers and industries alternative approaches that could replace existing chemical processes. Thus, a better understanding of African microbial/genetic diversity is crucial for the development of "greener" industries. A concerted drive to exploit the potential locked in biological resources has been previously seen with companies such as Diversa Incorporated and Verenium (Badische Anilin-und SodaFabrik-BASF) both building business models that pioneered the production of high-performance specialty enzymes for a variety of different industrial applications. The market potential and accompanying industry offerings have not been fully exploited in South Africa, nor in Africa at large. Utilization of the continent's indigenous microbial repositories could create long-lasting, sustainable growth in various production sectors, providing economic growth in resource-poor regions. By bolstering local manufacture of high-value bio-based products, scientific and engineering discoveries have the potential to generate new industries which in turn would provide employment avenues for many skilled and unskilled laborers. The positive implications of this could play a role in altering the face of business markets on the continent from costly import-driven markets to income-generating export markets. This review focuses on identifying microbially diverse areas located in South Africa while providing a profile for all associated microbial/genetically derived libraries in this country. A comprehensive list of all the relevant researchers and potential key players is presented, mapping out existing research networks for the facilitation of collaboration. The overall aim of this review is to facilitate a coordinated journey of exploration, one which will hopefully realize the value that South Africa's microbial diversity has to offer.

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.

Leveraging massively parallel reporter assays for evolutionary questions

A long-standing goal of evolutionary biology is to decode how gene regulation contributes to organismal diversity. Doing so is challenging because it is hard to predict function from non-coding sequence and to perform molecular research with non-model taxa. Massively parallel reporter assays (MPRAs) enable the testing of thousands to millions of sequences for regulatory activity simultaneously. Here, we discuss the execution, advantages, and limitations of MPRAs, with a focus on evolutionary questions. We propose solutions for extending MPRAs to rare taxa and those with limited genomic resources, and we underscore MPRA's broad potential for driving genome-scale, functional studies across organisms.

A Bioinformatic Guide to Identify Protein Effectors from Phytopathogens

Phytopathogenic fungi are a diverse and widespread group that has a significant detrimental impact on crops with an estimated annual average loss of 15% worldwide. Understanding the interaction between host plants and pathogenic fungi is critical to delineate underlying mechanisms of plant defense to mitigate agricultural losses. Fungal pathogens utilize suites of secreted molecules, called effectors, to modulate plant metabolism and immune response to overcome host defenses and promote colonization. Effectors come in many flavors including proteinaceous products, small RNAs, and metabolites such as mycotoxins. This review will focus on methods for identifying protein effectors from fungi. Excellent reviews have been published to identify secondary metabolites and small RNAs from fungi and therefore will not be part of this review.

Identifying Fungal Secondary Metabolites and Their Role in Plant Pathogenesis

Pathogenic fungi are the main infectious agents of plants. Secondary metabolites produced by these fungi, also recognized as natural products, are key mediators of plant-fungal interactions. Knowledge on the biosynthesis of these metabolites, the accessibility to fungal genome sequences, and the development of gene disruption techniques open up opportunities to identify many more of these metabolites both in vitro and in planta. This methodology chapter gives a detailed systematic approach aiming to discover new natural products from phytopathogenic fungi and characterize their role in triggering plant cell death and plant disease. This approach takes advantage of the global regulation of fungal secondary metabolite production by regulatory proteins reported in various fungal species.

Propagation, detection and correction of errors using the sequence database network

Nucleotide and protein sequences stored in public databases are the cornerstone of many bioinformatics analyses. The records containing these sequences are prone to a wide range of errors, including incorrect functional annotation, sequence contamination and taxonomic misclassification. One source of information that can help to detect errors are the strong interdependency between records. Novel sequences in one database draw their annotations from existing records, may generate new records in multiple other locations and will have varying degrees of similarity with existing records across a range of attributes. A network perspective of these relationships between sequence records, within and across databases, offers new opportunities to detect-or even correct-erroneous entries and more broadly to make inferences about record quality. Here, we describe this novel perspective of sequence database records as a rich network, which we call the sequence database network, and illustrate the opportunities this perspective offers for quantification of database quality and detection of spurious entries. We provide an overview of the relevant databases and describe how the interdependencies between sequence records across these databases can be exploited by network analyses. We review the process of sequence annotation and provide a classification of sources of error, highlighting propagation as a major source. We illustrate the value of a network perspective through three case studies that use network analysis to detect errors, and explore the quality and quantity of critical relationships that would inform such network analyses. This systematic description of a network perspective of sequence database records provides a novel direction to combat the proliferation of errors within these critical bioinformatics resources.

Development of a data science CURE in microbiology using publicly available microbiome datasets

Scientific and technological advances within the life sciences have enabled the generation of very large datasets that must be processed, stored, and managed computationally. Researchers increasingly require data science skills to work with these datasets at scale in order to convert information into actionable insights, and undergraduate educators have started to adapt pedagogies to fulfill this need. Course-based undergraduate research experiences (CUREs) have emerged as a leading model for providing large numbers of students with authentic research experiences including data science. Originally designed around wet-lab research experiences, CURE models have proliferated and diversified globally to accommodate a broad range of academic disciplines. Within microbiology, diversity metrics derived from microbiome sequence information have become standard data products in research. In some cases, researchers have deposited data in publicly accessible repositories, providing opportunities for reproducibility and comparative analysis. In 2020, with the onset of the COVID-19 pandemic and concomitant shift to remote learning, the University of British Columbia set out to develop an online data science CURE in microbiology. A team of faculty with collective domain expertise in microbiome research and CUREs developed and implemented a data science CURE in which teams of students learn to work with large publicly available datasets, develop and execute a novel scientific research project, and disseminate their findings in the online Undergraduate Journal of Experimental Microbiology and Immunology. Analysis of the resulting student-authored research articles, including comments from peer reviews conducted by subject matter experts, demonstrate high levels of learning effectiveness. Here, we describe core insights from course development and implementation based on a reverse course design model. Our approach to course design may be applicable to the development of other data science CUREs.

Microbes of traditional fermentation processes as synthetic biology chassis to tackle future food challenges

Microbial diversity is magnificent and essential to almost all life on Earth. Microbes are an essential part of every human, allowing us to utilize otherwise inaccessible resources. It is no surprise that humans started, initially unconsciously, domesticating microbes for food production: one may call this microbial domestication 1.0. Sourdough bread is just one of the miracles performed by microbial fermentation, allowing extraction of more nutrients from flour and at the same time creating a fluffy and delicious loaf. There are a broad range of products the production of which requires fermentation such as chocolate, cheese, coffee and vinegar. Eventually, with the rise of microscopy, humans became aware of microbial life. Today our knowledge and technological advances allow us to genetically engineer microbes - one may call this microbial domestication 2.0. Synthetic biology and microbial chassis adaptation allow us to tackle current and future food challenges. One of the most apparent challenges is the limited space on Earth available for agriculture and its major tolls on the environment through use of pesticides and the replacement of ecosystems with monocultures. Further challenges include transport and packaging, exacerbated by the 24/7 on-demand mentality of many customers. Synthetic biology already tackles multiple food challenges and will be able to tackle many future food challenges. In this perspective article, we highlight recent microbial synthetic biology research to address future food challenges. We further give a perspective on how synthetic biology tools may teach old microbes new tricks, and what standardized microbial domestication could look like.

Chloroplast Genome Annotation Tools: Prolegomena to the Identification of Inverted Repeats

The development of next-generation sequencing technology and the increasing amount of sequencing data have brought the bioinformatic tools used in genome assembly into focus. The final step of the process is genome annotation, which works on assembled genome sequences to identify the location of genome features. In the case of organelle genomes, specialized annotation tools are used to identify organelle genes and structural features. Numerous annotation tools target chloroplast sequences. Most chloroplast DNA genomes have a quadripartite structure caused by two copies of a large inverted repeat. We investigated the strategies of six annotation tools (Chloë, Chloroplot, GeSeq, ORG.Annotate, PGA, Plann) for identifying inverted repeats and analyzed their success using publicly available complete chloroplast sequences of taxa belonging to the asterid and rosid clades. The annotation tools use two different approaches to identify inverted repeats, using existing general search tools or implementing stand-alone solutions. The chloroplast sequences studied show that there are different types of imperfections in the assembled data and that each tool performs better on some sequences than the others.

Plant Proteome Dynamics

Proteins are intimately involved in executing and controlling virtually all cellular processes. To understand the molecular mechanisms that underlie plant phenotypes, it is essential to investigate protein expression, interactions, and modifications, to name a few. The proteome is highly dynamic in time and space, and a plethora of protein modifications, protein interactions, and network constellations are at play under specific conditions and developmental stages. Analysis of proteomes aims to characterize the entire protein complement of a particular cell type, tissue, or organism-a challenging task, given the dynamic nature of the proteome. Modern mass spectrometry-based proteomics technology can be used to address this complexity at a system-wide scale by the global identification and quantification of thousands of proteins. In this review, we present current methods and technologies employed in mass spectrometry-based proteomics and provide examples of dynamic changes in the plant proteome elucidated by proteomic approaches.

Applications of cell- and tissue-specific 'omics to improve plant productivity

The individual tissues and cell types of plants each have characteristic properties that contribute to the function of the plant as a whole. These are reflected by unique patterns of gene expression, protein and metabolite content, which enable cell-type-specific patterns of growth, development and physiology. Gene regulatory networks act within the cell types to govern the production and activity of these components. For the broader organism to grow and reproduce successfully, cell-type-specific activity must also function within the context of surrounding cell types, which is achieved by coordination of signalling pathways. We can investigate how gene regulatory networks are constructed and function using integrative 'omics technologies. Historically such experiments in plant biological research have been performed at the bulk tissue level, to organ resolution at best. In this review, we describe recent advances in cell- and tissue-specific 'omics technologies that allow investigation at much improved resolution. We discuss the advantages of these approaches for fundamental and translational plant biology, illustrated through the examples of specialised metabolism in medicinal plants and seed germination. We also discuss the challenges that must be overcome for such approaches to be adopted widely by the community.

The first genome for the Cape Primrose Streptocarpus rexii (Gesneriaceae), a model plant for studying meristem-driven shoot diversity

Cape Primroses (Streptocarpus, Gesneriaceae) are an ideal study system for investigating the genetics underlying species diversity in angiosperms. Streptocarpus rexii has served as a model species for plant developmental research for over five decades due to its unusual extended meristem activity present in the leaves. In this study, we sequenced and assembled the complete nuclear, chloroplast, and mitochondrial genomes of S. rexii using Oxford Nanopore Technologies long read sequencing. Two flow cells of PromethION sequencing resulted in 32 billion reads and were sufficient to generate a draft assembly including the chloroplast, mitochondrial and nuclear genomes, spanning 776 Mbp. The final nuclear genome assembly contained 5,855 contigs, spanning 766 Mbp of the 929-Mbp haploid genome with an N50 of 3.7 Mbp and an L50 of 57 contigs. Over 70% of the draft genome was identified as repeats. A genome repeat library of Gesneriaceae was generated and used for genome annotation, with a total of 45,045 genes annotated in the S. rexii genome. Ks plots of the paranomes suggested a recent whole genome duplication event, shared between S. rexii and Primulina huaijiensis. A new chloroplast and mitochondrial genome assembly method, based on contig coverage and identification, was developed, and successfully used to assemble both organellar genomes of S. rexii. This method was developed into a pipeline and proved widely applicable. The nuclear genome of S. rexii and other datasets generated and reported here will be invaluable resources for further research to aid in the identification of genes involved in morphological variation underpinning plant diversification.

Zooplankton diversity monitoring strategy for the urban coastal region using metabarcoding analysis

Marine ecosystems in urban coastal areas are exposed to many risks due to human activity. Thus, long-term and continuous monitoring of zooplankton diversity is necessary. High-throughput DNA metabarcoding has gained recognition as an efficient and highly sensitive approach to accurately describing the species diversity of marine zooplankton assemblages. In this study, we collected 30 zooplankton samples at about 2-week intervals for 1 year. Zooplankton diversity showing a typical four season pattern. Of the “total” and “common” zooplankton, we assigned 267 and 64 taxa. The cluster structure and seasonal diversity pattern were rough when only the “common” zooplankton was used. Our study examined how to maximize the benefits of metabarcoding for monitoring zooplankton diversity in urban coastal areas. The results suggest that to take full advantage of metabarcoding when monitoring a zooplankton community, it is necessary to carefully investigate potential ecosystem threats (non-indigenous species) through sufficient curation rather than disregarding low-abundance operational taxonomic units.

CStone: A de novo transcriptome assembler for short-read data that identifies non-chimeric contigs based on underlying graph structure

With the exponential growth of sequence information stored over the last decade, including that of de novo assembled contigs from RNA-Seq experiments, quantification of chimeric sequences has become essential when assembling read data. In transcriptomics, de novo assembled chimeras can closely resemble underlying transcripts, but patterns such as those seen between co-evolving sites, or mapped read counts, become obscured. We have created a de Bruijn based de novo assembler for RNA-Seq data that utilizes a classification system to describe the complexity of underlying graphs from which contigs are created. Each contig is labelled with one of three levels, indicating whether or not ambiguous paths exist. A by-product of this is information on the range of complexity of the underlying gene families present. As a demonstration of CStones ability to assemble high-quality contigs, and to label them in this manner, both simulated and real data were used. For simulated data, ten million read pairs were generated from cDNA libraries representing four species, Drosophila melanogaster, Panthera pardus, Rattus norvegicus and Serinus canaria. These were assembled using CStone, Trinity and rnaSPAdes; the latter two being high-quality, well established, de novo assembers. For real data, two RNA-Seq datasets, each consisting of ≈30 million read pairs, representing two adult D. melanogaster whole-body samples were used. The contigs that CStone produced were comparable in quality to those of Trinity and rnaSPAdes in terms of length, sequence identity of aligned regions and the range of cDNA transcripts represented, whilst providing additional information on chimerism. Here we describe the details of CStones assembly and classification process, and propose that similar classification systems can be incorporated into other de novo assembly tools. Within a related side study, we explore the effects that chimera’s within reference sets have on the identification of differentially expression genes. CStone is available at: https://sourceforge.net/projects/cstone/.

Within transcriptome reference sets, non-chimeric sequences are representations of transcribed genes, while artificially generated chimeric ones are mosaics of two or more pieces of DNA incorrectly pieced together. One area where such sets are utilized is in the quantification of gene expression patterns; where RNA-Seq reads are mapped to the sequences within, and subsequent count values reflect expression levels. Artificial chimeras can have a negative impact on count values by erroneously increasing variation in relation to the reads being mapped. Reference sets can be created from de novo assembled contigs, but chimeras can be introduced during the assembly process via the required traversal of graphs, representing gene families, constructed from the RNA-Seq data. Graph complexity determines how likely chimeras will arise. We have created CStone, a de novo assembler that utilizes a classification system to describe such complexity. Contigs created by CStone are labelled in a manner that indicates whether or not they are non-chimeric. This encourages contig dependent results to be presented with increased objectivity by maintaining the context of ambiguity associated with the assembly process. CStone has been tested extensively. Additionally, we have quantified the relationship between chimeras within reference sets and the identification of differentially expressed genes.

Importance of Thai macaque bioresources for biological research and human health

During the past century, macaque bioresources have provided remarkable scientific and biomedical discoveries related to the understanding of human physiology, neuroanatomy, reproduction, development, cognition, and pathology. Considerable progress has been made, and an urgent need has arisen to develop infrastructure and viable settings to meet the current global demand in research models during the so-called new normal after COVID-19 era. This review highlights the critical need for macaque bioresources and proposes the establishment of a designated primate research center to integrate research in primate laboratories for the rescue and rehabilitation of wild macaques. Key areas where macaque models have been and continue to be essential for advancing fundamental knowledge in biomedical and biological research are outlined. Detailed genetic studies on macaque bioresources of Thai origin can further facilitate the rapid pace of vaccine discovery.

Whole-Genome Sequencing Analysis to Identify Infection with Multiple Species of Nontuberculous Mycobacteria

Mixed infection with multiple species of nontuberculous mycobacteria (NTM) is difficult to identify and to treat. Current conventional molecular-based methods for identifying mixed infections are limited due to low specificity. Here, we evaluated the utility of whole-genome sequencing (WGS) analysis to detect and identify mixed NTM infections. Analytical tools used included PubMLST, MetaPhlAn3, Kraken2, Mykrobe-Predictor and analysis of heterozygous SNP frequencies. The ability of each to identify mixed infections of NTM species was compared. Sensitivity was tested using 101 samples (sequence sets) including 100 in-silico simulated mixed samples with various proportions of known NTM species and one sample of known mixed NTM species from a public database. Single-species NTM control samples (155 WGS samples from public databases and 15 samples from simulated reads) were tested for specificity. Kraken2 exhibited 100% sensitivity and 98.23% specificity for detection and identification of mixed NTM species with accurate estimation of relative abundance of each species in the mixture. PubMLST (99% and 96.47%) and MetaPhlAn3 (95.04% and 83.52%) had slightly lower sensitivity and specificity. Mykrobe-Predictor had the lowest sensitivity (57.42%). Analysis of read frequencies supporting single nucleotide polymorphisms (SNPs) could not detect mixed NTM samples. Clinical NTM samples (n = 16), suspected on the basis of a 16S-23S rRNA gene sequence-based line-probe assay (LPA) to contain more than one NTM species, were investigated using WGS-analysis tools. This identified only a small proportion (37.5%, 6/16 samples) of the samples as mixed infections and exhibited only partial agreement with LPA results. LPAs seem to be inadequate for detecting mixed NTM species infection. This study demonstrated that WGS-analysis tools can be used for diagnosis of mixed infections with different species of NTM.