Hidden Biodiversity in an Ecologically Important Freshwater Amphipod: Differences in Genetic Structure between Two Cryptic Species

Cryptic species, i.e. species that are morphologically hard to distinguish, have been detected repeatedly in various taxa and ecosystems. In order to evaluate the importance of this finding, we have to know in how far cryptic species differ in various aspects of their biology. The amphipod Gammarus fossarum is a key invertebrate in freshwater streams and contains several cryptic species. We examined the population genetic structure, genetic diversity and demographic history of two of them (type A and type B) using microsatellite markers and asked whether they show significant differences. We present results of population genetic analyses based on a total of 37 populations from the headwaters of two major European drainages, Rhine and Rhone. We found that, in both species, genetic diversity was geographically structured among and within drainages. For type A in the Rhine and type B in the Rhone, we detected significant patterns of isolation by distance. The increase of genetic differentiation with geographical distance, however, was much higher in type A than in type B. This result indicates substantial interspecific differences in population history and/or the extent of current gene flow between populations. In the Rhine, type B does not show evidence of isolation by distance, and population differentiation is relatively low across hundreds of kilometres. The majority of these populations also show signatures of recent bottlenecks. These patterns are consistent with a recent expansion of type B into the Rhine drainage. In summary, our results suggest considerable and previously unrecognized interspecific differences in the genetic structure of these cryptic keystone species.

Parentage analysis with genetic markers in natural populations. I. The expected proportion of offspring with unambiguous paternity.

Recent studies indicate that polymorphic genetic markers are potentially helpful in resolving genealogical relationships among individuals in a natural population. Genetic data provide opportunities for paternity exclusion when genotypic incompatibilities are observed among individuals, and the present investigation examines the resolving power of genetic markers in unambiguous positive determination of paternity. Under the assumption that the mother for each offspring in a population is unambiguously known, an analytical expression for the fraction of males excluded from paternity is derived for the case where males and females may be derived from two different gene pools. This theoretical formulation can also be used to predict the fraction of births for each of which all but one male can be excluded from paternity. We show that even when the average probability of exclusion approaches unity, a substantial fraction of births yield equivocal mother-father-offspring determinations. The number of loci needed to increase the frequency of unambiguous determinations to a high level is beyond the scope of current electrophoretic studies in most species. Applications of this theory to electrophoretic data on Chamaelirium luteum (L.) shows that in 2255 offspring derived from 273 males and 70 females, only 57 triplets could be unequivocally determined with eight polymorphic protein loci, even though the average combined exclusionary power of these loci was 73%. The distribution of potentially compatible male parents, based on multilocus genotypes, was reasonably well predicted from the allele frequency data available for these loci. We demonstrate that genetic paternity analysis in natural populations cannot be reliably based on exclusionary principles alone. In order to measure the reproductive contributions of individuals in natural populations, more elaborate likelihood principles must be deployed.

Population genetics of VNTR loci and their applications in evolutionary studies

Variable number of tandem repeats (VNTR) are genetic loci at which short sequence motifs are found repeated different numbers of times among chromosomes. To explore the potential utility of VNTR loci in evolutionary studies, I have conducted a series of studies to address the following questions: (1) What are the population genetic properties of these loci? (2) What are the mutational mechanisms of repeat number change at these loci? (3) Can DNA profiles be used to measure the relatedness between a pair of individuals? (4) Can DNA fingerprint be used to measure the relatedness between populations in evolutionary studies? (5) Can microsatellite and short tandem repeat (STR) loci which mutate stepwisely be used in evolutionary analyses?^ A large number of VNTR loci typed in many populations were studied by means of statistical methods developed recently. The results of this work indicate that there is no significant departure from Hardy-Weinberg expectation (HWE) at VNTR loci in most of the human populations examined, and the departure from HWE in some VNTR loci are not solely caused by the presence of population sub-structure.^ A statistical procedure is developed to investigate the mutational mechanisms of VNTR loci by studying the allele frequency distributions of these loci. Comparisons of frequency distribution data on several hundreds VNTR loci with the predictions of two mutation models demonstrated that there are differences among VNTR loci grouped by repeat unit sizes.^ By extending the ITO method, I derived the distribution of the number of shared bands between individuals with any kinship relationship. A maximum likelihood estimation procedure is proposed to estimate the relatedness between individuals from the observed number of shared bands between them.^ It was believed that classical measures of genetic distance are not applicable to analysis of DNA fingerprints which reveal many minisatellite loci simultaneously in the genome, because the information regarding underlying alleles and loci is not available. I proposed a new measure of genetic distance based on band sharing between individuals that is applicable to DNA fingerprint data.^ To address the concern that microsatellite and STR loci may not be useful for evolutionary studies because of the convergent nature of their mutation mechanisms, by a theoretical study as well as by computer simulation, I conclude that the possible bias caused by the convergent mutations can be corrected, and a novel measure of genetic distance that makes the correction is suggested. In summary, I conclude that hypervariable VNTR loci are useful in evolutionary studies of closely related populations or species, especially in the study of human evolution and the history of geographic dispersal of Homo sapiens. (Abstract shortened by UMI.) ^

Genetic diversity, phenotypic variation and local adaptation in the alpine landscape: case studies with alpine plant species

Plant survival in alpine landscapes is constantly challenged by the harsh and often unpredictable environmental conditions. Steep environmental gradients and patchy distribution of habitats lead to small size and spatial isolation of populations and restrict gene flow. Agricultural land use has further increased the diversity of habitats below and above the treeline. We studied the consequences of the highly structured alpine landscape for evolutionary processes in four study plants: Epilobium fleischeri, Geum reptans, Campanula thyrsoides and Poa alpina. The main questions were: (1) How is genetic diversity distributed within and among populations and is it affected by altitude, population size or land use? (2) Do reproductive traits such as allocation to sexual or vegetative reproduction vary with altitude or land use? Furthermore, we studied if seed weight increases with altitude. Within-population genetic diversity of the four species was high and mostly not related to altitude and population size. Nevertheless, genetic differentiation among populations was pronounced and strongly increasing with distance. In Poa alpina genetic diversity was affected by land use. Results suggest considerable genetic drift among populations of alpine plants. Reproductive allocation was affected by altitude and land use in Poa alpina and by succession in Geum reptans. Seed weight was usually higher in alpine species than in related lowland species. We conclude that the evolutionary potential to respond to global change is mostly intact in alpine plants, even at high altitude. Phenotypic variability is shaped by adaptive as well as by random evolutionary processes; moreover plastic responses to growth conditions seem to be crucial for survival of plants in the alpine landscape.

Global diversity, population stratification and selection of human copy number variation

In order to explore the diversity and selective signatures of duplication and deletion human copy number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single nucleotide variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load.

Genetic composition, genetic diversity, and small-scale environmental variation matter for the experimental reintroduction of a rare plant

Aims Reintroduction has become an important tool for the management of endangered plant species. We tested the little-explored effects of small-scale environmental variation, genotypic composition (i.e. identity of genotypes), and genotypic diversity on the population survival of the regionally rare clonal plant Ranunculus reptans. For this species of periodically inundated lakeshores genetic differentiation had been reported between populations and between short-flooded and long-flooded microsites within populations.Methods We established 306 experimental test populations at a previously unoccupied lake shore, comprising either monocultures of 32 genotypes, mixtures of genotypes within populations or mixtures of genotypes between populations. In 2000, three years after planting out at the experimental site, a long-lasting flood caused the death of half of the experimental populations. In 2003, an extreme drought resulted in the lowest summer water levels ever measured.Important findings Despite these climatic extremes, 27 of the established populations survived until the end of the experiment in December 2003. The success of experimental populations largely differed between microsites. Moreover, the success of genotype monocultures depended on genotype and source population. Genetic differentiation between microsites played a minor role for the success of reintroduction. After the flood, populations planted with genotypes from different source populations increased in abundance, whereas populations with genotypes from single source populations and genotype monocultures decreased. We conclude that sources for reintroductions need to be selected carefully. Moreover, mixtures of plants from different populations appear to be the best choice for successful reintroduction, at least in unpredictably varying environments.

Low but contrasting neutral genetic differentiation shaped by winter temperature in European great tits

Gene flow is usually thought to reduce genetic divergence and impede local adaptation by homogenising gene pools between populations. However, evidence for local adaptation and phenotypic differentiation in highly mobile species, experiencing high levels of gene flow, is emerging. Assessing population genetic structure at different spatial scales is thus a crucial step towards understanding mechanisms underlying intraspecific differentiation and diversification. Here, we studied the population genetic structure of a highly mobile species – the great tit Parus major – at different spatial scales. We analysed 884 individuals from 30 sites across Europe including 10 close-by sites (< 50 km), using 22 microsatellite markers. Overall we found a low but significant genetic differentiation among sites (FST = 0.008). Genetic differentiation was higher, and genetic diversity lower, in south-western Europe. These regional differences were statistically best explained by winter temperature. Overall, our results suggest that great tits form a single patchy metapopulation across Europe, in which genetic differentiation is independent of geographical distance and gene flow may be regulated by environmental factors via movements related to winter severity. This might have important implications for the evolutionary trajectories of sub-populations, especially in the context of climate change, and calls for future investigations of local differences in costs and benefits of philopatry at large scales.

Bayesian models for the analysis of genetic structure when populations are correlated

Motivation: Population allele frequencies are correlated when populations have a shared history or when they exchange genes. Unfortunately, most models for allele frequency and inference about population structure ignore this correlation. Recent analytical results show that among populations, correlations can be very high, which could affect estimates of population genetic structure. In this study, we propose a mixture beta model to characterize the allele frequency distribution among populations. This formulation incorporates the correlation among populations as well as extending the model to data with different clusters of populations. Results: Using simulated data, we show that in general, the mixture model provides a good approximation of the among-population allele frequency distribution and a good estimate of correlation among populations. Results from fitting the mixture model to a dataset of genotypes at 377 autosomal microsatellite loci from human populations indicate high correlation among populations, which may not be appropriate to neglect. Traditional measures of population structure tend to over-estimate the amount of genetic differentiation when correlation is neglected. Inference is performed in a Bayesian framework.

Multilocus analysis of extracellular putative virulence proteins made by group A Streptococcus: Population genetics, human serologic response, and gene transcription

Species of pathogenic microbes are composed of an array of evolutionarily distinct chromosomal genotypes characterized by diversity in gene content and sequence (allelic variation). The occurrence of substantial genetic diversity has hindered progress in developing a comprehensive understanding of the molecular basis of virulence and new therapeutics such as vaccines. To provide new information that bears on these issues, 11 genes encoding extracellular proteins in the human bacterial pathogen group A Streptococcus identified by analysis of four genomes were studied. Eight of the 11 genes encode proteins with a LPXTG(L) motif that covalently links Gram-positive virulence factors to the bacterial cell surface. Sequence analysis of the 11 genes in 37 geographically and phylogenetically diverse group A Streptococcus strains cultured from patients with different infection types found that recent horizontal gene transfer has contributed substantially to chromosomal diversity. Regions of the inferred proteins likely to interact with the host were identified by molecular population genetic analysis, and Western immunoblot analysis with sera from infected patients confirmed that they were antigenic. Real-time reverse transcriptase–PCR (TaqMan) assays found that transcription of six of the 11 genes was substantially up-regulated in the stationary phase. In addition, transcription of many genes was influenced by the covR and mga trans-acting gene regulatory loci. Multilocus investigation of putative virulence genes by the integrated approach described herein provides an important strategy to aid microbial pathogenesis research and rapidly identify new targets for therapeutics research.

Molecular markers reveal that population structure of the human pathogen Candida albicans exhibits both clonality and recombination.

The life history of Candida albicans presents an enigma: this species is thought to be exclusively asexual, yet strains show extensive phenotypic variation. To address the population genetics of C. albicans, we developed a genetic typing method for codominant single-locus markers by screening randomly amplified DNA for single-strand conformation polymorphisms. DNA fragments amplified by arbitrary primers were initially screened for single-strand conformation polymorphisms and later sequenced using locus-specific primers. A total of 12 single base mutations and insertions were detected from six out of eight PCR fragments. Patterns of sequence-level polymorphism observed for individual strains detected considerable heterozygosity at the DNA sequence level, supporting the view that most C. albicans strains are diploid. Population genetic analyses of 52 natural isolates from Duke University Medical Center provide evidence for both clonality and recombination in C. albicans. Evidence for clonality is supported by the presence of several overrepresented genotypes, as well as by deviation of genotypic frequencies from random (Hardy-Weinberg) expectations. However, tests for nonrandom association of alleles across loci reveal less evidence for linkage disequilibrium than expected for strictly clonal populations. Although C. albicans populations are primarily clonal, evidence for recombination suggests that sexual reproduction or some other form of genetic exchange occurs in this species.

Mitochondrial control region diversity of the houbara bustard Chlamydotis undulata complex and genetic structure along the Atlantic seaboard of North Africa

The houbara bustard, Chlamydotis undulata, is a declining cryptic desert bird whose range extends from North Africa to Central Asia. Three subspecies are currently recognized by geographical distribution and morphology: C.u.fuertaventurae, C.u.undulata and C.u.macqueenii. We have sequenced 854 bp of mitochondrial control region from 73 birds to describe their population genetic structure with a particular sampling focus on the connectivity between C.u.fuertaventurae and C.u.undulata along the Atlantic seaboard of North Africa. Nucleotide and haplotypic diversity varied among the subspecies being highest in C.u.undulata, lowest in C.u.fuertaventurae and intermediate in C.u.macqueenii. C.u.fuertaventurae and C.u.undulata are paraphyletic and an average nucleotide divergence of 2.08% splits the later from C.u.macqueenii. We estimate that C.u.fuertaventurae and C.u.undulata split from C.u.macqueenii approximately 430 000 years ago. C.u.fuertaventurae and C.u.undulata are weakly differentiated (F-ST = 0.27, N-m = 1.3), indicative of a recent shared history. Archaeological evidence indicates that houbara bustards have been present on the Canary Islands for 130-170 000 years. However, our genetic data point to a more recent separation of C.u.fuertaventurae and C.u.undulata at around 20-25 000 years. Concordant archaeological, climatic opportunities for colonization and genetic data point to a scenario of: (i) initial colonization of the Canary Islands about 130 000 years ago; (ii) a period of secondary contact 19-30 000 years ago homogenizing any pre-existing genetic structure followed by; (iii) a period of relative isolation that persists today.

The genetic structure of Australian populations of Mycosphaerella musicola suggests restricted gene flow at the continental scale

Mycosphaerello musicolo causes Sigatoka disease of banana and is endemic to Australia. The population genetic structure of M. musicola in Australia was examined by applying single-copy restriction fragment length polymorphism probes to hierarchically sampled populations collected along the Australian cast coast. The 363 isolates studied were from 16 plantations at 12 sites in four different regions, and comprised 11 populations. These populations displayed moderate levels of gene diversity (H = 0.142 to 0.369) and similar levels of genotypic richness and evenness. Populations were dominated by unique genotypes, but isolates sharing the same genotype (putative clones) were detected. Genotype distribution was highly localized within each population, and the majority of putative clones were detected for isolates sampled from different sporodochia in the same lesion or different lesions on a plant. Multilocus gametic disequilibrium tests provided further evidence of a degree of clonality within the populations at the plant scale. A complex pattern of population differentiation was detected for M. musicola in Australia. Populations sampled from plantations outside the two major production areas were genetically very different to all other populations. Differentiation was much lower between populations of the two major production areas, despite their geographic separation of over 1,000 km. These results suggest low gene flow at the continental scale due to limited spore dispersal and the movement of infected plant material.

Use of community genetic screening to prevent HFE-associated hereditary haemochromatosis

HFE-associated hereditary haemochromatosis is a recessive, iron-overload disorder that affects about one in 200 north Europeans and that can be easily prevented. However, genetic screening for this disease is controversial, and so we assessed whether such screening was suitable for communities. Cheek-brush screening for the Cys282Tyr HFE mutation was offered to individuals in the workplace. Outcomes were assessed by questionnaires before and after testing. 11307 individuals were screened. We recorded no increase in anxiety. in individuals who were homozygous for the Cys282Tyr mutation or non-homozygous. Self-reported tiredness before testing was significantly higher in homozygous participants than in non-homozygous participants (chi(2) test, p=0.029). Of the 47 homozygous individuals identified, 46 have taken steps to treat or prevent iron accumulation. Population genetic screening for HFE-associated hereditary haemochromatosis can be practicable and acceptable.

Contrasting quantitative traits and neutral genetic markers for genetic resource assessment of Mesoamerican Cedrela odorata

We compared within-population variability and degree of population differentiation for neutral genetic markers (RAPDS) and eight quantitative traits in Central American populations of the endangered tree, Cedrela odorata. Whilst population genetic diversity for neutral markers (Shannon index) and quantitative traits (heritability, coefficient of additive genetic variation) were uncorrelated, both marker types revealed strong differentiation between populations from the Atlantic coast of Costa Rica and the rest of the species' distribution. The degree of interpopulation differentiation was higher for RAPD markers (F-ST 0.67 for the sampled Mesoamerican range) than for quantitative traits (Q(ST) = 0.30). Hence, the divergence in quantitative traits was lower than could have been achieved by genetic drift alone, suggesting that balancing selection for similar phenotypes in different populations of this species. Nevertheless, a comparison of pair-wise estimates of population differentiation in neutral genetic markers and quantitative traits revealed a strong positive correlation (r = 0.66) suggesting that, for C. odorata, neutral marker divergence could be used as a surrogate for adaptive gene divergence for conservation planning. The utility of this finding and suggested further work are discussed.

Meta-analysis of the population and phenotypic expression of CYP2C9/2C19 polymorphism on drug metabolism in different ethnicities

The term "pharmacogenetics" has been defined as the scientific study of inherited factors that affect the human drug response. Many pharmacogenetie studies have been published since 1995 and have focussed on the principal enzyme family involved in drug metabolism, the cytochrome P450 family, particularly cytochrome P4502C9 and 2C19. In order to investigate the pharmacogenetic aspect of pharmacotherapy, the relevant studies describing the association of pharmacogenetic factor(s) in drug responses must be retrieved from existing literature using a systematic review approach. In addition, the estimation of variant allele prevalence for the gene under study between different ethnic populations is important for pharmacogenetic studies. In this thesis, the prevalence of CYP2C9/2C19 alleles between different ethnicities has been estimated through meta-analysis and the population genetic principle. The clinical outcome of CYP2C9/2C19 allelic variation on the pharmacotherapy of epilepsy has been investigated; although many new antiepileptic drugs have been launched into the market, carbamazepine, phenobarbital and phenytoin are still the major agents in the pharmacotherapy of epilepsy. Therefore, phenytoin was chosen as a model AED and the effect of CYP2C9/2C19 genetic polymorphism on phenytoin metabolism was further examined.An estimation of the allele prevalence was undertaken for three CYP2C9/2C19 alleles respectively using a meta-analysis of studies that fit the Hardy-Weinberg equilibrium. The prevalence of CYP2C9*1 is approximately 81%, 96%, 97% and 94% in Caucasian, Chinese, Japanese, African populations respectively; the pooled prevalence of CYP2C19*1 is about 86%, 57%, 58% and 85% in these ethnic populations respectively. However, the studies of association between CYP2C9/2C19 polymorphism and phenytoin metabolism failed to achieve any qualitative or quantitative conclusion. Therefore, mephenytoin metabolism was examined as a probe drug for association between CYP2C19 polymorphism and mephenytoin metabolic ratio. Similarly, analysis of association between CYP2C9 polymorphism and warfarin dose requirement was undertaken.It was confirmed that subjects carrying two mutated CYP2C19 alleles have higher S/R mephenytoin ratio due to deficient CYP2C19 enzyme activity. The studies of warfarin and CYP2C9 polymorphism did not provide a conclusive result due to poor comparability between studies.The genetic polymorphism of drug metabolism enzymes has been studied extensively, however other genetic factors, such as multiple drug resistance genes (MDR) and genes encoding ion channels, which may contribute to variability in function of drug transporters and targets, require more attention in future pharmacogenetic studies of antiepileptic drugs.