Genetic structure and relatedness of brown trout (Salmo trutta) populations in the drainage basin of the Ölfusá river, South-Western Iceland

Background

Lake Þingvallavatn in Iceland, a part of the river Ölfusá drainage basin, was presumably populated by brown trout soon after it formed at the end of the last Ice Age. The genetic relatedness of the brown trout in Þingvallavatn to other populations in the Ölfusá drainage basin is unknown. After the building of a dam at the outlet of the lake in 1959 brown trout catches declined, though numbers have now increased. The aim of this study was to assess effects of geographic isolation and potential downstream gene flow on the genetic structure and diversity in brown trout sampled in several locations in the western side of the watershed of River Ölfusá. We hypothesized that brown trout in Lake Þingvallavatn constituted several local spawning populations connected by occasional gene flow before the damming of the lake. We also estimated the effective population size (NE) of some of these populations and tested for signs of a recent population bottleneck in Lake Þingvallavatn.

Methods

We sampled brown trout inhabiting four lakes and 12 rivers within and near the watershed of River Ölfusá by means of electro- and net- fishing. After stringent data filtering, 2,597 polymorphic loci obtained from ddRADseq data from 317 individuals were ascertained as putative neutral markers.

Results

Overall, the genetic relatedness of brown trout in the Ölfusá watershed reflected the connectivity and topography of the waterways. Ancestry proportion analyses and a phylogenetic tree revealed seven distinct clusters, some of which corresponded to small populations with reduced genetic diversity. There was no evidence of downstream gene flow from Lake Þingvallavatn, although gene flow was observed from much smaller mountain populations. Most locations showed low NE values (i.e., ~14.6 on average) while the putative anadromous trout from River Sog and the spawning population from River Öxará, that flows into Lake Þingvallavatn, showed notably higher NE values (i.e., 71.2 and 56.5, respectively). No signals of recent population bottlenecks were detected in the brown trout of Lake Þingvallavatn.

Discussion

This is the first time that the genetic structure and diversity of brown trout in the watershed of River Ölfusá have been assessed. Our results point towards the presence of a metapopulation in the watershed of Lake Þingvallavatn, which has been influenced by restoration efforts and is now dominated by a genetic component originated in River Öxará. Many of the locations studied represent different populations. Those that are isolated in headwater streams and lakes are genetically distinct presenting low genetic diversity, yet they can be important in increasing the genetic variation in downstream populations. These populations should be considered for conservation and direct management.

Variation in egg size and offspring phenotype among and within seven Arctic charr morphs

Maternal effects have the potential to alter early developmental processes of offspring and contribute to adaptive diversification. Egg size is a major contributor to offspring phenotype, which can influence developmental trajectories and potential resource use. However, to what extent intraspecific variation in egg size facilitates evolution of resource polymorphism is poorly understood. We studied multiple resource morphs of Icelandic Arctic charr, ranging from an anadromous morph-with a phenotype similar to the proposed ancestral phenotype-to sympatric morphs that vary in their degree of phenotypic divergence from the ancestral anadromous morph. We characterized variation in egg size and tested whether egg size influenced offspring phenotype at early life stages (i.e., timing of- and size at- hatching and first feeding [FF]). We predicted that egg size would differ among morphs and be less variable as morphs diverge away from the ancestral anadromous phenotype. We also predicted that egg size would correlate with offspring size and developmental timing. We found morphs had different egg size, developmental timing, and size at hatching and FF. Egg size increased as phenotypic proximity to the ancestral anadromous phenotype decreased, with larger eggs generally giving rise to larger offspring, especially at FF, but egg size had no effect on developmental rate. The interaction between egg size and the environment may have a profound impact on offspring fitness, where the resulting differences in early life-history traits may act to initiate and/or maintain resource morphs diversification.

DNA methylation differences during development distinguish sympatric morphs of Arctic charr (Salvelinus alpinus)

Changes in DNA methylation in specific coding or non-coding regions can influence development and potentially divergence in traits within species and groups. While the impact of epigenetic variation on developmental pathways associated with evolutionary divergence is the focus of intense investigation, few studies have looked at recently diverged systems. Phenotypic diversity between closely related populations of Arctic charr (Salvelinus alpinus), which diverged within the last 10 000 years, offers an interesting ecological model to address such effects. Using bisulfite sequencing, we studied general DNA methylation patterns during development in the four sympatric morphs of Arctic charr from Lake Thingvallavatn. The data revealed strong differences between developmental timepoints and between morphs (mainly along the benthic - limnetic axis), both at single CpG sites and in 1,000bp-regions. Genes located close to differentially methylated CpG sites were involved in nucleosome assembly, regulation of osteoclast differentiation, and cell-matrix adhesion. Differentially methylated regions were enriched in tRNA and rRNA sequences, and half of them were located close to transcription start sites. The expression of 14 genes showing methylation differences over time or between morphs was further investigated by qPCR and nine of these were found to be differentially expressed between morphs. Four genes (ARHGEF37-like, H3-like, MPP3 and MEGF9) showed a correlation between methylation and expression. Lastly, histone gene clusters displayed interesting methylation differences between timepoints and morphs, as well as intragenic methylation variation. The results presented here provide a motivation for further studies on the contribution of epigenetic traits, such as DNA methylation, to phenotypic diversity and developmental mechanisms.

RAD sequencing of common whelk, Buccinum undatum, reveals fine-scale population structuring in Europe and cryptic speciation within the North Atlantic

Buccinum undatum is a subtidal gastropod that exhibits clear spatial variation in several phenotypic shell traits (color, shape, and thickness) across its North Atlantic distribution. Studies of spatial phenotypic variation exist for the species; however, population genetic studies have thus far relied on a limited set of mitochondrial and microsatellite markers. Here, we greatly expand on previous work by characterizing population genetic structure in B. undatum across the North Atlantic from SNP variation obtained by RAD sequencing. There was a high degree of genetic differentiation between Canadian and European populations (Iceland, Faroe Islands, and England) consistent with the divergence of populations in allopatry (F ST > 0.57 for all pairwise comparisons). In addition, B. undatum populations within Iceland, the Faroe Islands, and England are typified by weak but significant genetic structuring following an isolation-by-distance model. Finally, we established a significant correlation between genetic structuring in Iceland and two phenotypic traits: shell shape and color frequency. The works detailed here enhance our understanding of genetic structuring in B. undatum and establish the species as an intriguing model for future genome-wide association studies.

Differences among families in craniofacial shape at early life-stages of Arctic charr (Salvelinus alpinus)

BACKGROUND

Organismal fitness can be determined at early life-stages, but phenotypic variation at early life-stages is rarely considered in studies on evolutionary diversification. The trophic apparatus has been shown to contribute to sympatric resource-mediated divergence in several taxa. However, processes underlying diversification in trophic traits are poorly understood. Using phenotypically variable Icelandic Arctic charr (Salvelinus alpinus), we reared offspring from multiple families under standardized laboratory conditions and tested to what extent family (i.e. direct genetic and maternal effects) contributes to offspring morphology at hatching (H) and first feeding (FF). To understand the underlying mechanisms behind early life-stage variation in morphology, we examined how craniofacial shape varied according to family, offspring size, egg size and candidate gene expression.

RESULTS

Craniofacial shape (i.e. the Meckel's cartilage and hyoid arch) was more variable between families than within families both across and within developmental stages. Differences in craniofacial morphology between developmental stages correlated with offspring size, whilst within developmental stages only shape at FF correlated with offspring size, as well as female mean egg size. Larger offspring and offspring from females with larger eggs consistently had a wider hyoid arch and contracted Meckel's cartilage in comparison to smaller offspring.

CONCLUSIONS

This study provides evidence for family-level variation in early life-stage trophic morphology, indicating the potential for parental effects to facilitate resource polymorphism.

Extensive genetic differentiation between recently evolved sympatric Arctic charr morphs

The availability of diverse ecological niches can promote adaptation of trophic specializations and related traits, as has been repeatedly observed in evolutionary radiations of freshwater fish. The role of genetics, environment, and history in ecologically driven divergence and adaptation, can be studied on adaptive radiations or populations showing ecological polymorphism. Salmonids, especially the Salvelinus genus, are renowned for both phenotypic diversity and polymorphism. Arctic charr (Salvelinus alpinus) invaded Icelandic streams during the glacial retreat (about 10,000 years ago) and exhibits many instances of sympatric polymorphism. Particularly, well studied are the four morphs in Lake Þingvallavatn in Iceland. The small benthic (SB), large benthic (LB), planktivorous (PL), and piscivorous (PI) charr differ in many regards, including size, form, and life history traits. To investigate relatedness and genomic differentiation between morphs, we identified variable sites from RNA-sequencing data from three of those morphs and verified 22 variants in population samples. The data reveal genetic differences between the morphs, with the two benthic morphs being more similar and the PL-charr more genetically different. The markers with high differentiation map to all linkage groups, suggesting ancient and pervasive genetic separation of these three morphs. Furthermore, GO analyses suggest differences in collagen metabolism, odontogenesis, and sensory systems between PL-charr and the benthic morphs. Genotyping in population samples from all four morphs confirms the genetic separation and indicates that the PI-charr are less genetically distinct than the other three morphs. The genetic separation of the other three morphs indicates certain degree of reproductive isolation. The extent of gene flow between the morphs and the nature of reproductive barriers between them remain to be elucidated.

Revised phylogeography of the common whelk Buccinum undatum (Gastropoda: Buccinidae) across the North Atlantic

The purpose of this study was to revisit the mitochondrial genetic divergence of North Atlantic populations of the subtidal gastropod Buccinum undatum, because previous mitochondrial DNA (mtDNA) analysis by Pálsson et al. (2014) included 16S ribosomal RNA sequences that were incorrectly assigned to the species. In the present study, population mtDNA variation is now assessed using COI sequences obtained from previous research (Pálsson et al., 2014) and, to increase the geographical cover of the study, data from recently assembled transcriptomes of 96 Icelandic whelks (Jónsson et al., 2019) and GenBank are also used. To estimate the mtDNA divergence in B. undatum across the North Atlantic, two species screening indices were used: automatic barcode gap discovey (ABGD; Puillandre et al., 2012a) and species screening threshold index (SSTI; Witt et al., 2006). Phylogenetic reconstruction revealed monophyletic Eastern and Western North Atlantic whelk lineages, which diverged early in the Pleistocene glaciation (2.1 Mya), followed by a subsequent divergence event between Greenlandic and Canadian populations at 1.3 Mya. Species screening indices, ABGD and SSTI, indicated cryptic speciation or allopatric divergence. Genetic distances between populations from the two continents were similar to or greater than interspecific genetic distances across several North Pacific and North Atlantic Buccinum species.

Morphological variation in genetically divergent populations of the common whelk, Buccinum undatum (Gastropoda: Buccinidae), across the North Atlantic

The variation in shelled marine gastropod morphology across small spatial scales can reflect restricted population connectivity, resulting in evolution or plastic responses to environmental heterogeneity. The common whelk, Buccinum undatum, is a subtidal gastropod, ubiquitous in the North Atlantic, that exhibits considerable spatial variation in shell morphology and colour. Given that species delimitation in shelled marine gastropods is often based on shell characteristics, such morphological variation can lead to taxonomic confusion. Phylogeographical analysis based on mitochondrial DNA and microsatellites suggested cryptic species composed of Western and Eastern North Atlantic common whelk populations, the separation of which dates to the onset of the Pleistocene glaciation ~2.1 Mya. Divergence within the Eastern North Atlantic is more recent and characterized by isolation by distance. In the present study, phenotypic variation in shell morphology across the North Atlantic range is analysed and compared with molecular divergence. The morphological variation of B. undatum populations reflected the pattern observed for the molecular markers only for certain comparisons of populations and might, in other cases, reflect larger constraints on the morphological variation and, possibly, the impact of environmental influences.

Gene expression in the phenotypically plastic Arctic charr (Salvelinus alpinus): A focus on growth and ossification at early stages of development

Gene expression during development shapes the phenotypes of individuals. Although embryonic gene expression can have lasting effects on developmental trajectories, few studies consider the role of maternal effects, such as egg size, on gene expression. Using qPCR, we characterize relative expression of 14 growth and/or skeletal promoting genes across embryonic development in Arctic charr (Salvelinus alpinus). We test to what extent their relative expression is correlated with egg size and size at early life‐stages within the study population. We predict smaller individuals to have higher expression of growth and skeletal promoting genes, due to less maternal resources (i.e., yolk) and prioritization of energy toward ossification. We found expression levels to vary across developmental stages and only three genes (Mmp9, Star, and Sgk1) correlated with individual size at a given developmental stage. Contrary to our hypothesis, expression of Mmp9 and Star showed a non‐linear relationship with size (at post fertilization and hatching, respectively), whilst Sgk1 was higher in larger embryos at hatching. Interestingly, these genes are also associated with craniofacial divergence of Arctic charr morphs. Our results indicate that early life‐stage variation in gene expression, concomitant to maternal effects, can influence developmental plasticity and potentially the evolution of resource polymorphism in fishes.

Shell morphology and color of the subtidal whelk Buccinum undatum exhibit fine-scaled spatial patterns

Geographical patterns in morphology can be the result of divergence among populations due to neutral or selective changes and/or phenotypic plasticity in response to different environments. Marine gastropods are ideal subjects on which to explore these patterns, by virtue of the remarkable intraspecific variation in life‐history traits and morphology often observed across relatively small spatial scales. The ubiquitous N‐Atlantic common whelk (Buccinum undatum) is well known for spatial variation in life‐history traits and morphology. Previous studies on genetic population structure have revealed that it exhibits significant differentiation across geographic distances. Within Breiðafjörður Bay, a large and shallow bay in W‐Iceland, genetic differentiation was demonstrated between whelks from sites separated by just 20 km. Here, we extended our previous studies on the common whelk in Breiðafjörður Bay by quantifying phenotypic variation in shell morphology and color throughout the Bay. We sought to test whether trait differentiation is dependent on geographic distance and/or environmental variability. Whelk in Breiðafjörður Bay displayed fine‐scale patterns of spatial variation in shape, thickness, and color diversity. Differentiation increased with increasing distance between populations, indicating that population connectivity is limited. Both shape and color varied along a gradient from the inner part of the bay in the east to the outer part in the west. Whelk shells in the innermost part of Breiðafjörður Bay were thick with an elongate shell, round aperture, and low color diversity, whereas in the outer part of the bay the shells were thinner, rounder, with a more elongate aperture and richer color diversity. Significant site‐specific difference in shell traits of the common whelk in correlation with environmental variables indicates the presence of local ecotypes and limited demographic connectivity.

Inducers of salmon innate immunity: An in vitro and in vivo approach

Maintaining fish health is one of the most important aims in aquaculture. Prevention of fish diseases therefore is crucial and can be achieved by various different strategies, including most often a combination of different methods such as optimal feed and fish density, as well as strengthening the immune system. Understanding the fish innate immune system and developing methods to activate it, in an effort to prevent infections in the first place, has been a goal in recent years. In this study we choose different inducers of the innate immune system and examined their effects in vitro on the salmon cell line CHSE-214. We found that the butyrate derivatives 4-phenyl butyrate (PBA) and β-hydroxy-β-methyl butyrate (HMB) induce the expression of various innate immune genes differentially over 24-72 h. Similarly, lipids generated from fish oils were found to have an effect on the expression of the antimicrobial peptides cathelicidin and hepcidin, as well as iNOS and the viral receptor RIG-1. Interestingly we found that vitamin D3, similar as in mammals, was able to increase cathelicidin expression in fish cells. The observed induction of these different innate immune factors correlated with antibacterial activity against Aeromonas salmonicida and antiviral activity against IPNV and ISAV in vitro. To relate this data to the in vivo situation we examined cathelicidin expression in juvenile salmon and found that salmon families vary greatly in their basal cathelicidin levels. Examining cathelicidin levels in families known to be resistant to IPNV showed that these QTL-families had lower basal levels of cathelicidin in gills, than non QTL-families. Feeding fish with HMB caused a robust increase in cathelicidin expression in gills, but not skin and this was independent of the fish being resistant to IPNV. These findings support the use of fish cell lines as a tool to develop new inducers of the fish innate immune system, but also highlight the importance of the tissue studied in vivo. Understanding the response of the innate immune system in different tissues and what effect this might have on infections and downstream cellular pathways is an interesting research topic for the future.

Investigation of the effects of estrogen on skeletal gene expression during zebrafish larval head development

The development of craniofacial skeletal structures requires well-orchestrated tissue interactions controlled by distinct molecular signals. Disruptions in normal function of these molecular signals have been associated with a wide range of craniofacial malformations. A pathway mediated by estrogens is one of those molecular signals that plays role in formation of bone and cartilage including craniofacial skeletogenesis. Studies in zebrafish have shown that while higher concentrations of 17-β estradiol (E 2) cause severe craniofacial defects, treatment with lower concentrations result in subtle changes in head morphology characterized with shorter snouts and flatter faces. The molecular basis for these morphological changes, particularly the subtle skeletal effects mediated by lower E 2 concentrations, remains unexplored. In the present study we address these effects at a molecular level by quantitative expression analysis of sets of candidate genes in developing heads of zebrafish larvae treated with two different E 2 concentrations. To this end, we first validated three suitable reference genes, ppia2, rpl8 and tbp, to permit sensitive quantitative real-time PCR analysis. Next, we profiled the expression of 28 skeletogenesis-associated genes that potentially respond to estrogen signals and play role in craniofacial development. We found E 2 mediated differential expression of genes involved in extracellular matrix (ECM) remodelling, mmp2/9/13, sparc and timp2a, as well as components of skeletogenic pathways, bmp2a, erf, ptch1/2, rankl, rarab and sfrp1a. Furthermore, we identified a co-expressed network of genes, including cpn1, dnajc3, esr1, lman1, rrbp1a, ssr1 and tram1 with a stronger inductive response to a lower dose of E 2 during larval head development.

Differential expression of the aryl hydrocarbon receptor pathway associates with craniofacial polymorphism in sympatric Arctic charr

Background

The developmental basis of craniofacial morphology hinges on interactions of numerous signalling systems. Extensive craniofacial variation in the polymorphic Arctic charr, a member of the salmonid family, from Lake Thingvallavatn (Iceland), offers opportunities to find and study such signalling pathways and their key regulators, thereby shedding light on the developmental pathways, and the genetics of trophic divergence.

Results

To identify genes involved in the craniofacial differences between benthic and limnetic Arctic charr, we used transcriptome data from different morphs, spanning early development, together with data on craniofacial expression patterns and skeletogenesis in model vertebrate species. Out of 20 genes identified, 7 showed lower gene expression in benthic than in limnetic charr morphs. We had previously identified a conserved gene network involved in extracellular matrix (ECM) organization and skeletogenesis, showing higher expression in developing craniofacial elements of benthic than in limnetic Arctic charr morphs. The present study adds a second set of genes constituting an expanded gene network with strong, benthic–limnetic differential expression. To identify putative upstream regulators, we performed knowledge-based motif enrichment analyses on the regulatory sequences of the identified genes which yielded potential binding sites for a set of known transcription factors (TFs). Of the 8 TFs that we examined using qPCR, two (Ahr2b and Ap2) were found to be differentially expressed between benthic and limnetic charr. Expression analysis of several known AhR targets indicated higher activity of the AhR pathway during craniofacial development in benthic charr morphotypes.

Conclusion

These results suggest a key role of the aryl hydrocarbon receptor (AhR) pathway in the observed craniofacial differences between distinct charr morphotypes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0022-6) contains supplementary material, which is available to authorized users.

Transcriptional dynamics of a conserved gene expression network associated with craniofacial divergence in Arctic charr

Background

Understanding the molecular basis of craniofacial variation can provide insights into key developmental mechanisms of adaptive changes and their role in trophic divergence and speciation. Arctic charr (Salvelinus alpinus) is a polymorphic fish species, and, in Lake Thingvallavatn in Iceland, four sympatric morphs have evolved distinct craniofacial structures. We conducted a gene expression study on candidates from a conserved gene coexpression network, focusing on the development of craniofacial elements in embryos of two contrasting Arctic charr morphotypes (benthic and limnetic).

Results

Four Arctic charr morphs were studied: one limnetic and two benthic morphs from Lake Thingvallavatn and a limnetic reference aquaculture morph. The presence of morphological differences at developmental stages before the onset of feeding was verified by morphometric analysis. Following up on our previous findings that Mmp2 and Sparc were differentially expressed between morphotypes, we identified a network of genes with conserved coexpression across diverse vertebrate species. A comparative expression study of candidates from this network in developing heads of the four Arctic charr morphs verified the coexpression relationship of these genes and revealed distinct transcriptional dynamics strongly correlated with contrasting craniofacial morphologies (benthic versus limnetic). A literature review and Gene Ontology analysis indicated that a significant proportion of the network genes play a role in extracellular matrix organization and skeletogenesis, and motif enrichment analysis of conserved noncoding regions of network candidates predicted a handful of transcription factors, including Ap1 and Ets2, as potential regulators of the gene network. The expression of Ets2 itself was also found to associate with network gene expression. Genes linked to glucocorticoid signalling were also studied, as both Mmp2 and Sparc are responsive to this pathway. Among those, several transcriptional targets and upstream regulators showed differential expression between the contrasting morphotypes. Interestingly, although selected network genes showed overlapping expression patterns in situ and no morph differences, Timp2 expression patterns differed between morphs.

Conclusion

Our comparative study of transcriptional dynamics in divergent craniofacial morphologies of Arctic charr revealed a conserved network of coexpressed genes sharing functional roles in structural morphogenesis. We also implicate transcriptional regulators of the network as targets for future functional studies.

Electronic supplementary material

The online version of this article (doi:10.1186/2041-9139-5-40) contains supplementary material, which is available to authorized users.

Patterns of miRNA expression in Arctic charr development

Micro-RNAs (miRNAs) are now recognized as a major class of developmental regulators. Sequences of many miRNAs are highly conserved, yet they often exhibit temporal and spatial heterogeneity in expression among species and have been proposed as an important reservoir for adaptive evolution and divergence. With this in mind we studied miRNA expression during embryonic development of offspring from two contrasting morphs of the highly polymorphic salmonid Arctic charr (Salvelinus alpinus), a small benthic morph from Lake Thingvallavatn (SB) and an aquaculture stock (AC). These morphs differ extensively in morphology and adult body size. We established offspring groups of the two morphs and sampled at several time points during development. Four time points (3 embryonic and one just before first feeding) were selected for high-throughput small-RNA sequencing. We identified a total of 326 conserved and 427 novel miRNA candidates in Arctic charr, of which 51 conserved and 6 novel miRNA candidates were differentially expressed among developmental stages. Furthermore, 53 known and 19 novel miRNAs showed significantly different levels of expression in the two contrasting morphs. Hierarchical clustering of the 53 conserved miRNAs revealed that the expression differences are confined to the embryonic stages, where miRNAs such as sal-miR-130, 30, 451, 133, 26 and 199a were highly expressed in AC, whereas sal-miR-146, 183, 206 and 196a were highly expressed in SB embryos. The majority of these miRNAs have previously been found to be involved in key developmental processes in other species such as development of brain and sensory epithelia, skeletogenesis and myogenesis. Four of the novel miRNA candidates were only detected in either AC or SB. miRNA candidates identified in this study will be combined with available mRNA expression data to identify potential targets and involvement in developmental regulation.

Validation of reference genes for expression studies during craniofacial development in arctic charr

Arctic charr (Salvelinus alpinus) is a highly polymorphic species and in Lake Thingvallavatn, Iceland, four phenotypic morphs have evolved. These differences in morphology, especially in craniofacial structures are already apparent during embryonic development, indicating that genes important in the formation of the craniofacial features are expressed differentially between the morphs. In order to generate tools to examine these expression differences in Arctic charr, the aim of the present study was to identify reference genes for quantitative real-time PCR (qPCR). The specific aim was to select reference genes which are able to detect very small expression differences among different morphs. We selected twelve candidate reference genes from the literature, identified corresponding charr sequences using data derived from transcriptome sequencing (RNA-seq) and examined their expression using qPCR. Many of the candidate reference genes were found to be stably expressed, yet their quality-rank as reference genes varied considerably depending on the type of analysis used. In addition to commonly used software for reference gene validation, we used classical statistics to evaluate expression profiles avoiding a bias for reference genes with similar expression patterns (co-regulation). Based on these analyses we chose three reference genes, ACTB, UB2L3 and IF5A1 for further evaluation. Their consistency was assessed in an expression study of three known craniofacially expressed genes, sparc (or osteonectin), matrix metalloprotease 2 (mmp2) and sox9 (sex-determining region Y box 9 protein) using qPCR in embryo heads derived from four charr groups at three developmental time points. The three reference genes were found to be very suitable for studying expression differences between the morphotypes, enabling robust detection of small relative expression changes during charr development. Further, the results showed that sparc and mmp2 are differentially expressed in embryos of different Arctic charr morphotypes.

Mitochondrial genomes from the lichenized fungi Peltigera membranacea and Peltigera malacea: features and phylogeny

Mitochondrial genomes from the fungal partners of two terricolous foliose lichen symbioses, Peltigera membranacea and Peltigera malacea, have been determined using metagenomic approaches, including RNA-seq. The roughly 63 kb genomes show all the major features found in other Pezizomycotina, such as unidirectional transcription, 14 conserved protein genes, genes for the two subunit rRNAs and for a set of 26 tRNAs used in translating the 62 amino acid codons. In one of the tRNAs a CAU anticodon is proposed to be modified, via the action of the nuclear-encoded enzyme, tRNA Ile lysidine synthase, so that it recognizes the codon AUA (Ile) instead of AUG (Met). The overall arrangements and sequences of the two circular genomes are similar, the major difference being the inversion and deterioration of a gene encoding a type B DNA polymerase. Both genomes encode the RNA component of RNAse P, a feature seldom found in ascomycetes. The difference in genome size from the minimal ascomycete mitochondrial genomes is largely due to 17 and 20 group I introns, respectively, most associated with homing endonucleases and all found within protein-coding genes and the gene encoding the large subunit rRNA. One new intron insertion point was found, and an unusually small exon of seven nucleotides (nt) was identified and verified by RNA sequencing. Comparative analysis of mitochondrion-encoded proteins places the Peltigera spp., representatives of the class Lecanoromycetes, close to Leotiomycetes, Dothidiomycetes, and Sordariomycetes, in contrast to phylogenies found using nuclear genes.

Cloning and heterologous transcription of a polyketide synthase gene from the lichen Solorina crocea

Lichens and most ascomycete fungi produce polyketide secondary metabolites often with valuable biological activities. Their biosynthesis is primarily governed by large iterative multifunctional type I polyketide synthases. Although there has been good progress studying filamentous non-lichenized fungi, there is limited information on polyketide biosynthesis in lichens and their mycobionts, due to their slow growth, difficulties in establishing pure cultures, and the absence of methods for direct genetic manipulation. However, heterologous expression in a surrogate host offers an alternative approach for exploring lichen polyketide biosynthesis. Here, we report cloning of a type I polyketide synthase gene from the foliose lichen Solorina crocea and its heterologous transcription in the filamentous fungus Aspergillus oryzae, including processing of the transcript. No new polyketide product was detected. The lichen polyketide synthase showed greatest homology with uncharacterized genes from filamentous fungi and lower homology with proteins catalysing biosynthesis of the decaketide alternapyrone and the tetraketide side-chain of squalestatin. The technology platform utilized here presents a useful tool for functional characterization of fungal biosynthetic genes and provides a means for novel production of valuable compounds.

Rvb1p/Rvb2p recruit Arp5p and assemble a functional Ino80 chromatin remodeling complex

The Rvb1p and Rvb2p (or TIP48 and TIP49) nuclear ATP binding proteins are universally conserved in eukaryotes and essential for viability of yeasts. Rvbp associate with each other as a double hexamer, with YHR034c and with two complexes involved in chromatin remodeling, Ino80.com and Swr1.com. Loss of Rvb1p or Ino80p affects many yeast promoters similarly. Rvbp are not essential for the recruitment of Ino80p to promoters but are essential for the catalytic activity of Ino80.com. Loss of Rvbp leads to loss of the functionally critical Arp5p in Ino80.com. Rvb2p associates with Arp5p in vitro in a reaction dependent on the presence of ATP and Ino80p. Therefore, Rvbp are required for the structural and functional integrity of the Ino80 chromatin remodeling complex.

Small RNAs with imperfect match to endogenous mRNA repress translation. Implications for off-target activity of small inhibitory RNA in mammalian cells

A 21-base pair RNA duplex that perfectly matches an endogenous target mRNA selectively degrades the mRNA and suppresses gene expression in mammalian tissue culture cells. A single base mismatch with the target is believed to protect the mRNA from degradation, making this type of interference highly specific to the targeted gene. A short RNA with mismatches to a target sequence present in multiple copies in the 3'-untranslated region of an exogenously expressed gene can, however, silence it by translational repression. Here we report that a mismatched RNA, targeted to a single site in the coding sequence of an endogenous gene, can efficiently silence gene expression by repressing translation. The antisense strand of such a mismatched RNA requires a 5'-phosphate but not a 3'-hydroxyl group. G.U wobble base pairing is tolerated as a match for both RNA degradation and translation repression. Together, these findings suggest that a small inhibitory RNA duplex can suppress expression of off-target cellular proteins by RNA degradation or translation repression. Proper design of experimental small inhibitory RNAs or a search for targets of endogenous micro-RNAs must therefore take into account that these short RNAs can affect expression of cellular genes with as many as 3-4 base mismatches and additional G.U mismatches.

Rvb1p and Rvb2p are essential components of a chromatin remodeling complex that regulates transcription of over 5% of yeast genes

Eukaryotic Rvb1p and Rvb2p are two highly conserved proteins related to the helicase subset of the AAA+ family of ATPases. Conditional mutants in both genes show rapid changes in the transcription of over 5% of yeast genes, with a similar number of genes being repressed and activated. Both Rvb1p and Rvb2p are required for maintaining the induced state of many inducible promoters. ATP binding and hydrolysis by Rvb1p and Rvb2p is individually essential in vivo, and the two proteins are associated with each other in a high molecular weight complex that shows ATP-dependent chromatin remodeling activity in vitro. Our findings show that Rvb1p and Rvb2p are essential components of a chromatin remodeling complex and determine genes regulated by the complex.

In eukaryotic flap endonuclease 1, the C terminus is essential for substrate binding

Flap endonuclease 1 (Fen1) is a structure-specific metallonuclease with important functions in DNA replication and DNA repair. It interacts like many other proteins involved in DNA metabolic events with proliferating cell nuclear antigen (PCNA), and its enzymatic activity is stimulated by PCNA in vitro. The PCNA interaction site is located close to the C terminus of Fen1 and is flanked by a conserved basic region of 35-38 amino acids in eukaryotic species but not in archaea. We have constructed two deletion mutants of human Fen1 that lack either the PCNA interaction motif or a part of its adjacent C-terminal region and analyzed them in a variety of assays. Remarkably, deletion of the basic C-terminal region did not affect PCNA interaction but resulted in a protein with significantly reduced enzymatic activity. Electrophoretic mobility shift analysis revealed that this mutant displayed a severe defect in substrate binding. Our results suggest that the C terminus of eukaryotic Fen1 consists of two functionally distinct regions that together might form an important regulatory domain.

A coordinated interplay: proteins with multiple functions in DNA replication, DNA repair, cell cycle/checkpoint control, and transcription

In eukaryotic cells, DNA transactions such as replication, repair, and transcription require a large set of proteins. In all of these events, complexes of more than 30 polypetides appear to function in highly organized and structurally well-defined machines. We have learned in the past few years that the three essential macromolecular events, replication, repair, and transcription, have common functional entities and are coordinated by complex regulatory mechanisms. This can be documented for replication and repair, for replication and checkpoint control, and for replication and cell cycle control, as well as for replication and transcription. In this review we cover the three different protein classes: DNA polymerases, DNA polymerase accessory proteins, and selected transcription factors. The "common enzyme-different pathway strategy" is fascinating from several points of view: first, it might guarantee that these events are coordinated; second, it can be viewed from an evolutionary angle; and third, this strategy might provide cells with backup mechanisms for essential physiological tasks.

The PCNA from Thermococcus fumicolans functionally interacts with DNA polymerase delta

We have cloned the gene encoding proliferating cell nuclear antigen (PCNA) from the hyperthermophilic euryarchaeote Thermococcus fumicolans (Tfu). Tfu PCNA contains 250 amino acids with a calculated M(r) of 28,000 and is 26% identical to human PCNA. Next, Tfu PCNA was overexpressed in Escherichia coli and it showed an apparent molecular mass of 33.5 kDa. The purified Tfu PCNA was tested first with recombinant Tfu DNA polymerase I (Tfu pol) and second with calf thymus DNA polymerase delta (pol delta). When tested with the homologous Tfu pol on bacteriophage lambda DNA, large amounts of Tfu PCNA were required to obtain two- to threefold stimulation. Surprisingly, however, Tfu PCNA was much more efficient than human PCNA in stimulating calf thymus pol delta. Our data suggest that PCNA has been functionally conserved not only within eukaryotes but also from hyperthermophilic euryarchaeotes to mammals.

A direct interaction between proliferating cell nuclear antigen (PCNA) and Cdk2 targets PCNA-interacting proteins for phosphorylation

Proliferating cell nuclear antigen is best known as a DNA polymerase accessory protein but has more recently also been shown to have different functions in important cellular processes such as DNA replication, DNA repair, and cell cycle control. PCNA has been found in quaternary complexes with the cyclin kinase inhibitor p21 and several pairs of cyclin-dependent protein kinases and their regulatory partner, the cyclins. Here we show a direct interaction between PCNA and Cdk2. This interaction involves the regions of the PCNA trimer close to the C termini. We found that PCNA and Cdk2 form a complex together with cyclin A. This ternary PCNA-Cdk2-cyclin A complex was able to phosphorylate the PCNA binding region of the large subunit of replication factor C as well as DNA ligase I. Furthermore, PCNA appears to be a connector between Cdk2 and DNA ligase I and to stimulate phosphorylation of DNA ligase I. Based on our results, we propose the model that PCNA brings Cdk2 to proteins involved in DNA replication and possibly might act as an "adaptor" for Cdk2-cyclin A to PCNA-binding DNA replication proteins.

A DNA ligase from the psychrophile Pseudoalteromonas haloplanktis gives insights into the adaptation of proteins to low temperatures

The cloning, overexpression and characterization of a cold-adapted DNA ligase from the Antarctic sea water bacterium Pseudoalteromonas haloplanktis are described. Protein sequence analysis revealed that the cold-adapted Ph DNA ligase shows a significant level of sequence similarity to other NAD+-dependent DNA ligases and contains several previously described sequence motifs. Also, a decreased level of arginine and proline residues in Ph DNA ligase could be involved in the cold-adaptation strategy. Moreover, 3D modelling of the N-terminal domain of Ph DNA ligase clearly indicates that this domain is destabilized compared with its thermophilic homologue. The recombinant Ph DNA ligase was overexpressed in Escherichia coli and purified to homogeneity. Mass spectroscopy experiments indicated that the purified enzyme is mainly in an adenylated form with a molecular mass of 74 593 Da. Ph DNA ligase shows similar overall catalytic properties to other NAD+-dependent DNA ligases but is a cold-adapted enzyme as its catalytic efficiency (kcat/Km) at low and moderate temperatures is higher than that of its mesophilic counterpart E. coli DNA ligase. A kinetic comparison of three enzymes adapted to different temperatures (P. haloplanktis, E. coli and Thermus scotoductus DNA ligases) indicated that an increased kcat is the most important adaptive parameter for enzymatic activity at low temperatures, whereas a decreased Km for the nicked DNA substrate seems to allow T. scotoductus DNA ligase to work efficiently at high temperatures. Besides being useful for investigation of the adaptation of enzymes to extreme temperatures, P. haloplanktis DNA ligase, which is very efficient at low temperatures, offers a novel tool for biotechnology.

Optimised ligation of oligonucleotides by thermal ligases: comparison of Thermus scotoductus and Rhodothermus marinus DNA ligases to other thermophilic ligases

We describe the characterisation of four thermo-stable NAD(+)-dependent DNA ligases, from Thermus thermophilus (Tth), Thermus scotoductus (Ts), Rhodothermus marinus (Rm) and Thermus aquaticus (Taq), by an assay which measures ligation rate and mismatch discrimination. Complete libraries of octa-, nona- and decanucleotides were used as substrates. The assay comprised the polymerisation of oligo-nucleotides initiated from a 17 base 'primer', using M13mp18 ssDNA as template. Polymers of ligation products were analysed by polyacrylamide gel electro-phoresis. Under optimum conditions, the enzymes produced polymers ranging from 8 to 16 additions; there was variation between enzymes and the length of the oligonucleotides had a strong effect. The optimal total oligonucleotide concentration for each library was approximately 4 nmol. We compared the rates of ligation between the four ligases using an octanucleotide library as substrate. By this criterion, the Ts and Rm ligases are far more active compared to the more commonly available thermostable ligases.

A CAF-1-PCNA-mediated chromatin assembly pathway triggered by sensing DNA damage

Sensing DNA damage is crucial for the maintenance of genomic integrity and cell cycle progression. The participation of chromatin in these events is becoming of increasing interest. We show that the presence of single-strand breaks and gaps, formed either directly or during DNA damage processing, can trigger the propagation of nucleosomal arrays. This nucleosome assembly pathway involves the histone chaperone chromatin assembly factor 1 (CAF-1). The largest subunit (p150) of this factor interacts directly with proliferating cell nuclear antigen (PCNA), and critical regions for this interaction on both proteins have been mapped. To isolate proteins specifically recruited during DNA repair, damaged DNA linked to magnetic beads was used. The binding of both PCNA and CAF-1 to this damaged DNA was dependent on the number of DNA lesions and required ATP. Chromatin assembly linked to the repair of single-strand breaks was disrupted by depletion of PCNA from a cell-free system. This defect was rescued by complementation with recombinant PCNA, arguing for role of PCNA in mediating chromatin assembly linked to DNA repair. We discuss the importance of the PCNA-CAF-1 interaction in the context of DNA damage processing and checkpoint control.

Long patch base excision repair with purified human proteins. DNA ligase I as patch size mediator for DNA polymerases delta and epsilon

Among the different base excision repair pathways known, the long patch base excision repair of apurinic/apyrimidinic sites is an important mechanism that requires proliferating cell nuclear antigen. We have reconstituted this pathway using purified human proteins. Our data indicated that efficient repair is dependent on six components including AP endonuclease, replication factor C, proliferating cell nuclear antigen, DNA polymerases delta or epsilon, flap endonuclease 1, and DNA ligase I. Fine mapping of the nucleotide replacement events showed that repair patches extended up to a maximum of 10 nucleotides 3' to the lesion. However, almost 70% of the repair synthesis was confined to 2-4-nucleotide patches and DNA ligase I appeared to be responsible for limiting the repair patch length. Moreover, both proliferating cell nuclear antigen and flap endonuclease 1 are required for the production and ligation of long patch repair intermediates suggesting an important role of this complex in both excision and resynthesis steps.

Dual mode of interaction of DNA polymerase epsilon with proliferating cell nuclear antigen in primer binding and DNA synthesis

Proliferating cell nuclear antigen can interact with DNA polymerase epsilon on linear DNA templates, even in the absence of other auxiliary factors (replication factor C, replication protein A), and thereby stimulate its primer recognition and DNA synthesis. Using four characterized mutants of proliferating cell nuclear antigen containing three or four alanine residue substitutions on the C-terminal side and the back side of the trimer, we have tested the kinetics of primer binding and nucleotide incorporation by DNA polymerase epsilon in different assays. In contrast with what has been found in interaction studies between DNA polymerase delta and proliferating cell nuclear antigen, our data suggested that stimulation of DNA polymerase epsilon primer binding involves interactions with both the C-terminal side and the back side of proliferating cell nuclear antigen. However, for stimulation of DNA polymerase epsilon DNA synthesis, exclusively the C-terminal side appears to be sufficient. The significance of this dual interaction is discussed with reference to the physiological roles of DNA polymerase epsilon and its interaction with the clamp proliferating cell nuclear antigen.

Regulation of DNA replication and repair proteins through interaction with the front side of proliferating cell nuclear antigen

The DNA polymerase accessory factor proliferating cell nuclear antigen (PCNA) has been caught in interaction with an ever increasing number of proteins. To characterize the sites and functions of some of these interactions, we constructed four mutants of human PCNA and analysed them in a variety of assays. By targeting loops on the surface of the PCNA trimer and changing three or four residues at a time to alanine, we found that a region including part of the domain-connecting loop of PCNA and loops on one face of the trimer, close to the C-termini, is involved in binding to all of the following proteins: DNA polymerase delta, replication factor C, the flap endonuclease Fen1, the cyclin dependent kinase inhibitor p21 and DNA ligase I. An inhibition of DNA ligation caused by the interaction of PCNA with DNA ligase I was found, and we show that DNA ligase I and Fen1 can inhibit DNA synthesis by DNA polymerase delta/PCNA. We demonstrate that PCNA must be located below a 5' flap on a forked template to stimulate Fen1 activity, and considering the interacting region on PCNA for Fen1, this suggests an orientation for PCNA during DNA replication with the C-termini facing forwards, in the direction of DNA synthesis.

The solution structure of functionally active human proliferating cell nuclear antigen determined by small-angle neutron scattering

The function of proliferating cell nuclear antigen (PCNA) in DNA replication and repair is to form a sliding clamp with replication factor C (RF-C) tethering DNA polymerase delta or epsilon to DNA. In addition, PCNA has been found to interact directly with various proteins involved in cell cycle regulation. The crystal structure of yeast PCNA shows that the protein forms a homotrimeric ring lining a hole through which double-stranded DNA can thread, thus forming a moving platform for DNA synthesis. Human and yeast PCNA are highly conserved at a structural and functional level. We determined the solution structure of functionally active human PCNA by small-angle neutron scattering. Our measurements strongly support a trimeric ring-like structure of functionally active PCNA in solution, and the data are in good agreement with model calculations based on the crystal structure from yeast PCNA. The human PCNA used in the small-angle neutron scattering experiments was active before and after the measurements in a RF-C independent and a RF-C dependent assay suggesting that the trimeric structure is the in vivo functional form.

Proliferating cell nuclear antigen: more than a clamp for DNA polymerases

DNA metabolic events such as replication, repair and recombination require the concerted action of several enzymes and cofactors. Nature has provided a set of proteins that support DNA polymerases in performing processive, accurate and rapid DNA synthesis. Two of them, the proliferating cell nuclear antigen and its adapter protein replication factor C, cooperate to form a moving platform that was initially thought of only as an anchor point for DNA polymerases delta and epsilon. It now appears that proliferating cell nuclear antigen is also a communication point between a variety of important cellular processes including cell cycle control, DNA replication, nucleotide excision repair, post-replication mismatch repair, base excision repair and at least one apoptotic pathway. The dynamic movement of proliferating cell nuclear antigen on and off the DNA renders this protein an ideal communicator for a variety of proteins that are essential for DNA metabolic events in eukaryotic cells.

Replication factor C interacts with the C-terminal side of proliferating cell nuclear antigen

Replication factor C (RF-C) is a heteropentameric protein essential for DNA replication and repair. It is a molecular matchmaker required for loading of proliferating cell nuclear antigen (PCNA) onto double-stranded DNA and, thus, for PCNA-dependent DNA elongation by DNA polymerases delta and epsilon. To elucidate the mode of RF-C binding to the PCNA clamp, modified forms of human PCNA were used that could be 32P-labeled in vitro either at the C or the N terminus. Using a kinase protection assay, we show that the heteropentameric calf thymus RF-C was able to protect the C-terminal region but not the N-terminal region of human PCNA from phosphorylation, suggesting that RF-C interacts with the PCNA face at which the C termini are located (C-side). A similar protection profile was obtained with the recently identified PCNA binding region (residues 478-712), but not with the DNA binding region (residues 366-477), of the human RF-C large subunit (Fotedar, R., Mossi, R., Fitzgerald, P., Rousselle, T., Maga, G., Brickner, H., Messner, H., Khastilba, S., Hübscher, U., and Fotedar, A., (1996) EMBO J., 15, 4423-4433). Furthermore, we show that the RF-C 36 kDa subunit of human RF-C could interact independently with the C-side of PCNA. The RF-C large subunit from a third species, namely Drosophila melanogaster, interacted similarly with the modified human PCNA, indicating that the interaction between RF-C and PCNA is conserved through eukaryotic evolution.

Interaction studies between the p21Cip1/Waf1 cyclin-dependent kinase inhibitor and proliferating cell nuclear antigen (PCNA) by surface plasmon resonance

The cyclin-dependent kinase (CDK) inhibitor p21Cip1 consists of two domains that interact with CDKs and proliferating cell nuclear antigen (PCNA), respectively. We have investigated the interaction between p21Cip1 and PCNA using surface plasmon resonance (SPR) technology and compared the results with those obtained from other sources such as the yeast two-hybrid system. Whilst other methods are only semi-quantitative, the SPR technique allowed us to determine the kinetic parameters of the interaction. The apparent equilibrium constant KD calculated for these kinetic parameters was 3.2 x 10(-7) M. We further demonstrate the use of SPR to study the interaction between mutant proteins and to determine their actual KD. The interaction between p21Cip1/PCNA is shown to be dependent upon the trimeric conformation of PCNA since a point mutant that abolishes PCNA-PCNA interaction also abolishes PCNA's interaction with p21Cip1. Finally, we demonstrate that SPR can be used to characterise the interaction of p21Cip1 and PCNA in the presence of short competitive peptides.

Tyrosine 114 is essential for the trimeric structure and the functional activities of human proliferating cell nuclear antigen

In order to study the effect of trimerization of proliferating cell nuclear antigen (PCNA) on its interaction with DNA polymerase (pol) delta and its loading onto DNA by replication factor C (RF-C) we have mutated a single tyrosine residue located at the subunit interface (Tyr114) to alanine. This mutation (Y114A) had a profound effect on PCNA, since it completely abolished trimer formation as seen by glycerol gradient sedimentation and native gel electrophoresis. Furthermore, the mutant protein was unable to stimulate DNA synthesis by pol delta and did not compete effectively with wild-type PCNA for pol delta, although it was able to oligomerize and could to some extent interact with subunits of functionally active PCNA. We thus conclude that PCNA molecules that are not part of a circular trimeric complex cannot interact with the pol delta core. furthermore, the mutant protein could not be loaded onto DNA by RF-C and did not compete with wild-type PCNA for loading onto DNA, indicating that PCNA trimerization may also be a prerequisite for its recognition by RF-C. The adverse effects caused by this single mutation suggest that trimerization of PCNA is essential for the monomers to keep their overall structure and that the structural changes imposed by trimerization are important for interaction with other proteins.

Cloning and sequence analysis of the DNA ligase-encoding gene of Rhodothermus marinus, and overproduction, purification and characterization of two thermophilic DNA ligases

In this paper we describe the cloning and sequence analysis of a gene encoding DNA ligase (Lig; EC 6.5.1.2) from the thermophilic bacterium Rhodothermus marinus (Rm). We also describe the overexpression of the Lig-encoding genes of Rm and the thermophile, Thermus scotoductus (Ts), in Escherichia coli, and the purification and characterization of the overproduced Lig. The Rm lig gene encodes a protein of 712 amino acids (aa) with a calculated molecular mass of 79,487 Da. Comparison with published sequences of bacterial Lig revealed significant homology between the NAD(+)-utilizing Lig, and alignment of their aa sequences revealed several blocks of conserved residues. Both of the purified Lig exhibit nick-closing activity over a wide range of temperatures. Under our assay conditions the Rm Lig was active at 5-75 degrees C with apparent optimal activity above 55 degrees C. The Ts enzyme showed activity at 15-75 degrees C with optimal activity above 65 degrees C. The half-life of the Lig at 91 degrees C was estimated to be 7 min for the Rm Lig and 26 min for the Ts Lig.

Sequence of the DNA ligase-encoding gene from Thermus scotoductus and conserved motifs in DNA ligases

By dideoxynucleotide sequencing of a genomic clone, we have determined the complete nucleotide sequence of the gene encoding NAD(+)-dependent DNA ligase (EC 6.5.1.2) of the thermophilic bacterium Thermus scotoductus. The gene encodes a 674-amino-acid thermostable enzyme highly similar to other bacterial DNA ligases and to parts of the deduced gene product of Escherichia coli ORF f562, 5' to the spoR gene encoding 5' guanosyl kinase.