Identification and classification of papain-like cysteine proteinases.

Papain-like cysteine peptidases form a big and highly diverse superfamily of proteins involved in many important biological functions, such as protein turnover, deubiquitination, tissue remodeling, blood clotting, virulence, defense, and cell wall remodeling. High sequence and structure diversity observed within these proteins hinders their comprehensive classification as well as the identification of new representatives. Moreover, in general protein databases, many families already classified as papain like lack details regarding their mechanism of action or biological function. Here, we use transitive remote homology searches and 3D modeling to newly classify 21 families to the papain-like cysteine peptidase superfamily. We attempt to predict their biological function and provide structural characterization of 89 protein clusters defined based on sequence similarity altogether spanning 106 papain-like families. Moreover, we systematically discuss observed diversity in sequences, structures, and catalytic sites. Eventually, we expand the list of human papain-related proteins by seven representatives, including dopamine receptor-interacting protein 1 as potential deubiquitinase, and centriole duplication regulating CEP76 as retaining catalytically active peptidase-like domain. The presented results not only provide structure-based rationales to already existing peptidase databases but also may inspire further experimental research focused on peptidase-related biological processes.

Toward an understanding of the conformational plasticity of S100A8 and S100A9 Ca2+-binding proteins.

S100A8 and S100A9 are small, human, Ca2+-binding proteins with multiple intracellular and extracellular functions in signaling, regulation, and defense. The two proteins are not detected as monomers but form various noncovalent homo- or hetero-oligomers related to specific activities in human physiology. Because of their significant roles in numerous medical conditions, there has been intense research on the conformational properties of various S100A8 and S100A9 proteoforms as essential targets of drug discovery. NMR or crystal structures are currently available only for mutated or truncated protein complexes, mainly with bound metal ions, that may well reflect the proteins' properties outside cells but not in other biological contexts in which they perform. Here, we used structural mass spectrometry methods combined with molecular dynamics simulations to compare the conformations of wildtype full-length S100A8 and S100A9 subunits in biologically relevant homo- and heterodimers and in higher oligomers formed in the presence of calcium or zinc ions. We provide, first, rationales for their functional response to changing environmental conditions, by elucidating differences between proteoforms in flexible protein regions that may provide the plasticity of the binding sites for the multiple targets, and second, the key factors contributing to the variable stability of the oligomers. The described methods and a systematic view of the conformational properties of S100A8 and S100A9 complexes provide a basis for further research to characterize and modulate their functions for basic science and therapies.

Unique Properties of the Alpha-Helical DNA-Binding Protein KfrA Encoded by the IncU Incompatibility Group Plasmid RA3 and Its Host-Dependent Role in Plasmid Maintenance.

KfrA, encoded on the broad-host-range RA3 plasmid, is an alpha-helical DNA-binding protein that acts as a transcriptional autoregulator. The KfrARA3 operator site overlaps the kfrA promoter and is composed of five 9-bp direct repeats (DRs). Here, the biological properties of KfrA were studied using both in vivo and in vitro approaches. Localization of the DNA-binding helix-turn-helix motif (HTH) was mapped to the N29-R52 region by protein structure modeling and confirmed by alanine scanning. KfrA repressor ability depended on the number and orientation of DRs in the operator, as well as the ability of the protein to oligomerize. The long alpha-helical tail from residues 54 to 355 was shown to be involved in self-interactions, whereas the region from residue 54 to 177 was involved in heterodimerization with KfrC, another RA3-encoded alpha-helical protein. KfrA also interacted with the segrosome proteins IncC (ParA) and KorB (ParB), representatives of the class Ia active partition systems. Deletion of the kfr genes from the RA3 stability module decreased the plasmid retention in diverse hosts in a species-dependent manner. The specific interactions of KfrA with DNA are essential not only for the transcriptional regulatory function but also for the accessory role of KfrA in stable plasmid maintenance.IMPORTANCE Alpha-helical coiled-coil KfrA-type proteins are encoded by various broad-host-range low-copy-number conjugative plasmids. The DNA-binding protein KfrA encoded on the RA3 plasmid, a member of the IncU incompatibility group, oligomerizes, forms a complex with another plasmid-encoded, alpha-helical protein, KfrC, and interacts with the segrosome proteins IncC and KorB. The unique mode of KfrA dimer binding to the repetitive operator is required for a KfrA role in the stable maintenance of RA3 plasmid in distinct hosts.

Searching for Biological Function of the Mysterious PA2504 Protein from Pseudomonas aeruginosa.

For nearly half of the proteome of an important pathogen, Pseudomonas aeruginosa, the function has not yet been recognised. Here, we characterise one such mysterious protein PA2504, originally isolated by us as a sole partner of the RppH RNA hydrolase involved in transcription regulation of multiple genes. This study aims at elucidating details of PA2504 function and discussing its implications for bacterial biology. We show that PA2504 forms homodimers and is evenly distributed in the cytoplasm of bacterial cells. Molecular modelling identified the presence of a Tudor-like domain in PA2504. Transcriptomic analysis of a ΔPA2504 mutant showed that 42 transcripts, mainly coding for proteins involved in sulphur metabolism, were affected by the lack of PA2504. In vivo crosslinking of cellular proteins in the exponential and stationary phase of growth revealed several polypeptides that bound to PA2504 exclusively in the stationary phase. Mass spectrometry analysis identified them as the 30S ribosomal protein S4, the translation elongation factor TufA, and the global response regulator GacA. These results indicate that PA2504 may function as a tether for several important cellular factors.

Genomic characterization of Parengyodontium americanum sp. nov.

Modern genome analysis and phylogenomic methods have increased the number of fungal species, as well as enhanced appreciation of the degree of diversity within the fungal kingdom. In this context, we describe a new Parengyodontium species, P. americanum, which is phylogenetically related to the opportunistic human fungal pathogen P. album. Five unusual fungal isolates were recovered from five unique and confirmed coccidioidomycosis patients, and these isolates were subsequently submitted to detailed molecular and morphological identification procedures to determine identity. Molecular and morphological diagnostic analyses showed that the isolates belong to the Cordycipitaceae. Subsequently, three representative genomes were sequenced and annotated, and a new species, P. americanum, was identified. Using various genomic analyses, gene family expansions related to novel compounds and potential for ability to grow in diverse habitats are predicted. A general description of the genomic composition of this newly described species and comparison of genome content with Beauveria bassiana, Isaria fumosorosea and Cordyceps militaris shows a shared core genome of 6371 genes, and 148 genes that appear to be specific for P. americanum. This work provides the framework for future investigations of this interesting fungal species.

Comprehensive classification of the PIN domain-like superfamily.

PIN-like domains constitute a widespread superfamily of nucleases, diverse in terms of the reaction mechanism, substrate specificity, biological function and taxonomic distribution. Proteins with PIN-like domains are involved in central cellular processes, such as DNA replication and repair, mRNA degradation, transcription regulation and ncRNA maturation. In this work, we identify and classify the most complete set of PIN-like domains to provide the first comprehensive analysis of sequence-structure-function relationships within the whole PIN domain-like superfamily. Transitive sequence searches using highly sensitive methods for remote homology detection led to the identification of several new families, including representatives of Pfam (DUF1308, DUF4935) and CDD (COG2454), and 23 other families not classified in the public domain databases. Further sequence clustering revealed relationships between individual sequence clusters and showed heterogeneity within some families, suggesting a possible functional divergence. With five structural groups, 70 defined clusters, over 100,000 proteins, and broad biological functions, the PIN domain-like superfamily constitutes one of the largest and most diverse nuclease superfamilies. Detailed analyses of sequences and structures, domain architectures, and genomic contexts allowed us to predict biological function of several new families, including new toxin-antitoxin components, proteins involved in tRNA/rRNA maturation and transcription/translation regulation.

Systematic classification of the His-Me finger superfamily.

The His-Me finger endonucleases, also known as HNH or ßßα-metal endonucleases, form a large and diverse protein superfamily. The His-Me finger domain can be found in proteins that play an essential role in cells, including genome maintenance, intron homing, host defense and target offense. Its overall structural compactness and non-specificity make it a perfectly-tailored pathogenic module that participates on both sides of inter- and intra-organismal competition. An extremely low sequence similarity across the superfamily makes it difficult to identify and classify new His-Me fingers. Using state-of-the-art distant homology detection methods, we provide an updated and systematic classification of His-Me finger proteins. In this work, we identified over 100 000 proteins and clustered them into 38 groups, of which three groups are new and cannot be found in any existing public domain database of protein families. Based on an analysis of sequences, structures, domain architectures, and genomic contexts, we provide a careful functional annotation of the poorly characterized members of this superfamily. Our results may inspire further experimental investigations that should address the predicted activity and clarify the potential substrates, to provide more detailed insights into the fundamental biological roles of these proteins.

Poly-Saturated Dolichols from Filamentous Fungi Modulate Activity of Dolichol-Dependent Glycosyltransferase and Physical Properties of Membranes.

 Mono-saturated polyprenols (dolichols) have been found in almost all Eukaryotic cells, however, dolichols containing additional saturated bonds at the ω-end, have been identified in A. fumigatus and A. niger. Here we confirm using an LC-ESI-QTOF-MS analysis, that poly-saturated dolichols are abundant in other filamentous fungi, Trichoderma reesei, A. nidulans and Neurospora crassa, while the yeast Saccharomyces cerevisiae only contains the typical mono-saturated dolichols. We also show, using differential scanning calorimetry (DSC) and fluorescence anisotropy of 1,6-diphenyl-l,3,5-hexatriene (DPH) that the structure of dolichols modulates the properties of membranes and affects the functioning of dolichyl diphosphate mannose synthase (DPMS). The activity of this enzyme from T. reesei and S. cerevisiae was strongly affected by the structure of dolichols. Additionally, the structure of phosphatidylcholine (PC) and phosphatidylethanolamine (PE) model membranes was more strongly disturbed by the poly-saturated dolichols from Trichoderma than by the mono-saturated dolichols from yeast. By comparing the lipidome of filamentous fungi with that from S. cerevisiae, we revealed significant differences in the PC/PE ratio and fatty acids composition. Filamentous fungi differ from S. cerevisiae in the lipid composition of their membranes and the structure of dolichols. The structure of dolichols profoundly affects the functioning of dolichol-dependent enzyme, DPMS.

FAM46 proteins are novel eukaryotic non-canonical poly(A) polymerases.

FAM46 proteins, encoded in all known animal genomes, belong to the nucleotidyltransferase (NTase) fold superfamily. All four human FAM46 paralogs (FAM46A, FAM46B, FAM46C, FAM46D) are thought to be involved in several diseases, with FAM46C reported as a causal driver of multiple myeloma; however, their exact functions remain unknown. By using a combination of various bioinformatics analyses (e.g. domain architecture, cellular localization) and exhaustive literature and database searches (e.g. expression profiles, protein interactors), we classified FAM46 proteins as active non-canonical poly(A) polymerases, which modify cytosolic and/or nuclear RNA 3' ends. These proteins may thus regulate gene expression and probably play a critical role during cell differentiation. A detailed analysis of sequence and structure diversity of known NTases possessing PAP/OAS1 SBD domain, combined with state-of-the-art comparative modelling, allowed us to identify potential active site residues responsible for catalysis and substrate binding. We also explored the role of single point mutations found in human cancers and propose that FAM46 genes may be involved in the development of other major malignancies including lung, colorectal, hepatocellular, head and neck, urothelial, endometrial and renal papillary carcinomas and melanoma. Identification of these novel enzymes taking part in RNA metabolism in eukaryotes may guide their further functional studies.

The RNase H-like superfamily: new members, comparative structural analysis and evolutionary classification.

Ribonuclease H-like (RNHL) superfamily, also called the retroviral integrase superfamily, groups together numerous enzymes involved in nucleic acid metabolism and implicated in many biological processes, including replication, homologous recombination, DNA repair, transposition and RNA interference. The RNHL superfamily proteins show extensive divergence of sequences and structures. We conducted database searches to identify members of the RNHL superfamily (including those previously unknown), yielding >60 000 unique domain sequences. Our analysis led to the identification of new RNHL superfamily members, such as RRXRR (PF14239), DUF460 (PF04312, COG2433), DUF3010 (PF11215), DUF429 (PF04250 and COG2410, COG4328, COG4923), DUF1092 (PF06485), COG5558, OrfB_IS605 (PF01385, COG0675) and Peptidase_A17 (PF05380). Based on the clustering analysis we grouped all identified RNHL domain sequences into 152 families. Phylogenetic studies revealed relationships between these families, and suggested a possible history of the evolution of RNHL fold and its active site. Our results revealed clear division of the RNHL superfamily into exonucleases and endonucleases. Structural analyses of features characteristic for particular groups revealed a correlation between the orientation of the C-terminal helix with the exonuclease/endonuclease function and the architecture of the active site. Our analysis provides a comprehensive picture of sequence-structure-function relationships in the RNHL superfamily that may guide functional studies of the previously uncharacterized protein families.

Identification of a novel human mitochondrial endo-/exonuclease Ddk1/c20orf72 necessary for maintenance of proper 7S DNA levels.

Although the human mitochondrial genome has been investigated for several decades, the proteins responsible for its replication and expression, especially nucleolytic enzymes, are poorly described. Here, we characterized a novel putative PD-(D/E)XK nuclease encoded by the human C20orf72 gene named Ddk1 for its predicted catalytic residues. We show that Ddk1 is a mitochondrially localized metal-dependent DNase lacking detectable ribonuclease activity. Ddk1 degrades DNA mainly in a 3'-5' direction with a strong preference for single-stranded DNA. Interestingly, Ddk1 requires free ends for its activity and does not degrade circular substrates. In addition, when a chimeric RNA-DNA substrate is provided, Ddk1 can slide over the RNA fragment and digest DNA endonucleolytically. Although the levels of the mitochondrial DNA are unchanged on RNAi-mediated depletion of Ddk1, the mitochondrial single-stranded DNA molecule (7S DNA) accumulates. On the other hand, overexperssion of Ddk1 decreases the levels of 7S DNA, suggesting an important role of the protein in 7S DNA regulation. We propose a structural model of Ddk1 and discuss its similarity to other PD-(D/E)XK superfamily members.

Sequence, structure and functional diversity of PD-(D/E)XK phosphodiesterase superfamily.

Proteins belonging to PD-(D/E)XK phosphodiesterases constitute a functionally diverse superfamily with representatives involved in replication, restriction, DNA repair and tRNA-intron splicing. Their malfunction in humans triggers severe diseases, such as Fanconi anemia and Xeroderma pigmentosum. To date there have been several attempts to identify and classify new PD-(D/E)KK phosphodiesterases using remote homology detection methods. Such efforts are complicated, because the superfamily exhibits extreme sequence and structural divergence. Using advanced homology detection methods supported with superfamily-wide domain architecture and horizontal gene transfer analyses, we provide a comprehensive reclassification of proteins containing a PD-(D/E)XK domain. The PD-(D/E)XK phosphodiesterases span over 21,900 proteins, which can be classified into 121 groups of various families. Eleven of them, including DUF4420, DUF3883, DUF4263, COG5482, COG1395, Tsp45I, HaeII, Eco47II, ScaI, HpaII and Replic_Relax, are newly assigned to the PD-(D/E)XK superfamily. Some groups of PD-(D/E)XK proteins are present in all domains of life, whereas others occur within small numbers of organisms. We observed multiple horizontal gene transfers even between human pathogenic bacteria or from Prokaryota to Eukaryota. Uncommon domain arrangements greatly elaborate the PD-(D/E)XK world. These include domain architectures suggesting regulatory roles in Eukaryotes, like stress sensing and cell-cycle regulation. Our results may inspire further experimental studies aimed at identification of exact biological functions, specific substrates and molecular mechanisms of reactions performed by these highly diverse proteins.

Human telomerase model shows the role of the TEN domain in advancing the double helix for the next polymerization step.

Telomerases constitute a group of specialized ribonucleoprotein enzymes that remediate chromosomal shrinkage resulting from the "end-replication" problem. Defects in telomere length regulation are associated with several diseases as well as with aging and cancer. Despite significant progress in understanding the roles of telomerase, the complete structure of the human telomerase enzyme bound to telomeric DNA remains elusive, with the detailed molecular mechanism of telomere elongation still unknown. By application of computational methods for distant homology detection, comparative modeling, and molecular docking, guided by available experimental data, we have generated a three-dimensional structural model of a partial telomerase elongation complex composed of three essential protein domains bound to a single-stranded telomeric DNA sequence in the form of a heteroduplex with the template region of the human RNA subunit, TER. This model provides a structural mechanism for the processivity of telomerase and offers new insights into elongation. We conclude that the RNADNA heteroduplex is constrained by the telomerase TEN domain through repeated extension cycles and that the TEN domain controls the process by moving the template ahead one base at a time by translation and rotation of the double helix. The RNA region directly following the template can bind complementarily to the newly synthesized telomeric DNA, while the template itself is reused in the telomerase active site during the next reaction cycle. This first structural model of the human telomerase enzyme provides many details of the molecular mechanism of telomerase and immediately provides an important target for rational drug design.

An Adenosine Triphosphate- Dependent 5'-3' DNA Helicase From sk1-Like Lactococcus lactis F13 Phage.

Here, we describe functional characterization of an early gene (gp46) product of a virulent Lactococcus lactis sk1-like phage, vB_Llc_bIBBF13 (abbr. F13). The GP46 F13 protein carries a catalytically active RecA-like domain belonging to the P-loop NTPase superfamily. It also retains features characteristic for ATPases forming oligomers. In order to elucidate its detailed molecular function, we cloned and overexpressed the gp46 gene in Escherichia coli. Purified GP46 F13 protein binds to DNA and exhibits DNA unwinding activity on branched substrates in the presence of adenosine triphosphate (ATP). Size exclusion chromatography with multi-angle light scattering (SEC-MALS) experiments demonstrate that GP46 F13 forms oligomers, and further pull-down assays show that GP46 F13 interacts with host proteins involved in replication (i.e., DnaK, DnaJ, topoisomerase I, and single-strand binding protein). Taking together the localization of the gene and the obtained results, GP46 F13 is the first protein encoded in the early-expressed gene region with helicase activity that has been identified among lytic L. lactis phages up to date.

The Genome of Rhyzopertha dominica (Fab.) (Coleoptera: Bostrichidae): Adaptation for Success.

The lesser grain borer, Rhyzopertha dominica (F.) (Coleoptera: Bostrichidae), is a major global pest of cereal grains. Infestations are difficult to control as larvae feed inside grain kernels, and many populations are resistant to both contact insecticides and fumigants. We sequenced the genome of R. dominica to identify genes responsible for important biological functions and develop more targeted and efficacious management strategies. The genome was assembled from long read sequencing and long-range scaffolding technologies. The genome assembly is 479.1 Mb, close to the predicted genome size of 480.4 Mb by flow cytometry. This assembly is among the most contiguous beetle assemblies published to date, with 139 scaffolds, an N50 of 53.6 Mb, and L50 of 4, indicating chromosome-scale scaffolds. Predicted genes from biologically relevant groups were manually annotated using transcriptome data from adults and different larval tissues to guide annotation. The expansion of carbohydrase and serine peptidase genes suggest that they combine to enable efficient digestion of cereal proteins. A reduction in the copy number of several detoxification gene families relative to other coleopterans may reflect the low selective pressure on these genes in an insect that spends most of its life feeding internally. Chemoreceptor genes contain elevated numbers of pseudogenes for odorant receptors that also may be related to the recent ontogenetic shift of R. dominica to a diet consisting primarily of stored grains. Analysis of repetitive sequences will further define the evolution of bostrichid beetles compared to other species. The data overall contribute significantly to coleopteran genetic research.

Utilization of cobalamin is ubiquitous in early-branching fungal phyla.

Cobalamin is a cofactor present in essential metabolic pathways in animals and one of the water-soluble vitamins. It is a complex compound synthesized solely by prokaryotes. Cobalamin dependence is scattered across the tree of life. In particular, fungi and plants were deemed devoid of cobalamin. We demonstrate that cobalamin is utilized by all non-Dikarya fungi lineages. This observation is supported by the genomic presence of both B12-dependent enzymes and cobalamin modifying enzymes. Fungal cobalamin-dependent enzymes are highly similar to their animal homologs. Phylogenetic analyses support a scenario of vertical inheritance of the cobalamin usage with several losses. Cobalamin usage was probably lost in Mucorinae and at the base of Dikarya which groups most of the model organisms and which hindered B12-dependent metabolism discovery in fungi. Our results indicate that cobalamin dependence was a widely distributed trait at least in Opisthokonta, across diverse microbial eukaryotes and was likely present in the LECA.

Expanding Diversity of Firmicutes Single-Strand Annealing Proteins: A Putative Role of Bacteriophage-Host Arms Race.

Bacteriophage-encoded single strand annealing proteins (SSAPs) are recombinases which can substitute the classical, bacterial RecA and manage the DNA metabolism at different steps of phage propagation. SSAPs have been shown to efficiently promote recombination between short and rather divergent DNA sequences and were exploited for in vivo genetic engineering mainly in Gram-negative bacteria. In opposition to the conserved and almost universal bacterial RecA protein, SSAPs display great sequence diversity. The importance for SSAPs in phage biology and phage-bacteria evolution is underlined by their role as key players in events of horizontal gene transfer (HGT). All of the above provoke a constant interest for the identification and study of new phage recombinase proteins in vivo, in vitro as well as in silico. Despite this, a huge body of putative ssap genes escapes conventional classification, as they are not properly annotated. In this work, we performed a wide-scale identification, classification and analysis of SSAPs encoded by the Firmicutes bacteria and their phages. By using sequence similarity network and gene context analyses, we created a new high quality dataset of phage-related SSAPs, substantially increasing the number of annotated SSAPs. We classified the identified SSAPs into seven distinct families, namely RecA, Gp2.5, RecT/Redß, Erf, Rad52/22, Sak3, and Sak4, organized into three superfamilies. Analysis of the relationships between the revealed protein clusters led us to recognize Sak3-like proteins as a new distinct SSAP family. Our analysis showed an irregular phylogenetic distribution of ssap genes among different bacterial phyla and specific phages, which can be explained by the high rates of ssap HGT. We propose that the evolution of phage recombinases could be tightly linked to the dissemination of bacterial phage-resistance mechanisms (e.g., abortive infection and CRISPR/Cas systems) targeting ssap genes and be a part of the constant phage-bacteria arms race.

A model of full-length RAGE in complex with S100B.

The receptor for advanced glycation end products (RAGE) is an immunoglobulin-type multiligand transmembrane protein expressed in numerous cell types, including the central nervous system cells. RAGE interaction with S100B, released during brain tissue damage, leads to RAGE upregulation and initialization of a spiral proinflammatory associated with different neural disorders. Here, we present the structural characterization of the hetero-oligomeric complex of the full-length RAGE with S100B, obtained by a combination of mass spectrometry-based methods and molecular modeling. We predict that RAGE functions as a tightly packed tetramer exposing a positively charged surface formed by V domains for S100B binding. Based on HDX results we demonstrate an allosteric coupling of the distal extracellular V domains and the transmembrane region, indicating a possible mechanism of signal transmission by RAGE across the membrane. Our model provides an insight into RAGE-ligand interactions, providing a basis for the rational design of the therapeutic modifiers of its activity.

Metabolic Potential, Ecology and Presence of Associated Bacteria Is Reflected in Genomic Diversity of Mucoromycotina.

Mucoromycotina are often considered mainly in pathogenic context but their biology remains understudied. We describe the genomes of six Mucoromycotina fungi representing distant saprotrophic lineages within the subphylum (i.e., Umbelopsidales and Mucorales). We selected two Umbelopsis isolates from soil (i.e., U. isabellina, U. vinacea), two soil-derived Mucor isolates (i.e., M. circinatus, M. plumbeus), and two Mucorales representatives with extended proteolytic activity (i.e., Thamnidium elegans and Mucor saturninus). We complement computational genome annotation with experimental characteristics of their digestive capabilities, cell wall carbohydrate composition, and extensive total lipid profiles. These traits inferred from genome composition, e.g., in terms of identified encoded enzymes, are in accordance with experimental results. Finally, we link the presence of associated bacteria with observed characteristics. Thamnidium elegans genome harbors an additional, complete genome of an associated bacterium classified to Paenibacillus sp. This fungus displays multiple altered traits compared to the remaining isolates, regardless of their evolutionary distance. For instance, it has expanded carbon assimilation capabilities, e.g., efficiently degrades carboxylic acids, and has a higher diacylglycerol:triacylglycerol ratio and skewed phospholipid composition which suggests a more rigid cellular membrane. The bacterium can complement the host enzymatic capabilities, alter the fungal metabolism, cell membrane composition but does not change the composition of the cell wall of the fungus. Comparison of early-diverging Umbelopsidales with evolutionary younger Mucorales points at several subtle differences particularly in their carbon source preferences and encoded carbohydrate repertoire. Nevertheless, all tested Mucoromycotina share features including the ability to produce 18:3 gamma-linoleic acid, use TAG as the storage lipid and have fucose as a cell wall component.

Transposable elements contribute to fungal genes and impact fungal lifestyle.

The last decade brought a still growing experimental evidence of mobilome impact on host's gene expression. We systematically analysed genomic location of transposable elements (TEs) in 625 publicly available fungal genomes from the NCBI database in order to explore their potential roles in genome evolution and correlation with species' lifestyle. We found that non-autonomous TEs and remnant copies are evenly distributed across genomes. In consequence, they also massively overlap with regions annotated as genes, which suggests a great contribution of TE-derived sequences to host's coding genome. Younger and potentially active TEs cluster with one another away from genic regions. This non-randomness is a sign of either selection against insertion of TEs in gene proximity or target site preference among some types of TEs. Proteins encoded by genes with old transposable elements insertions have significantly less repeat and protein-protein interaction motifs but are richer in enzymatic domains. However, genes only proximal to TEs do not display any functional enrichment. Our findings show that adaptive cases of TE insertion remain a marginal phenomenon, and the overwhelming majority of TEs are evolving neutrally. Eventually, animal-related and pathogenic fungi have more TEs inserted into genes than fungi with other lifestyles. This is the first systematic, kingdom-wide study concerning mobile elements and their genomic neighbourhood. The obtained results should inspire further research concerning the roles TEs played in evolution and how they shape the life we know today.