Explore the vast range of titles available.

Assembled genome sequences are being generated at an exponential rate. Here we present FCS-GX, part of NCBI’s Foreign Contamination Screen (FCS) tool suite, optimized to identify and remove contaminant sequences in new genomes. FCS-GX screens most genomes in 0.1–10 min. Testing FCS-GX on artificially fragmented genomes demonstrates high sensitivity and specificity for diverse contaminant species. We used FCS-GX to screen 1.6 million GenBank assemblies and identified 36.8 Gbp of contamination, comprising 0.16% of total bases, with half from 161 assemblies. We updated assemblies in NCBI RefSeq to reduce detected contamination to 0.01% of bases. FCS-GX is available at https://github.com/ncbi/fcs/ or https://doi.org/10.5281/zenodo.10651084 .

Background Urea amidolyase breaks down urea into ammonia and carbon dioxide in a two-step process, while another enzyme, urease, does this in a one step-process. Urea amidolyase has been found only in some fungal species among eukaryotes. It contains two major domains: the amidase and urea carboxylase domains. A shorter form of urea amidolyase is known as urea carboxylase and has no amidase domain. Eukaryotic urea carboxylase has been found only in several fungal species and green algae. In order to elucidate the evolutionary origin of urea amidolyase and urea carboxylase, we studied the distribution of urea amidolyase, urea carboxylase, as well as other proteins including urease, across kingdoms. Results Among the 64 fungal species we examined, only those in two Ascomycota classes (Sordariomycetes and Saccharomycetes) had the urea amidolyase sequences. Urea carboxylase was found in many but not all of the species in the phylum Basidiomycota and in the subphylum Pezizomycotina (phylum Ascomycota). It was completely absent from the class Saccharomycetes (phylum Ascomycota; subphylum Saccharomycotina). Four Sordariomycetes species we examined had both the urea carboxylase and the urea amidolyase sequences. Phylogenetic analysis showed that these two enzymes appeared to have gone through independent evolution since their bacterial origin. The amidase domain and the urea carboxylase domain sequences from fungal urea amidolyases clustered strongly together with the amidase and urea carboxylase sequences, respectively, from a small number of beta- and gammaproteobacteria. On the other hand, fungal urea carboxylase proteins clustered together with another copy of urea carboxylases distributed broadly among bacteria. The urease proteins were found in all the fungal species examined except for those of the subphylum Saccharomycotina. Conclusions We conclude that the urea amidolyase genes currently found only in fungi are the results of a horizontal gene transfer event from beta-, gamma-, or related species of proteobacteria. The event took place before the divergence of the subphyla Pezizomycotina and Saccharomycotina but after the divergence of the subphylum Taphrinomycotina. Urea carboxylase genes currently found in fungi and other limited organisms were also likely derived from another ancestral gene in bacteria. Our study presented another important example showing plastic and opportunistic genome evolution in bacteria and fungi and their evolutionary interplay.

We determined that extrachromosomal 2μ plasmid was present in 67 of the Saccharomyces cerevisiae 100-genome strains; in addition to variation in the size and copy number of 2μ, we identified three distinct classes of 2μ. We identified 2μ presence/absence and class associations with populations, clinical origin and nuclear genotypes. We also screened genome sequences of S. paradoxus, S. kudriavzevii, S. uvarum, S. eubayanus, S. mikatae, S. arboricolus and S. bayanus strains for both integrated and extrachromosomal 2μ. Similar to S. cerevisiae, we found no integrated 2μ sequences in any S. paradoxus strains. However, we identified part of 2μ integrated into the genomes of some S. uvarum, S. kudriavzevii, S. mikatae and S. bayanus strains, which were distinct from each other and from all extrachromosomal 2μ. We identified extrachromosomal 2μ in one S. paradoxus, one S. eubayanus, two S. bayanus and 13 S. uvarum strains. The extrachromosomal 2μ in S. paradoxus, S. eubayanus and S. cerevisiae were distinct from each other. In contrast, the extrachromosomal 2μ in S. bayanus and S. uvarum strains were identical with each other and with one of the three classes of S. cerevisiae 2μ, consistent with interspecific transfer. We show extensive 2μ sequence variation in Saccharomyces cerevisiae and in other species of Saccharomyces and 2μ associations with S. cerevisiae populations, clinical origin and genotypes.

We characterized previously identified RNA viruses (L-A, L-BC, 20S, and 23S), L-A–dependent M satellites (M1, M2, M28, and Mlus), and M satellite–dependent killer phenotypes in the Saccharomyces cerevisiae 100-genomes genetic resource population. L-BC was present in all strains, albeit in 2 distinct levels, L-BChi and L-BClo; the L-BC level is associated with the L-BC genotype. L-BChi, L-A, 20S, 23S, M1, M2, and Mlus (M28 was absent) were in fewer strains than the similarly inherited 2µ plasmid. Novel L-A–dependent phenotypes were identified. Ten M+ strains exhibited M satellite–dependent killing (K+) of at least 1 of the naturally M0 and cured M0 derivatives of the 100-genomes strains; in these M0 strains, sensitivities to K1+, K2+, and K28+ strains varied. Finally, to complement our M satellite–encoded killer toxin analysis, we assembled the chromosomal KHS1 and KHR1 killer genes and used naturally M0 and cured M0 derivatives of the 100-genomes strains to assess and characterize the chromosomal killer phenotypes.

Mitochondrial genome variation and its effects on phenotypes have been widely analyzed in higher eukaryotes but less so in the model eukaryote Saccharomyces cerevisiae. Here, we describe mitochondrial genome variation in 96 diverse S. cerevisiae strains and assess associations between mitochondrial genotype and phenotypes as well as nuclear-mitochondrial epistasis. We associate sensitivity to the ATP synthase inhibitor oligomycin with SNPs in the mitochondrially encoded ATP6 gene. We describe the use of iso-nuclear F1 pairs, the mitochondrial genome equivalent of reciprocal hemizygosity analysis, to identify and analyze mitochondrial genotype-dependent phenotypes. Using iso-nuclear F1 pairs, we analyze the oligomycin phenotype-ATP6 association and find extensive nuclear-mitochondrial epistasis. Similarly, in iso-nuclear F1 pairs, we identify many additional mitochondrial genotype-dependent respiration phenotypes, for which there was no association in the 96 strains, and again find extensive nuclear-mitochondrial epistasis that likely contributes to the lack of association in the 96 strains. Finally, in iso-nuclear F1 pairs, we identify novel mitochondrial genotype-dependent nonrespiration phenotypes: resistance to cycloheximide, ketoconazole, and copper. We discuss potential mechanisms and the implications of mitochondrial genotype and of nuclear-mitochondrial epistasis effects on respiratory and nonrespiratory quantitative traits.

An important issue in genome evolution is the mechanism by which tandem duplications are generated from single-copy genes. In the yeast Saccharomyces cerevisiae, most strains contain tandemly duplicated copies of CUP1, a gene that encodes a copper-binding metallothionein. By screening 101 natural isolates of S. cerevisiae, we identified five different types of CUP1-containing repeats, as well as strains that only had one copy of CUP1. A comparison of the DNA sequences of these strains indicates that the CUP1 tandem arrays were generated by unequal nonhomologous recombination events from strains that had one CUP1 gene.

Assembled genome sequences are being generated at an exponential rate. Here we present FCS-GX, part of NCBI's Foreign Contamination Screen (FCS) tool suite, optimized to identify and remove contaminant sequences in new genomes. FCS-GX screens most genomes in 0.1-10 minutes. Testing FCS-GX on artificially fragmented genomes demonstrates sensitivity >95% for diverse contaminant species and specificity >99.93%. We used FCS-GX to screen 1.6 million GenBank assemblies and identified 36.8 Gbp of contamination (0.16% of total bases), with half from 161 assemblies. We updated assemblies in NCBI RefSeq to reduce detected contamination to 0.01% of bases. FCS-GX is available at https://github.com/ncbi/fcs/.

Transmembrane proteins and multi-domain proteins together make up more than 80% of the total proteins in any eukaryotic proteome. Therefore accurately classifying such proteins into functional classes is an important task. Furthermore, understanding the molecular evolution of multi-domain proteins is important because it shows how various domains fuse to form more complex proteins, and acquire new functions possibly affecting the organismal level of evolution. In this thesis, I first investigated the performance of several protein classifiers using one of the most divergent transmembrane protein families, the G-protein-coupled receptor (GPCR) superfamily, as an example. Alignment-free classifiers based on support vector machines using simple amino acid compositions were effective in remote-similarity detection even from short fragmented sequences. While a support vector machine using local pairwise-alignment scores showed very well-balanced performance, profile hidden Markov models were generally highly specific and well suited for classifying well-established protein family members. We suggested that different types of protein classifiers should be applied to gain the optimal mining power. Including some of these methods, combinations of multiple protein classification methods were applied to identify especially divergent plant GPCRs (or seven-transmembrane receptors) from the Arabidopsis thaliana genome. We identified 394 proteins as the candidates and provided a prioritized list including 54 proteins for further investigation. For multi-domain protein families, the distribution of urea amidolyase, urea carboxylase, and sterol-sensing domain (SSD) proteins across kingdoms was investigated. Molecular evolutionary analysis showed that the urea amidolyase genes currently found only in fungi among eukaryotes are the results of a horizontal gene transfer event from proteobacteria. Urea carboxylase genes currently found in fungi and other limited organisms were also likely derived from another ancestral gene in bacteria. Finally we showed the possibility of the bacterial origin of the eukaryotic SSD-containing proteins and that these ancestral sequences evolved into four different SSD-containing proteins acquiring specific functions. Two groups of SSD-containing proteins seemed to have been formed before the divergence of fungal and metazoan lineages by domain acquisition.