Whole-Genome Duplication and Yeast’s Fruitful Way of Life

The immediate metabolomic effects of whole-genome duplication in the greater duckweed, Spirodela polyrhiza

PREMISE

In plants, whole-genome duplication (WGD) is a common mutation with profound evolutionary potential. Given the costs associated with a superfluous genome copy, polyploid establishment is enigmatic. However, in the right environment, immediate phenotypic changes following WGD can facilitate establishment. Metabolite abundances are the direct output of the cell's regulatory network and determine much of the impact of environmental and genetic change on the phenotype. While it is well known that an increase in the bulk amount of genetic material can increase cell size, the impact of gene dosage multiplication on the metabolome remains largely unknown.

METHODS

We used untargeted metabolomics on four genetically distinct diploid-neoautotetraploid pairs of the greater duckweed, Spirodela polyrhiza, to investigate how WGD affects metabolite abundances per cell and per biomass.

RESULTS

Autopolyploidy increased metabolite levels per cell, but the response of individual metabolites varied considerably. However, the impact on metabolite level per biomass was restricted because the increased cell size reduced the metabolite concentration per cell. Nevertheless, we detected both quantitative and qualitative effects of WGD on the metabolome. Many effects were strain-specific, but some were shared by all four strains.

CONCLUSIONS

The nature and impact of metabolic changes after WGD depended strongly on the genotype. Dosage effects have the potential to alter the plant metabolome qualitatively and quantitatively, but were largely balanced out by the reduction in metabolite concentration due to an increase in cell size in this species.

The genetic and molecular basis of haploinsufficiency in flowering plants

In diploid organisms, haploinsufficiency can be defined as the requirement for more than one fully functional copy of a gene. In contrast to most genes, whose loss-of-function alleles are recessive, loss-of-function alleles of haploinsufficient genes are dominant. However, forward and reverse genetic screens are biased toward obtaining recessive, loss-of-function mutations, and therefore, dominant mutations of all types are underrepresented in mutant collections. Despite this underrepresentation, haploinsufficient loci have intriguing implications for studies of genome evolution, gene dosage, stability of protein complexes, genetic redundancy, and gene expression. Here we review examples of haploinsufficiency in flowering plants and describe the underlying molecular mechanisms and evolutionary forces driving haploinsufficiency. Finally, we discuss the masking of haploinsufficiency by genetic redundancy, a widespread phenomenon among angiosperms.

Proteome dynamics during transition from exponential to stationary phase under aerobic and anaerobic conditions in yeast

The yeast Saccharomyces cerevisiae is a widely used eukaryotic model organism and a promising cell factory for industry. However, despite decades of research, the regulation of its metabolism is not yet fully understood, and its complexity represents a major challenge for engineering and optimising biosynthetic routes. Recent studies have demonstrated the potential of resource and proteomic allocation data in enhancing models for metabolic processes. However, comprehensive and accurate proteome dynamics data that can be used for such approaches are still very limited. Therefore, we performed a quantitative proteome dynamics study to comprehensively cover the transition from exponential to stationary phase for both aerobically and anaerobically grown yeast cells. The combination of highly controlled reactor experiments, biological replicates and standardised sample preparation procedures ensured reproducibility and accuracy. Additionally, we selected the CEN.PK lineage for our experiments because of its relevance for both fundamental and applied research. Together with the prototrophic, standard haploid strain CEN.PK113-7D, we also investigated an engineered strain with genetic minimisation of the glycolytic pathway, resulting in the quantitative assessment of 54 proteomes. The anaerobic cultures showed remarkably less proteome-level changes compared to the aerobic cultures, during transition from the exponential to the stationary phase as a consequence of the lack of the diauxic shift in the absence of oxygen. These results support the notion that anaerobically growing cells lack resources to adequately adapt to starvation. This proteome dynamics study constitutes an important step towards better understanding of the impact of glucose exhaustion and oxygen on the complex proteome allocation process in yeast. Finally, the established proteome dynamics data provide a valuable resource for the development of resource allocation models as well as for metabolic engineering efforts.

Top-Down, Knowledge-Based Genetic Reduction of Yeast Central Carbon Metabolism

Saccharomyces cerevisiae, whose evolutionary past includes a whole-genome duplication event, is characterized by a mosaic genome configuration with substantial apparent genetic redundancy. This apparent redundancy raises questions about the evolutionary driving force for genomic fixation of “minor” paralogs and complicates modular and combinatorial metabolic engineering strategies. While isoenzymes might be important in specific environments, they could be dispensable in controlled laboratory or industrial contexts. The present study explores the extent to which the genetic complexity of the central carbon metabolism (CCM) in S. cerevisiae, here defined as the combination of glycolysis, the pentose phosphate pathway, the tricarboxylic acid cycle, and a limited number of related pathways and reactions, can be reduced by elimination of (iso)enzymes without major negative impacts on strain physiology. Cas9-mediated, groupwise deletion of 35 of the 111 genes yielded a “minimal CCM” strain which, despite the elimination of 32% of CCM-related proteins, showed only a minimal change in phenotype on glucose-containing synthetic medium in controlled bioreactor cultures relative to a congenic reference strain. Analysis under a wide range of other growth and stress conditions revealed remarkably few phenotypic changes from the reduction of genetic complexity. Still, a well-documented context-dependent role of GPD1 in osmotolerance was confirmed. The minimal CCM strain provides a model system for further research into genetic redundancy of yeast genes and a platform for strategies aimed at large-scale, combinatorial remodeling of yeast CCM.

Mitogenomics and mitochondrial gene phylogeny decipher the evolution of Saccharomycotina yeasts

Saccharomycotina yeasts belong to diverse clades within the kingdom of fungi and are important to human everyday life. This work investigates the evolutionary relationships among these yeasts from a mitochondrial (mt) genomic perspective. A comparative study of 155 yeast mt genomes representing all major phylogenetic lineages of Saccharomycotina was performed, including genome size and content variability, intron and intergenic regions’ diversity, genetic code alterations, and syntenic variation. Findings from this study suggest that mt genome size diversity is the result of a ceaseless random process, mainly based on genetic recombination and intron mobility. Gene order analysis revealed conserved syntenic units and many occurring rearrangements, which can be correlated with major evolutionary events as shown by the phylogenetic analysis of the concatenated mt protein matrix. For the first time, molecular dating indicated a slower mt genome divergence rate in the early stages of yeast evolution, in contrast with a faster rate in the late evolutionary stages, compared to their nuclear time divergence. Genetic code reassignments of mt genomes are a perpetual process happening in many different parallel evolutionary steps throughout the evolution of Saccharomycotina. Overall, this work shows that phylogenetic studies based on the mt genome of yeasts highlight major evolutionary events.

Two alpha isopropylmalate synthase isozymes with similar kinetic properties are extant in the yeast Lachancea kluyveri

The first committed step in the leucine biosynthetic pathway is catalyzed by α-isopropylmalate synthase (α-IPMS, EC 2.3.3.13), which in the Saccaromycotina subphylum of Ascomycete yeasts is frequently encoded by duplicated genes. Following a gene duplication event, the two copies may be preserved presumably because the encoded proteins diverge in either functional properties and/or cellular localization. The genome of the petite-negative budding yeast Lachancea kluyveri includes two SAKL0E10472 (LkLEU4) and SAKL0F05170g (LKLEU4BIS) paralogous genes, which are homologous to other yeast α-IPMS sequences. Here, we investigate whether these paralogous genes encode functional α-IPMS isozymes and whether their functions have diverged. Molecular phylogeny suggested that the LkLeu4 isozyme is located in the mitochondria and LkLeu4BIS in the cytosol. Comparison of growth rates, leucine intracellular pools and mRNA levels, indicate that the LkLeu4 isozyme is the predominant α-IPMS enzyme during growth on glucose as carbon source. Determination of the kinetic parameters indicates that the isozymes have similar affinities for the substrates and for the feedback inhibitor leucine. Thus, the diversification of the physiological roles of the genes LkLEU4 and LKLEU4BIS involves preferential transcription of the LkLEU4 gene during growth on glucose and different subcellular localization, although ligand interactions have not diverged.

Engineering and functional analysis of yeast with a monotypic 40S ribosome subunit

Ribosomes are not monolithic but dynamic machines composed of heterogeneous ribosomal protein (RP) paralogs with elusive functions. Isolation and characterization of monotypic ribosomes with homogeneous RP paralog compositions represent ideal approaches to understand the role of pervasive RP paralogs in customizing translation abilities but are largely hurdled by the complexity of the cellular ribosome pool (e.g., in Saccharomyces cerevisiae, 59 RP paralog pairs allow >10¹⁷ potential RP combinations). Here, we engineered a yeast with monotypic 40S ribosomes, including both defined and homogenous RP paralogs, and further functional studies revealed that duplicated RP paralogs impart robustness and phenotypic plasticity (such as paromomycin tolerance) through both gene dose amplification and paralog-specific regulation, paving a way for the study of monotypic ribosomes.

Emerging evidence reveals that ribosomes are not monolithic but dynamic machines with heterogeneous protein compositions that can reshape ribosomal translational abilities and cellular adaptation to environmental changes. Duplications of ribosomal protein (RP) genes are ubiquitous among organisms and are believed to affect cell function through paralog-specific regulation (e.g., by generating heterogeneous ribosomes) and/or gene dose amplification. However, direct evaluations of their impacts on cell function remain elusive due to the highly heterogeneous cellular RP pool. Here, we engineered a yeast with homogeneous 40S RP paralog compositions, designated homo-40S, by deleting the entire set of alternative duplicated genes encoding yeast 40S RP paralogs. Homo-40S displayed mild growth defects along with high sensitivity to the translation inhibitor paromomycin and a significantly increased stop codon readthrough. Moreover, doubling of the remaining RP paralogous genes in homo-40S rescued these phenotypes markedly, although not fully, compared to the wild-type phenotype, indicating that the dose of 40S RP genes together with the heterogeneity of the contents was vital for maintaining normal translational functionalities and growth robustness. Additional experiments revealed that homo-40S improved paromomycin tolerance via acquisition of bypass mutations or evolved to be diploid to generate fast-growing derivatives, highlighting the mutational robustness of engineered yeast to accommodate environmental and genetic changes. In summary, our work demonstrated that duplicated RP paralogs impart robustness and phenotypic plasticity through both gene dose amplification and paralog-specific regulation, paving the way for the direct study of ribosome biology through monotypic ribosomes with a homogeneous composition of specific RP paralogs.

Identification of Antibacterial Peptide Candidates Encrypted in Stress-Related and Metabolic Saccharomyces cerevisiae Proteins

The protein-rich nature of Saccharomyces cerevisiae has led this yeast to the spotlight concerning the search for antimicrobial peptides. Herein, a <10 kDa peptide-rich extract displaying antibacterial activity was obtained through the autolysis of yeast biomass under mild thermal treatment with self-proteolysis by endogenous peptidases. Estimated IC₅₀ for the peptide pools obtained by FPLC gel filtration indicated improved antibacterial activities against foodborne bacteria and bacteria of clinical interest. Similarly, the estimated cytotoxicity concentrations against healthy human fibroblasts, alongside selective indices ≥10, indicates the fractions are safe, at least in a mixture format, for human tissues. Nano-LC-MS/MS analysis revealed that the peptides in FPLC fractions could be derived from both induced-proteolysis and proteasome activity in abundant proteins, up-regulated under stress conditions during S. cerevisiae biomass manufacturing, including those coded by TDH1/2/3, HSP12, SSA1/2, ADH1/2, CDC19, PGK1, PPI1, PDC1, and GMP1, as well as by other non-abundant proteins. Fifty-eight AMP candidate sequences were predicted following an in silico analysis using four independent algorithms, indicating their possible contribution to the bacterial inactivation observed in the peptides pool, which deserve special attention for further validation of individual functionality. S. cerevisiae-biomass peptides, an unconventional but abundant source of pharmaceuticals, may be promissory adjuvants to treat infectious diseases that are poorly sensitive to conventional antibiotics.

Modular, synthetic chromosomes as new tools for large scale engineering of metabolism

The construction of powerful cell factories requires intensive genetic engineering for the addition of new functionalities and the remodeling of native pathways and processes. The present study demonstrates the feasibility of extensive genome reprogramming using modular, specialized de novo-assembled neochromosomes in yeast. The in vivo assembly of linear and circular neochromosomes, carrying 20 native and 21 heterologous genes, enabled the first de novo production in a microbial cell factory of anthocyanins, plant compounds with a broad range pharmacological properties. Turned into exclusive expression platforms for heterologous and essential metabolic routes, the neochromosomes mimic native chromosomes regarding mitotic and genetic stability, copy number, harmlessness for the host and editability by CRISPR/Cas9. This study paves the way for future microbial cell factories with modular genomes in which core metabolic networks, localized on satellite, specialized neochromosomes can be swapped for alternative configurations and serve as landing pads for the addition of functionalities.

Saccharomyces cerevisiae biomass as a source of next-generation food preservatives: Evaluating potential proteins as a source of antimicrobial peptides

Saccharomyces cerevisiae is the main biotechnological tool for the production of Baker's or Brewer's biomasses, largely applied in beverage and fermented-food production. Through its gene expression reprogramming and production of compounds that inactivate the growth of other microorganisms, S. cerevisiae is able to grow in adverse environments and in complex microbial consortia, as in fruit pulps and root flour fermentations. The distinct set of up-regulated genes throughout yeast biomass propagation includes those involved in sugar fermentation, ethanol metabolization, and in protective responses against abiotic stresses. These high abundant proteins are precursors of several peptides with promising health-beneficial activities such as antihypertensive, antioxidant, antimicrobial, immunomodulatory, anti-obesity, antidiabetes, and mitogenic properties. An in silico investigation of these S. cerevisiae derived peptides produced during yeast biomass propagation or induced by physicochemical treatments were performed using four algorithms to predict antimicrobial candidates encrypted in abundantly expressed stress-related proteins encoded by different genes like AHP1, TSA1, HSP26, SOD1, HSP10, and UTR2, or metabolic enzymes involved in carbon source utilization, like ENO1/2, TDH1/2/3, ADH1/2, FBA1, and PDC1. Glyceraldehyde-3-phosphate dehydrogenase and enolase II are noteworthy precursor proteins, since they exhibited the highest scores concerning the release of antimicrobial peptide candidates. Considering the set of genes upregulated during biomass propagation, we conclude that S. cerevisiae biomass, a food-grade product consumed and marketed worldwide, should be considered a safe and nonseasonal source for designing next-generation bioactive agents, especially protein encrypting antimicrobial peptides that display broad spectra activity and could reduce the emergence of microbial resistance while also avoiding cytotoxicity.

Duplication and Functional Divergence of Branched-Chain Amino Acid Biosynthesis Genes in Aspergillus nidulans

Fungi, bacteria, and plants, but not animals, synthesize the branched-chain amino acids: leucine, isoleucine, and valine. While branched-chain amino acid (BCAA) biosynthesis has been well characterized in the yeast Saccharomyces cerevisiae, it is incompletely understood in filamentous fungi. The three BCAAs share several early biosynthesis steps before divergence into specific pathways. In Aspergillus nidulans, the genes for the first two dedicated steps in leucine biosynthesis have been characterized, but the final two have not. We used sequence searches of the A. nidulans genome to identify two genes encoding β-isopropylmalate dehydrogenase, which catalyzes the penultimate step of leucine biosynthesis, and six genes encoding BCAA aminotransferase, which catalyzes the final step in biosynthesis of all three BCAA. We have used combinations of gene knockouts to determine the relative contribution of each of these genes to BCAA biosynthesis. While both β-isopropylmalate dehydrogenase genes act in leucine biosynthesis, the two most highly expressed BCAA aminotransferases are responsible for BCAA biosynthesis. We have also characterized the expression of leucine biosynthesis genes using reverse transcriptase-quantitative PCR and found regulation in response to leucine availability is mediated through the Zn(II)₂Cys₆ transcription factor LeuB. IMPORTANCE Branched-chain amino acid (BCAA) biosynthesis is important for pathogenic fungi to successfully cause disease in human and plant hosts. The enzymes for their production are absent from humans and, therefore, provide potential antifungal targets. While BCAA biosynthesis is well characterized in yeasts, it is poorly understood in filamentous fungal pathogens. Developing a thorough understanding of both the genes encoding the metabolic enzymes for BCAA biosynthesis and how their expression is regulated will inform target selection for antifungal drug development.

Eukaryotic Elongation Factor 3 Protects Saccharomyces cerevisiae Yeast from Oxidative Stress

Translation is a core process of cellular protein homeostasis and, thus, needs to be tightly regulated. The production of newly synthesized proteins adapts to the current needs of the cell, including the response to conditions of oxidative stress. Overall protein synthesis decreases upon oxidative stress. However, the selective production of proteins is initiated to help neutralize stress conditions. In contrast to higher eukaryotes, fungi require three translation elongation factors, eEF1, eEF2, and eEF3, for protein synthesis. eEF1 and eEF2 are evolutionarily conserved, but they alone are insufficient for the translation elongation process. eEF3 is encoded by two paralogous genes, YEF3 and HEF3. However, only YEF3 is essential in yeast, whereas the function of HEF3 remains unknown. To elucidate the cellular function of Hef3p, we used cells that were depleted of HEF3 and treated with H₂O₂ and analyzed the growth of yeast, global protein production, and protein levels. We found that HEF3 is necessary to withstand oxidative stress conditions, suggesting that Hef3p is involved in the selective production of proteins that are necessary for defense against reactive oxygen species.

Genetic Engineering and Synthetic Genomics in Yeast to Understand Life and Boost Biotechnology

The field of genetic engineering was born in 1973 with the "construction of biologically functional bacterial plasmids in vitro". Since then, a vast number of technologies have been developed allowing large-scale reading and writing of DNA, as well as tools for complex modifications and alterations of the genetic code. Natural genomes can be seen as software version 1.0; synthetic genomics aims to rewrite this software with "build to understand" and "build to apply" philosophies. One of the predominant model organisms is the baker's yeast Saccharomyces cerevisiae. Its importance ranges from ancient biotechnologies such as baking and brewing, to high-end valuable compound synthesis on industrial scales. This tiny sugar fungus contributed greatly to enabling humankind to reach its current development status. This review discusses recent developments in the field of genetic engineering for budding yeast S. cerevisiae, and its application in biotechnology. The article highlights advances from Sc1.0 to the developments in synthetic genomics paving the way towards Sc2.0. With the synthetic genome of Sc2.0 nearing completion, the article also aims to propose perspectives for potential Sc3.0 and subsequent versions as well as its implications for basic and applied research.

Evaluating the probability of CRISPR-based gene drive contaminating another species

The probability D that a given clustered regularly interspaced short palindromic repeats (CRISPR)-based gene drive element contaminates another, nontarget species can be estimated by the following Drive Risk Assessment Quantitative Estimate (DRAQUE) Equation: D = h y b + t r a n s f × e x p r e s s × c u t × f l a n k × i m m u n e × n o n e x t i n c t with hyb = probability of hybridization between the target species and a nontarget species; transf = probability of horizontal transfer of a piece of DNA containing the gene drive cassette from the target species to a nontarget species (with no hybridization); express = probability that the Cas9 and guide RNA genes are expressed; cut = probability that the CRISPR-guide RNA recognizes and cuts at a DNA site in the new host; flank = probability that the gene drive cassette inserts at the cut site; immune = probability that the immune system does not reject Cas9-expressing cells; nonextinct = probability of invasion of the drive within the population. We discuss and estimate each of the seven parameters of the equation, with particular emphasis on possible transfers within insects, and between rodents and humans. We conclude from current data that the probability of a gene drive cassette to contaminate another species is not insignificant. We propose strategies to reduce this risk and call for more work on estimating all the parameters of the formula.

Impact of genome duplication on secondary metabolite composition in non-cultivated species: a systematic meta-analysis

BACKGROUND AND AIMS

Whole-genome duplication is known to influence ecological interactions and plant physiology; however, despite abundant case studies, much is still unknown about the typical impact of genome duplication on plant secondary metabolites (PSMs). In this study, we assessed the impact of polyploidy events on PSM characteristics in non-cultivated plants.

METHODS

We conducted a systematic review and meta-analysis to compare composition and concentration of PSMs among closely related plant species or species complexes differing in ploidy level.

KEY RESULTS

We assessed 53 studies that focus on PSMs among multiple cytotypes, of which only 14 studies compared concentration quantitatively among cytotypes. We found that whole-genome duplication can have a significant effect on PSM concentration; however, these effects are highly inconsistent.

CONCLUSION

Overall, there was no consistent effect of whole-genome duplication on PSM concentrations or profiles.

Transcriptional Rewiring, Adaptation, and the Role of Gene Duplication in the Metabolism of Ethanol of Saccharomyces cerevisiae

Ethanol is the main by-product of yeast sugar fermentation that affects microbial growth parameters, being considered a dual molecule, a nutrient and a stressor. Previous works demonstrated that the budding yeast arose after an ancient hybridization process resulted in a tier of duplicated genes within its genome, many of them with implications in this ethanol "produce-accumulate-consume" strategy. The evolutionary link between ethanol production, consumption, and tolerance versus ploidy and stability of the hybrids is an ongoing debatable issue. The implication of ancestral duplicates in this metabolic rewiring, and how these duplicates differ transcriptionally, remains unsolved. Here, we study the transcriptomic adaptive signatures to ethanol as a nonfermentative carbon source to sustain clonal yeast growth by experimental evolution, emphasizing the role of duplicated genes in the adaptive process. As expected, ethanol was able to sustain growth but at a lower rate than glucose. Our results demonstrate that in asexual populations a complete transcriptomic rewiring was produced, strikingly by downregulation of duplicated genes, mainly whole-genome duplicates, whereas small-scale duplicates exhibited significant transcriptional divergence between copies. Overall, this study contributes to the understanding of evolution after gene duplication, linking transcriptional divergence with duplicates' fate in a multigene trait as ethanol tolerance.IMPORTANCE Gene duplication events have been related with increasing biological complexity through the tree of life, but also with illnesses, including cancer. Early evolutionary theories indicated that duplicated genes could explore alternative functions due to relaxation of selective constraints in one of the copies, as the other remains as ancestral-function backup. In unicellular eukaryotes like yeasts, it has been demonstrated that the fate and persistence of duplicates depend on duplication mechanism (whole-genome or small-scale events), shaping their actual genomes. Although it has been shown that small-scale duplicates tend to innovate and whole-genome duplicates specialize in ancestral functions, the implication of duplicates' transcriptional plasticity and transcriptional divergence on environmental and metabolic responses remains largely obscure. Here, by experimental adaptive evolution, we show that Saccharomyces cerevisiae is able to respond to metabolic stress (ethanol as nonfermentative carbon source) due to the persistence of duplicated genes. These duplicates respond by transcriptional rewiring, depending on their transcriptional background. Our results shed light on the mechanisms that determine the role of duplicates, and on their evolvability.

Expansion of fish CCL20_like chemokines by genome and local gene duplication: Characterisation and expression analysis of 10 CCL20_like chemokines in rainbow trout (Oncorhynchus mykiss)

Mammalian CCL20, or macrophage inflammatory protein-3α, can function as a homeostatic and inflammatory chemokine. In relation to the latter, it is responsible for the chemoattraction of lymphocytes and dendritic cells to mucosal immune sites under inflammatory and pathological conditions. CK1, CK8A and CK8B are rainbow trout (Oncorhynchus mykiss) CC chemokines that were reported previously to be phylogenetically related to mammalian CCL20. In the current study, an additional seven CCL20_L paralogues in rainbow trout are reported, that are divided into three subgroups and have been designated here as: CCL20_L1a (also referred to as CK1), CCL20_L1b1-2, CCL20_L2a (CK8A), CCL20_L2b (CK8B), CCL20_L3a, and CCL20_L3b1-4. Multiple CCL20_L genes were also identified in other salmonids that arose from both whole genome duplication and local gene duplication. Phylogenetic tree, homology and synteny analysis support that CCL20_L1-3 found in salmonids are also present in most teleosts arose from the 3 R whole genome duplication and in some species, local gene duplication. Like mammalian CCL20, rainbow trout CCL20_L molecules possess a high positive net charge with a pI of 9.34-10.16, that is reported to be important for antimicrobial activity. Rainbow trout CCL20_L paralogues are differentially expressed and in general highly expressed in mucosal tissues, such as gills, thymus and intestine. The expression levels of rainbow trout CCL20_L paralogues are increased during development and following PAMP/cytokine stimulation. For example, in RTS-11 cells CCL20_L3b1 and CCL20_L3b2 are highly up-regulated by LPS, Poly I:C, recombinant(r) IFNa and rIL-1β. Trout CCL20_L paralogues are also increased after Yersinia ruckeri infection or Poly I:C stimulation in vivo, with CCL20_L3b1 and CCL20_L3b2 again highly up-regulated. Overall, this is the first report of the complete CCL20 chemokine subfamily in rainbow trout, and the analysis of their expression and modulation in vitro and in vivo. These results suggest that teleosts possess divergent CCL20_L molecules that may have important roles in anti-viral/anti-bacterial defence and in mucosal immunity.