Comparative molecular population genetics of phycoerythrin locus in Prochlorococcus

As the only remainder type of phycobiliproteins in Prochlorococcus, the actual role of phycoerythrin still remains unknown. Previous studies revealed that two different forms of phycoerythrin gene were found in two ecotypes of Prochlorococcus that are specifically adapted to either high light (HL) or low light (LL) conditions. Here we analyze patterns of phycoerythrin nucleotide variation in the HL- and LL-Prochlorococcus populations. Our analyses reveal a significantly greater number of non-synonymous fixed substitutions in peB and peA than expected based on interspecific comparisons. This pattern of excess non-synonymous fixed substitutions is not seen in other five phycoerythrin-related genes (peZ/V/Y/T/S). Several neutrality statistical tests indicate an excess of rare frequency polymorphisms in the LL-Prochlorococcus data, but an excess of intermediate frequency polymorphisms in the HL-Prochlorococcus data. Distributions of the positively selected sites identified using the likelihood ratio test, when mapped onto the phycoerythrin tertiary structure, reveal that HL- and LL-phycoerythrin should be under different selective patterns. These findings may provide insights into the likely role of selection at the phycoerythrin locus and motivate further research to unveil the function of phycoerythrin in Prochlorococcus.

Molecular estimation of dispersal for ecology and population genetics

The dispersal process, by which individuals or other dispersing agents such as gametes or seeds move from birthplace to a new settlement locality, has important consequences for the dynamics of genes, individuals, and species. Many of the questions addressed by ecology and evolutionary biology require a good understanding of species' dispersal patterns. Much effort has thus been devoted to overcoming the difficulties associated with dispersal measurement. In this context, genetic tools have long been the focus of intensive research, providing a great variety of potential solutions to measuring dispersal. This methodological diversity is reviewed here to help (molecular) ecologists find their way toward dispersal inference and interpretation and to stimulate further developments.

Bacillus anthracis: Molecular taxonomy, population genetics, phylogeny and patho-evolution

Bacillus anthracis, the etiological agent of anthrax, manifests a particular bimodal lifestyle. This bacterial species alternates between short replication phases of 20-40 generations that strictly require infection of the host, normally causing death, interrupted by relatively long, mostly dormant phases as spores in the environment. Hence, the B. anthracis genome is highly homogeneous. This feature and the fact that strains from nearly all parts of the world have been analysed for canonical single nucleotide polymorphisms (canSNPs) and variable number tandem repeats (VNTRs) has allowed the development of molecular epidemiological and molecular clock models to estimate the age of major diversifications in the evolution of B. anthracis and to trace the global spread of this pathogen, which was mostly promoted by movement of domestic cattle with settlers and by international trade of contaminated animal products. From a taxonomic and phylogenetic point of view, B. anthracis is a member of the Bacillus cereus group. The differentiation of B. anthracis from B. cereus sensu strict, solely based on chromosomal markers, is difficult. However, differences in pathogenicity clearly differentiate B. anthracis from B. cereus and are marked by the strict presence of virulence genes located on the two virulence plasmids pXO1 and pXO2, which both are required by the bacterium to cause anthrax. Conversely, anthrax-like symptoms can also be caused by organisms with chromosomal features that are more closely related to B. cereus, but which carry these virulence genes on two plasmids that largely resemble the B. anthracis virulence plasmids. (C) 2011 Elsevier B.V. All rights reserved.

Living on the edge : Population genetics of Finno-Ugric-speaking humans in North Eurasia

The major aim of this thesis was to examine the origins and distribution of uniparental and autosomal genetic variation among the Finno-Ugric-speaking human populations living in Boreal and Arctic regions of North Eurasia. In more detail, I aimed to disentangle the underlying molecular and population genetic factors which have produced the patterns of uniparental and autosomal genetic diversity in these populations. Among Finno-Ugrics the genetic amalgamation and clinal distribution of West and East Eurasian gene pools were observed within uniparental markers. This admixture indicates that North Eurasia was colonized through Central Asia/ South Siberia by human groups already carrying both West and East Eurasian lineages. The complex combination of founder effects, gene flow and genetic drift underlying the genetic diversity of the Finno-Ugric- speaking populations were emphasized by low haplotype diversity within and among uniparental and biparental markers. A high prevalence of lactase persistence allele among the North Eurasian Finno- Ugric agriculturalist populations was also shown indicating a local adaptation to subsistence change with lactose rich diet. Moreover, the haplotype background of lactase persistence allele among the Finno- Ugric-speakers strongly suggested that the lactase persistence T-13910 mutation was introduced independently more than once to the North Eurasian gene pool. A significant difference in genetic diversity, haplotype structure and LD distribution within the cytochrome P450 CYP2C and CYP2D regions revealed the unique gene pool of the Finno-Ugric Saami created mainly by population genetic processes compared to other Europeans and sub-Saharan Mandenka population. From all studied populations the Saami showed also significantly the highest allele frequency of a CYP2C19 gene mutation causing variable drug reactions. The diversity patterns observed within CYP2C and CYP2D regions emphasize the strong effect of demographic history shaping genetic diversity and LD especially among such small and constant size populations as the Finno-Ugric-speaking Saami. Moreover, the increased LD in Saami due to genetic drift and/or admixture was shown to offer an advantage for further attempts to identify alleles associated to common complex pharmacogenetic traits.

Bayesian Inference for Retrospective Population Genetics Models Using Markov Chain Monte Carlo Methods

Genetics, the science of heredity and variation in living organisms, has a central role in medicine, in breeding crops and livestock, and in studying fundamental topics of biological sciences such as evolution and cell functioning. Currently the field of genetics is under a rapid development because of the recent advances in technologies by which molecular data can be obtained from living organisms. In order that most information from such data can be extracted, the analyses need to be carried out using statistical models that are tailored to take account of the particular genetic processes. In this thesis we formulate and analyze Bayesian models for genetic marker data of contemporary individuals. The major focus is on the modeling of the unobserved recent ancestry of the sampled individuals (say, for tens of generations or so), which is carried out by using explicit probabilistic reconstructions of the pedigree structures accompanied by the gene flows at the marker loci. For such a recent history, the recombination process is the major genetic force that shapes the genomes of the individuals, and it is included in the model by assuming that the recombination fractions between the adjacent markers are known. The posterior distribution of the unobserved history of the individuals is studied conditionally on the observed marker data by using a Markov chain Monte Carlo algorithm (MCMC). The example analyses consider estimation of the population structure, relatedness structure (both at the level of whole genomes as well as at each marker separately), and haplotype configurations. For situations where the pedigree structure is partially known, an algorithm to create an initial state for the MCMC algorithm is given. Furthermore, the thesis includes an extension of the model for the recent genetic history to situations where also a quantitative phenotype has been measured from the contemporary individuals. In that case the goal is to identify positions on the genome that affect the observed phenotypic values. This task is carried out within the Bayesian framework, where the number and the relative effects of the quantitative trait loci are treated as random variables whose posterior distribution is studied conditionally on the observed genetic and phenotypic data. In addition, the thesis contains an extension of a widely-used haplotyping method, the PHASE algorithm, to settings where genetic material from several individuals has been pooled together, and the allele frequencies of each pool are determined in a single genotyping.

Working at the interface of phylogenetics and population genetics: a biogeographical analysis of Triaenops spp. (Chiroptera: Hipposideridae).

New applications of genetic data to questions of historical biogeography have revolutionized our understanding of how organisms have come to occupy their present distributions. Phylogenetic methods in combination with divergence time estimation can reveal biogeographical centres of origin, differentiate between hypotheses of vicariance and dispersal, and reveal the directionality of dispersal events. Despite their power, however, phylogenetic methods can sometimes yield patterns that are compatible with multiple, equally well-supported biogeographical hypotheses. In such cases, additional approaches must be integrated to differentiate among conflicting dispersal hypotheses. Here, we use a synthetic approach that draws upon the analytical strengths of coalescent and population genetic methods to augment phylogenetic analyses in order to assess the biogeographical history of Madagascar's Triaenops bats (Chiroptera: Hipposideridae). Phylogenetic analyses of mitochondrial DNA sequence data for Malagasy and east African Triaenops reveal a pattern that equally supports two competing hypotheses. While the phylogeny cannot determine whether Africa or Madagascar was the centre of origin for the species investigated, it serves as the essential backbone for the application of coalescent and population genetic methods. From the application of these methods, we conclude that a hypothesis of two independent but unidirectional dispersal events from Africa to Madagascar is best supported by the data.

Application of SNPs for population genetics of nonmodel organisms:new opportunities and challenges

Recent improvements in the speed, cost and accuracy of next generation sequencing are revolutionizing the discovery of single nucleotide polymorphisms (SNPs). SNPs are increasingly being used as an addition to the molecular ecology toolkit in nonmodel organisms, but their efficient use remains challenging. Here, we discuss common issues when employing SNP markers, including the high numbers of markers typically employed, the effects of ascertainment bias and the inclusion of nonneutral loci in a marker panel. We provide a critique of considerations specifically associated with the application and population genetic analysis of SNPs in nonmodel taxa, focusing specifically on some of the most commonly applied methods.

Population genetics of C. Maenas : oceanography and larval dispersal

Decifrar a complexa interacção entre os ciclos de vida de espécies marinhas e a oceanografia revela-se fundamental para a compreensão do fluxo genético e da conectividade no meio marinho. Nas espécies marinhas com desenvolvimento indirecto o fluxo de genes entre populações depende da distância que separa as populações, bem como da interacção entre a duração do desenvolvimento larvar, do comportamento das larvas e dos padrões de circulação oceânica. A conectividade larvar influencia uma variedade de processos como a dinâmica de stocks e de populações, a distribuição e limites geográficos das espécies, a estrutura genética das populações e a dispersão de espécies invasivas e reveste-se consequentemente de uma importância fundamental na identificação das unidades populacionais evolucionariamente relevantes e para a gestão e conservação marinhas. Os marcadores genéticos e os Modelos Individuais Acoplados a Modelos Físico-Biológicos (“ICPBMs”) são actualmente ferramentas fundamentais para o estudo dos padrões de dispersão larvar e para avaliar o nível de conectividade populacional. A presente tese respeita à avaliação das escalas espaciais de conectividade de populações de uma espécie costeira, o caranguejo Carcinus maenas, e utiliza conjuntamente informação de marcadores genéticos, análise de séries temporais de fornecimento de larvas e um modelo numérico de circulação oceânica. O primeiro capítulo introduz a temática da conectividade em espécies marinhas e inclui algumas referências aos métodos moleculares, analíticos e de modelação seguidos ao longo da tese. Através da utilização de múltiplas ferramentas – avaliação da estrutura genética geográfica de C. maenas na sua distribuição nativa com recurso a marcadores de DNA (microssatélites) (Capítulo 2), avaliação da estrutura genética temporal das larvas que formam os eventos de fornecimento larvar à Ria de Aveiro, NW Portugal (Capítulo 3), descrição da variabilidade inter-anual do fornecimento larvar à Ria de Aveiro, NW Portugal (Capítulo 4) e validação de um modelo ICPBM que descreve os padrões observados de fornecimento (Capítulo 5) – esta tese espera poder contribuir para uma melhor compreensão dos mecanismos que regulam o fluxo de genes e a conectividade entre populações de organismos marinhos. No Capítulo 6 são apresentadas as principais conclusões da investigação. A análise genética com recurso a microssatélites indicou que as populações de C. maenas são geneticamente homogéneas ao longo de várias centenas de km, dentro da distribuição nativa da espécie. Paralelamente, não foram encontrados indícios da existência de reprodução por “sweepstakes” em C. maenas de populações da costa oeste da Península Ibérica, visto que não se obtiveram diferenças genéticas significativas entre os eventos larvares. Também não se encontrou qualquer estrutura familiar entre as larvas que formam cada episódio de fornecimento, e não houve nenhuma redução significativa da variabilidade genética das larvas quando comparada com a de caranguejos adultos. A análise de séries temporais de suprimento de larvas na Ria de Aveiro em cinco anos estudados indica que este é um fenómeno episódico e variável, sendo os maiores episódios de fornecimento coincidentes com as marés vivas e acentuados por fortes ventos de sul. O modelo ICPBM foi validado com sucesso e parece fornecer uma estimativa realística das escalas espaciais e temporais de dispersão larvar, de acordo com as observações da estrutura genética e da ausência de reprodução por “sweepstake” em C. maenas da costa oeste da Península Ibérica

Epidemiology and population genetics of Podosphaera fusca and Golovinomyces orontii, causal agents of cucurbit powdery mildew

In 2010, 2011 and 2012 growing seasons, the occurrence of the ascomycetes Podosphaera fusca and Golovinomyces orontii, causal agents of powdery mildew disease, was monitored on cultivated cucurbits located in Bologna and Mantua provinces to determine the epidemiology of the species. To identify the pathogens, both morphological and molecular identifications were performed on infected leaf samples and a Multiplex-PCR was performed to identify the mating type genes of P. fusca isolates. The investigations indicated a temporal succession of the two species with the earlier infections caused by G. orontii, that seems to be the predominant species till the middle of July when it progressively disappears and P. fusca becomes the main species infecting cucurbits till the end of October. The temporal variation is likely due to the different overwintering strategies of the two species instead of climatic conditions. Only chasmothecia of P. fusca were recorded and mating type alleles ratio tended to be 1:1. Considering that only chasmothecia of P. fusca were found, molecular-genetic analysis were carried out to find some evidence of recombination within this species by MLST and AFLP methods. Surprisingly, no variations were observed within isolates for the 8 MLST markers used. According to this result, AFLP analysis showed a high similarity within isolates, with SM similarity coefficient ranging between 0.91-1.00 and also, sequencing of 12 polymorphic bands revealed identity to some gene involved in mutation and selection. The results suggest that populations of P. fusca are likely to be a clonal population with some differences among isolates probably due to agricultural practices such as fungicides treatments and cultivated hosts. Therefore, asexual reproduction, producing a lot of fungal biomass that can be easily transported by wind, is the most common and useful way to the spread and colonization of the pathogen.

Impact of anthropogenic stressors on the population genetics of the common bioindicator Dreissena polymorpha (Pallas, 1771)

Toxicant inputs from agriculture, industry and human settlements have been shown to severely affect freshwater ecosystems. Pollution can lead to changes in population genetic patterns through various genetic and stochastic processes. In my thesis, I investigated the impact of anthropogenic stressors on the population genetics of the zebra mussel Dreissena polymorpha. In order to analyze the genetics of zebra mussel populations, I isolated five new highly polymorphic microsatellite loci. Out of those and other already existing microsatellite markers for this species, I established a robust marker set of six microsatellite loci for D. polymorpha. rnMonitoring the biogeographical background is an important requirement when integrating population genetic measures into ecotoxicological studies. I analyzed the biogeographical background of eleven populations in a section of the River Danube (in Hungary and Croatia) and some of its tributaries, and another population in the River Rhine as genetic outgroup. Moreover, I measured abiotic water parameters at the sampling sites and analyzed if they were correlated with the genetic parameters of the populations. The genetic differentiation was basically consistent with the overall biogeographical history of the populations in the study region. However, the genetic diversity of the populations was not influenced by the geographical distance between the populations, but by the environmental factors oxygen and temperature and also by other unidentified factors. I found strong evidence that genetic adaptation of zebra mussel populations to local habitat conditions had influenced the genetic constitution of the populations. Moreover, by establishing the biogeographical baseline of molecular variance in the study area, I laid the foundation for interpreting population genetic results in ecotoxicological experiments in this region.rnIn a cooperation project with the Department of Zoology of the University of Zagreb, I elaborated an integrated approach in biomonitoring with D. polymorpha by combining the analysis techniques of microsatellite analysis, Comet assay and micronucleus test (MNT). This approach was applied in a case study on freshwater contamination by an effluent of a wastewater treatment plant (WWTP) in the River Drava (Croatia) and a complementary laboratory experiment. I assessed and compared the genetic status of two zebra mussel populations from a contaminated and a reference site. Microsatellite analysis suggested that the contaminated population had undergone a genetic bottleneck, caused by random genetic drift and selection, whereas a bottleneck was not detected in the reference population. The Comet assay did not indicate any difference in DNA damage between the two populations, but MNT revealed that the contaminated population had an increased percentage of micronuclei in hemocytes in comparison to the reference population. The laboratory experiment with mussels exposed to municipal wastewater revealed that mussels from the contaminated site had a lower percentage of tail DNA and a higher percentage of micronuclei than the reference population. These differences between populations were probably caused by an overall decreased fitness of mussels from the contaminated site due to genetic drift and by an enhanced DNA repair mechanism due to adaptation to pollution in the source habitat. Overall, the combination of the three biomarkers provided sufficient information on the impact of both treated and non-treated municipal wastewater on the genetics of zebra mussels at different levels of biological organization.rnIn my thesis, I could show that the newly established marker set of six microsatellite loci provided reliable and informative data for population genetic analyses of D. polymorpha. The adaptation of the analyzed zebra mussel populations to the local conditions of their habitat had a strong influence on their genetic constitution. We found evidence that the different genetic constitutions of two populations had influenced the outcome of our ecotoxicological experiment. Overall, the integrated approach in biomonitoring gave comprehensive information about the impact of both treated and non-treated municipal wastewater on the genetics of zebra mussels at different levels of biological organization and was well practicable in a first case study.

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.