Search

MitoFish, MitoAnnotator, and MiFish Pipeline are comprehensive databases of fish mitochondrial genomes (mitogenomes), accurate annotation software of fish mitogenomes, and a web platform for metabarcoding analysis of fish mitochondrial environmental DNA (eDNA), respectively. The MitoFish Suite currently receives over 48,000 visits worldwide every year; however, the performance and usefulness of the online platforms can still be improved. Here, we present essential updates on these platforms, including an enrichment of the reference data sets, an enhanced searching function, substantially faster genome annotation and eDNA analysis with the denoising of sequencing errors, and a multisample comparative analysis function. These updates have made our platform more intuitive, effective, and reliable. These updated platforms are freely available at http://mitofish.aori.u-tokyo.ac.jp/. [ABSTRACT FROM AUTHOR]

Copyright of Food Science & Technology & Economy is the property of Grain Science & Technology & Economy Editorial Office and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Here we report the draft genome sequence of the Latilactobacillus curvatus PN39MY strain. The strain was isolated from Lithuanian traditionally (homemade) fermented cucumber. The genome consisted of 83 contigs with a total size of 1,899,018 bp, an N50 of 40562 and a GC% of 42.1. After sequence trimming, 83 contigs were annotated and 1910 genes were coding sequences. The average nucleotide identity (ANI) between PN39MY and Latilactobacillus curvatus_ZJUNIT8 was 99.45% identifying the strain as Latilactobacillus curvatus. No genes related to antimicrobial resistance or virulence factors were found. The data presented here can be used in comparative genomics to identify antimicrobial resistant genes, plasmids and/or virulence factors that may be present in related Latilactobacillus species. The draft genome sequence data was deposited at NCBI under Bioproject with the accession number PRJNA941180.

Abstract 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.

Abstract Since the entry into genome‐enabled biology several decades ago, much progress has been made in determining, describing, and disseminating the functions of genes and their products. Yet, this information is still difficult to access for many scientists and for most genomes. To provide easy access and a graphical summary of the status of genome function annotation for model organisms and bioenergy and food crop species, we created a web application (https://genomeannotation.rheelab.org) to visualize, search, and download genome annotation data for 28 species. The summary graphics and data tables will be updated semi‐annually, and snapshots will be archived to provide a historical record of the progress of genome function annotation efforts. Clear and simple visualization of up‐to‐date genome function annotation status, including the extent of what is unknown, will help address the grand challenge of elucidating the functions of all genes in organisms.

Compared to nuclear genomes, mitochondrial genomes (mitogenomes) are small and usually code for only a few dozen genes. Still, identifying genes and their structure can be challenging and time-consuming. Even automated tools for mitochondrial genome annotation often require manual analysis and curation by skilled experts. The most difficult steps are (i) the structural modelling of intron-containing genes; (ii) the identification and delineation of Group I and II introns; and (iii) the identification of moderately conserved, non-coding RNA (ncRNA) genes specifying 5S rRNAs, tmRNAs and RNase P RNAs. Additional challenges arise through genetic code evolution which can redefine the translational identity of both start and stop codons, thus obscuring protein-coding genes. Further, RNA editing can render gene identification difficult, if not impossible, without additional RNA sequence data. Current automated mito- and plastid-genome annotators are limited as they are typically tailored to specific eukaryotic groups. The MFannot annotator we developed is unique in its applicability to a broad taxonomic scope, its accuracy in gene model inference, and its capabilities in intron identification and classification. The pipeline leverages curated profile Hidden Markov Models (HMMs), covariance (CMs) and ERPIN models to better capture evolutionarily conserved signatures in the primary sequence (HMMs and CMs) as well as secondary structure (CMs and ERPIN). Here we formally describe MFannot, which has been available as a web-accessible service (https://megasun.bch.umontreal.ca/apps/mfannot/) to the research community for nearly 16 years. Further, we report its performance on particularly intron-rich mitogenomes and describe ongoing and future developments.

ABSTRACT Lactiplantibacillus plantarum (previously known as Lactobacillus plantarum) is increasingly used as a probiotic to treat human diseases, but its phages in the human gut remain unexplored. Here, we report its first gut phage, Gut-P1, which we systematically screened using metagenomic sequencing, virus-like particle (VLP) sequencing, and enrichment culture from 35 fecal samples. Gut-P1 is virulent, belongs to the Douglaswolinvirus genus, and is highly prevalent in the gut (~11% prevalence); it has a genome of 79,928 bp consisting of 125 protein coding genes and displaying low sequence similarities to public L. plantarum phages. Physiochemical characterization shows that it has a short latent period and adapts to broad ranges of temperatures and pHs. Furthermore, Gut-P1 strongly inhibits the growth of L. plantarum strains at a multiplicity of infection (MOI) of 1e−6. Together, these results indicate that Gut-P1 can greatly impede the application of L. plantarum in humans. Strikingly, Gut-P1 was identified only in the enrichment culture, not in our metagenomic or VLP sequencing data nor in any public human phage databases, indicating the inefficiency of bulk sequencing in recovering low-abundance but highly prevalent phages and pointing to the unexplored hidden diversity of the human gut virome despite recent large-scale sequencing and bioinformatics efforts. IMPORTANCE As Lactiplantibacillus plantarum (previously known as Lactobacillus plantarum) is increasingly used as a probiotic to treat human gut-related diseases, its bacteriophages may pose a certain threat to their further application and should be identified and characterized more often from the human intestine. Here, we isolated and identified the first gut L. plantarum phage that is prevalent in a Chinese population. This phage, Gut-P1, is virulent and can strongly inhibit the growth of multiple L. plantarum strains at low MOIs. Our results also show that bulk sequencing is inefficient at recovering low-abundance but highly prevalent phages such as Gut-P1, suggesting that the hidden diversity of human enteroviruses has not yet been explored. Our results call for innovative approaches to isolate and identify intestinal phages from the human gut and to rethink our current understanding of the enterovirus, particularly its underestimated diversity and overestimated individual specificity.

Two novel mycobacteriophages (Prann and LeoAvram) belonging to the family Siphoviridae were isolated from soil samples of Northern India. Genomic DNA of both the phages was extracted, and further sequenced using Illumina technology. Complete genome annotation of both the isolates was performed using DNA Master. Prann and LeoAvram had linear genomes of 68398bp and 47079bp, respectively, with G+C contents of 60-70%. A total of 99 and 75 ORFs were predicted in Prann and LeoAvram, respectively. Based on sequence similarity to known phage proteins, functions were assigned to 44 and 53 genes, respectively. These proteins could be classified into five major groups, viz., phage structural proteins, proteins for recombination, lytic enzymes, proteins involved in DNA / RNA metabolism, and in regulation. Mycobacterium smegmatis was used in this work as a safe surrogate for Mycobacterium tuberculosis, the causative agent for tuberculosis, a major infectious disease worldwide with developing antibiotic resistance. This is the first report of M. smegmatis phages from Northern India.

A potential strain, Paenibacillus sp. JNUCC32, was isolated and subjected to whole-genome sequencing. Genome functional annotation revealed its active metabolic capabilities. This study aimed to investigate the pivotal secondary metabolites in the biological system. Fermentation and extraction were performed, resulting in the isolation of seven known compounds: tryptophol (1), 3-(4-hydroxyphenyl)propionic acid (2), ferulic acid (3), maculosin (4), brevianamide F (5), indole-3-acetic acid (6), and butyric acid (7). Tryptophol exhibited favorable pharmacokinetic properties and demonstrated certain tyrosinase inhibitory activity (IC50 = 999 μM). For further analysis of its inhibition mechanism through molecular docking and molecular dynamics (MD) simulation, tryptophol formed three hydrogen bonds and a pro-Michaelis complex with tyrosinase (binding energy = −5.3 kcal/mol). The MD simulation indicated favorable stability for the tryptophol–mushroom tyrosinase complex, primarily governed by hydrogen bond interactions. The crucial residues VAL-283 and HIS-263 in the docking were also validated. This study suggests tryptophol as a potential candidate for antibrowning agents and dermatological research.

Objectives: Citrus yellow vein clearing virus (CYVCV) is an emerging disease that poses a significant threat to the citrus industry in California. In this study, the viral genomic RNA was isolated from Eureka lemon plants in the greenhouse exhibiting CYVCV symptoms. Subsequently, the corresponding DNA genome amplicon was sequenced and annotated. These efforts expand the genotype database of CYVCV, which aims to enhance detection assays, promote understanding of the virus's genetics and evolution, and support the management of this disease. Data description: In this report, we present the complete genome sequence of the CYVCV California isolate (CA1). The genome was found to be 7,530 bp in length, with a G + C content of 51.7%. The 5′ and 3′ termini were determined using 5′ and 3′ termini rapid amplification of cDNA ends (RACE) systems. Furthermore, our analysis revealed the presence of six open reading frames (ORFs) potentially encoding proteins. All sequence data and annotation have been deposited in GenBank under the accession number OR037276.1. [ABSTRACT FROM AUTHOR]

Accurate annotation of plant genomes remains complex due to the presence of many pseudogenes arising from whole-genome duplication-generated redundancy or the capture and movement of gene fragments by transposable elements. Machine learning on genome-wide epigenetic marks, informed by transcriptomic and proteomic training data, could be used to improve annotations through classification of all putative protein-coding genes as either constitutively silent or able to be expressed. Expressed genes were subclassified as able to express both mRNAs and proteins or only RNAs, and CG gene body methylation was associated only with the former subclass. More than 60,000 protein-coding genes have been annotated in the reference genome of maize inbred B73. About two-thirds of these genes are transcribed and are designated the filtered gene set (FGS). Classification of genes by our trained random forest algorithm was accurate and relied only on histone modifications or DNA methylation patterns within the gene body; promoter methylation was unimportant. Other inbred lines are known to transcribe significantly different sets of genes, indicating that the FGS is specific to B73. We accurately classified the sets of transcribed genes in additional inbred lines, arising from inbred-specific DNA methylation patterns. This approach highlights the potential of using chromatin information to improve annotations of functional genes.

A high-quality genome annotation greatly facilitates successful cell line engineering. Standard draft genome annotation pipelines are based largely on de novo gene prediction, homology, and RNA-Seq data. However, draft annotations can suffer from incorrect predictions of translated sequence, inaccurate splice isoforms, and missing genes. Here, we generated a draft annotation for the newly assembled Chinese hamster genome and used RNA-Seq, proteomics, and Ribo-Seq to experimentally annotate the genome. We identified 3529 new proteins compared to the hamster RefSeq protein annotation and 2256 novel translational events (e.g., alternative splices, mutations, and novel splices). Finally, we used this pipeline to identify the source of translated retroviruses contaminating recombinant products from Chinese hamster ovary (CHO) cell lines, including 119 type-C retroviruses, thus enabling future efforts to eliminate retroviruses to reduce the costs incurred with retroviral particle clearance. In summary, the improved annotation provides a more accurate resource for CHO cell line engineering, by facilitating the interpretation of omics data, defining of cellular pathways, and engineering of complex phenotypes.

Cowpea (Vigna unguiculata [L.] Walp.) is a major crop for worldwide food and nutritional security, especially in sub-Saharan Africa, that is resilient to hot and drought-prone environments. An assembly of the single-haplotype inbred genome of cowpea IT97K-499-35 was developed by exploiting the synergies between single-molecule real-time sequencing, optical and genetic mapping, and an assembly reconciliation algorithm. A total of 519 Mb is included in the assembled sequences. Nearly half of the assembled sequence is composed of repetitive elements, which are enriched within recombination-poor pericentromeric regions. A comparative analysis of these elements suggests that genome size differences between Vigna species are mainly attributable to changes in the amount of Gypsy retrotransposons. Conversely, genes are more abundant in more distal, high-recombination regions of the chromosomes; there appears to be more duplication of genes within the NBS-LRR and the SAUR-like auxin superfamilies compared with other warm-season legumes that have been sequenced. A surprising outcome is the identification of an inversion of 4.2 Mb among landraces and cultivars, which includes a gene that has been associated in other plants with interactions with the parasitic weed Striga gesnerioides. The genome sequence facilitated the identification of a putative syntelog for multiple organ gigantism in legumes. A revised numbering system has been adopted for cowpea chromosomes based on synteny with common bean (Phaseolus vulgaris). An estimate of nuclear genome size of 640.6 Mbp based on cytometry is presented.

Objective: To study phage applications in the control of bacteria contaminations, Pseudomonas aeruginosa phage K4 was characterized for its biological characteristics and inhibition effects of bacterial growth. Methods: One-step growth curve, phage stability, genome sequencing, comparative genomics analysis, and growth inhibition curve were conducted in this work. Results: One-step growth curve showed that phage K4 had a latent phase of 15 min and a burst size of about 85 PFU/infection center, displaying a strong infectivity. Under different MOIs (multiplicity of infection), phage K4 efficiently inhibited the growth of the host cells. Genomic analysis revealed that phage K4 had a genome of 50358 bp with 77 coding sequences and 1 tRNA-Arg gene. Comparative genomics showed that the genome of phage K4 shared an identity of 95.08% with Pseudomonas virus PA11 which belongs to the Paundecimvirus genus of the Zobellviridae family, suggesting that phage K4 might be a new member of the Paundecimvirus genus. In application assays, phage K4 could significantly control the growth of the host cells in sterilized milk and ready-to-eat beef samples along with obvious reproductions of phage K4. In conclusion, the data suggested that phage K4 was a lytic virus that replicated rapidly with relatively high bactericidal activity. Conclusion: The genome of phage K4 did not carry genes for integration, antibiotic resistance, and virulence, appropriately to be selected as a biocontrol agent in the control and prevention of P. aeruginosa contaminations in foods.

Abstract Background Red clover (Trifolium pratense L.) is a diploid perennial temperate legume with 14 chromosomes (2n = 14) native to Europe and West Asia, with high nutritional and economic value. It is a very important forage grass and is widely grown in marine climates, such as the United States and Sweden. Genetic research and molecular breeding are limited by the lack of high-quality reference genomes. In this study, we used Illumina, PacBio HiFi, and Hi-C to obtain a high-quality chromosome-scale red clover genome and used genome annotation results to analyze evolutionary relationships among related species. Results The red clover genome obtained by PacBio HiFi assembly sequencing was 423 M. The assembly quality was the highest among legume genome assemblies published to date. The contig N50 was 13 Mb, scaffold N50 was 55 Mb, and BUSCO completeness was 97.9%, accounting for 92.8% of the predicted genome. Genome annotation revealed 44,588 gene models with high confidence and 52.81% repetitive elements in red clover genome. Based on a comparison of genome annotation results, red clover was closely related to Trifolium medium and distantly related to Glycine max, Vigna radiata, Medicago truncatula, and Cicer arietinum among legumes. Analyses of gene family expansions and contractions and forward gene selection revealed gene families and genes related to environmental stress resistance and energy metabolism. Conclusions We report a high-quality de novo genome assembly for the red clover at the chromosome level, with a substantial improvement in assembly quality over those of previously published red clover genomes. These annotated gene models can provide an important resource for molecular genetic breeding and legume evolution studies. Furthermore, we analyzed the evolutionary relationships among red clover and closely related species, providing a basis for evolutionary studies of clover leaf and legumes, genomics analyses of forage grass, the improvement of agronomic traits.

Abstract Background An optimal starting point for relating genome function to organismal biology is a high-quality nuclear genome assembly, and long-read sequencing is revolutionizing the production of this genomic resource in insects. Despite this, nuclear genome assemblies have been under-represented for agricultural insect pests, particularly from the order Coleoptera. Here we present a de novo genome assembly and structural annotation for the coconut rhinoceros beetle, Oryctes rhinoceros (Coleoptera: Scarabaeidae), based on Oxford Nanopore Technologies (ONT) long-read data generated from a wild-caught female, as well as the assembly process that also led to the recovery of the complete circular genome assemblies of the beetle’s mitochondrial genome and that of the biocontrol agent, Oryctes rhinoceros nudivirus (OrNV). As an invasive pest of palm trees, O. rhinoceros is undergoing an expansion in its range across the Pacific Islands, requiring new approaches to management that may include strategies facilitated by genome assembly and annotation. Results High-quality DNA isolated from an adult female was used to create four ONT libraries that were sequenced using four MinION flow cells, producing a total of 27.2 Gb of high-quality long-read sequences. We employed an iterative assembly process and polishing with one lane of high-accuracy Illumina reads, obtaining a final size of the assembly of 377.36 Mb that had high contiguity (fragment N50 length = 12 Mb) and accuracy, as evidenced by the exceptionally high completeness of the benchmarked set of conserved single-copy orthologous genes (BUSCO completeness = 99.1%). These quality metrics place our assembly ahead of the published Coleopteran genomes, including that of an insect model, the red flour beetle (Tribolium castaneum). The structural annotation of the nuclear genome assembly contained a highly-accurate set of 16,371 protein-coding genes, with only 2.8% missing BUSCOs, and the expected number of non-coding RNAs. The number and structure of paralogous genes in a gene family like Sigma GST is lower than in another scarab beetle (Onthophagus taurus), but higher than in the red flour beetle (Tribolium castaneum), which suggests expansion of this GST class in Scarabaeidae. The quality of our gene models was also confirmed with the correct placement of O. rhinoceros among other members of the rhinoceros beetles (subfamily Dynastinae) in a phylogeny based on the sequences of 95 protein-coding genes in 373 beetle species from all major lineages of Coleoptera. Finally, we provide a list of 30 candidate dsRNA targets whose orthologs have been experimentally validated as highly effective targets for RNAi-based control of several beetles. Conclusions The genomic resources produced in this study form a foundation for further functional genetic research and management programs that may inform the control and surveillance of O. rhinoceros populations, and we demonstrate the efficacy of de novo genome assembly using long-read ONT data from a single field-caught insect.

ABSTRACT: Pasteurella multocida (P. multocida) is a zoonotic bacterium that can cause diseases in a variety of animals. It was divided into 5 serogroups, and serogroup A is mainly prevalent in avian hosts. We isolated a virulent and multidrug-resistant P. multocida strain from Guangdong duck liver and named it PMWSG-4 (GenBank accession no. CP077723.1). To understand the pathogenicity of this strain, the pathogenicity test was carried out with mice and ducks. The results showed that PMSWG-4 was highly pathogenic to ducks and mice, and the LD50 is 4.5 and 73 CFU, respectively. In order to study its genetic characteristics, pathogenicity, and relationship with the host, we performed a whole genome sequencing. The genome size of the isolated PMWSG-4 was 2.38 Mbp, with a G+C content of 40.3%, and coding 2,313 Coding DNA Sequence (CDS). The genome carries 162 potential virulence-associated genes, 32 different drug resistance phenotypes, 102 genes possibly involved in pathogen–host interaction, 2 gene island groups, and 4 prophages. In addition, we also found a new drug-resistant plasmid from strain PMWSG-4, named pXL001 (GenBank accession no. CP077724.1). After verified, the plasmid is a new plasmid carrying the floR florfenicol resistance gene. The whole genome is of great significance for further studying the pathogenesis and genetic characteristics of duck-derived P. multocida.

Edwardsiella tarda is a Gram-negative, facultative anaerobic rod-shaped bacterium and the causative agent of the systemic disease “Edwardsiellosis”. It is commonly prevalent in aquatic organisms with subsequent economic loss and hence has attracted increasing attention from researchers. In this study, we investigated the complete genome sequence of a highly virulent isolate Edwardsiella tarda SC002 isolated from hatchlings of the Siamese crocodile. The genome of SC002 consisted of one circular chromosome of length 3,662,469 bp with a 57.29% G+C content and four novel plasmids. A total of 3,734 protein-coding genes, 12 genomic islands (GIs), 7 prophages, 48 interspersed repeat sequences, 248 tandem repeat sequences, a CRISPR component with a total length of 175 bp, and 171 ncRNAs (tRNA = 106, sRNA = 37, and rRNA = 28) were predicted. In addition, the coding genes of assembled genome were successfully annotated against eight general databases (NR = 3,618/3,734, COG = 2,947/3,734, KEGG = 3,485/3,734, SWISS-PROT = 2,787/3,734, GO = 2,648/3,734, Pfam = 2,648/3,734, CAZy = 130/3,734, and TCDB = 637/3,734) and four pathogenicity-related databases (ARDB = 11/3,734, CARD = 142/3,734, PHI = 538/3,734, and VFDB = 315/3,734). Pan-genome and comparative genome analyses of the complete sequenced genomes confirmed their evolutionary relationships. The present study confirmed that E. tarda SC002 is a potential pathogen bearing a bulk amount of antibiotic resistance, virulence, and pathogenic genes and its open pan-genome may enhance its host range in the future.

Simple Summary: We identified the neuropeptides and their genomic loci on the draft genome sequences of Gryllus bimaculatus. These annotations were additionally assigned to the draft genome annotation. This addition to the draft genome annotation improved the convenience of research by consolidating the knowledge of neuropeptides, such as the sequence information and functional annotation, which could not be obtained without searching each article, into the draft genome annotation. This contributed to the infrastructure for facilitating genome-wide research using high-throughput sequencing technology. Genome annotation is critically important data that can support research. Draft genome annotations cover representative genes; however, they often do not include genes that are expressed only in limited tissues and stages, or genes with low expression levels. Neuropeptides are responsible for regulation of various physiological and biological processes. A recent study disclosed the genome draft of the two-spotted cricket Gryllus bimaculatus, which was utilized to understand the intriguing physiology and biology of crickets. Thus far, only two of the nine reported neuropeptides in G. bimaculatus were annotated in the draft genome. Even though de novo assembly using transcriptomic analyses can comprehensively identify neuropeptides, this method does not follow those annotations on the genome locus. In this study, we performed the annotations based on the reference mapping, de novo transcriptome assembly, and manual curation. Consequently, we identified 41 neuropeptides out of 43 neuropeptides, which were reported in the insects. Further, 32 of the identified neuropeptides on the genomic loci in G. bimaculatus were annotated. The present annotation methods can be applicable for the neuropeptide annotation of other insects. Furthermore, the methods will help to generate useful infrastructures for studies relevant to neuropeptides. [ABSTRACT FROM AUTHOR]

Douglas-fir (Pseudotsuga menziesii) is native to western North America. It grows in a wide range of environmental conditions and is an important timber tree. Although there are several studies on the gene expression responses of Douglas-fir to abiotic cues, the absence of highquality transcriptome and genome data is a barrier to further investigation. Like for most conifers, the available transcriptome and genome reference dataset for Douglas-fir remains fragmented and requires refinement. We aimed to generate a highly accurate, and complete reference transcriptome and genome annotation. We deep-sequenced the transcriptome of Douglas-fir needles from seedlings that were grown under nonstress control conditions or a combination of heat and drought stress conditions using long-read (LR) and short-read (SR) sequencing platforms. We used 2 computational approaches, namely de novo and genome-guided LR transcriptome assembly. Using the LR de novo assembly, we identified 1.3X more high-quality transcripts, 1.85X more "complete" genes, and 2.7X more functionally annotated genes compared to the genome-guided assembly approach. We predicted 666 long noncoding RNAs and 12,778 unique protein-coding transcripts including 2,016 putative transcription factors. We leveraged the LR de novo assembled transcriptome with pairedend SR and a published single-end SR transcriptome to generate an improved genome annotation. This was conducted with BRAKER2 and refined based on functional annotation, repetitive content, and transcriptome alignment. This high-quality genome annotation has 51,419 unique gene models derived from 322,631 initial predictions. Overall, our informatics approach provides a new reference Douglas-fir transcriptome assembly and genome annotation with considerably improved completeness and functional annotation. [ABSTRACT FROM AUTHOR]

A contiguous assembly of the inbred 'EL10' sugar beet (Beta vulgaris ssp. vulgaris) genome was constructed using PacBio long-read sequencing, BioNano optical mapping, Hi-C scaffolding, and Illumina short-read error correction. The EL10.1 assembly was 540 Mb, of which 96.2% was contained in nine chromosome-sized pseudomolecules with lengths from 52 to 65 Mb, and 31 contigs with a median size of 282 kb that remained unassembled. Gene annotation incorporating RNA-seq data and curated sequences via the MAKER annotation pipeline generated 24,255 gene models. Results indicated that the EL10.1 genome assembly is a contiguous genome assembly highly congruent with the published sugar beet reference genome. Gross duplicate gene analyses of EL10.1 revealed little large-scale intra-genome duplication. Reduced gene copy number for well-annotated gene families relative to other core eudicots was observed, especially for transcription factors. Variation in genome size in B. vulgaris was investigated by flow cytometry among 50 individuals producing estimates from 633 to 875 Mb/1C. Read-depth mapping with short-read whole-genome sequences from other sugar beet germplasm suggested that relatively few regions of the sugar beet genome appeared associated with high-copy number variation. [ABSTRACT FROM AUTHOR]

Background: Melaleuca alternifolia is a commercially important medicinal tea tree native to Australia. Tea tree oil, the essential oil distilled from its branches and leaves, has broad-spectrum germicidal activity and is highly valued in the pharmaceutical and cosmetic industries. Thus, the study of genome, which can provide reference for the investigation of genes involved in terpinen-4-ol biosynthesis, is quite crucial for improving the productivity of Tea tree oil. Methods and results: In our study, the next-generation sequencing was used to investigate the whole genome of Melaleuca alternifolia. About 114 Gb high quality sequence data were obtained and assembled into 1,838,159 scafolds with an N50 length of 1021 bp. The assembled genome size is about 595 Mb, twice of that predicted by flow cytometer (300 Mb) and k-mer analysis (345 Mb). Benchmarking Universal Single-Copy Orthologs analyses indicated that only 11.3% of the conserved single-copy genes were miss. Repetitive regions cover over 40.43% of the genome. A total of 44,369 protein-coding genes were predicted and annotated against Nr, Swissprot, Refseq, COG, KOG, and KEGG database. Among these genes, 32,909 and 16,241 genes were functionally annotated in Nr and KEGG, respectively. Moreover, 29,411 and 14,435 genes were functionally annotated in COG and KOG. Additionally, 457,661 simple sequence repeats and 1109 transcription factors (TFs) form 67 TF families were identified in the assembled genome. Conclusion: Our findings provide a draft genome sequencing of M. alternifolia which can act as a reference for the deep sequencing strategies, and are useful for future functional and comparative genomics analyses. [ABSTRACT FROM AUTHOR]

Background: Massive parallel sequencing technologies have enabled the elucidation of plant phylogenetic relationships from chloroplast genomes at a high pace. These include members of the family Rhamnaceae. The current Rhamnaceae phylogenetic tree is from 13 out of 24 Rhamnaceae chloroplast genomes, and only one chloroplast genome of the genus Ventilago is available. Hence, the phylogenetic relationships in Rhamnaceae remain incomplete, and more representative species are needed. Results: The complete chloroplast genome of Ventilago harmandiana Pierre was outlined using a hybrid assembly of long- and short-read technologies. The accuracy and validity of the final genome were confirmed with PCR amplifications and investigation of coverage depth. Sanger sequencing was used to correct for differences in lengths and nucleotide bases between inverted repeats because of the homopolymers. The phylogenetic trees reconstructed using prevalent methods for phylogenetic inference were topologically similar. The clustering based on codon usage was congruent with the molecular phylogenetic tree. The groups of genera in each tribe were in accordance with tribal classification based on molecular markers. We resolved the phylogenetic relationships among six Hovenia species, three Rhamnus species, and two Ventilago species. Our reconstructed tree provides the most complete and reliable low-level taxonomy to date for the family Rhamnaceae. Similar to other higher plants, the RNA editing mostly resulted in converting serine to leucine. Besides, most genes were subjected to purifying selection. Annotation anomalies, including indel calling errors, unaligned open reading frames of the same gene, inconsistent prediction of intergenic regions, and misannotated genes, were identified in the published chloroplast genomes used in this study. These could be a result of the usual imperfections in computational tools, and/or existing errors in reference genomes. Importantly, these are points of concern with regards to utilizing published chloroplast genomes for comparative genomic analysis. Conclusions: In summary, we successfully demonstrated the use of comprehensive genomic data, including DNA and amino acid sequences, to build a reliable and high-resolution phylogenetic tree for the family Rhamnaceae. Additionally, our study indicates that the revision of genome annotation before comparative genomic analyses is necessary to prevent the propagation of errors and complications in downstream analysis and interpretation. [ABSTRACT FROM AUTHOR]

SUMMARY: The king protea (Protea cynaroides), an early‐diverging eudicot, is the most iconic species from the Megadiverse Cape Floristic Region, and the national flower of South Africa. Perhaps best known for its iconic flower head, Protea is a key genus for the South African horticulture industry and cut‐flower market. Ecologically, the genus and the family Proteaceae are important models for radiation and adaptation, particularly to soils with limited phosphorus bio‐availability. Here, we present a high‐quality chromosome‐scale assembly of the P. cynaroides genome as the first representative of the fynbos biome. We reveal an ancestral whole‐genome duplication event that occurred in the Proteaceae around the late Cretaceous that preceded the divergence of all crown groups within the family and its extant diversity in all Southern continents. The relatively stable genome structure of P. cynaroides is invaluable for comparative studies and for unveiling paleopolyploidy in other groups, such as the distantly related sister group Ranunculales. Comparative genomics in sequenced genomes of the Proteales shows loss of key arbuscular mycorrhizal symbiosis genes likely ancestral to the family, and possibly the order. The P. cynaroides genome empowers new research in plant diversification, horticulture and adaptation, particularly to nutrient‐poor soils. Significance Statement: Our chromosome‐level Protea cynaroides genome assembly is the first endemic South African plant genome to be published, and can signify a new wave of research focusing on indigenous fynbos and endemic African plant species. We foresee this serving as increasing visibility of plant genomics research from the African continent and the Global South, and enriching global biodiversity, sustainability, evolution and conservation research. [ABSTRACT FROM AUTHOR]

Poa annua L. is a globally distributed grass with economic and horticultural significance as a weed and as a turfgrass. This dual significance, and its phenotypic plasticity and ecological adaptation, have made P. annua an intriguing plant for genetic and evolutionary studies. Because of the lack of genomic resources and its allotetraploid (2 n = 4 x = 28) nature, a reference genome sequence would be a valuable asset to better understand the significance and polyploid origin of P. annua. Here we report a genome assembly with scaffolds representing the 14 haploid chromosomes that are 1.78 Gb in length with an N50 of 112 Mb and 96.7% of BUSCO orthologs. Seventy percent of the genome was identified as repetitive elements, 91.0% of which were Copia- or Gypsy-like long-terminal repeats. The genome was annotated with 76,420 genes spanning 13.3% of the 14 chromosomes. The two subgenomes originating from Poa infirma (Knuth) and Poa supina (Schrad) were sufficiently divergent to be distinguishable but syntenic in sequence and annotation with repetitive elements contributing to the expansion of the P. infirma subgenome. [ABSTRACT FROM AUTHOR]

Eukaryotic genomes contain exons and introns, and it is necessary to accurately identify exon-intron boundaries, i.e., splice sites, to annotate genomes. To address this problem, many previous works have proposed annotation methods/tools based on RNA-seq evidence. Many recent works exploit neural networks (NNs) as their prediction models, but only a few can be used to generate new genome annotation in practice. In this study, we propose AtLASS, a fully automated method for predicting splice sites from genomic and RNA-seq data using attention-based Bi-LSTM (Bidirectional Long Short-Term Memory). We exploit two-stage training on RNA-seq data to address the problem of biased label problem, thereby reducing the false positives. The experiments on the genomes of three species show that the performance of the proposed method itself is comparable to that of existing methods, but we can achieve better performance by combining the outputs of the proposed method and the existing method. The proposed method is the first program specialized in end-to-end splice site prediction using NNs. [ABSTRACT FROM AUTHOR]

Apple scab, caused by Venturia inaequalis, is one of the world's most commercially significant apple diseases. The fungi have a catastrophic impact on apples, causing considerable losses in fruit quality and productivity in many apple-growing locations despite numerous control agents. Fungi secrete various effectors and other virulence-associated proteins that suppress or alter the host's immune system, and several such proteins were discovered in this work. Using state-of-the-art bioinformatics techniques, we examined the V. inaequalis reference genome (EU-B04), resulting in the identification of 647 secreted proteins, of which 328 were classified as small secreted proteins (SSPs), with 76.52% of SSPs identified as anticipated effector proteins. The more prevalent CAZyme proteins were the enzymes engaged in plant cell wall disintegration (targeting pectin and xylanase), adhesion and penetration (Cutinases/acetyl xylan esterase), and reactive oxygen species formation (multicopper oxidases). Furthermore, members of the S9 prolyl oligopeptidase family were identified as the most abundant host defense peptidases. Several known effector proteins were discovered to be expressed during the V. inaequalis infection process on apple leaves. The present study provides valuable data that can be used to develop new strategies for controlling apple scab. [ABSTRACT FROM AUTHOR]

Recently developed methods to predict three-dimensional protein structure with high accuracy have opened new avenues for genome and proteome research. We explore a new hypothesis in genome annotation, namely whether computationally predicted structures can help to identify which of multiple possible gene isoforms represents a functional protein product. Guided by protein structure predictions, we evaluated over 230,000 isoforms of human protein-coding genes assembled from over 10,000 RNA sequencing experiments across many human tissues. From this set of assembled transcripts, we identified hundreds of isoforms with more confidently predicted structure and potentially superior function in comparison to canonical isoforms in the latest human gene database. We illustrate our new method with examples where structure provides a guide to function in combination with expression and evolutionary evidence. Additionally, we provide the complete set of structures as a resource to better understand the function of human genes and their isoforms. These results demonstrate the promise of protein structure prediction as a genome annotation tool, allowing us to refine even the most highly curated catalog of human proteins. More generally we demonstrate a practical, structure-guided approach that can be used to enhance the annotation of any genome. [ABSTRACT FROM AUTHOR]

Plant, animal and protist genomes often contain endogenous viral elements (EVEs), which correspond to partial and sometimes entire viral genomes that have been captured in the genome of their host organism through a variety of integration mechanisms. While the number of sequenced eukaryotic genomes is rapidly increasing, the annotation and characterization of EVEs remains largely overlooked. EVEs that derive from members of the family Caulimoviridae are widespread across tracheophyte plants, and sometimes they occur in very high copy numbers. However, existing programs for annotating repetitive DNA elements in plant genomes are poor at identifying and then classifying these EVEs. Other than accurately annotating plant genomes, there is intrinsic value in a tool that could identify caulimovirid EVEs as they testify to recent or ancient host-virus interactions and provide valuable insights into virus evolution. In response to this research need, we have developed CAULIFINDER, an automated and sensitive annotation software package. CAULIFINDER consists of two complementary workflows, one to reconstruct, annotate and group caulimovirid EVEs in a given plant genome and the second to classify these genetic elements into officially recognized or tentative genera in the Caulimoviridae. We have benchmarked the CAULIFINDER package using the Vitis vinifera reference genome, which contains a rich assortment of caulimovirid EVEs that have previously been characterized using manual methods. The CAULIFINDER package is distributed in the form of a Docker image. [ABSTRACT FROM AUTHOR]

Simple Summary: The coral reef-dwelling grouper Epinephelus cyanopodus has huge economic and ecological value. Due to its special reproductive strategy, complex social structure, and classification controversy, this species is a good model to study the coral reef ecosystem and the classification and speciation of groupers. The lack of genomic resources has hampered research into the genetic basis of their biological traits and adaptive evolution. Therefore, we have assembled a high-quality genome of E. cyanopodus and provided insights into the genetic basis of its adaptive evolution and rapid differentiation at the genomic level, as well as a foundation for subsequent studies on mechanisms of speciation, resistance breeding and genetic conservation for this species. Epinephelus cyanopodus is a coral reef-dwelling grouper with important economic and ecological value and is widely distributed in the western Pacific Ocean. The lack of genomic resources for E. cyanopodus hinders its adaptive evolution and phylogeny research. We constructed the first high-quality genome of E. cyanopodus based on DNBSEQ, PacBio, and Hic sequencing technologies, with a genome size of 998.82 Mb, contig N50 of 5.855 Mb, and scaffold N50 of 41.98 Mb. More than 99.7% of contigs were anchored to 24 pseudochromosomes, and 94.2% of BUSCO genes were found in the E. cyanopodus genome, indicating a high genome assembly completeness. A total of 26,337 protein-coding genes were predicted, of which 98.77% were functionally annotated. Phylogenetic analysis showed that E. cyanopodus separated from its closely related species Epinephelus akaara about 11.5–26.5 million years ago, and the uplift of the Indo-Australian archipelago may have provided an opportunity for its rapid radiation. Moreover, several gene families associated with innate and adaptive immunity were significantly expanded in speckled blue grouper compared to other teleost genomes. Additionally, we identified several genes associated with immunity, growth and reproduction that are under positive selection in E. cyanopodus compared to other groupers, suggesting that E. cyanopodus has evolved broad adaptability in response to complex survival environment, which may provide the genetic basis for its rapid radiation. In brief, the high-quality reference genome of the speckled blue grouper provides a foundation for research on its biological traits and adaptive evolution and will be an important genetic tool to guide aquaculture and resolve its taxonomic controversies in future studies. [ABSTRACT FROM AUTHOR]

Background: Chrysanthemum seticuspe has emerged as a model plant species of cultivated chrysanthemums, especially for studies involving diploid and self-compatible pure lines (Gojo-0). Its genome was sequenced and assembled into chromosomes. However, the genome annotation of C. seticuspe still needs to be improved to elucidate the complex regulatory networks in this species. Results: In addition to the 74,259 mRNAs annotated in the C. seticuspe genome, we identified 18,265 novel mRNAs, 51,425 novel lncRNAs, 501 novel miRNAs and 22,065 novel siRNAs. Two C-class genes and YABBY family genes were highly expressed in disc florets, while B-class genes were highly expressed in ray florets. A WGCNA was performed to identify the hub lncRNAs and mRNAs in ray floret- and disc floret-specific modules, and CDM19, BBX22, HTH, HSP70 and several lncRNAs were identified. ceRNA and lncNAT networks related to flower development were also constructed, and we found a latent functional lncNAT–mRNA combination, LXLOC_026470 and MIF2. Conclusions: The annotations of mRNAs, lncRNAs and small RNAs in the C. seticuspe genome have been improved. The expression profiles of flower development-related genes, ceRNA networks and lncNAT networks were identified, laying a foundation for elucidating the regulatory mechanisms underlying disc floret and ray floret formation. [ABSTRACT FROM AUTHOR]

Background: The annotation of genomes is a crucial step regarding the analysis of new genomic data and resulting insights, and this especially for emerging organisms which allow researchers to access unexplored lineages, so as to expand our knowledge of poorly represented taxonomic groups. Complete pipelines for eukaryotic genome annotation have been proposed for more than a decade, but the issue is still challenging. One of the most widely used tools in the field is MAKER2, an annotation pipeline using experimental evidence (mRNA-seq and proteins) and combining different gene prediction tools. MAKER2 enables individual laboratories and small-scale projects to annotate non-model organisms for which pre-existing gene models are not available. The optimal use of MAKER2 requires gathering evidence data (by searching and assembling transcripts, and/or collecting homologous proteins from related organisms), elaborating the best annotation strategy (training of gene models) and efficiently orchestrating the different steps of the software in a grid computing environment, which is tedious, time-consuming and requires a great deal of bioinformatic skills. Methods: To address these issues, we present AMAW (Automated MAKER2 Annotation Wrapper), a wrapper pipeline for MAKER2 that automates the above-mentioned tasks. Importantly, AMAW also exists as a Singularity container recipe easy to deploy on a grid computer, thereby overcoming the tricky installation of MAKER2. Use case: The performance of AMAW is illustrated through the annotation of a selection of 32 protist genomes, for which we compared its annotations with those produced with gene models directly available in AUGUSTUS. Conclusions: Importantly, AMAW also exists as a Singularity container recipe easy to deploy on a grid computer, thereby overcoming the tricky installation of MAKER2

Abstract Background Gongronella butleri is a fungus with many industrial applications including the composting of solid biowaste. Kerala Agricultural University, India, has developed a microbial consortium of which GbKAU strain of G. butleri is a major component. Even with great industrial significance, genome of this fungus is not published, and the genes and pathways contributing to the applications are not understood. This study had the objective to demonstrate the solid biowaste decomposing capability of the strain, to sequence and annotate the genome, and to reveal the genes and pathways contributing to its biodegradation potential. Results Strain GbKAU of G. butleri isolated and purified from the organic compost was found to produce higher levels of laccase and amylase, compared to Bacillus subtilis which is being widely used in biosolid waste management. Both were shown to be equally efficient in the in vivo composting capabilities. Whole genome sequencing has given ~11 million paired-end good quality reads. De novo assembly using dual-fold approach has yielded 44,639 scaffolds with draft genome size of 29.8 Mb. A total of 11,428 genes were predicted and classified into 359 groups involved in diverse pathways, of which 14 belonged to the enzymes involved in the degradation of macromolecules. Seven previously sequenced strains of the fungus were assembled and annotated. A direct comparison showed that the number of genes present in those strains was comparable to our strain, while all the important biodegrading genes were conserved across the genomes. Gene Ontology analysis had classified the genes according to their molecular function, biological process, and cellular component. A total of 104,718 SSRs were mined and classified to mono- to hexa-nucleotide repeats. The variant analysis in comparison with the closely related genus Cunninghamella has revealed 1156 variants. Conclusions Apart from demonstrating the biodegradation capabilities of the GbKAU strain of G. butleri, the genome of this industrially important fungus was sequenced, de novo assembled, and annotated. GO analysis has classified the genes based on their functions, and the genes involved in biodegradation were revealed. Biodegradation potential, genome features in comparison with other strains, and the functions of the identified genes are discussed.

Abstract Background Paspalum notatum ‘Flugge’ is a diploid with 20 chromosomes (2n = 20) multi-purpose subtropical herb native to South America and has a high ecological significance. It is currently widely planted in tropical and subtropical regions. Despite the gene pool of P. notatum ‘Flugge’ being unearthed to a large extent in the past decade, no details about the genomic information of relevant species in Paspalum have been reported. In this study, the complete genome information of P. notatum was established and annotated through sequencing and de novo assembly of its genome. Results The latest PacBio third-generation HiFi assembly and sequencing revealed that the genome size of P. notatum ‘Flugge’ is 541 M. The assembly result is the higher index among the genomes of the gramineous family published so far, with a contig N50 = 52Mbp, scaffold N50 = 49Mbp, and BUSCOs = 98.1%, accounting for 98.5% of the estimated genome. Genome annotation revealed 36,511 high-confidence gene models, thus providing an important resource for future molecular breeding and evolutionary research. A comparison of the genome annotation results of P. notatum ‘Flugge’ with other closely related species revealed that it had a close relationship with Zea mays but not close compared to Brachypodium distachyon, Setaria viridis, Oryza sativa, Puccinellia tenuiflora, Echinochloa crusgalli. An analysis of the expansion and contraction of gene families suggested that P. notatum ‘Flugge’ contains gene families associated with environmental resistance, increased reproductive ability, and molecular evolution, which explained its excellent agronomic traits. Conclusion This study is the first to report the high-quality chromosome-scale-based genome of P. notatum ‘Flugge’ assembled using the latest PacBio third-generation HiFi sequencing reads. The study provides an excellent genetic resource bank for gramineous crops and invaluable perspectives regarding the evolution of gramineous plants.

Abstract Objective Elizabethkingia meningoseptica is a multidrug resistance strain which primarily causes meningitis in neonates and immunocompromised patients. Being a nosocomial infection causing agent, less information is available in literature, specifically, about its genomic makeup and associated features. An attempt is made to study them through bioinformatics tools with respect to compositions, embedded periodicities, open reading frames, origin of replication, phylogeny, orthologous gene clusters analysis and pathways. Results Complete DNA and protein sequence pertaining to E. meningoseptica were thoroughly analyzed as part of the study. E. meningoseptica G4076 genome showed 7593 ORFs it is GC rich. Fourier based analysis showed the presence of typical three base periodicity at the genome level. Putative origin of replication has been identified. Phylogenetically, E. meningoseptica is relatively closer to E. anophelis compared to other Elizabethkingia species. A total of 2606 COGs were shared by all five Elizabethkingia species. Out of 3391 annotated proteins, we could identify 18 unique ones involved in metabolic pathway of E. meningoseptica and this can be an initiation point for drug designing and development. Our study is novel in the aspect in characterizing and analyzing the whole genome data of E. meningoseptica.

Microbial cell factories (MCFs) have been leveraged to construct sustainable platforms for value-added compound production. To optimize metabolism and reach optimal productivity, synthetic biology has developed various genetic devices to engineer microbial systems by gene editing, high-throughput protein engineering, and dynamic regulation. However, current synthetic biology methodologies still rely heavily on manual design, laborious testing, and exhaustive analysis. The emerging interdisciplinary field of artificial intelligence (AI) and biology has become pivotal in addressing the remaining challenges. AI-aided microbial production harnesses the power of processing, learning, and predicting vast amounts of biological data within seconds, providing outputs with high probability. With well-trained AI models, the conventional Design-Build-Test (DBT) cycle has been transformed into a multidimensional Design-Build-Test-Learn-Predict (DBTLP) workflow, leading to significantly improved operational efficiency and reduced labor consumption. Here, we comprehensively review the main components and recent advances in AI-aided microbial production, focusing on genome annotation, AI-aided protein engineering, artificial functional protein design, and AI-enabled pathway prediction. Finally, we discuss the challenges of integrating novel AI techniques into biology and propose the potential of large language models (LLMs) in advancing microbial production. • The integration of artificial intelligence (AI) and biology accelerates the progress of microbial production. • Machine learning have been used in genome mining, protein discovery, protein design, and biosynthetic pathway design. • Large language models (LLMs) represent a novel and efficient approach to interpreting biological language. • Significant challenges remain, but the future of AI-aided microbial production is foreseeable. [ABSTRACT FROM AUTHOR]

The identification of oleaginous yeast species capable of simultaneously utilizing xylose and glucose as substrates to generate value-added biological products is an area of key economic interest. We have previously demonstrated that the Cutaneotrichosporon dermatis NICC30027 yeast strain is capable of simultaneously assimilating both xylose and glucose, resulting in considerable lipid accumulation. However, as no high-quality genome sequencing data or associated annotations for this strain are available at present, it remains challenging to study the metabolic mechanisms underlying this phenotype. Herein, we report a 39,305,439 bp draft genome assembly for C. dermatis NICC30027 comprised of 37 scaffolds, with 60.15% GC content. Within this genome, we identified 524 tRNAs, 142 sRNAs, 53 miRNAs, 28 snRNAs, and eight rRNA clusters. Moreover, repeat sequences totaling 1,032,129 bp in length were identified (2.63% of the genome), as were 14,238 unigenes that were 1,789.35 bp in length on average (64.82% of the genome). The NCBI non-redundant protein sequences (NR) database was employed to successfully annotate 11,795 of these unigenes, while 3,621 and 11,902 were annotated with the Swiss-Prot and TrEMBL databases, respectively. Unigenes were additionally subjected to pathway enrichment analyses using the Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Cluster of Orthologous Groups of proteins (COG), Clusters of orthologous groups for eukaryotic complete genomes (KOG), and Non-supervised Orthologous Groups (eggNOG) databases. Together, these results provide a foundation for future studies aimed at clarifying the mechanistic basis for the ability of C. dermatis NICC30027 to simultaneously utilize glucose and xylose to synthesize lipids.

Fragaria vesca, a wild diploid strawberry, serves as a fundamental research model for cultivated strawberry. The current reference genomes available are limited to two closely-related accessions, Hawaii 4 and CFRA2339. The widely-used model accession 'Yellow Wonder' does not yet have its reference genome. In this study, the genome of a 7th generation inbred 'Yellow Wonder' was assembled using a combination of Oxford Nanopore long reads and Illumina short reads. The de novo chromosome-scale assembly of this 220 megabase genome possesses 34,007 genes which were annotated through lift over from the Hawaii 4 genome annotation. Genome comparisons show that the 'Yellow Wonder' genome is relatively distinct from the two previously published F. vesca accessions, Hawaii 4 and CFRA2339. The availability of a 'Yellow Wonder' reference genome adds another important genomic resource to Fragaria vesca and enables rapid research progress in strawberry.

Dollar spot is one of the most damaging diseases in turfgrass, reducing its quality and playability. Two species, Clarireedia monteithiana and C. jacksonii (formerly Sclerotinia homoeocarpa) have been reported so far in the United States To study the Clarireedia genome, two isolates H2 and H3, sampled from seashore paspalum in Hawaii in 2019 were sequenced via Illumina paired-end sequencing by synthesis technology and PacBio SMRT sequencing. Both isolates were identified as C. aff. paspali, a novel species in the United States Using short and long reads, C. aff. paspali H3 contained 193 contigs with 48.6 Mbp and presented the most completed assembly and annotation among Clarireedia species. Out of the 13,428 protein models from AUGUSTUS, 349 cytoplasmic effectors and 13 apoplastic effectors were identified by EffectorP. To further decipher Clarireedia pathogenicity, C. aff. paspali genomes (H2 and H3), as well as available C. jacksonii (LWC-10 and HRI11), C. monteithiana (DRR09 and RB-19) genomes were screened for fifty-four pathogenesis determinants, previously identified in S. sclerotiorum. Seventeen orthologs of pathogenicity genes have been identified in Clarireedia species involved in oxalic acid production (pac1, nox1), mitogen-activated protein kinase cascade (pka1, smk3, ste12), appressorium formation (caf1, pks13, ams2, rgb1, rhs1) and glycolytic pathway (gpd). Within these genes, 366 species-specific SNPs were recorded between Clarireedia species; twenty-eight were non-synonymous and non-conservative. The predicted protein structure of six of these genes showed superimposition of the models among Clarireedia spp. The genomic variations revealed here could potentially lead to differences in pathogenesis and other physiological functions among Clarireedia species.

Recently developed methods to predict three-dimensional protein structure with high accuracy have opened new avenues for genome and proteome research. We explore a new hypothesis in genome annotation, namely whether computationally predicted structures can help to identify which of multiple possible gene isoforms represents a functional protein product. Guided by protein structure predictions, we evaluated over 230,000 isoforms of human protein-coding genes assembled from over 10,000 RNA sequencing experiments across many human tissues. From this set of assembled transcripts, we identified hundreds of isoforms with more confidently predicted structure and potentially superior function in comparison to canonical isoforms in the latest human gene database. We illustrate our new method with examples where structure provides a guide to function in combination with expression and evolutionary evidence. Additionally, we provide the complete set of structures as a resource to better understand the function of human genes and their isoforms. These results demonstrate the promise of protein structure prediction as a genome annotation tool, allowing us to refine even the most highly curated catalog of human proteins. More generally we demonstrate a practical, structure-guided approach that can be used to enhance the annotation of any genome.

Assignment of gene function has been a crucial, laborious, and time-consuming step in genomics. Due to a variety of sequencing platforms that generates increasing amounts of data, manual annotation is no longer feasible. Thus, the need for an integrated, automated pipeline allowing the use of experimental data towards validation of in silico prediction of gene function is of utmost relevance. Here, we present a computational workflow named AnnotaPipeline that integrates distinct software and data types on a proteogenomic approach to annotate and validate predicted features in genomic sequences. Based on FASTA (i) nucleotide or (ii) protein sequences or (iii) structural annotation files (GFF3), users can input FASTQ RNA-seq data, MS/MS data from mzXML or similar formats, as the pipeline uses both transcriptomic and proteomic information to corroborate annotations and validate gene prediction, providing transcription and expression evidence for functional annotation. Reannotation of the available Arabidopsis thaliana, Caenorhabditis elegans, Candida albicans, Trypanosoma cruzi, and Trypanosoma rangeli genomes was performed using the AnnotaPipeline, resulting in a higher proportion of annotated proteins and a reduced proportion of hypothetical proteins when compared to the annotations publicly available for these organisms. AnnotaPipeline is a Unix-based pipeline developed using Python and is available at: https://github.com/bioinformatics-ufsc/AnnotaPipeline. [ABSTRACT FROM AUTHOR]

Advances in medicine and biotechnology rely on a deep understanding of biological processes. Despite the increasingly available types and amounts of omics data, significant knowledge gaps remain, with current approaches to identify and curate missing annotations being limited to a set of already known reactions. Here, we introduce Network Integrated Computational Explorer for Gap Annotation of Metabolism (NICEgame), a workflow to identify and curate nonannotated metabolic functions in genomes using the ATLAS of Biochemistry and genome-scale metabolic models (GEMs). To resolve gaps in GEMs, NICEgame provides alternative sets of known and hypothetical reactions, assesses their thermodynamic feasibility, and suggests candidate genes to catalyze these reactions. We identified metabolic gaps and applied NICEgame in the latest GEM of Escherichia coli, iML1515, and enhanced the E. coli genome annotation by resolving 47% of these gaps. NICEgame, applicable to any GEM and functioning from open-source software, should thus enhance all GEM-based predictions and subsequent biotechnological and biomedical applications. [ABSTRACT FROM AUTHOR]

Background: Identification of deleterious genetic variants using DNA sequencing data relies on increasingly detailed filtering strategies to isolate the small subset of variants that are more likely to underlie a disease phenotype. Datasets reflecting population allele frequencies of different types of variants serve as powerful filtering tools, especially in the context of rare disease analysis. While such population-scale allele frequency datasets now exist for structural variants (SVs), it remains a challenge to match SV calls between multiple datasets, thereby complicating estimates of a putative SV's population allele frequency. Results: We introduce SVAFotate, a software tool that enables the annotation of SVs with variant allele frequency and related information from existing SV datasets. As a result, VCF files annotated by SVAFotate offer a variety of metrics to aid in the stratification of SVs as common or rare in the broader human population. Conclusions: Here we demonstrate the use of SVAFotate in the classification of SVs with regards to their population frequency and illustrate how SVAFotate's annotations can be used to filter and prioritize SVs. Lastly, we detail how best to utilize these SV annotations in the analysis of genetic variation in studies of rare disease. [ABSTRACT FROM AUTHOR]

The European starling, Sturnus vulgaris, is an ecologically significant, globally invasive avian species that is also suffering from a major decline in its native range. Here, we present the genome assembly and long‐read transcriptome of an Australian‐sourced European starling (S. vulgaris vAU), and a second, North American, short‐read genome assembly (S. vulgaris vNA), as complementary reference genomes for population genetic and evolutionary characterization. S. vulgaris vAU combined 10× genomics linked‐reads, low‐coverage Nanopore sequencing, and PacBio Iso‐Seq full‐length transcript scaffolding to generate a 1050 Mb assembly on 6222 scaffolds (7.6 Mb scaffold N50, 94.6% busco completeness). Further scaffolding against the high‐quality zebra finch (Taeniopygia guttata) genome assigned 98.6% of the assembly to 32 putative nuclear chromosome scaffolds. Species‐specific transcript mapping and gene annotation revealed good gene‐level assembly and high functional completeness. Using S. vulgaris vAU, we demonstrate how the multifunctional use of PacBio Iso‐Seq transcript data and complementary homology‐based annotation of sequential assembly steps (assessed using a new tool, saaga) can be used to assess, inform, and validate assembly workflow decisions. We also highlight some counterintuitive behaviour in traditional busco metrics, and present buscomp, a complementary tool for assembly comparison designed to be robust to differences in assembly size and base‐calling quality. This work expands our knowledge of avian genomes and the available toolkit for assessing and improving genome quality. The new genomic resources presented will facilitate further global genomic and transcriptomic analysis on this ecologically important species. [ABSTRACT FROM AUTHOR]

A major limiting factor in target discovery for both basic research and therapeutic intervention is the identification of structural and/or functional RNA elements in genomes and transcriptomes. This was the impetus for the original ScanFold algorithm, which provides maps of local RNA structural stability, evidence of sequence-ordered (potentially evolved) structure, and unique model structures comprised of recurring base pairs with the greatest structural bias. A key step in quantifying this propensity for ordered structure is the prediction of secondary structural stability for randomized sequences which, in the original implementation of ScanFold, is explicitly evaluated. This slow process has limited the rapid identification of ordered structures in large genomes/transcriptomes, which we seek to overcome in this current work introducing ScanFold 2.0. In this revised version of ScanFold, we no longer explicitly evaluate randomized sequence folding energy, but rather estimate it using a machine learning approach. For high randomization numbers, this can increase prediction speeds over 100-fold compared to ScanFold 1.0, allowing for the analysis of large sequences, as well as the use of additional folding algorithms that may be computationally expensive. In the testing of ScanFold 2.0, we re-evaluate the Zika, HIV, and SARS-CoV-2 genomes and compare both the consistency of results and the time of each run to ScanFold 1.0. We also re-evaluate the SARS-CoV-2 genome to assess the quality of ScanFold 2.0 predictions vs several biochemical structure probing datasets and compare the results to those of the original ScanFold program. [ABSTRACT FROM AUTHOR]

RNA sequencing (RNA-Seq) and mass-spectrometry-based proteomics data are often integrated in proteogenomic studies to assist in the prediction of eukaryote genome features, such as genes, splicing, single-nucleotide (SNVs), and single-amino-acid variants (SAAVs). Most genomes of parasite nematodes are draft versions that lack transcript- and protein-level information and whose gene annotations rely only on computational predictions. Angiostrongylus costaricensis is a roundworm species that causes an intestinal inflammatory disease, known as abdominal angiostrongyliasis (AA). Currently, there is no drug available that acts directly on this parasite, mostly due to the sparse understanding of its molecular characteristics. The available genome of A. costaricensis, specific to the Costa Rica strain, is a draft version that is not supported by transcript- or protein-level evidence. This study used RNA-Seq and MS/MS data to perform an in-depth annotation of the A. costaricensis genome. Our prediction improved the reference annotation with (a) novel coding and non-coding genes; (b) pieces of evidence of alternative splicing generating new proteoforms; and (c) a list of SNVs between the Brazilian (Crissiumal) and the Costa Rica strain. To the best of our knowledge, this is the first time that a multi-omics approach has been used to improve the genome annotation of A. costaricensis. We hope this improved genome annotation can assist in the future development of drugs, kits, and vaccines to treat, diagnose, and prevent AA caused by either the Brazil strain (Crissiumal) or the Costa Rica strain. [ABSTRACT FROM AUTHOR]

Burkholderia glumae is an important rice pathogen, thus the genomic and evolutionary history may be helpful to control this notorious pathogen. Here, we present two complete genomes of the B. glumae strains HN1 and HN2, which were isolated from diseased rice seed in China. Average nucleotide identity (ANI) analysis shows greater than 99% similarity of the strains HN1 and HN2 with other published B. glumae genomes. Genomic annotation revealed that the genome of strain HN1 consists of five replicons (6,680,415 bp) with an overall G + C content of 68.06%, whereas the genome of strain HN2 comprises of three replicons (6,560,085 bp) with an overall G + C content of 68.34%. The genome of HN1 contains 5434 protein-coding genes, 351 pseudogenes, and 1 CRISPR, whereas the genome of HN2 encodes 5278 protein-coding genes, 357 pseudogenes, and 2 CRISPR. Both strains encode many pathogenic-associated genes (143 genes in HN1 vs. 141 genes in HN2). Moreover, comparative genomic analysis shows the extreme plasticity of B. glumae, which may contribute to its pathogenicity. In total, 259 single-copy genes were affected by positive selection. These genes may contribute to the adaption to different environments. Notably, six genes were characterized as virulence factors which may be an additional way to assist the pathogenicity of B. glumae. [ABSTRACT FROM AUTHOR]

The European earwig Forficula auricularia is an important model for studies of maternal care, sexual selection, sociality, and host--parasite interactions. However, detailed genetic investigations of this species are hindered by a lack of genomic resources. Here, we present a high-quality hybrid genome assembly for Forficula auricularia using Nanopore long-reads and 10x linked-reads. The final assembly is 1.06Gb in length with 31.03% GC content. It consists of 919 scaffolds with an N50 of 12.55Mb. Half of the genome is present in only 20 scaffolds. Benchmarking Universal Single-Copy Orthologs scores are ~90% from 3 sets of single-copy orthologs (eukaryotic, insect, and arthropod). The total repeat elements in the genome are 64.62%. The MAKER2 pipeline annotated 12,876 protein-coding genes and 21,031 mRNAs. Phylogenetic analysis revealed the assembled genome as that of species B, one of the 2 known genetic subspecies of Forficula auricularia. The genome assembly, annotation, and associated resources will be of high value to a large and diverse group of researchers working on dermapterans. [ABSTRACT FROM AUTHOR]

As a newly described Sillaginidae species, Chinese sillago (Sillago sinica) needs a better understanding of gene annotation information. In this study, we reported the first full-length transcriptome data of S. sinica using the PacBio isoform sequencing Iso-seq and a description of transcriptome structure analysis. A total of 454,979 high-quality full-length transcripts were obtained by single-molecule real-time (SMRT) sequencing, which was corrected by Illumina sequencing data. After that, 66,948 non-redundant full-length transcripts were generated after mapping to the reference genome of S. sinica, including 49 fusion isoforms and 9,250 novel isoforms. 63,459 isoforms were successfully annotated by one of the Nr, Nt, SwissProt, Pfam, KOG, GO, and KEGG databases. Additionally, 30,987 alternative polyadenylation (APA) sites, 451,867 alternative splicing (AS) events, 21,928 long non-coding RNAs (lncRNAs) and 12,911 transcription factors (TFs) were identified. The full-length transcripts of S. sinica would provide a precious resource for characterizing the transcriptome of S. sinica and for the further study of gene function and regulatory mechanism of this species. [ABSTRACT FROM AUTHOR]