The TPLATE adaptor complex drives clathrin-mediated endocytosis in plants.

Clathrin-mediated endocytosis is the major mechanism for eukaryotic plasma membrane-based proteome turn-over. In plants, clathrin-mediated endocytosis is essential for physiology and development, but the identification and organization of the machinery operating this process remains largely obscure. Here, we identified an eight-core-component protein complex, the TPLATE complex, essential for plant growth via its role as major adaptor module for clathrin-mediated endocytosis. This complex consists of evolutionarily unique proteins that associate closely with core endocytic elements. The TPLATE complex is recruited as dynamic foci at the plasma membrane preceding recruitment of adaptor protein complex 2, clathrin, and dynamin-related proteins. Reduced function of different complex components severely impaired internalization of assorted endocytic cargoes, demonstrating its pivotal role in clathrin-mediated endocytosis. Taken together, the TPLATE complex is an early endocytic module representing a unique evolutionary plant adaptation of the canonical eukaryotic pathway for clathrin-mediated endocytosis.

Closing the gap: Oxford Nanopore Technologies R10 sequencing allows comparable results to Illumina sequencing for SNP-based outbreak investigation of bacterial pathogens.

Whole-genome sequencing has become the method of choice for bacterial outbreak investigation, with most clinical and public health laboratories currently routinely using short-read Illumina sequencing. Recently, long-read Oxford Nanopore Technologies (ONT) sequencing has gained prominence and may offer advantages over short-read sequencing, particularly with the recent introduction of the R10 chemistry, which promises much lower error rates than the R9 chemistry. However, limited information is available on its performance for bacterial single-nucleotide polymorphism (SNP)-based outbreak investigation. We present an open-source workflow, Prokaryotic Awesome variant Calling Utility (PACU) (https://github.com/BioinformaticsPlatformWIV-ISP/PACU), for constructing SNP phylogenies using Illumina and/or ONT R9/R10 sequencing data. The workflow was evaluated using outbreak data sets of Shiga toxin-producing Escherichia coli and Listeria monocytogenes by comparing ONT R9 and R10 with Illumina data. The performance of each sequencing technology was evaluated not only separately but also by integrating samples sequenced by different technologies/chemistries into the same phylogenomic analysis. Additionally, the minimum sequencing time required to obtain accurate phylogenetic results using nanopore sequencing was evaluated. PACU allowed accurate identification of outbreak clusters for both species using all technologies/chemistries, but ONT R9 results deviated slightly more from the Illumina results. ONT R10 results showed trends very similar to Illumina, and we found that integrating data sets sequenced by either Illumina or ONT R10 for different isolates into the same analysis produced stable and highly accurate phylogenomic results. The resulting phylogenies for these two outbreaks stabilized after ~20 hours of sequencing for ONT R9 and ~8 hours for ONT R10. This study provides a proof of concept for using ONT R10, either in isolation or in combination with Illumina, for rapid and accurate bacterial SNP-based outbreak investigation.

Comparison of 6 DNA extraction methods for isolation of high yield of high molecular weight DNA suitable for shotgun metagenomics Nanopore sequencing to detect bacteria.

BACKGROUND: Oxford Nanopore Technologies (ONT) offers an accessible platform for long-read sequencing, which improves the reconstruction of genomes and helps to resolve complex genomic contexts, especially in the case of metagenome analysis. To take the best advantage of long-read sequencing, DNA extraction methods must be able to isolate pure high molecular weight (HMW) DNA from complex metagenomics samples, without introducing any bias. New methods released on the market, and protocols developed at the research level, were specifically designed for this application and need to be assessed. RESULTS: In this study, with different bacterial cocktail mixes, analyzed as pure or spiked in a synthetic fecal matrix, we evaluated the performances of 6 DNA extraction methods using various cells lysis and purification techniques, from quick and easy, to more time-consuming and gentle protocols, including a portable method for on-site application. In addition to the comparison of the quality, quantity and purity of the extracted DNA, the performance obtained when doing Nanopore sequencing on a MinION flow cell was also tested. From the obtained results, the Quick-DNA HMW MagBead Kit (Zymo Research) was selected as producing the best yield of pure HMW DNA. Furthermore, this kit allowed an accurate detection, by Nanopore sequencing, of almost all the bacterial species present in a complex mock community. CONCLUSION: Amongst the 6 tested methods, the Quick-DNA HMW MagBead Kit (Zymo Research) was considered as the most suitable for Nanopore sequencing and would be recommended for bacterial metagenomics studies using this technology.

The genome of the seagrass Zostera marina reveals angiosperm adaptation to the sea.

Seagrasses colonized the sea on at least three independent occasions to form the basis of one of the most productive and widespread coastal ecosystems on the planet. Here we report the genome of Zostera marina (L.), the first, to our knowledge, marine angiosperm to be fully sequenced. This reveals unique insights into the genomic losses and gains involved in achieving the structural and physiological adaptations required for its marine lifestyle, arguably the most severe habitat shift ever accomplished by flowering plants. Key angiosperm innovations that were lost include the entire repertoire of stomatal genes, genes involved in the synthesis of terpenoids and ethylene signalling, and genes for ultraviolet protection and phytochromes for far-red sensing. Seagrasses have also regained functions enabling them to adjust to full salinity. Their cell walls contain all of the polysaccharides typical of land plants, but also contain polyanionic, low-methylated pectins and sulfated galactans, a feature shared with the cell walls of all macroalgae and that is important for ion homoeostasis, nutrient uptake and O2/CO2 exchange through leaf epidermal cells. The Z. marina genome resource will markedly advance a wide range of functional ecological studies from adaptation of marine ecosystems under climate warming, to unravelling the mechanisms of osmoregulation under high salinities that may further inform our understanding of the evolution of salt tolerance in crop plants.

A Bioinformatics Whole-Genome Sequencing Workflow for Clinical Mycobacterium tuberculosis Complex Isolate Analysis, Validated Using a Reference Collection Extensively Characterized with Conventional Methods and In Silico Approaches.

The use of whole-genome sequencing (WGS) for routine typing of bacterial isolates has increased substantially in recent years. For Mycobacterium tuberculosis (MTB), in particular, WGS has the benefit of drastically reducing the time required to generate results compared to most conventional phenotypic methods. Consequently, a multitude of solutions for analyzing WGS MTB data have been developed, but their successful integration in clinical and national reference laboratories is hindered by the requirement for their validation, for which a consensus framework is still largely absent. We developed a bioinformatics workflow for (Illumina) WGS-based routine typing of MTB complex (MTBC) member isolates allowing complete characterization, including (sub)species confirmation and identification (16S, csb/RD, hsp65), single nucleotide polymorphism (SNP)-based antimicrobial resistance (AMR) prediction, and pathogen typing (spoligotyping, SNP barcoding, and core genome multilocus sequence typing). Workflow performance was validated on a per-assay basis using a collection of 238 in-house-sequenced MTBC isolates, extensively characterized with conventional molecular biology-based approaches supplemented with public data. For SNP-based AMR prediction, results from molecular genotyping methods were supplemented with in silico modified data sets, allowing us to greatly increase the set of evaluated mutations. The workflow demonstrated very high performance with performance metrics of >99% for all assays, except for spoligotyping, where sensitivity dropped to ∼90%. The validation framework for our WGS-based bioinformatics workflow can aid in the standardization of bioinformatics tools by the MTB community and other SNP-based applications regardless of the targeted pathogen(s). The bioinformatics workflow is available for academic and nonprofit use through the Galaxy instance of our institute at https://galaxy.sciensano.be.

Large diversity of linezolid-resistant isolates discovered in food-producing animals through linezolid selective monitoring in Belgium in 2019.

BACKGROUND: Linezolid is a critically important antibiotic used to treat human infections caused by MRSA and VRE. While linezolid is not licensed for food-producing animals, linezolid-resistant (LR) isolates have been reported in European countries, including Belgium. OBJECTIVES: To: (i) assess LR occurrence in staphylococci and enterococci isolated from different Belgian food-producing animals in 2019 through selective monitoring; and (ii) investigate the genomes and relatedness of these isolates. METHODS: Faecal samples (n = 1325) and nasal swab samples (n = 148) were analysed with a protocol designed to select LR bacteria, including a 44-48 h incubation period. The presence of LR chromosomal mutations, transferable LR genes and their genetic organizations and other resistance genes, as well as LR isolate relatedness (from this study and the NCBI database) were assessed through WGS. RESULTS: The LR rate differed widely between animal host species, with the highest rates occurring in nasal samples from pigs and sows (25.7% and 20.5%, respectively) and faecal samples from veal calves (16.4%). WGS results showed that LR determinants are present in a large diversity of isolates circulating in the agricultural sector, with some isolates closely related to human isolates, posing a human health risk. CONCLUSIONS: LR dedicated monitoring with WGS analysis could help to better understand the spread of LR. Cross-selection of LR transferable genes through other antibiotic use should be considered in future action plans aimed at combatting antimicrobial resistance and in future objectives for the rational use of antibiotics in a One Health perspective.

Evaluation of the added value of viral genomic information for predicting severity of influenza infection.

BACKGROUND: The severity of an influenza infection is influenced by both host and viral characteristics. This study aims to assess the relevance of viral genomic data for the prediction of severe influenza A(H3N2) infections among patients hospitalized for severe acute respiratory infection (SARI), in view of risk assessment and patient management. METHODS: 160 A(H3N2) influenza positive samples from the 2016-2017 season originating from the Belgian SARI surveillance were selected for whole genome sequencing. Predictor variables for severity were selected using a penalized elastic net logistic regression model from a combined host and genomic dataset, including patient information and nucleotide mutations identified in the viral genome. The goodness-of-fit of the model combining host and genomic data was compared using a likelihood-ratio test with the model including host data only. Internal validation of model discrimination was conducted by calculating the optimism-adjusted area under the Receiver Operating Characteristic curve (AUC) for both models. RESULTS: The model including viral mutations in addition to the host characteristics had an improved fit ([Formula: see text]=12.03, df = 3, p = 0.007). The optimism-adjusted AUC increased from 0.671 to 0.732. CONCLUSIONS: Adding genomic data (selected season-specific mutations in the viral genome) to the model containing host characteristics improved the prediction of severe influenza infection among hospitalized SARI patients, thereby offering the potential for translation into a prospective strategy to perform early season risk assessment or to guide individual patient management.

First detection of a plasmid-encoded New-Delhi metallo-beta-lactamase-1 (NDM-1) producing Acinetobacter baumannii using whole genome sequencing, isolated in a clinical setting in Benin.

BACKGROUND: Carbapenem-resistant Acinetobacter baumannii is considered a top priority pathogen by the World Health Organization for combatting increasing antibiotic resistance and development of new drugs. Since it was originally reported in Klebsiella pneumoniae in 2009, the quick spread of the blaNDM-1 gene encoding a New-Delhi metallo-beta-lactamase-1 (NDM-1) is increasingly recognized as a serious threat. This gene is usually carried by large plasmids and has already been documented in diverse bacterial species, including A. baumannii. Here, we report the first detection of a NDM-1-producing A. baumannii strain isolated in Benin. CASE PRESENTATION: A 31-year-old woman was admitted to a surgical unit with a diagnosis of post-cesarean hematoma. An extensively-drug resistant A. baumannii strain solely susceptible to amikacin, colistin and ciprofloxacin, and resistant to several other antibiotics including ceftazidime, imipenem, meropenem, gentamicin, tobramycin, ceftazidime/avibactam, and sulfamethoxazole-trimethoprim, was isolated from the wound. Production of NDM-1 was demonstrated by immunochromatographic testing. Whole genome sequencing of the isolate confirmed the presence of blaNDM-1, but also antibiotic resistance genes against multiple beta-lactamases and other classes of antibiotics, in addition to several virulence genes. Moreover, the blaNDM-1 gene was found to be present in a Tn125 transposon integrated on a plasmid. CONCLUSIONS: The discovery of this extensively-drug resistant A. baumannii strain carrying blaNDM-1 in Benin is worrying, especially because of its high potential risk of horizontal gene transfer due to being integrated into a transposon located on a plasmid. Strict control and prevention measures should be taken, once NDM-1 positive A. baumannii has been identified to prevent transfer of this resistance gene to other Enterobacterales. Capacity building is required by governmental agencies to provide suitable antibiotic treatment options and strategies, in combination with strengthening laboratory services for detection and surveillance of this pathogen.

Isolation of Drug-Resistant Gallibacterium anatis from Calves with Unresponsive Bronchopneumonia, Belgium.

Gallibacterium anatis is an opportunistic pathogen, previously associated with deaths in poultry, domestic birds, and occasionally humans. We obtained G. anatis isolates from bronchoalveolar lavage samples of 10 calves with bronchopneumonia unresponsive to antimicrobial therapy. Collected isolates were multidrug-resistant to extensively drug-resistant, exhibiting resistance against 5-7 classes of antimicrobial drugs. Whole-genome sequencing revealed 24 different antimicrobial-resistance determinants, including genes not previously described in the Gallibacterium genus or even the Pasteurellaceae family, such as aadA23, blaCARB-8, tet(Y), and qnrD1. Some resistance genes were closely linked in resistance gene cassettes with either transposases in close proximity or situated on putative mobile elements or predicted plasmids. Single-nucleotide polymorphism genotyping revealed large genetic variation between the G. anatis isolates, including isolates retrieved from the same farm. G. anatis might play a hitherto unrecognized role as a respiratory pathogen and resistance gene reservoir in cattle and has unknown zoonotic potential.

Increased viral read counts and metagenomic full genome characterization of porcine astrovirus 4 and Posavirus 1 in sows in a swine farm with unexplained neonatal piglet diarrhea.

Neonatal diarrhea in piglets may cause major losses in affected pig herds. The present study used random high-throughput RNA sequencing (metagenomic next generation sequencing, mNGS) to investigate the virome of sows from a farm with persistent neonatal piglet diarrhea in comparison to two control farms without diarrhea problems. A variety of known swine gastrointestinal viruses was detected in the control farms as well as in the problem farm (Mamastrovirus, Enterovirus, Picobirnavirus, Posavirus 1, Kobuvirus, Proprismacovirus). A substantial increase in normalized viral read counts was observed in the affected farm compared to the control farms. The increase was attributable to a single viral species in each of the sampled sows (porcine astrovirus 4 and Posavirus 1). The complete genomes of a porcine astrovirus 4 and two co-infecting Posavirus 1 were de novo assembled and characterized. The 6734 nt single-stranded RNA genome of porcine astrovirus 4 (PoAstV-4) strain Belgium/2019 contains three overlapping open reading frames (nonstructural protein 1ab, nonstructural protein 1a, capsid protein). Posavirus 1 strains Belgium/01/2019 and Belgium/02/2019 have a 9814 nt single-stranded positive-sense RNA genome encoding a single open reading frame (polyprotein precursor) containing the five expected Picornavirales-conserved protein domains. The study highlights the potential of mNGS workflows to study unexplained neonatal diarrhea in piglets and contributes to the scarce availability of both PoAstV-4 and Posavirus-1 whole genome sequences from Western Europe.

The genome of Eucalyptus grandis.

Eucalypts are the world's most widely planted hardwood trees. Their outstanding diversity, adaptability and growth have made them a global renewable resource of fibre and energy. We sequenced and assembled >94% of the 640-megabase genome of Eucalyptus grandis. Of 36,376 predicted protein-coding genes, 34% occur in tandem duplications, the largest proportion thus far in plant genomes. Eucalyptus also shows the highest diversity of genes for specialized metabolites such as terpenes that act as chemical defence and provide unique pharmaceutical oils. Genome sequencing of the E. grandis sister species E. globulus and a set of inbred E. grandis tree genomes reveals dynamic genome evolution and hotspots of inbreeding depression. The E. grandis genome is the first reference for the eudicot order Myrtales and is placed here sister to the eurosids. This resource expands our understanding of the unique biology of large woody perennials and provides a powerful tool to accelerate comparative biology, breeding and biotechnology.

Use of Whole Genome Sequencing Data for a First in Silico Specificity Evaluation of the RT-qPCR Assays Used for SARS-CoV-2 Detection.

The current COronaVIrus Disease 2019 (COVID-19) pandemic started in December 2019. COVID-19 cases are confirmed by the detection of SARS-CoV-2 RNA in biological samples by RT-qPCR. However, limited numbers of SARS-CoV-2 genomes were available when the first RT-qPCR methods were developed in January 2020 for initial in silico specificity evaluation and to verify whether the targeted loci are highly conserved. Now that more whole genome data have become available, we used the bioinformatics tool SCREENED and a total of 4755 publicly available SARS-CoV-2 genomes, downloaded at two different time points, to evaluate the specificity of 12 RT-qPCR tests (consisting of a total of 30 primers and probe sets) used for SARS-CoV-2 detection and the impact of the virus' genetic evolution on four of them. The exclusivity of these methods was also assessed using the human reference genome and 2624 closely related other respiratory viral genomes. The specificity of the assays was generally good and stable over time. An exception is the first method developed by the China Center for Disease Control and prevention (CDC), which exhibits three primer mismatches present in 358 SARS-CoV-2 genomes sequenced mainly in Europe from February 2020 onwards. The best results were obtained for the assay of Chan et al. (2020) targeting the gene coding for the spiking protein (S). This demonstrates that our user-friendly strategy can be used for a first in silico specificity evaluation of future RT-qPCR tests, as well as verifying that the former methods are still capable of detecting circulating SARS-CoV-2 variants.

The Norway spruce genome sequence and conifer genome evolution.

Conifers have dominated forests for more than 200 million years and are of huge ecological and economic importance. Here we present the draft assembly of the 20-gigabase genome of Norway spruce (Picea abies), the first available for any gymnosperm. The number of well-supported genes (28,354) is similar to the >100 times smaller genome of Arabidopsis thaliana, and there is no evidence of a recent whole-genome duplication in the gymnosperm lineage. Instead, the large genome size seems to result from the slow and steady accumulation of a diverse set of long-terminal repeat transposable elements, possibly owing to the lack of an efficient elimination mechanism. Comparative sequencing of Pinus sylvestris, Abies sibirica, Juniperus communis, Taxus baccata and Gnetum gnemon reveals that the transposable element diversity is shared among extant conifers. Expression of 24-nucleotide small RNAs, previously implicated in transposable element silencing, is tissue-specific and much lower than in other plants. We further identify numerous long (>10,000 base pairs) introns, gene-like fragments, uncharacterized long non-coding RNAs and short RNAs. This opens up new genomic avenues for conifer forestry and breeding.

TARGETING THE 16S RRNA GENE FOR BACTERIAL IDENTIFICATION IN COMPLEX MIXED SAMPLES: COMPARATIVE EVALUATION OF SECOND (ILLUMINA) AND THIRD (OXFORD NANOPORE TECHNOLOGIES) GENERATION SEQUENCING TECHNOLOGIES.

Rapid, accurate bacterial identification in biological samples is an important task for microbiology laboratories, for which 16S~rRNA gene Sanger sequencing of cultured isolates is frequently used. In contrast, next-generation sequencing does not require intermediate culturing steps and can be directly applied on communities, but its performance has not been extensively evaluated. We present a comparative evaluation of second (Illumina) and third (Oxford Nanopore Technologies (ONT)) generation sequencing technologies for 16S targeted genomics using a well-characterized reference sample. Different 16S gene regions were amplified and sequenced using the Illumina MiSeq, and analyzed with Mothur. Correct classification was variable, depending on the region amplified. Using a majority vote over all regions, most false positives could be eliminated at the genus level but not the species level. Alternatively, the entire 16S gene was amplified and sequenced using the ONT MinION, and analyzed with Mothur, EPI2ME, and GraphMap. Although >99\% of reads were correctly classified at the genus level, up to $\approx$40\% were misclassified at the species level. Both~technologies, therefore, allow reliable identification of bacterial genera, but can potentially misguide identification of bacterial species, and constitute viable alternatives to Sanger sequencing for rapid analysis of mixed samples without requiring any culturing steps.

Application of whole genome data for in silico evaluation of primers and probes routinely employed for the detection of viral species by RT-qPCR using dengue virus as a case study.

BACKGROUND: Viral infection by dengue virus is a major public health problem in tropical countries. Early diagnosis and detection are increasingly based on quantitative reverse transcriptase real-time polymerase chain reaction (RT-qPCR) directed against genomic regions conserved between different isolates. Genetic variation can however result in mismatches of primers and probes with their targeted nucleic acid regions. Whole genome sequencing allows to characterize and track such changes, which in turn enables to evaluate, optimize, and (re-)design novel and existing RT-qPCR methods. The immense amount of available sequence data renders this however a labour-intensive and complex task. RESULTS: We present a bioinformatics approach that enables in silico evaluation of primers and probes intended for routinely employed RT-qPCR methods. This approach is based on analysing large amounts of publically available whole genome data, by first employing BLASTN to mine the genomic regions targeted by the RT-qPCR method(s), and afterwards using BLASTN-SHORT to evaluate whether primers and probes will anneal based on a set of simple in silico criteria. Using dengue virus as a case study, we evaluated 18 published RT-qPCR methods using more than 3000 publically available genomes in the NCBI Virus Variation Resource, and provide a systematic overview of method performance based on in silico sensitivity and specificity. CONCLUSIONS: We provide a comprehensive overview of dengue virus RT-qPCR method performance that will aid appropriate method selection allowing to take specific measures that aim to contain and prevent viral spread in afflicted regions. Notably, we find that primer-template mismatches at their 3' end may represent a general issue for dengue virus RT-qPCR detection methods that merits more attention in their development process. Our approach is also available as a public tool, and demonstrates how utilizing genomic data can provide meaningful insights in an applied public health setting such as the detection of viral species in human diagnostics.

Isolation of Burkholderia pseudomallei from a Pet Green Iguana, Belgium.

We isolated Burkholderia pseudomallei, the causative agent of melioidosis, from liver granulomas of a pet green iguana (Iguana iguana) in Belgium. This case highlights a risk for imported green iguanas acting as a reservoir for introduction of this high-threat, zoonotic pathogen into nonendemic regions.

Horsetails Are Ancient Polyploids: Evidence from Equisetum giganteum.

Horsetails represent an enigmatic clade within the land plants. Despite consisting only of one genus (Equisetum) that contains 15 species, they are thought to represent the oldest extant genus within the vascular plants dating back possibly as far as the Triassic. Horsetails have retained several ancient features and are also characterized by a particularly high chromosome count (n = 108). Whole-genome duplications (WGDs) have been uncovered in many angiosperm clades and have been associated with the success of angiosperms, both in terms of species richness and biomass dominance, but remain understudied in nonangiosperm clades. Here, we report unambiguous evidence of an ancient WGD in the fern lineage, based on sequencing and de novo assembly of an expressed gene catalog (transcriptome) from the giant horsetail (Equisetum giganteum). We demonstrate that horsetails underwent an independent paleopolyploidy during the Late Cretaceous prior to the diversification of the genus but did not experience any recent polyploidizations that could account for their high chromosome number. We also discuss the specific retention of genes following the WGD and how this may be linked to their long-term survival.

Detection of Plasmid-Mediated Colistin Resistance, mcr-1 and mcr-2 Genes, in Salmonella spp. Isolated from Food at Retail in Belgium from 2012 to 2015.

A collection of 105 colistin-resistant Salmonella isolates collected from 2012 to 2015 in the national surveillance program in Belgium was screened by PCR for the presence of genes mcr-1 and mcr-2. Of these, 1.90% (2/105) and 0.95% (1/105) tested positive for mcr-1 and mcr-2, respectively. The presence of the mcr-1 or mcr-2 determinant has been confirmed by whole genome sequencing and allowed the localization of these two genes on IncX4 type plasmids. We report here the presence of mcr-1 and the first mcr-2 gene in Salmonella ever isolated in the Belgian food chain. Although present at retail since 2012, the occurrence is low and sporadic.

Analysis of 41 plant genomes supports a wave of successful genome duplications in association with the Cretaceous-Paleogene boundary.

Ancient whole-genome duplications (WGDs), also referred to as paleopolyploidizations, have been reported in most evolutionary lineages. Their attributed role remains a major topic of discussion, ranging from an evolutionary dead end to a road toward evolutionary success, with evidence supporting both fates. Previously, based on dating WGDs in a limited number of plant species, we found a clustering of angiosperm paleopolyploidizations around the Cretaceous-Paleogene (K-Pg) extinction event about 66 million years ago. Here we revisit this finding, which has proven controversial, by combining genome sequence information for many more plant lineages and using more sophisticated analyses. We include 38 full genome sequences and three transcriptome assemblies in a Bayesian evolutionary analysis framework that incorporates uncorrelated relaxed clock methods and fossil uncertainty. In accordance with earlier findings, we demonstrate a strongly nonrandom pattern of genome duplications over time with many WGDs clustering around the K-Pg boundary. We interpret these results in the context of recent studies on invasive polyploid plant species, and suggest that polyploid establishment is promoted during times of environmental stress. We argue that considering the evolutionary potential of polyploids in light of the environmental and ecological conditions present around the time of polyploidization could mitigate the stark contrast in the proposed evolutionary fates of polyploids.

Reconstruction of ancestral metabolic enzymes reveals molecular mechanisms underlying evolutionary innovation through gene duplication.

Gene duplications are believed to facilitate evolutionary innovation. However, the mechanisms shaping the fate of duplicated genes remain heavily debated because the molecular processes and evolutionary forces involved are difficult to reconstruct. Here, we study a large family of fungal glucosidase genes that underwent several duplication events. We reconstruct all key ancestral enzymes and show that the very first preduplication enzyme was primarily active on maltose-like substrates, with trace activity for isomaltose-like sugars. Structural analysis and activity measurements on resurrected and present-day enzymes suggest that both activities cannot be fully optimized in a single enzyme. However, gene duplications repeatedly spawned daughter genes in which mutations optimized either isomaltase or maltase activity. Interestingly, similar shifts in enzyme activity were reached multiple times via different evolutionary routes. Together, our results provide a detailed picture of the molecular mechanisms that drove divergence of these duplicated enzymes and show that whereas the classic models of dosage, sub-, and neofunctionalization are helpful to conceptualize the implications of gene duplication, the three mechanisms co-occur and intertwine.