Effects of natural selection on patterns of DNA sequence variation at the transferrin, somatolactin, and p53 genes within and among chinook salmon (Oncorhynchus tshawytscha) populations

Evolutionary conservation of transferrin genomic organization and expression characterization in seven freshwater turtles

Serum transferrin (tf), encoding an iron-binding glycoprotein, has been revealed to play important roles in iron transportation and immune response, and it also has been demonstrated to be valuable for phylogenetic analysis in vertebrates. However, the evolutionary conservation, expression profiles and positive selection of transferrin genes among freshwater turtle species remain largely unclear. Here, the genomic DNA and coding sequences of transferrin genes were cloned and characterized in seven freshwater turtles including Mauremys mutica, Mauremys sinensis, Cyclemys dentate, Mauremyssi reevesi, Heosemys grandis, Trachemys scripta and Chrysemys picta. The isolated coding sequences of turtles' tf genes were 2118 bp or 2121 bp, encoding 706 or 707 amino acids. The predicted Tf proteins of turtles share high identities with M. mutica Tf, up to 91%-98% and the M. mutica Tf has the highest identity (91%) in amino acid with the Chelomia mydas Tf, the moderate with other reptiles' Tfs (65%-59%), chicken (58%), and Human Tf (∼55%), and the lowest with zebrafish Tf (41%). Additionally, tf genes were consistently composed of 17 exons and 16 introns with the same splicing sites in introns in all the turtles examined. Moreover, 12 positive selected sites were detected in these turtles' Tf and mainly distributed on the surface of transferrin protein. Importantly, it was found that transferrin genes in all turtles examined were predominantly expressed in adult liver via real-time quantitative PCR. The molecular characterizations and expression profiles of transferrin would shed new insights into understanding the conversations and divergences of transferrin genes in turtles, even in vertebrates.

Polymorphism, selection and tandem duplication of transferrin genes in Atlantic cod (Gadus morhua)--conserved synteny between fish monolobal and tetrapod bilobal transferrin loci

BACKGROUND

The two homologous iron-binding lobes of transferrins are thought to have evolved by gene duplication of an ancestral monolobal form, but any conserved synteny between bilobal and monolobal transferrin loci remains unexplored. The important role played by transferrin in the resistance to invading pathogens makes this polymorphic gene a highly valuable candidate for studying adaptive divergence among local populations.

RESULTS

The Atlantic cod genome was shown to harbour two tandem duplicated serum transferrin genes (Tf1, Tf2), a melanotransferrin gene (MTf), and a monolobal transferrin gene (Omp). Whereas Tf1 and Tf2 were differentially expressed in liver and brain, the Omp transcript was restricted to the otoliths. Fish, chicken and mammals showed highly conserved syntenic regions in which monolobal and bilobal transferrins reside, but contrasting with tetrapods, the fish transferrin genes are positioned on three different linkage groups. Sequence alignment of cod Tf1 cDNAs from Northeast (NE) and Northwest (NW) Atlantic populations revealed 22 single nucleotide polymorphisms (SNP) causing the replacement of 16 amino acids, including eight surface residues revealed by the modelled 3D-structures, that might influence the binding of pathogens for removal of iron. SNP analysis of a total of 375 individuals from 14 trans-Atlantic populations showed that the Tf1-NE variant was almost fixed in the Baltic cod and predominated in the other NE Atlantic populations, whereas the NW Atlantic populations were more heterozygous and showed high frequencies of the Tf-NW SNP alleles.

CONCLUSIONS

The highly conserved synteny between fish and tetrapod transferrin loci infers that the fusion of tandem duplicated Omp-like genes gave rise to the modern transferrins. The multiple nonsynonymous substitutions in cod Tf1 with putative structural effects, together with highly divergent allele frequencies among different cod populations, strongly suggest evidence for positive selection and local adaptation in trans-Atlantic cod populations.

Maximum likelihood analysis of mammalian p53 indicates the presence of positively selected sites and higher tumorigenic mutations in purifying sites

The tumor suppressor gene TP53 (p53) maintains genome stability. Mutation or loss of p53 is found in most cancers. Analysis of evolutionary constrains and p53 mutations reveal important sites for concomitant functional studies. In this study, phylogenetic analyses of the coding sequences of p53 from 26 mammals were carried out by applying a maximum likelihood method. The results display two branches under adaptive evolution in mammals. Moreover, each codon of p53 was analyzed by the PAML method for presence of positively selected sites. PAML identified several statistically significant amino acids that undergo positive selection. The data indicates that amino acids responsible for the core functions of p53 are highly conserved, while positively selected sites are predominantly located in the N- and C-terminus of p53. Further analysis of evolutionary pressure and mutations showed the occurrence of more frequent tumorigenic mutations in purifying sites of p53.

Extent and scale of local adaptation in salmonid fishes: review and meta-analysis

What is the extent and scale of local adaptation (LA)? How quickly does LA arise? And what is its underlying molecular basis? Our review and meta-analysis on salmonid fishes estimates the frequency of LA to be ∼55-70%, with local populations having a 1.2 times average fitness advantage relative to foreign populations or to their performance in new environments. Salmonid LA is evident at a variety of spatial scales (for example, few km to>1000 km) and can manifest itself quickly (6-30 generations). As the geographic scale between populations increases, LA is generally more frequent and stronger. Yet the extent of LA in salmonids does not appear to differ from that in other assessed taxa. Moreover, the frequency with which foreign salmonid populations outperform local populations (∼23-35%) suggests that drift, gene flow and plasticity often limit or mediate LA. The relatively few studies based on candidate gene and genomewide analyses have identified footprints of selection at both small and large geographical scales, likely reflecting the specific functional properties of loci and the associated selection regimes (for example, local niche partitioning, pathogens, parasites, photoperiodicity and seasonal timing). The molecular basis of LA in salmonids is still largely unknown, but differential expression at the same few genes is implicated in the convergent evolution of certain phenotypes. Collectively, future research will benefit from an integration of classical and molecular approaches to understand: (i) species differences and how they originate, (ii) variation in adaptation across scales, life stages, population sizes and environmental gradients, and (iii) evolutionary responses to human activities.

Analysis of single nucleotide polymorphisms in three chromosomes of European sea bass Dicentrarchus labrax

Single nucleotide polymorphisms (SNPs) are believed to contain relevant information and have been therefore extensively used as genetic markers in population and conservation genetics, and molecular ecology studies. This study reports on the identification of potential SNPs in a diploid European sea bass Dicentrarchus labrax genome by using reference sequences from three assembled chromosomes and mapping all WGS datasets onto them (3× Sanger, 3× 454 and 20× SOLEXA). A total of 20,779 SNPs were identified over the 1469 gene loci and intergenic space analysed. Within chromosomes the occurrence of SNPs was the lowest in exons and higher in introns and intergenic regions, which may be explained by the fact, that coding regions are under strong selective pressure to maintain their biological function. The ratio of nonsynonymous to synonymous mutations was smaller than one for all the chromosomes, suggesting that most of deleterious nonsynonymous mutations were eliminated by negative selection. SNPs were not uniformly distributed over the chromosomes. Two chromosomes exhibited large regions with extremely low SNP density, which might represent homozygous regions in the diploid genome. The results of this study show how SNP detection can take profit from sequencing a single diploid individual, but also uncover the limits of such an approach. SNPs that have been identified will support marker development for genetic linkage mapping, population genetics and aquaculture related questions in general.

Evidence for two transferrin loci in the Salmo trutta genome

To determine the organization of transferrin (TF) locus in the Salmo trutta genome, partial DNA and cDNA sequencing, fluorescent in situ hybridization (FISH) and Salmo salar BAC analysis were performed. TF expression levels and copy number prediction were assessed using real-time PCR. In addition to two previously reported DNA TF variant sequences of S. trutta and Salmo marmoratus (TF1), two novel variant sequences (TF2) were revealed in both species. Variant-specific sequence tags, characterizing two variants for each TF type (TF1 and TF2), were identified in genomic clones from each of the F1 hybrids between S. trutta and S. marmoratus. These clearly documented double heterozygote status at the TF loci. The real-time PCR data showed that each of the two TF types (TF1 and TF2) existed in one copy only and that the transcription of TF2 was considerably lower compared with TF1. Using FISH, hybridization signals were observed on two medium-sized acrocentric chromosomes of S. trutta karyotype. A TF type-specific PCR followed by a restriction analysis revealed the presence of two TF loci in the majority of analysed BAC clones. It was concluded that the TF gene is duplicated in the genome of S. trutta, and that the two TF loci are located adjacent to one another on the same chromosome. The differing transcription levels of TF1 and TF2 appear to depend on the corresponding promoter activity, which at least for TF2 seems to vary between different Salmo congeners.

Variation at two flowering time genes within and among populations of Arabidopsis thaliana: comparison with markers and traits

Flowering Locus C (FLC) and Frigida are two interacting genes controlling flowering time variation in Arabidopsis thaliana. Variation at these genes was surveyed in 12 A. thaliana populations sampled in France. These populations were also screened for variation at molecular markers [12 microsatellites and 19 cleaved amplified polymorphic sequence (CAPS) markers] and at seven quantitative traits measured with and without vernalization. Seven populations were highly polymorphic at markers (H(S) = 0.57 at microsatellites, 0.24 at CAPS) and showed heritable variation for bolting time and some other traits. Five populations were genetically fixed or nearly fixed. Q(ST) for bolting time without vernalization was significantly higher than F(ST), suggesting local divergent selection. One of the two haplotype groups at FLC (FLC(A)) was very predominant (frequency of 99%). The first exon of Frigida showed elevated nonsynonymous variation, and nine loss-of-function mutations were found throughout the gene. The association between loss-of-function and earlier bolting was confirmed. Overall, 18 Frigida haplotypes were detected. The pattern of variation at Frigida was largely similar to that found at markers and traits, with the same populations being fixed or highly diverse. Metapopulation dynamics is thus probably the main factor shaping genetic variation in A. thaliana. However, F(ST) for functional (FRI) vs. nonfunctional (FRI(Delta)) haplotypes was significantly higher than F(ST) at markers. This suggested that loss-of-function at Frigida is under local selection for flowering time.

Evolution of duplications in the transferrin family of proteins

The transferrin family is a group of proteins, defined by conserved amino acid motifs and putative function, found in both vertebrates and invertebrates. Included in this group are molecules known to bind iron, including serum transferrin, ovotransferrin, lactotransferrin, and melanotransferrin (MTF). Additional members of this family include inhibitor of carbonic anhydrase (ICA; mammals), major yolk protein (sea urchins), saxiphilin (frog), pacifastin (crayfish), and TTF-1 (algae). Most family members contain two lobes (N and C) of around 340 amino acids, the result of an ancient duplication event. In this article, we review the known functions of these proteins and speculate as to when the different homologs arose. From multiple-sequence alignments and neighbor-joining trees using 71 transferrin family sequences from 51 different species, including several novel sequences found in the Takifugu and Ciona genome databases, we conclude that melanotransferrins are much older (>670 MY) and more pervasive than previously thought, and the serum transferrin/melanotransferrin split may have occurred not long after lobe duplication. All subsequent duplication events diverged from the serum transferrin gene. The creation of such a large multiple-sequence alignment provides important information and could, in the future, highlight the role of specific residues in protein function.

Distribution of genetic variation in the growth hormone 1 gene in Atlantic salmon (Salmo salar) populations from Europe and North America

The level and hierarchical distribution of genetic variation in complete sequences of the Atlantic salmon (Salmo salar) growth hormone (GH1) gene were investigated in populations from Europe and North America with a view to inferring the major evolutionary forces affecting genetic variation at this locus. Seventeen polymorphic sites were identified in complete sequences from nine populations, with levels of noncoding (intron and untranslated region sequences) nucleotide diversity being similar to those observed in other species. No variation, however, was observed in exonic sequences, indicating that nucleotide diversity in the Atlantic salmon GH1 gene is three and 25 times less than that estimated for human and Drosophila coding sequences, respectively. This suggests that purifying selection is the predominant contemporary force controlling the molecular evolution of GH1 coding sequences. Comparison of haplotype relationships within and between populations indicated that differentiation between populations from Europe and North America was greater than within-continent comparisons. However, several haplotypes observed in the northernmost European populations were more similar to those observed in North American than to any other haplotypes observed in Europe. This is most likely to be a result of historical, rather than contemporary, gene flow. Neutrality test statistics, such as Tajima's D, were significantly positive in the European populations in which North American-like haplotypes were observed. Although a positive Tajima's D is commonly interpreted as the signal of balancing selection, a more likely explanation in this case is that either historical migration or ascertainment bias, rather than within population local adaptation, has given rise to an excess of intermediate frequency alleles.

Genomic structure and functional analysis of promoter region of somatolactin gene of sea bream (Sparus aurata)

Somatolactin (SL) is a pituitary hormone belonging to the growth hormone-prolactin family and is produced in the intermediate lobe of teleosts. The SL gene was isolated from a sea bream genomic library and found to be composed of 5 exons distributed within a 9-kb length of DNA. Sequence analysis of the proximal promoter region showed the presence of a classical TATA box located 59 bp upstream from the initial start ATG codon, 5 consensus sequences corresponding to the Pit-1 binding element, and a putative CREB site. In CHO cells cotransfected with the DNA from 2 plasmids, one encoding sea bream Pit-1 under Rous sarcoma virus long terminal repeat regulation and one encoding the SL promoter driving the expression of luciferase, Pit-1 was found to enhance the expression of luciferase. Only one Pit-1 binding site was necessary for enhancement. Analysis by immunoblots of in vitro culture of pituitaries of Sparus aurata showed that several agents, including estradiol, verapamil, and phorbol myristate acetate, had different inhibitory effects on SL and growth hormone released to the culture medium.

Genetic variation and population structure of moose (Alces alces) at neutral and functional DNA loci

Genetic variation was examined for moose (Alces alces) from Riding Mountain, Isle Royale, and Pukaskwa national parks; northwestern, Nipigon, northeastern, and central Ontario; New Brunswick; and Newfoundland. The national parks were identified as maintaining potentially different local selection pressures due to the absence of hunting and the presence or absence of the parasite Parelaphostrongylus tenuis. Genetic variation was estimated using neutral DNA markers, assessed by multilocus DNA fingerprinting and five microsatellite loci, and the functional antigen binding region (ARS) (exon 2) of the major histocompatibility complex (MHC) gene DRB. There was discordance in the allelic diversity observed at the neutral loci compared with the MHC DRB locus in a number of populations. Ontario populations demonstrated higher levels of variability at the neutral loci and relatively low levels at the DRB locus. Conversely, the Isle Royale population has the lowest genetic variability, consistent with a historic small founding event, at the neutral DNA markers and relatively high variability at the MHC gene. Relatively high levels of genetic variation at the DRB locus were observed in protected park populations concomitant with the absence of white-tailed deer (Odocoileus virginianus) or the parasite P. tenuis and an absence of hunting. Gene flow was observed among the neighboring geographic regions within Ontario, including Pukaskwa National Park, with evidence of isolation-by-distance among more distant regions within Ontario. The discordant patterns between DNA markers suggest that neutral DNA markers may not accurately reflect adaptive variation present at functional loci.

Natural selection and the genetic differentiation of coastal and Arctic populations of the Atlantic cod in northern Norway: a test involving nucleotide sequence variation at the pantophysin (PanI) locus

To examine the role of contemporary selection in maintaining significant allele frequency differences at the pantophysin (PanI) locus among populations of the Atlantic cod, Gadus morhua, in northern Norway, we sequenced 127 PanIA alleles sampled from six coastal and two Barents Sea populations. The distributions of variable sites segregating within the PanIA allelic class were then compared among the populations. Significant differences were detected in the overall frequencies of PanIA alleles among populations within coastal and Arctic regions that was similar in magnitude to heterogeneity in the distributions of polymorphic sites segregating within the PanIA allelic class. The differentiation observed at silent sites in the PanIA allelic class contradicts the predicted effects of widescale gene flow and suggests that postsettlement selection acting on cohorts cannot be responsible for the genetic differences described between coastal and Arctic populations. Our results suggest that the marked differences observed between coastal and Arctic populations of G. morhua in northern Norway at the PanI locus reflect the action of recent diversifying selection and that populations throughout the region may be more independent than suggested by previous studies.

The role of nuclear genes in intraspecific evolutionary inference: genealogy of the transferrin gene in the brown trout

Technical and biological hurdles have precluded the retrieval of nuclear gene genealogies within most species. Among these obstacles, the possibility of intragenic recombination is one of the most demanding challenges. We studied the utility of nuclear genes for intraspecific evolutionary inferences by selecting a nuclear gene that exhibits patterns of considerable geographic differentiation in the brown trout (Salmo trutta) species complex. Haplotype variation from a nucleotide sequence of approximately 3.7 kb encompassing a portion of the transferrin (TF) gene was surveyed in 31 brown trout individuals collected across the native Eurasian range. Statistically significant recombination and gene conversion events were detected. However, we showed that the substantial cladistic structure was not disrupted by recombination or gene conversion events and the additional structure was estimated to have emerged after those events. Because loci with unusually high levels of variation might indicate the presence of selection, we tested the hypothesis of neutrality and found some evidence for directional selection. The strong geographic signal observed in the TF genealogy, coupled with the current spatial distribution of electromorphs, gave us the ability to draw empirical phylogeographic inferences. We delineated the composition of current brown trout populations on the basis of 3,625 individuals electrophoretically scored for the TF locus. We hypothesized scenarios of historical radiation and dispersal events, thus providing new insights refining previous allozyme and mtDNA inferences. We infer that the most ancestral brown trout populations inhabited tributaries from the Black, Caspian, and Aral Sea drainages. An early radiation of the species occurred throughout the Mediterranean, followed by independent dispersal events from the Adriatic to the southernmost Iberian Atlantic and, more recently, a rapid expansion throughout most of the Atlantic drainages.

Applications of selective neutrality tests to molecular ecology

This paper reviews how statistical tests of neutrality have been used to address questions in molecular ecology are reviewed. The work consists of four major parts: a brief review of the current status of the neutral theory; a review of several particularly interesting examples of how statistical tests of neutrality have led to insight into ecological problems; a brief discussion of the pitfalls of assuming a strictly neutral model if it is false; and a discussion of some of the opportunities and problems that molecular ecologists face when using neutrality tests to study natural selection.

Complex signatures of natural selection at the Duffy blood group locus

The Duffy blood group locus (FY) has long been considered a likely target of natural selection, because of the extreme pattern of geographic differentiation of its three major alleles (FY*B, FY*A, and FY*O). In the present study, we resequenced the FY region in samples of Hausa from Cameroon (fixed for FY*O), Han Chinese (fixed for FY*A), Italians, and Pakistanis. Our goals were to characterize the signature of directional selection on FY*O in sub-Saharan Africa and to understand the extent to which natural selection has also played a role in the extreme geographic differentiation of the other derived allele at this locus, FY*A. The data from the FY region are compared with the patterns of variation observed at 10 unlinked, putatively neutral loci from the same populations, as well as with theoretical expectations from the neutral-equilibrium model. The FY region in the Hausa shows evidence of directional selection in two independent properties of the data (i.e., level of sequence variation and frequency spectrum), observations that are consistent with the FY*O mutation being the target. The Italian and Chinese FY data show patterns of variation that are very unusual, particularly with regard to frequency spectrum and linkage disequilibrium, but do not fit the predictions of any simple model of selection. These patterns may represent a more complex and previously unrecognized signature of positive selection.

Adaptive hitchhiking effects on genome variability

The continuing deluge of nucleotide polymorphism data is providing insights into the role of adaptation in shaping genome-wide patterns of variability and molecular evolution. Population genetic models in which linkage and selection interact (i.e. hitchhiking) predict that selection can leave 'footprints' in closely linked genomic regions. New analytical approaches show promise for distinguishing the signature of adaptation from that of several non-adaptive alternatives. Accounting for the effects of population structure and history poses a challenge for future investigations.