Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Molecular genetics of natural populations

Significant progress in evolutionary genetics has been made by studying, on the one hand, patterns of DNA sequence polymorphism and, on the other, genetic architecture of complex adaptive traits. However, connections between nucleotide variants under selection and adaptively relevant phenotypes are missing. Such connections can be established using precise gene replacement. We review the recent successful introduction of this technique to the analysis of two evolutionarily interesting loci--Odysseus and desaturase2. Both genes have subtle phenotypes that nevertheless could be identified using gene replacement, demonstrating that effects of naturally occurring alleles can be measured in the laboratory. This is an important first step in connecting statistical signatures of selection with adaptation in nature. More candidate genes involved in adaptation, for example, through cloning of genes responsible for reproductive isolation, now need to be identified. Molecular genetic manipulation, DNA polymorphism analysis, and field studies then have to be integrated to provide fresh insights into the mechanisms of evolutionary change.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Improved consensus network techniques for genome-scale phylogeny

Although recent studies indicate that estimating phylogenies from alignments of concatenated genes greatly reduces the stochastic error, the potential for systematic error still remains, heightening the need for reliable methods to analyze multigene data sets. Consensus methods provide an alternative, more inclusive, approach for analyzing collections of trees arising from multiple genes. We extend a previously described consensus network method for genome-scale phylogeny (Holland, B. R., K. T. Huber, V. Moulton, and P. J. Lockhart. 2004. Using consensus networks to visualize contradictory evidence for species phylogeny. Mol. Biol. Evol. 21:1459-1461) to incorporate additional information. This additional information could come from bootstrap analysis, Bayesian analysis, or various methods to find confidence sets of trees. The new methods can be extended to include edge weights representing genetic distance. We use three data sets to illustrate the approach: 61 genes from 14 angiosperm taxa and one gymnosperm, 106 genes from eight yeast taxa, and 46 members of a gene family from 15 vertebrate taxa.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Insight into the diversity and evolution of the cryptomonad nucleomorph genome

The cryptomonads are an enigmatic group of marine and freshwater unicellular algae that acquired their plastids through the engulfment and retention of a eukaryotic (secondary) endosymbiont. Together with the chlorarachniophyte algae, the cryptomonads are unusual in that they have retained the nucleus of their endosymbiont in a miniaturized form called a nucleomorph. The nucleomorph genome of the cryptomonad Guillardia theta has been completely sequenced and with only three chromosomes and a total size of 551 kb, is a model of nuclear genome compaction. Using this genome as a reference, we have investigated the structure and content of nucleomorph genomes in a wide range of cryptomonad algae. In this study, we have sequenced nine new cryptomonad nucleomorph 18S ribosomal DNA (rDNA) genes and four heat shock protein 90 (hsp90) gene fragments, and using pulsed-field gel electrophoresis and Southern hybridizations, have obtained nucleomorph genome size estimates for nine different species. We also used long-range polymerase chain reaction to obtain nucleomorph genomic fragments from Hanusia phi CCMP325 and Proteomonas sulcata CCMP704 that are syntenic with the subtelomeric region of nucleomorph chromosome I in G. theta. Our results indicate that (1) the presence of three chromosomes is a common feature of the nucleomorph genomes of these organisms, (2) nucleomorph genome size varies dramatically in the cryptomonads examined, (3) unidentified cryptomonad species CCMP1178 has the largest nucleomorph genome identified to date at approximately 845 kb, (4) nucleomorph genome size reductions appear to have occurred multiple times independently during cryptomonad evolution, (5) the relative positions of the 18S rDNA, ubc4, and hsp90 genes are conserved in three different cryptomonad genera, and (6) interchromosomal recombination appears to be rapidly changing the size and sequence of a repetitive subtelomeric region of the nucleomorph genome between the 18S rDNA and ubc4 loci. These results provide a glimpse into the genetic diversity of nucleomorph genomes in cryptomonads and set the stage for more comprehensive sequence-based studies in closely and distantly related taxa.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Accurate inference and estimation in population genomics

Both intra- and interspecific genomic comparisons have revealed local similarities in the level and frequency of mutational variation, as well as in patterns of gene expression. This autocorrelation between measurements leads to violations of assumptions of independence in many statistical methods, resulting in misleading and incorrect inferences. Here I show that autocorrelation can be due to many factors and is present across the genome. Using a one-dimensional spatial stochastic model, I further show how previous results can be employed to correct for autocorrelation along chromosomes in population and comparative genomics research. When multiple hypothesis tests are autocorrelated, I demonstrate that a simple correction can lead to increased power in statistical inference. I present a preliminary analysis of population genomic data from Drosophila simulans to show the ubiquity of autocorrelation and applicability of the methods proposed here.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Mutation rate and the cost of complexity

Two recent theoretical studies of adaptation suggest that more complex organisms tend to adapt more slowly. Specifically, in Fisher's "geometric" model of a finite population where multiple traits are under optimizing selection, the average progress ensuing from a single mutation decreases as the number of traits increases--the "cost of complexity." Here, I draw on molecular and histological data to assess the extent to which on a large phylogenetic scale, this predicted decrease in the rate of adaptation per mutation is mitigated by an increase in the number of mutations per generation as complexity increases. As an index of complexity for multicellular organisms, I use the number of visibly distinct types of cell in the body. Mutation rate is the product of mutational target size and population mutation rate per unit target. Despite much scatter, genome size appears to be positively correlated with complexity (as indexed by cell-type number), which along with other considerations suggests that mutational target size tends to increase with complexity. In contrast, effective population mutation rate per unit target appears to be negatively correlated with complexity. The net result is that mutation rate probably does tend to increase with complexity, although probably not fast enough to eliminate the cost of complexity.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Investigating the intron recognition mechanism in eukaryotes

Recent studies indicate that many introns, as well as the complex spliceosomal mechanism to remove them, were present early in eukaryotic evolution. This study examines intron and exon characteristics from annotations of whole genomes to investigate the intron recognition mechanism. Exon definition uses the exon as the unit of recognition, placing length constraints on the exon but not on the intron (allowing it a greater range of lengths). In contrast, intron definition uses the intron itself as the unit of recognition and thus removes constraints on internal exon length forced by the use of an exon definition mechanism. Thus, intron and exon lengths within a genome can reflect the constraints imposed by its splicing. This study shows that it is possible firstly to recover valid intron and exon information from genome annotation. We then compare internal intron and exon information from a range of eukaryotic genomes and investigate possible evolutionary length constraints on introns and exons and how they can impact on the intron recognition mechanism. Results indicate that exon definition-based mechanisms may predominate in vertebrates although the exact system in fish is expected to show some differences with the better characterized system from mammals. We also raise the possibility that the last common ancestor of plants and animals contained some type of exon definition and that this mechanism was replaced in some genes and lineages by intron definition, possibly as a result of intron loss and/or intron shortening.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Unusual evolution of interspersed repeat sequences in the Drosophila ananassae subgroup

New repeat sequences were found in the Drosophila ananassae genome sequence. They accounted for approximately 1.2% of the D. ananassae genome and were estimated to be more abundant in genomes of its closely related species belonging to the Drosophila bipectinata complex, whereas it was entirely absent in the Drosophila melanogaster genome. They were interspersed throughout euchromatic regions of the genome, usually as short tandem arrays of unit sequences, which were mostly 175-200 bp long with two distinct peaks at 180 and 189 bp in the length distribution. The nucleotide differences among unit sequences within the same array (locus) were much smaller than those between separate loci, suggesting within-locus concerted evolution. The phylogenetic tree of the repeat sequences from different loci showed that divergences between sequences from different chromosome arms occurred only at earlier stages of evolution, while those within the same chromosome arm occurred thereafter, resulting in the increase in copy number. We found RNA polymerase III promoter sequences (A box and B box), which play a critical role in retroposition of short interspersed elements. We also found conserved stem-loop structures, which are possibly associated with certain DNA rearrangements responsible for the increase in copy number within a chromosome arm. Such an atypical combination of characteristics (i.e., wide dispersal and tandem repetition) may have been generated by these different transposition mechanisms during the course of evolution.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Genome-wide associations between hybrid sterility QTL and marker transmission ratio distortion

Marker transmission ratio distortion (TRD) in genetic mapping populations is frequently ascribed to selection against allelic combinations that cause hybrid incompatibility. Accordingly, genomic regions of TRD should be nonrandomly associated (colocated) with loci that underlie hybrid incompatibility. To directly test this hypothesis, we evaluated the genome-wide qualitative and quantitative agreement between chromosomal regions exhibiting marker TRD and those known to contain hybrid incompatibility quantitative trait locus (QTL). Incompatibility data came from a near-isogenic line (NIL) analysis of pollen and seed sterility in a cross between two Solanum (formerly Lycopersicon) species. We assessed (1) whether these incompatibility loci are colocated with markers that show significant TRD in two earlier generations preceding these introgression lines and (2) whether the magnitude of marker distortion quantitatively matches the estimated strength of selection against each incompatibility locus. We found evidence that TRD regions are chromosomally colocated with hybrid incompatibility loci more frequently than is expected by chance: pollen sterility QTLs were most closely associated with distorted heterozygote frequencies in later-generation backcrosses. Nonetheless, there was no evidence for an association between TRD and seed sterility and little evidence of a quantitative association between the magnitude of marker TRD and the fitness effects of heterospecific alleles at each chromosomal location. We propose and test a model (the "dance partner" model) to explain several cases where regions of TRD are not associated with hybrid incompatibility loci. Under this model, some NILs containing greater than one heterospecific introgression may not express hybrid incompatibility phenotypes because they carry both appropriate genetic dance partners required for a fully functional interaction. Accordingly, negative interactions expressed in earlier backcross generations are masked in these double-introgression NILs. Based on this model, we identify the location of several new putative pairwise interactors underlying hybrid incompatibility in this species cross.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Coalescent-based estimation of population parameters when the number of demes changes over time

We expand a coalescent-based method that uses serially sampled genetic data from a subdivided population to incorporate changes to the number of demes and patterns of colonization. Often, when estimating population parameters or other parameters of interest from genetic data, the demographic structure and parameters are not constant over evolutionary time. In this paper, we develop a Bayesian Markov chain Monte Carlo method that allows for step changes in mutation, migration, and population sizes, as well as changing numbers of demes, where the times of these changes are also estimated. We show that in parameter ranges of interest, reliable estimates can often be obtained, including the historical times of parameter changes. However, posterior densities of migration rates can be quite diffuse and estimators somewhat biased, as reported by other authors.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Reconstructing the origins and dispersal of the Polynesian bottle gourd (Lagenaria siceraria)

The origin of the Polynesian bottle gourd (Lagenaria siceraria), an important crop species in prehistoric Polynesia, has remained elusive. Most recently, a South American origin has been favored as the bottle gourd could have been introduced from this continent with the sweet potato by Polynesian voyagers around A.D. 1,000. To test the hypothesis of an American origin for the Polynesian bottle gourd, we developed seven markers specific to bottle gourd (two chloroplast and five nuclear). The nuclear markers were developed using a new technique where polymorphic inter simple sequence repeat (ISSR) markers are converted into single-locus polymerase chain reaction and sequencing markers--an approach that will be useful for developing markers in other taxa. All seven markers were sequenced in 36 cultivars of bottle gourd from Asia, the Americas, and Polynesia. The results support a dual origin for the Polynesian bottle gourd: the chloroplast markers are exclusively of Asian origin, but the nuclear markers show alleles originating in both the Americas and Asia. Because hybridization of Polynesian bottle gourds with post-European introductions cannot be excluded, ancient DNA from archaeological material will be useful for further elucidating the prehistoric movements of this species in Polynesia. This work has implications not only for the dispersal of the Polynesian bottle gourd but also for the domestication and dispersal of the species as a whole.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Baleen whale phylogeny and a past extensive radiation event revealed by SINE insertion analysis

Baleen whales (suborder Mysticeti) comprise 11 extant species that are classified into four families. Although several phylogenetic hypotheses about these taxa have been proposed, their phylogenetic relationships remain confused. We addressed this problem using short interspersed repetitive element (SINE) insertion data, which now are regarded as almost ideal shared, derived characters at the molecular level. We reconstructed the phylogenetic relationships of baleen whales by characterizing 36 informative SINE loci. One of the intriguing conclusions is that balaenopterids and eschrichtiids radiated very rapidly during a very short evolutionary period. During this period, speciation occurred in balaenopterids and eschrichtiids while newly inserted SINE loci remains polymorphic. Later on, these SINEs were sorted incompletely into each lineage. Thus, there are now inconsistencies among species regarding the presence or absence of a given SINE. This is in sharp contrast to the phylogeny of toothed whales, for which no SINE inconsistencies have been found. Furthermore, we found monophyletic groupings between humpback and fin whales as well as between (sei+Bryde's) whales and blue whales, both of which have not previously been recognized. The comprehensive SINE insertion data, together with the mitochondrial DNA phylogeny that was recently completed (Sasaki, T., M. Nikaido, H. Healy et al. 2005. Mitochondrial phylogenetics and evolution of mysticete whales. Syst. Biol. 56:77-90; Rychel, A. L., T. W. Reeder, and A. Berta. 2004. Phylogeny of mysticete whales based on mitochondrial and nuclear data. Mol. Phylogenet. Evol. 32:892-901), provide a nearly complete picture of the evolutionary history of baleen whales.     

The human bitter taste receptor hTAS2R39 is the primary receptor for the bitterness of theaflavins

We purified several hundred mgs of four major theaflavins (theaflavin, theaflavin-3-O-gallate, theaflavin-3′-O-gallate, and theaflavin-3,3′-O-digallate). Among the 25 hTAS2Rs expressed in HEK293T cells, hTAS2R39 and hTAS2R14 were activated by theaflavins. Both hTAS2R39 and hTAS2R14 responded to theaflavin-3′-O-gallate. In addition, hTAS2R39 was activated by theaflavin and theaflavin-3,3′-O-gallate, but not by theaflavin-3-O-gallate. In contrast, hTAS2R14 responded to theaflavin-3-O-gallate.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Control of the false discovery rate applied to the detection of positively selected amino acid sites

In this article, we consider the probabilistic identification of amino acid positions that evolve under positive selection as a multiple hypothesis testing problem. The null hypothesis "H0,s: site s evolves under a negative selection or under a neutral process of evolution" is tested at each codon site of the alignment of homologous coding sequences. Standard hypothesis testing is based on the control of the expected proportion of falsely rejected null hypotheses or type-I error rate. As the number of tests increases, however, the power of an individual test may become unacceptably low. Recent advances in statistics have shown that the false discovery rate--in this case, the expected proportion of sites that do not evolve under positive selection among those that are estimated to evolve under this selection regime--is a quantity that can be controlled. Keeping the proportion of false positives low among the significant results generally leads to an increase in power. In this article, we show that controlling the false detection rate is relevant when searching for positively selected sites. We also compare this new approach to traditional methods using extensive simulations.     

Functional diversity of bitter taste receptor TAS2R16 in primates

In mammals, bitter taste is mediated by TAS2R genes, which belong to the large family of seven transmembrane G protein-coupled receptors. Because TAS2Rs are directly involved in the interaction between mammals and their dietary sources, it is likely that these genes evolved to reflect species-specific diets during mammalian evolution. Here, we investigated the sensitivities of TAS2R16s of various primates by using a cultured cell expression system, and found that the sensitivity of each primate species varied according to the ligand. Especially, the sensitivity of TAS2R16 of Japanese macaques to salicin was much lower than that of human TAS2R16, which was supported by behavioural tests. These results suggest the possibility that bitter-taste sensitivities evolved independently by replacing specific amino acid residues of TAS2Rs in different primate species to adapt to food items they use.     

Mapping of the porcine oestrogen receptor 2 gene and association study with litter size in Iberian pigs

Genes encoding the oestrogen receptors (ESR) are considered candidate genes for prolificacy traits due to the key role these molecules play in the regulation of reproductive physiology. In this paper, we report the assignment of the pig ESR2 gene to porcine chromosome 1 by radiation hybrid mapping. Most of the ESR2 cDNA was sequenced from Iberian pig ovarian RNA samples and one A/G single nucleotide polymorphism (SNP) was found at exon 5, being associated with a Met/Val substitution at position 949. This SNP was genotyped using a PCR-RFLP (Hsp92II) protocol and its potential effect on litter size was evaluated in two Iberian pig populations. However, no statistically significant association between the ESR2 polymorphism and litter size was found.     

Association of polymorphisms in the dopamine D4 receptor gene and the activity-impulsivity endophenotype in dogs

A variable number of tandem repeats (VNTR) polymorphism in exon 3 of the human dopamine D4 receptor gene ( DRD4 ) has been associated with attention deficit hyperactivity disorder (ADHD). Rodents possess no analogous repeat sequence, whereas a similar tandem repeat polymorphism of the DRD4 gene was identified in dogs, horses and chimpanzees. Here, we present a genetic association study of the DRD4 VNTR and the activity-impulsivity dimension of the recently validated dog-ADHD Rating Scale. To avoid false positives arising from population stratification, a single breed of dogs (German shepherd) was studied. Two DRD4 alleles (referred to as 2 and 3a ) were detected in this breed, and genotype frequencies were in Hardy–Weinberg equilibrium. For modelling distinct environmental conditions, 'pet' and 'police' German shepherds were characterized. Police German shepherds possessing at least one 3a allele showed significantly higher scores in the activity-impulsivity dimension of the dog-ADHD Rating Scale than dogs without this allele ( P  = 0.0180). This difference was not significant in pet German shepherds. To the best of our knowledge, this is the first report of an association between a candidate gene and a behaviour trait in dogs, and it reinforces the functional role of DRD4 exon 3 polymorphism.     

Polymorphism of the melatonin receptor MT1 gene and its relationship with seasonal reproductive activity in the Sarda sheep breed

The aim was to study the polymorphisms of the melatonin receptor 1A gene (MTNR1A) and its relationship with seasonal reproduction in the Sarda sheep breed. Four-thousand multiparous ewes reared under natural photoperiod were randomly chosen. Genomic DNA was extracted and subjected to PCR for the amplification of the main part of exon II of the ovine MTNR1A gene (GenBank U14109). PCR products were subjected to restriction enzymes MnlI and RsaI and placed into +/+, +/− or −/− group for MnlI and C/C, C/T or T/T group for RsaI. Samples were cloned and sequenced. The sequences were aligned with the U14109 sequence of GenBank. Data were subjected to allelic frequency analysis and to the χ² test in order to evaluate the link between genotype and reproductive activity. After MnlI digestion, allelic frequency was 0.78 for allele +and 0.22 for allele −; genotype frequency of the +/+ homozygote was 68%, 20.5% for +/− and 11.5% for −/−. After RsaI, allelic frequency was 0.66 for allele C and 0.34 for allele T; genotype frequency of the C/C homozygote was 53.5%, 26% for C/T and 20.5% for T/T. The population was in Hardy-Weinberg disequilibrium both for the MnlI and RsaI. Lambing frequency of +/+ genotype ewes was higher in the period September–December while for −/− genotype in January–April (P < 0.01). Lambing of C/C genotype ewes showed a higher frequency in September–December while for T/T genotype in January–April (P < 0.01). Results confirmed that the polymorphism of the MTNR1A locus was also present in the Sarda with a higher incidence of the +/+ and C/C genotypes. The animals that carried one of these two gene isoforms showed a not seasonal reproductive activity with the lambing period in September–December.