Genomic evidence for the widespread presence of GH45 cellulases among soil invertebrates.

Lignocellulose is a major component of vascular plant biomass. Its decomposition is crucial for the terrestrial carbon cycle. Microorganisms are considered primary decomposers, but evidence increases that some invertebrates may also decompose lignocellulose. We investigated the taxonomic distribution and evolutionary origins of GH45 hydrolases, important enzymes for the decomposition of cellulose and hemicellulose, in a collection of soil invertebrate genomes. We found that these genes are common in springtails and oribatid mites. Phylogenetic analysis revealed that cellulase genes were acquired early in the evolutionary history of these groups. Domain architectures and predicted 3D enzyme structures indicate that these cellulases are functional. Patterns of presence and absence of these genes across different lineages prompt further investigation into their evolutionary and ecological benefits. The ubiquity of cellulase genes suggests that soil invertebrates may play a role in lignocellulose decomposition, independently or in synergy with microorganisms. Understanding the ecological and evolutionary implications might be crucial for understanding soil food webs and the carbon cycle.

The MetaInvert soil invertebrate genome resource provides insights into below-ground biodiversity and evolution.

Soil invertebrates are among the least understood metazoans on Earth. Thus far, the lack of taxonomically broad and dense genomic resources has made it hard to thoroughly investigate their evolution and ecology. With MetaInvert we provide draft genome assemblies for 232 soil invertebrate species, representing 14 common groups and 94 families. We show that this data substantially extends the taxonomic scope of DNA- or RNA-based taxonomic identification. Moreover, we confirm that theories of genome evolution cannot be generalised across evolutionarily distinct invertebrate groups. The soil invertebrate genomes presented here will support the management of soil biodiversity through molecular monitoring of community composition and function, and the discovery of evolutionary adaptations to the challenges of soil conditions.

Sequencing the Genomes of the First Terrestrial Fungal Lineages: What Have We Learned?

The first genome sequenced of a eukaryotic organism was for Saccharomyces cerevisiae, as reported in 1996, but it was more than 10 years before any of the zygomycete fungi, which are the early-diverging terrestrial fungi currently placed in the phyla Mucoromycota and Zoopagomycota, were sequenced. The genome for Rhizopus delemar was completed in 2008; currently, more than 1000 zygomycete genomes have been sequenced. Genomic data from these early-diverging terrestrial fungi revealed deep phylogenetic separation of the two major clades-primarily plant-associated saprotrophic and mycorrhizal Mucoromycota versus the primarily mycoparasitic or animal-associated parasites and commensals in the Zoopagomycota. Genomic studies provide many valuable insights into how these fungi evolved in response to the challenges of living on land, including adaptations to sensing light and gravity, development of hyphal growth, and co-existence with the first terrestrial plants. Genome sequence data have facilitated studies of genome architecture, including a history of genome duplications and horizontal gene transfer events, distribution and organization of mating type loci, rDNA genes and transposable elements, methylation processes, and genes useful for various industrial applications. Pathogenicity genes and specialized secondary metabolites have also been detected in soil saprobes and pathogenic fungi. Novel endosymbiotic bacteria and viruses have been discovered during several zygomycete genome projects. Overall, genomic information has helped to resolve a plethora of research questions, from the placement of zygomycetes on the evolutionary tree of life and in natural ecosystems, to the applied biotechnological and medical questions.

FUBP1 is a general splicing factor facilitating 3' splice site recognition and splicing of long introns.

Splicing of pre-mRNAs critically contributes to gene regulation and proteome expansion in eukaryotes, but our understanding of the recognition and pairing of splice sites during spliceosome assembly lacks detail. Here, we identify the multidomain RNA-binding protein FUBP1 as a key splicing factor that binds to a hitherto unknown cis-regulatory motif. By collecting NMR, structural, and in vivo interaction data, we demonstrate that FUBP1 stabilizes U2AF2 and SF1, key components at the 3' splice site, through multivalent binding interfaces located within its disordered regions. Transcriptional profiling and kinetic modeling reveal that FUBP1 is required for efficient splicing of long introns, which is impaired in cancer patients harboring FUBP1 mutations. Notably, FUBP1 interacts with numerous U1 snRNP-associated proteins, suggesting a unique role for FUBP1 in splice site bridging for long introns. We propose a compelling model for 3' splice site recognition of long introns, which represent 80% of all human introns.

Feature architecture aware phylogenetic profiling indicates a functional diversification of type IVa pili in the nosocomial pathogen Acinetobacter baumannii.

The Gram-negative bacterial pathogen Acinetobacter baumannii is a major cause of hospital-acquired opportunistic infections. The increasing spread of pan-drug resistant strains makes A. baumannii top-ranking among the ESKAPE pathogens for which novel routes of treatment are urgently needed. Comparative genomics approaches have successfully identified genetic changes coinciding with the emergence of pathogenicity in Acinetobacter. Genes that are prevalent both in pathogenic and a-pathogenic Acinetobacter species were not considered ignoring that virulence factors may emerge by the modification of evolutionarily old and widespread proteins. Here, we increased the resolution of comparative genomics analyses to also include lineage-specific changes in protein feature architectures. Using type IVa pili (T4aP) as an example, we show that three pilus components, among them the pilus tip adhesin ComC, vary in their Pfam domain annotation within the genus Acinetobacter. In most pathogenic Acinetobacter isolates, ComC displays a von Willebrand Factor type A domain harboring a finger-like protrusion, and we provide experimental evidence that this finger conveys virulence-related functions in A. baumannii. All three genes are part of an evolutionary cassette, which has been replaced at least twice during A. baumannii diversification. The resulting strain-specific differences in T4aP layout suggests differences in the way how individual strains interact with their host. Our study underpins the hypothesis that A. baumannii uses T4aP for host infection as it was shown previously for other pathogens. It also indicates that many more functional complexes may exist whose precise functions have been adjusted by modifying individual components on the domain level.

ncOrtho: efficient and reliable identification of miRNA orthologs.

MicroRNAs (miRNAs) are post-transcriptional regulators that finetune gene expression via translational repression or degradation of their target mRNAs. Despite their functional relevance, frameworks for the scalable and accurate detection of miRNA orthologs are missing. Consequently, there is still no comprehensive picture of how miRNAs and their associated regulatory networks have evolved. Here we present ncOrtho, a synteny informed pipeline for the targeted search of miRNA orthologs in unannotated genome sequences. ncOrtho matches miRNA annotations from multi-tissue transcriptomes in precision, while scaling to the analysis of hundreds of custom-selected species. The presence-absence pattern of orthologs to 266 human miRNA families across 402 vertebrate species reveals four bursts of miRNA acquisition, of which the most recent event occurred in the last common ancestor of higher primates. miRNA families are rarely modified or lost, but notable exceptions for both events exist. miRNA co-ortholog numbers faithfully indicate lineage-specific whole genome duplications, and miRNAs are powerful markers for phylogenomic analyses. Their exceptionally low genetic diversity makes them suitable to resolve clades where the phylogenetic signal is blurred by incomplete lineage sorting of ancestral alleles. In summary, ncOrtho allows to routinely consider miRNAs in evolutionary analyses that were thus far reserved to protein-coding genes.

FAS: assessing the similarity between proteins using multi-layered feature architectures.

MOTIVATION: Protein sequence comparison is a fundamental element in the bioinformatics toolkit. When sequences are annotated with features such as functional domains, transmembrane domains, low complexity regions or secondary structure elements, the resulting feature architectures allow better informed comparisons. However, many existing schemes for scoring architecture similarities cannot cope with features arising from multiple annotation sources. Those that do fall short in the resolution of overlapping and redundant feature annotations. RESULTS: Here, we introduce FAS, a scoring method that integrates features from multiple annotation sources in a directed acyclic architecture graph. Redundancies are resolved as part of the architecture comparison by finding the paths through the graphs that maximize the pair-wise architecture similarity. In a large-scale evaluation on more than 10 000 human-yeast ortholog pairs, architecture similarities assessed with FAS are consistently more plausible than those obtained using e-values to resolve overlaps or leaving overlaps unresolved. Three case studies demonstrate the utility of FAS on architecture comparison tasks: benchmarking of orthology assignment software, identification of functionally diverged orthologs, and diagnosing protein architecture changes stemming from faulty gene predictions. With the help of FAS, feature architecture comparisons can now be routinely integrated into these and many other applications. AVAILABILITY AND IMPLEMENTATION: FAS is available as python package: https://pypi.org/project/greedyFAS/.

Endothelial versus pronephron fate decision is modulated by the transcription factors Cloche/Npas4l, Tal1, and Lmo2.

Endothelial specification is a key event during embryogenesis; however, when, and how, endothelial cells separate from other lineages is poorly understood. In zebrafish, Npas4l is indispensable for endothelial specification by inducing the expression of the transcription factor genes etsrp, tal1, and lmo2. We generated a knock-in reporter in zebrafish npas4l to visualize endothelial progenitors and their derivatives in wild-type and mutant embryos. Unexpectedly, we find that in npas4l mutants, npas4l reporter-expressing cells contribute to the pronephron tubules. Single-cell transcriptomics and live imaging of the early lateral plate mesoderm in wild-type embryos indeed reveals coexpression of endothelial and pronephron markers, a finding confirmed by creERT2-based lineage tracing. Increased contribution of npas4l reporter-expressing cells to pronephron tubules is also observed in tal1 and lmo2 mutants and is reversed in npas4l mutants injected with tal1 mRNA. Together, these data reveal that Npas4l/Tal1/Lmo2 regulate the fate decision between the endothelial and pronephron lineages.

Evolutionarily stable gene clusters shed light on the common grounds of pathogenicity in the Acinetobacter calcoaceticus-baumannii complex.

Nosocomial pathogens of the Acinetobacter calcoaceticus-baumannii (ACB) complex are a cautionary example for the world-wide spread of multi- and pan-drug resistant bacteria. Aiding the urgent demand for novel therapeutic targets, comparative genomics studies between pathogens and their apathogenic relatives shed light on the genetic basis of human-pathogen interaction. Yet, existing studies are limited in taxonomic scope, sensing of the phylogenetic signal, and resolution by largely analyzing genes independent of their organization in functional gene clusters. Here, we explored more than 3,000 Acinetobacter genomes in a phylogenomic framework integrating orthology-based phylogenetic profiling and microsynteny conservation analyses. We delineate gene clusters in the type strain A. baumannii ATCC 19606 whose evolutionary conservation indicates a functional integration of the subsumed genes. These evolutionarily stable gene clusters (ESGCs) reveal metabolic pathways, transcriptional regulators residing next to their targets but also tie together sub-clusters with distinct functions to form higher-order functional modules. We shortlisted 150 ESGCs that either co-emerged with the pathogenic ACB clade or are preferentially found therein. They provide a high-resolution picture of genetic and functional changes that coincide with the manifestation of the pathogenic phenotype in the ACB clade. Key innovations are the remodeling of the regulatory-effector cascade connecting LuxR/LuxI quorum sensing via an intermediate messenger to biofilm formation, the extension of micronutrient scavenging systems, and the increase of metabolic flexibility by exploiting carbon sources that are provided by the human host. We could show experimentally that only members of the ACB clade use kynurenine as a sole carbon and energy source, a substance produced by humans to fine-tune the antimicrobial innate immune response. In summary, this study provides a rich and unbiased set of novel testable hypotheses on how pathogenic Acinetobacter interact with and ultimately infect their human host. It is a comprehensive resource for future research into novel therapeutic strategies.

InvL, an Invasin-Like Adhesin, Is a Type II Secretion System Substrate Required for Acinetobacter baumannii Uropathogenesis.

Acinetobacter baumannii is an opportunistic pathogen of growing concern, as isolates are commonly multidrug resistant. While A. baumannii is most frequently associated with pulmonary infections, a significant proportion of clinical isolates come from urinary sources, highlighting its uropathogenic potential. The type II secretion system (T2SS) of commonly used model Acinetobacter strains is important for virulence in various animal models, but the potential role of the T2SS in urinary tract infection (UTI) remains unknown. Here, we used a catheter-associated UTI (CAUTI) model to demonstrate that a modern urinary isolate, UPAB1, requires the T2SS for full virulence. A proteomic screen to identify putative UPAB1 T2SS effectors revealed an uncharacterized lipoprotein with structural similarity to the intimin-invasin family, which serve as type V secretion system (T5SS) adhesins required for the pathogenesis of several bacteria. This protein, designated InvL, lacked the ß-barrel domain associated with T5SSs but was confirmed to require the T2SS for both surface localization and secretion. This makes InvL the first identified T2SS effector belonging to the intimin-invasin family. InvL was confirmed to be an adhesin, as the protein bound to extracellular matrix components and mediated adhesion to urinary tract cell lines in vitro. Additionally, the invL mutant was attenuated in the CAUTI model, indicating a role in Acinetobacter uropathogenesis. Finally, bioinformatic analyses revealed that InvL is present in nearly all clinical isolates belonging to international clone 2, a lineage of significant clinical importance. In all, we conclude that the T2SS substrate InvL is an adhesin required for A. baumannii uropathogenesis. IMPORTANCE While pathogenic Acinetobacter can cause various infections, we recently found that 20% of clinical isolates come from urinary sources. Despite the clinical relevance of Acinetobacter as a uropathogen, few virulence factors involved in urinary tract colonization have been defined. Here, we identify a novel type II secretion system effector, InvL, which is required for full uropathogenesis by a modern urinary isolate. Although InvL has predicted structural similarity to the intimin-invasin family of autotransporter adhesins, InvL is predicted to be anchored to the membrane as a lipoprotein. Similar to other invasin homologs, however, we demonstrate that InvL is a bona fide adhesin capable of binding extracellular matrix components and mediating adhesion to urinary tract cell lines. In all, this work establishes InvL as an adhesin important for Acinetobacter's urinary tract virulence and represents the first report of a type II secretion system effector belonging to the intimin-invasin family.

Assembly and Functional Role of PACE Transporter PA2880 from Pseudomonas aeruginosa.

The recently identified proteobacterial antimicrobial compound efflux (PACE) transporters are multidrug transporters energized by the electrochemical gradient of protons. Here, we present the results of phylogenetic and functional studies on the PACE family transporter PA2880 from Pseudomonas aeruginosa. A phylogenetic analysis of the PACE family revealed that PA2880 and AceI from Acinetobacter baumannii are classified into evolutionarily distinct clades, although they both transport chlorhexidine. We demonstrate that PA2880 mainly exists as a dimer in solution, which is independent of pH, and its dimeric state is essential for its proper function. Electrogenicity studies revealed that the chlorhexidine/H+ antiport process is electrogenic. The function of several highly conserved residues was investigated. These findings provide further insights into the functional features of PACE family transporters, facilitating studies on their transport mechanisms. IMPORTANCE Pseudomonas aeruginosa is a pathogen that causes hospital-acquired (nosocomial) infections, such as ventilator-associated pneumonia and sepsis syndromes. Chlorhexidine diacetate is a disinfectant used for bacterial control in various environments potentially harboring P. aeruginosa. Therefore, investigation of the mechanism of the efflux of chlorhexidine mediated by PA2880, a PACE family transporter from P. aeruginosa, is of significance to combat bacterial infections. This study improves our understanding of the transport mechanism of PACE family transporters and will facilitate the effective utilization of chlorhexidine for P. aeruginosa control.

The structure of the Aquifex aeolicus MATE family multidrug resistance transporter and sequence comparisons suggest the existence of a new subfamily.

Multidrug and toxic compound extrusion (MATE) transporters are widespread in all domains of life. Bacterial MATE transporters confer multidrug resistance by utilizing an electrochemical gradient of H+ or Na+ to export xenobiotics across the membrane. Despite the availability of X-ray structures of several MATE transporters, a detailed understanding of the transport mechanism has remained elusive. Here we report the crystal structure of a MATE transporter from Aquifex aeolicus at 2.0-Å resolution. In light of its phylogenetic placement outside of the diversity of hitherto-described MATE transporters and the lack of conserved acidic residues, this protein may represent a subfamily of prokaryotic MATE transporters, which was proven by phylogenetic analysis. Furthermore, the crystal structure and substrate docking results indicate that the substrate binding site is located in the N bundle. The importance of residues surrounding this binding site was demonstrated by structure-based site-directed mutagenesis. We suggest that Aq_128 is functionally similar but structurally diverse from DinF subfamily transporters. Our results provide structural insights into the MATE transporter, which further advances our global understanding of this important transporter family.

Tracing Eukaryotic Ribosome Biogenesis Factors Into the Archaeal Domain Sheds Light on the Evolution of Functional Complexity.

Ribosome assembly is an essential and carefully choreographed cellular process. In eukaryotes, several 100 proteins, distributed across the nucleolus, nucleus, and cytoplasm, co-ordinate the step-wise assembly of four ribosomal RNAs (rRNAs) and approximately 80 ribosomal proteins (RPs) into the mature ribosomal subunits. Due to the inherent complexity of the assembly process, functional studies identifying ribosome biogenesis factors and, more importantly, their precise functions and interplay are confined to a few and very well-established model organisms. Although best characterized in yeast (Saccharomyces cerevisiae), emerging links to disease and the discovery of additional layers of regulation have recently encouraged deeper analysis of the pathway in human cells. In archaea, ribosome biogenesis is less well-understood. However, their simpler sub-cellular structure should allow a less elaborated assembly procedure, potentially providing insights into the functional essentials of ribosome biogenesis that evolved long before the diversification of archaea and eukaryotes. Here, we use a comprehensive phylogenetic profiling setup, integrating targeted ortholog searches with automated scoring of protein domain architecture similarities and an assessment of when search sensitivity becomes limiting, to trace 301 curated eukaryotic ribosome biogenesis factors across 982 taxa spanning the tree of life and including 727 archaea. We show that both factor loss and lineage-specific modifications of factor function modulate ribosome biogenesis, and we highlight that limited sensitivity of the ortholog search can confound evolutionary conclusions. Projecting into the archaeal domain, we find that only few factors are consistently present across the analyzed taxa, and lineage-specific loss is common. While members of the Asgard group are not special with respect to their inventory of ribosome biogenesis factors (RBFs), they unite the highest number of orthologs to eukaryotic RBFs in one taxon. Using large ribosomal subunit maturation as an example, we demonstrate that archaea pursue a simplified version of the corresponding steps in eukaryotes. Much of the complexity of this process evolved on the eukaryotic lineage by the duplication of ribosomal proteins and their subsequent functional diversification into ribosome biogenesis factors. This highlights that studying ribosome biogenesis in archaea provides fundamental information also for understanding the process in eukaryotes.

Minimized combinatorial CRISPR screens identify genetic interactions in autophagy.

Combinatorial CRISPR-Cas screens have advanced the mapping of genetic interactions, but their experimental scale limits the number of targetable gene combinations. Here, we describe 3Cs multiplexing, a rapid and scalable method to generate highly diverse and uniformly distributed combinatorial CRISPR libraries. We demonstrate that the library distribution skew is the critical determinant of its required screening coverage. By circumventing iterative cloning of PCR-amplified oligonucleotides, 3Cs multiplexing facilitates the generation of combinatorial CRISPR libraries with low distribution skews. We show that combinatorial 3Cs libraries can be screened with minimal coverages, reducing associated efforts and costs at least 10-fold. We apply a 3Cs multiplexing library targeting 12,736 autophagy gene combinations with 247,032 paired gRNAs in viability and reporter-based enrichment screens. In the viability screen, we identify, among others, the synthetic lethal WDR45B-PIK3R4 and the proliferation-enhancing ATG7-KEAP1 genetic interactions. In the reporter-based screen, we identify over 1,570 essential genetic interactions for autophagy flux, including interactions among paralogous genes, namely ATG2A-ATG2B, GABARAP-MAP1LC3B and GABARAP-GABARAPL2. However, we only observe few genetic interactions within paralogous gene families of more than two members, indicating functional compensation between them. This work establishes 3Cs multiplexing as a platform for genetic interaction screens at scale.

Ten Years of Collaborative Progress in the Quest for Orthologs.

Accurate determination of the evolutionary relationships between genes is a foundational challenge in biology. Homology-evolutionary relatedness-is in many cases readily determined based on sequence similarity analysis. By contrast, whether or not two genes directly descended from a common ancestor by a speciation event (orthologs) or duplication event (paralogs) is more challenging, yet provides critical information on the history of a gene. Since 2009, this task has been the focus of the Quest for Orthologs (QFO) Consortium. The sixth QFO meeting took place in Okazaki, Japan in conjunction with the 67th National Institute for Basic Biology conference. Here, we report recent advances, applications, and oncoming challenges that were discussed during the conference. Steady progress has been made toward standardization and scalability of new and existing tools. A feature of the conference was the presentation of a panel of accessible tools for phylogenetic profiling and several developments to bring orthology beyond the gene unit-from domains to networks. This meeting brought into light several challenges to come: leveraging orthology computations to get the most of the incoming avalanche of genomic data, integrating orthology from domain to biological network levels, building better gene models, and adapting orthology approaches to the broad evolutionary and genomic diversity recognized in different forms of life and viruses.

Tek (Tie2) is not required for cardiovascular development in zebrafish.

Angiopoietin/TIE signalling plays a major role in blood and lymphatic vessel development. In mouse, Tek (previously known as Tie2) mutants die prenatally due to a severely underdeveloped cardiovascular system. In contrast, in zebrafish, previous studies have reported that although embryos injected with tek morpholinos (MOs) exhibit severe vascular defects, tek mutants display no obvious vascular malformations. To further investigate the function of zebrafish Tek, we generated a panel of loss-of-function tek mutants, including RNA-less alleles, an allele lacking the MO-binding site, an in-frame deletion allele and a premature termination codon-containing allele. Our data show that all these mutants survive to adulthood with no obvious cardiovascular defects. MO injections into tek mutants lacking the MO-binding site or the entire tek locus cause similar vascular defects to those observed in MO-injected +/+ siblings, indicating off-target effects of the MOs. Surprisingly, comprehensive phylogenetic profiling and synteny analyses reveal that Tek was lost in the largest teleost clade, suggesting a lineage-specific shift in the function of TEK during vertebrate evolution. Altogether, these data show that Tek is dispensable for zebrafish development, and probably dispensable in most teleost species.

Peptidoglycan editing provides immunity to Acinetobacter baumannii during bacterial warfare.

Peptidoglycan (PG) is essential in most bacteria. Thus, it is often targeted by various assaults, including interbacterial attacks via the type VI secretion system (T6SS). Here, we report that the Gram-negative bacterium Acinetobacter baumannii strain ATCC 17978 produces, secretes, and incorporates the noncanonical d-amino acid d-lysine into its PG during stationary phase. We show that PG editing increases the competitiveness of A. baumannii during bacterial warfare by providing immunity against peptidoglycan-targeting T6SS effectors from various bacterial competitors. In contrast, we found that d-Lys production is detrimental to pathogenesis due, at least in part, to the activity of the human enzyme d-amino acid oxidase (DAO), which degrades d-Lys producing H2O2 toxic to bacteria. Phylogenetic analyses indicate that the last common ancestor of A. baumannii had the ability to produce d-Lys. However, this trait was independently lost multiple times, likely reflecting the evolution of A. baumannii as a human pathogen.

Adjustment of photosynthetic activity to drought and fluctuating light in wheat.

Drought is a major cause of losses in crop yield. Under field conditions, plants exposed to drought are usually also experiencing rapid changes in light intensity. Accordingly, plants need to acclimate to both, drought and light stress. Two crucial mechanisms in plant acclimation to changes in light conditions comprise thylakoid protein phosphorylation and dissipation of light energy as heat by non-photochemical quenching (NPQ). Here, we analyzed the acclimation efficacy of two different wheat varieties, by applying fluctuating light for analysis of plants, which had been subjected to a slowly developing drought stress as it usually occurs in the field. This novel approach allowed us to distinguish four drought phases, which are critical for grain yield, and to discover acclimatory responses which are independent of photodamage. In short-term, under fluctuating light, the slowdown of NPQ relaxation adjusts the photosynthetic activity to the reduced metabolic capacity. In long-term, the photosynthetic machinery acquires a drought-specific configuration by changing the PSII-LHCII phosphorylation pattern together with protein stoichiometry. Therefore, the fine-tuning of NPQ relaxation and PSII-LHCII phosphorylation pattern represent promising traits for future crop breeding strategies.

What Is in Umbilicaria pustulata? A Metagenomic Approach to Reconstruct the Holo-Genome of a Lichen.

Lichens are valuable models in symbiosis research and promising sources of biosynthetic genes for biotechnological applications. Most lichenized fungi grow slowly, resist aposymbiotic cultivation, and are poor candidates for experimentation. Obtaining contiguous, high-quality genomes for such symbiotic communities is technically challenging. Here, we present the first assembly of a lichen holo-genome from metagenomic whole-genome shotgun data comprising both PacBio long reads and Illumina short reads. The nuclear genomes of the two primary components of the lichen symbiosis-the fungus Umbilicaria pustulata (33 Mb) and the green alga Trebouxia sp. (53 Mb)-were assembled at contiguities comparable to single-species assemblies. The analysis of the read coverage pattern revealed a relative abundance of fungal to algal nuclei of ∼20:1. Gap-free, circular sequences for all organellar genomes were obtained. The bacterial community is dominated by Acidobacteriaceae and encompasses strains closely related to bacteria isolated from other lichens. Gene set analyses showed no evidence of horizontal gene transfer from algae or bacteria into the fungal genome. Our data suggest a lineage-specific loss of a putative gibberellin-20-oxidase in the fungus, a gene fusion in the fungal mitochondrion, and a relocation of an algal chloroplast gene to the algal nucleus. Major technical obstacles during reconstruction of the holo-genome were coverage differences among individual genomes surpassing three orders of magnitude. Moreover, we show that GC-rich inverted repeats paired with nonrandom sequencing error in PacBio data can result in missing gene predictions. This likely poses a general problem for genome assemblies based on long reads.