The absence of canonical respiratory complex I subunits in male-type mitogenomes of three Donax species.

Bivalves are an extraordinary class of animals in which species with a doubly uniparental inheritance (DUI) of mitochondrial DNA have been described. DUI is characterized as a mitochondrial homoplasmy of females and heteroplasmy of male individuals where F-type mitogenomes are passed to the progeny with mother egg cells and divergent M-type mitogenomes are inherited with fathers sperm cells. However, in most cases only male individuals retain divergent mitogenome inherited with spermatozoa. Additionally, in many of bivalves, unique mitochondrial features, like additional genes, gene duplication, gene extensions, mitochondrial introns, and recombination, were observed. In this study, we sequenced and assembled male-type mitogenomes of three Donax species. Comparative analysis of mitochondrial sequences revealed a lack of all seven NADH dehydrogenase subunits as well as the presence of three long additional open reading frames lacking identifiable homology to any of the existing genes.

Human parainfluenza virus 3 vaccine candidates attenuated by codon-pair deoptimization are immunogenic and protective in hamsters.

Human parainfluenza virus type 3 (HPIV3) is a major pediatric respiratory pathogen lacking available vaccines or antiviral drugs. We generated live-attenuated HPIV3 vaccine candidates by codon-pair deoptimization (CPD). HPIV3 open reading frames (ORFs) encoding the nucleoprotein (N), phosphoprotein (P), matrix (M), fusion (F), hemagglutinin-neuraminidase (HN), and polymerase (L) were modified singly or in combination to generate 12 viruses designated Min-N, Min-P, Min-M, Min-FHN, Min-L, Min-NP, Min-NPM, Min-NPL, Min-PM, Min-PFHN, Min-MFHN, and Min-PMFHN. CPD of N or L severely reduced growth in vitro and was not further evaluated. CPD of P or M was associated with increased and decreased interferon (IFN) response in vitro, respectively, but had little effect on virus replication. In Vero cells, CPD of F and HN delayed virus replication, but final titers were comparable to wild-type (wt) HPIV3. In human lung epithelial A549 cells, CPD F and HN induced a stronger IFN response, viral titers were reduced 100-fold, and the expression of F and HN proteins was significantly reduced without affecting N or P or the relative packaging of proteins into virions. Following intranasal infection in hamsters, replication in the nasal turbinates and lungs tended to be the most reduced for viruses bearing CPD F and HN, with maximum reductions of approximately 10-fold. Despite decreased in vivo replication (and lower expression of CPD F and HN in vitro), all viruses induced titers of serum HPIV3-neutralizing antibodies similar to wt and provided complete protection against HPIV3 challenge. In summary, CPD of HPIV3 yielded promising vaccine candidates suitable for further development.

Phosphorylation of phospholamban promotes SERCA2a activation by dwarf open reading frame (DWORF).

In cardiac myocytes, the type 2a sarco/endoplasmic reticulum Ca-ATPase (SERCA2a) plays a key role in intracellular Ca regulation. Due to its critical role in heart function, SERCA2a activity is tightly regulated by different mechanisms, including micropeptides. While phospholamban (PLB) is a well-known SERCA2a inhibitor, dwarf open reading frame (DWORF) is a recently identified SERCA2a activator. Since PLB phosphorylation is the most recognized mechanism of SERCA2a activation during adrenergic stress, we studied whether PLB phosphorylation also affects SERCA2a regulation by DWORF. By using confocal Ca imaging in a HEK293 expressing cell system, we analyzed the effect of the co-expression of PLB and DWORF using a bicistronic construct on SERCA2a-mediated Ca uptake. Under these conditions of matched expression of PLB and DWORF, we found that SERCA2a inhibition by non-phosphorylated PLB prevails over DWORF activating effect. However, when PLB is phosphorylated at PKA and CaMKII sites, not only PLB's inhibitory effect was relieved, but SERCA2a was effectively activated by DWORF. Förster resonance energy transfer (FRET) analysis between SERCA2a and DWORF showed that DWORF has a higher relative affinity for SERCA2a when PLB is phosphorylated. Thus, SERCA2a regulation by DWORF responds to the PLB phosphorylation status, suggesting that DWORF might contribute to SERCA2a activation during conditions of adrenergic stress.

SATIN: a micro and mini satellite mining tool of total genome and coding regions with analysis of perfect repeats polymorphism in coding regions.

BACKGROUND: Tandem repeats are specific sequences in genomic DNA repeated in tandem that are present in all organisms. Among the subcategories of TRs we have Satellite repeats, that is divided into macrosatellites, minisatellites, and microsatellites, being the last two of specific interest because they can identify polymorphisms between organisms due to their instability. Currently, most mining tools focus on Simple Sequence Repeats (SSR) mining, and only a few can identify SSRs in the coding regions. RESULTS: We developed a microsatellite mining software called SATIN (Micro and Mini SATellite IdentificatioN tool) based on a new sliding window algorithm written in C and Python. It represents a new approach to SSR mining by addressing the limitations of existing tools, particularly in coding region SSR mining. SATIN is available at https://github.com/labgm/SATIN.git . It was shown to be the second fastest for perfect and compound SSR mining. It can identify SSRs from coding regions plus SSRs with motif sizes bigger than 6. Besides the SSR mining, SATIN can also analyze SSRs polymorphism on coding-regions from pre-determined groups, and identify SSRs differentially abundant among them on a per-gene basis. To validate, we analyzed SSRs from two groups of Escherichia coli (K12 and O157) and compared the results with 5 known SSRs from coding regions. SATIN identified all 5 SSRs from 237 genes with at least one SSR on it. CONCLUSIONS: The SATIN is a novel microsatellite search software that utilizes an innovative sliding window technique based on a numerical list for repeat region search to identify perfect, and composite SSRs while generating comprehensible and analyzable outputs. It is a tool capable of using files in fasta or GenBank format as input for microsatellite mining, also being able to identify SSRs present in coding regions for GenBank files. In conclusion, we expect SATIN to help identify potential SSRs to be used as genetic markers.

Complete genome of a novel mycobacteriophage WXIN isolated in Wuhan, China.

OBJECTIVES: The rising of antibiotic resistance has sparked a renewed interest in mycobacteriophage as alternative therapeutic strategies against mycobacterial infections. So far, the vast majority of mycobacteriophages have been isolated using the model species Mycobacterium smegmatis, implying an overwhelming majority of mycobacteriophages in the environment remain uncultured, unclassified, and their specific hosts and infection strategies are still unknown. This study was undertaken to isolate and characterize novel mycobacteriophages targeting Mycobacterium septicum. DATA DESCRIPTION: Here a novel mycobacteriophage WXIN against M. septicum was isolated from soil samples in Wuhan, China. Whole genome analysis indicates that the phage genome consists of 115,158 bp with a GC content of 61.9%. Of the 260 putative open reading frames, 46 may be associated with phage packaging, structure, lysis, lysogeny, genome modification/replication, and other functional roles. The limited genome-wide similarity, along with phylogenetic trees constructed based on viral proteome and orthologous genes show that phage WXIN represents a novel cluster distantly related to cluster J mycobacteriophages (genus Omegavirus). Overall, these results provide novel insights into the genomic properties of mycobacteriophages, highlighting the great genetic diversity of mycobacteriophages in relation to their hosts.

CRISPR-Based Modular Assembly for High-Throughput Construction of a UAS-cDNA/ORF Plasmid Library.

Functional genomics screening offers a powerful approach to probe gene function and relies on the construction of genome-wide plasmid libraries. Conventional approaches for plasmid library construction are time-consuming and laborious. Therefore, we recently developed a simple and efficient method, CRISPR-based modular assembly (CRISPRmass), for high-throughput construction of a genome-wide upstream activating sequence-complementary DNA/open reading frame (UAS-cDNA/ORF) plasmid library. Here, we present a protocol for CRISPRmass, taking as an example the construction of a GAL4/UAS-based UAS-cDNA/ORF plasmid library. The protocol includes massively parallel two-step test tube reactions followed by bacterial transformation. The first step is to linearize the existing complementary DNA (cDNA) or open reading frame (ORF) cDNA or ORF library plasmids by cutting the shared upstream vector sequences adjacent to the 5' end of cDNAs or ORFs using CRISPR/Cas9 together with single guide RNA (sgRNA), and the second step is to insert a UAS module into the linearized cDNA or ORF plasmids using a single step reaction. CRISPRmass allows the simple, fast, efficient, and cost-effective construction of various plasmid libraries. The UAS-cDNA/ORF plasmid library can be utilized for gain-of-function screening in cultured cells and for constructing a genome-wide transgenic UAS-cDNA/ORF library in Drosophila.

Probing of plant transcriptomes reveals the hidden genetic diversity of the family Secoviridae.

Secoviruses are single-stranded RNA viruses that infect plants. In the present study, we identified 61 putative novel secoviral genomes in various plant species by mining publicly available plant transcriptome data. These viral sequences represent the genomes of 13 monopartite and 48 bipartite secovirids. The genome sequences of 52 secovirids were coding-complete, and nine were partial. Except for small open reading frames (ORFs) determined in waikaviral genomes and RNA2 of torradoviruses, all of the recovered genomes/genome segments contained a large ORF encoding a polyprotein. Based on genome organization and phylogeny, all but three of the novel secoviruses were assigned to different genera. The genome organization of two identified waika-like viruses resembled that of the recently identified waika-like virus Triticum aestivum secovirus. Phylogenetic analysis revealed a pattern of host-virus co-evolution in a few waika- and waika-like viruses and increased phylogenetic diversity of nepoviruses. The study provides a basis for further investigation of the biological properties of these novel secoviruses.

GNN Codon Adjacency Tunes Protein Translation.

The central dogma treats the ribosome as a molecular machine that reads one mRNA codon at a time as it adds each amino acid to its growing peptide chain. However, this and previous studies suggest that ribosomes actually perceive pairs of adjacent codons as they take three-nucleotide steps along the mRNA. We examined GNN codons, which we find are surprisingly overrepresented in eukaryote protein-coding open reading frames (ORFs), especially immediately after NNU codons. Ribosome profiling experiments in yeast revealed that ribosomes with NNU at their aminoacyl (A) site have particularly elevated densities when NNU is immediately followed (3') by a GNN codon, indicating slower mRNA threading of the NNU codon from the ribosome's A to peptidyl (P) sites. Moreover, if the assessment was limited to ribosomes that have only recently arrived at the next codon, by examining 21-nucleotide ribosome footprints (21-nt RFPs), elevated densities were observed for multiple codon classes when followed by GNN. This striking translation slowdown at adjacent 5'-NNN GNN codon pairs is likely mediated, in part, by the ribosome's CAR surface, which acts as an extension of the A-site tRNA anticodon during ribosome translocation and interacts through hydrogen bonding and pi stacking with the GNN codon. The functional consequences of 5'-NNN GNN codon adjacency are expected to influence the evolution of protein coding sequences.

Complete coding region sequence analyses and antigenic characterization of emerging lineage G-IX of foot- and-mouth disease virus serotype Asia1.

Foot-and-mouth disease Virus (FMDV) serotype Asia1 is prevalent in the Indian subcontinent, with only G-III and G-VIII reported in India until 2020. However, in 2019, a novel genetic group within serotype Asia1, designated as G-IX, emerged in Bangladesh, followed by its detection in India in 2020. This report presents analyses of the complete coding region sequences of the G-IX lineage isolates. The length of the open reading frame (ORF) of the two G-IX isolates was 6990 nucleotides without any deletion or insertion. The G-IX isolates showed the highest sequence similarity with an isolate of G-III at the ORF, L, P2, and P3 regions, and with an isolate of G-VIII at the P1 region. Phylogenetic analysis based on the capsid region (P1) supports the hypothesis that G-VIII and G-IX originated from a common ancestor, as speculated earlier. Further, VP1 region-based phylogenetic analyses revealed the re-emergence of G-VIII after a gap of 3 years. One isolate of G-VIII collected during 2023 revealed a codon insertion in the G-H loop of VP1. The vaccine matching studies support the suitability of the currently used Indian vaccine strain IND63/1972 to contain outbreaks due to viruses belonging to G-IX.

A Targeted Deep Sequencing Method to Quantify Endogenous Retrovirus Gag Sequence Variants and Open Reading Frames Expressed in Nonobese Diabetic Mice.

Endogenous retroviruses (ERVs) are involved in autoimmune diseases such as type 1 diabetes (T1D). ERV gene products homologous to murine leukemia retroviruses are expressed in the pancreatic islets of NOD mice, a model of T1D. One ERV gene, Gag, with partial or complete open reading frames (ORFs), is detected in the islets, and it contains many sequence variants. An amplicon deep sequencing analysis was established by targeting a conserved region within the Gag gene to compare NOD with T1D-resistant mice or different ages of prediabetic NOD mice. We observed that the numbers of different Gag variants and ORFs are linked to T1D susceptibility. More importantly, these numbers change during the course of diabetes development and can be quantified to calculate the levels of disease progression. Sequence alignment analysis led to identification of additional markers, including nucleotide mismatching and amino acid consensus at specific positions that can distinguish the early and late stages, before diabetes onset. Therefore, the expression of sequence variants and ORFs of ERV genes, particularly Gag, can be quantified as biomarkers to estimate T1D susceptibility and disease progression.

A new point mutation in the HC-Pro of potato virus Y is involved in tobacco vein necrosis.

Tobacco vein necrosis (TVN) is a complex phenomenon regulated by different genetic determinants mapped in the HC-Pro protein (amino acids N330, K391 and E410) and in two regions of potato virus Y (PVY) genome, corresponding to the cytoplasmic inclusion (CI) protein and the nuclear inclusion protein a-protease (NIa-Pro), respectively. A new determinant of TVN was discovered in the MK isolate of PVY which, although carried the HC-Pro determinants associated to TVN, did not induce TVN. The HC-Pro open reading frame (ORF) of the necrotic infectious clone PVY N605 was replaced with that of the non-necrotic MK isolate, which differed only by one amino acid at position 392 (T392 instead of I392). The cDNA clone N605_MKHCPro inoculated in tobacco induced only weak mosaics at the systemic level, demostrating that the amino acid at position 392 is a new determinant for TVN. No significant difference in accumulation in tobacco was observed between N605 and N605_MKHCPro. Since phylogenetic analyses showed that the loss of necrosis in tobacco has occurred several times independently during PVY evolution, these repeated evolutions strongly suggest that tobacco necrosis is a costly trait in PVY.

Complete characterization of the yak testicular development using accurate full-length transcriptome sequencing.

Alternative splicing is a prevalent phenomenon in testicular tissues. Due to the low assembly accuracy of short-read RNA sequencing technology in analyzing post-transcriptional regulatory events, full-length (FL) transcript sequencing is highly demanded to accurately determine FL splicing variants. In this study, we performed FL transcriptome sequencing of testicular tissues from 0.5, 1.5, 2.5, and 4-year-old yaks and 4-year-old cattle-yaks using Oxford Nanopore Technologies. The obtained sequencing data were predicted to have 47,185 open reading frames (ORFs), including 26,630 complete ORFs, detected 7645 fusion transcripts, 15,355 alternative splicing events, 25,798 simple sequence repeats, 7628 transcription factors, and 35,503 long non-coding RNAs. A total of 40,038 novel transcripts were obtained from the sequencing data, and the proportion was almost close to the number of known transcripts identified. Structural analysis and functional annotation of these novel transcripts resulted in the successful annotation of 9568 transcripts, with the highest and lowest annotation numbers in the Nr and KOG databases, respectively. Weighted gene co-expression network analysis revealed the key regulatory pathways and hub genes at various stages of yak testicular development. Our findings enhance our comprehension of transcriptome complexity, contribute to genome annotation refinement, and provide foundational data for further investigations into male sterility in cattle-yaks.

Analysis of a new negevirus-like sequence from Bemisia tabaci unveils a potential new taxon linking nelorpi- and centiviruses.

This study presents the complete genome sequence of a novel nege-like virus identified in whiteflies (Bemisia tabaci MEAM1), provisionally designated as whitefly negevirus 1 (WfNgV1). The virus possesses a single-stranded RNA genome comprising 11,848 nucleotides, organized into four open reading frames (ORFs). These ORFs encode the putative RNA-dependent-RNA-polymerase (RdRp, ORF 1), a glycoprotein (ORF 2), a structural protein with homology to those in the SP24 family, (ORF 3), and a protein of unknown function (ORF 4). Phylogenetic analysis focusing on RdRp and SP24 amino acid sequences revealed a close relationship between WfNgV1 and Bemisia tabaci negevirus 1, a negevirus sequence recently discovered in whiteflies from Israel. Both viruses form a clade sharing a most recent common ancestor with the proposed nelorpivirus and centivirus taxa. The putative glycoprotein from ORF 2 and SP24 (ORF 3) of WfNgV1 exhibit the characteristic topologies previously reported for negevirus counterparts. This marks the first reported negevirus-like sequence from whiteflies in the Americas.

Characterization of microsatellite markers in the coding regions of the Penaeus vannamei genome.

In this study, an extensive analysis of microsatellite markers (Single Tandem Repeats-STRs) in Penaeus vannamei was conducted at an advanced level. The markers were thoroughly examined, characterized, and specific markers located within coding regions were identified. Out of a total of 306 STRs, 117 were classified as perfect markers based on their single repeat motif. Among these perfect markers, 62 were found to be associated with predicted coding genes (mRNA), which were involved in various functions such as binding, catalytic activity, ATP-dependent activity, transcription, structural and molecular regulation. To validate the accuracy of the findings, a sample of nine markers was subjected to in vitro testing, which confirmed the presence of polymorphisms within the population. These results suggest the existence of different protein isoforms within the population, indicating the potential of these markers for application in both population and phenotype-genotype association studies. This innovative approach opens up new possibilities for investigating the impact of genomic plasticity in populations of P. vannamei.

UPF1 regulates mRNA stability by sensing poorly translated coding sequences.

Post-transcriptional mRNA regulation shapes gene expression, yet how cis-elements and mRNA translation interface to regulate mRNA stability is poorly understood. We find that the strength of translation initiation, upstream open reading frame (uORF) content, codon optimality, AU-rich elements, microRNA binding sites, and open reading frame (ORF) length function combinatorially to regulate mRNA stability. Machine-learning analysis identifies ORF length as the most important conserved feature regulating mRNA decay. We find that Upf1 binds poorly translated and untranslated ORFs, which are associated with a higher decay rate, including mRNAs with uORFs and those with exposed ORFs after stop codons. Our study emphasizes Upf1's converging role in surveilling mRNAs with exposed ORFs that are poorly translated, such as mRNAs with long ORFs, ORF-like 3' UTRs, and mRNAs containing uORFs. We propose that Upf1 regulation of poorly/untranslated ORFs provides a unifying mechanism of surveillance in regulating mRNA stability and homeostasis in an exon-junction complex (EJC)-independent nonsense-mediated decay (NMD) pathway that we term ORF-mediated decay (OMD).

An Inner Mitochondrial Membrane Microprotein from the SLC35A4 Upstream ORF Regulates Cellular Metabolism.

Upstream open reading frames (uORFs) are cis-acting elements that can dynamically regulate the translation of downstream ORFs by suppressing downstream translation under basal conditions and, in some cases, increasing downstream translation under stress conditions. Computational and empirical methods have identified uORFs in the 5'-UTRs of approximately half of all mouse and human transcripts, making uORFs one of the largest regulatory elements known. Because the prevailing dogma was that eukaryotic mRNAs produce a single functional protein, the peptides and small proteins, or microproteins, encoded by uORFs were rarely studied. We hypothesized that a uORF in the SLC35A4 mRNA is producing a functional microprotein (SLC35A4-MP) because of its conserved amino acid sequence. Through a series of biochemical and cellular experiments, we find that the 103-amino acid SLC35A4-MP is a single-pass transmembrane inner mitochondrial membrane (IMM) microprotein. The IMM contains the protein machinery crucial for cellular respiration and ATP generation, and loss of function studies with SLC35A4-MP significantly diminish maximal cellular respiration, indicating a vital role for this microprotein in cellular metabolism. The findings add SLC35A4-MP to the growing list of functional microproteins and, more generally, indicate that uORFs that encode conserved microproteins are an untapped reservoir of functional microproteins.

Evolution of a Human-Specific De Novo Open Reading Frame and Its Linked Transcriptional Silencer.

In the human genome, two short open reading frames (ORFs) separated by a transcriptional silencer and a small intervening sequence stem from the gene SMIM45. The two ORFs show different translational characteristics, and they also show divergent patterns of evolutionary development. The studies presented here describe the evolution of the components of SMIM45. One ORF consists of an ultra-conserved 68 amino acid (aa) sequence, whose origins can be traced beyond the evolutionary age of divergence of the elephant shark, ~462 MYA. The silencer also has ancient origins, but it has a complex and divergent pattern of evolutionary formation, as it overlaps both at the 68 aa ORF and the intervening sequence. The other ORF consists of 107 aa. It develops during primate evolution but is found to originate de novo from an ancestral non-coding genomic region with root origins within the Afrothere clade of placental mammals, whose evolutionary age of divergence is ~99 MYA. The formation of the complete 107 aa ORF during primate evolution is outlined, whereby sequence development is found to occur through biased mutations, with disruptive random mutations that also occur but lead to a dead-end. The 107 aa ORF is of particular significance, as there is evidence to suggest it is a protein that may function in human brain development. Its evolutionary formation presents a view of a human-specific ORF and its linked silencer that were predetermined in non-primate ancestral species. The genomic position of the silencer offers interesting possibilities for the regulation of transcription of the 107 aa ORF. A hypothesis is presented with respect to possible spatiotemporal expression of the 107 aa ORF in embryonic tissues.

Neuronal activity rapidly reprograms dendritic translation via eIF4G2:uORF binding.

Learning and memory require activity-induced changes in dendritic translation, but which mRNAs are involved and how they are regulated are unclear. In this study, to monitor how depolarization impacts local dendritic biology, we employed a dendritically targeted proximity labeling approach followed by crosslinking immunoprecipitation, ribosome profiling and mass spectrometry. Depolarization of primary cortical neurons with KCl or the glutamate agonist DHPG caused rapid reprogramming of dendritic protein expression, where changes in dendritic mRNAs and proteins are weakly correlated. For a subset of pre-localized messages, depolarization increased the translation of upstream open reading frames (uORFs) and their downstream coding sequences, enabling localized production of proteins involved in long-term potentiation, cell signaling and energy metabolism. This activity-dependent translation was accompanied by the phosphorylation and recruitment of the non-canonical translation initiation factor eIF4G2, and the translated uORFs were sufficient to confer depolarization-induced, eIF4G2-dependent translational control. These studies uncovered an unanticipated mechanism by which activity-dependent uORF translational control by eIF4G2 couples activity to local dendritic remodeling.

Arabidopsis Fhit-like tumor suppressor resumes early terminated constitutive triple response1-10 mRNA translation.

The Arabidopsis (Arabidopsis thaliana) constitutive triple response1-10 (ctr1-10) mutant produces a reduced level of CTR1 protein and exhibits a weak ctr1 mutant phenotype. Sequence analysis revealed highly active translation of the upstream open reading frame (uORF) at the extended 5'-UTR of the ctr1-10 mRNA, resulting from T-DNA insertion. Enhancer screening for ctr1-10 isolated the fragile histidine triad-1 (fhit-1) mutation. The fhit-1 ctr1-10 mutant phenotypically resembled strong ctr1 mutants and barely produced CTR1, and the fhit-1 mutation reduced the translation efficiency of ctr1-10 but not that of CTR1 mRNA. The human (Homo sapiens) Fhit that involves tumorigenesis and genome instability has the in vitro dinucleotide 5',5'″-P1, P3-triphosphate hydrolase activity, and expression of the human HsFHIT or the hydrolase-defective HsFHITH96N transgene reversed the fhit-1 ctr1-10 mutant phenotype and restored CTR1 levels. Genetic editing that in situ disrupts individual upstream ATG codons proximal to the ctr1-10 mORF elevated CTR1 levels in ctr1-10 plants independent of FHIT. EUKARYOTIC INITIATION FACTOR3G (eIF3G), which is involved in translation and reinitiation, interacted with FHIT, and both were associated with the polysome. We propose that FHIT resumes early terminated ctr1-10 mORF translation in the face of active and complex uORF translation. Our study unveils a niche that may lead to investigations on the molecular mechanism of Fhit-like proteins in translation reinitiation. The biological significance of FHIT-regulated translation is discussed.

Challenges in computational discovery of bioactive peptides in 'omics data.

Peptides have a plethora of activities in biological systems that can potentially be exploited biotechnologically. Several peptides are used clinically, as well as in industry and agriculture. The increase in available 'omics data has recently provided a large opportunity for mining novel enzymes, biosynthetic gene clusters, and molecules. While these data primarily consist of DNA sequences, other types of data provide important complementary information. Due to their size, the approaches proven successful at discovering novel proteins of canonical size cannot be naïvely applied to the discovery of peptides. Peptides can be encoded directly in the genome as short open reading frames (smORFs), or they can be derived from larger proteins by proteolysis. Both of these peptide classes pose challenges as simple methods for their prediction result in large numbers of false positives. Similarly, functional annotation of larger proteins, traditionally based on sequence similarity to infer orthology and then transferring functions between characterized proteins and uncharacterized ones, cannot be applied for short sequences. The use of these techniques is much more limited and alternative approaches based on machine learning are used instead. Here, we review the limitations of traditional methods as well as the alternative methods that have recently been developed for discovering novel bioactive peptides with a focus on prokaryotic genomes and metagenomes.