LOW GENETIC DIVERSITY IN Wolbachia-INFECTED Culex quinquefasciatus (DIPTERA: CULICIDAE) FROM BRAZIL AND ARGENTINA

Culex quinquefasciatus is a vector of human pathogens, including filarial nematodes and several viruses. Although its epidemiological relevance is known to vary across geographical regions, an understanding of its population genetic structure is still incipient. In light of this, we evaluated the genetic diversity of Cx. quinquefasciatus and Cx. pipiens x Cx. quinquefasciatus hybrids collected from nine localities in Brazil and one site in Argentina. We used mitochondrial genes cox1 and nd4, along with the coxA and wsp genes of the maternally-inherited Wolbachia endosymbiont. The nd4 fragment was invariant between samples, whilst cox1 exhibited four haplotypes that separated two types of Cx. quinquefasciatus, one clustered in southern Brazil. Low sequence diversity was generally observed, being discussed. Both Brazilian and Argentinian mosquitoes were infected with a single Wolbachia strain. As reported in previous studies with these populations, cox1 and nd4 diversity is not congruent with the population structure revealed by nuclear markers or alar morphology. Future Cx. quinquefasciatus research should, if possible, evaluate mtDNA diversity in light of other markers.

Mannose-binding lectin and MBL-associated serine protease-2 gene polymorphisms in a Brazilian population from Rio de Janeiro

Mannose-binding lectin (MBL) is a protein able to bind to carbohydrate patterns on pathogen membranes; upon MBL binding, its associated serine protease MBL-associated serine protease type 2 (MASP2) is autoactivated, promoting the activation of complement via the lectin pathway. For both MBL2 and MASP2 genes, the frequencies of polymorphisms are extremely variable between different ethnicities, and this aspect has to be carefully considered when performing genetic studies. While polymorphisms in the MBL-encoding gene (MBL2) have been associated, depending upon ethnicity, with several diseases in different populations, little is known about the distribution of MASP2 gene polymorphisms in human populations. The aim of our study was thus to determine the frequencies of MBL2 (exon 1 and promoter) and MASP2 (p.D371Y) polymorphisms in a Brazilian population from Rio de Janeiro. A total of 294 blood donor samples were genotyped for 27 polymorphisms in the MBL2 gene by direct sequencing of a region spanning from the promoter polymorphism H/L rs11003125 to the rs1800451 polymorphism (at codon 57 in the first exon of the gene). Genotyping for MASP2 p.D371Y was carried out using fluorogenic probes. To our knowledge, this is the first study reporting the prevalence of the MASP2 p.D371Y polymorphism in a Brazilian population. The C allele frequency 39% is something intermediate between the reported 14% in Europeans and 90% in Sub-Saharan Africans. MBL2 polymorphisms frequencies were quite comparable to those previously reported for admixed Brazilians. Both MBL2 and MASP2 polymorphisms frequencies reported in our study for the admixed Brazilian population are somehow intermediate between those reported in Europeans and Africans, reflecting the ethnic composition of the southern Brazilian population, estimated to derive from an admixture of Caucasian (31%), African (34%) and Native American (33%) populations. In conclusion, our population genetic study describes the frequencies of MBL2 and MASP2 functional SNPs in a population from Rio de Janeiro, with the aim of adding new information concerning the distribution of these SNPs in a previously unanalysed Brazilian population, thus providing a new genetic tool for the evaluation of the association of MBL2 and MASP2 functional SNPs with diseases in Brazil, with particular emphasis on the state of Rio de Janeiro.

Low genetic diversity in Wolbachia-Infected Culex quinquefasciatus (Diptera: Culicidae) from Brazil and Argentina

Culex quinquefasciatus is a vector of human pathogens, including filarial nematodes and several viruses. Although its epidemiological relevance is known to vary across geographical regions, an understanding of its population genetic structure is still incipient. In light of this, we evaluated the genetic diversity of Cx. quinquefasciatus and Cx. pipiens x Cx. quinquefasciatus hybrids collected from nine localities in Brazil and one site in Argentina. We used mitochondrial genes cox1 and nd4, along with the coxA and wsp genes of the maternally-inherited Wolbachia endosymbiont. The nd4 fragment was invariant between samples, whilst cox1 exhibited four haplotypes that separated two types of Cx. quinquefasciatus, one clustered in southern Brazil. Low sequence diversity was generally observed, being discussed. Both Brazilian and Argentinian mosquitoes were infected with a single Wolbachia strain. As reported in previous studies with these populations, cox1 and nd4 diversity is not congruent with the population structure revealed by nuclear markers or alar morphology. Future Cx. quinquefasciatus research should, if possible, evaluate mtDNA diversity in light of other markers.

Genetic structure and connectivity between populations of two common Mediterranean sessile invertebrates

Population genetic and phylogeography of two common mediterranean species were studied in 10 localities located on the coasts of Toscana, Puglia and Calabria. The aim of the study was to verify the extent of genetic breaks, in areas recognized as boundaries between Mediterranean biogeographic sectors. From about 100 sequences obtained from the mitochondrial Cytochrome Oxidase subunit I (COI) gene of Halocynthia papillosa and Hexaplex trunculus genetic diversity, genetic structure at small and large distances and demographic history of both specieswere analyzed. No evidences of genetic breaks were found for the two species in Toscana and Puglia. The genetic structure of H. trunculus evidences the extent of a barrier to gene flow localized in Calabria, which could be represented by the Siculo-Tunisian Strait and the Strait of Messina. The observed patterns showed similar level of gene flow at small distances in both species, although the two species have different larval ecology. These results suggest that other factors, such as currents, local dynamics and seasonal temperatures, influence the connectivity along the Italian peninsula. The geographic distribution of the haplotypes shows that H. papillosacould represent a single genetic pool in expansion, whereas H. trunculus has two distinct genetic pools in expansion. The demographic pattern of the two species suggests that Pleistocene sea level oscillations, in particular of the LGM, may have played a key role in shaping genetic structure of the two species. This knowledge provides basic information, useful for the definition of management plans, or for the design of a network of marine protected areas along the Italian peninsula.

Study of variability and genetic structure of European populations of Myotis emarginatus and Myotis capaccinii (Chiroptera, Vespertilionidae)

The Geoffroy’s bat Myotis emarginatus is mainly present in southern, south-eastern and central Europe (Červerný, 1999) and is often recorded from northern Spain (Quetglas, 2002; Flaquer et al., 2004). It has demonstrated the species’ preference for forest. Myotis capaccinii, confined to the Mediterranean (Guille´n, 1999), is classified as ‘vulnerable’ on a global scale (Hutson, Mickleburgh & Racey, 2001). In general, the species preferred calm waters bordered by well-developed riparian vegetation and large (> 5 m) inter-bank distances (Biscardi et al. 2007). In this study we present the first results about population genetic structure of these two species of genus Myotis. We used two methods of sampling: invasive and non-invasive techniques. A total of 323 invasive samples and a total of 107 non-invasive samples were collected and analyzed. For Myotis emarginatus we have individuated for the first time a set of 7 microsatellites, which can work on this species, started from a set developed on Myotis myotis (Castella et al. 2000). We developed also a method for analysis of non-invasive samples, that given a good percentage of positive analyzed samples. The results have highlighted for the species Myotis emarginatus the presence on the European territory of two big groups, discovered by using the microsatellites tracers. On this species, 33 haplotypes of Dloop have been identified, some of them are presented only in some colonies. We identified respectively 33 haplotypes of Dloop and 10 of cytB for Myotis emarginatus and 25 of dloop and 15 of cytB for Myotis capaccinii. Myotis emarginatus’ results, both microsatellites and mtDNA, show that there is a strong genetic flow between different colonies across Europe. The results achieved on Myotis capaccinii are very interesting, in this case either for the microsatellites or the mitochondrial DNA sequences, and it has been highlighted a big difference between different colonies.

Novel tools for conservation genetics in marine fish: population structure and evolution of the European hake (Merluccius merluccius) inferred by SNP variation and applications to traceability

The research presented in my PhD thesis is part of a wider European project, FishPopTrace, focused on traceability of fish populations and products. My work was aimed at developing and analyzing novel genetic tools for a widely distributed marine fish species, the European hake (Merluccius merluccius), in order to investigate population genetic structure and explore potential applications to traceability scenarios. A total of 395 SNPs (Single Nucleotide Polymorphisms) were discovered from a massive collection of Expressed Sequence Tags, obtained by high-throughput sequencing, and validated on 19 geographic samples from Atlantic and Mediterranean. Genome-scan approaches were applied to identify polymorphisms on genes potentially under divergent selection (outlier SNPs), showing higher genetic differentiation among populations respect to the average observed across loci. Comparative analysis on population structure were carried out on putative neutral and outlier loci at wide (Atlantic and Mediterranean samples) and regional (samples within each basin) spatial scales, to disentangle the effects of demographic and adaptive evolutionary forces on European hake populations genetic structure. Results demonstrated the potential of outlier loci to unveil fine scale genetic structure, possibly identifying locally adapted populations, despite the weak signal showed from putative neutral SNPs. The application of outlier SNPs within the framework of fishery resources management was also explored. A minimum panel of SNP markers showing maximum discriminatory power was selected and applied to a traceability scenario aiming at identifying the basin (and hence the stock) of origin, Atlantic or Mediterranean, of individual fish. This case study illustrates how molecular analytical technologies have operational potential in real-world contexts, and more specifically, potential to support fisheries control and enforcement and fish and fish product traceability.

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.

Rapid parallel adaptive radiations from a single hybridogenic ancestral population

The Alpine lake whitefish (Coregonus lavaretus) species complex is a classic example of a recent radiation, associated with colonization of the Alpine lakes following the glacial retreat (less than 15 kyr BP). They have formed a unique array of endemic lake flocks, each with one to six described sympatric species differing in morphology, diet and reproductive ecology. Here, we present a genomic investigation of the relationships between and within the lake flocks. Comparing the signal between over 1000 AFLP loci and mitochondrial control region sequence data, we use phylogenetic tree-based and population genetic methods to reconstruct the phylogenetic history of the group and to delineate the principal centres of genetic diversity within the radiation. We find significant cytonuclear discordance showing that the genomically monophyletic Alpine whitefish clade arose from a hybrid swarm of at least two glacial refugial lineages. Within this radiation, we find seven extant genetic clusters centred on seven lake systems. Most interestingly, we find evidence of sympatric speciation within and parallel evolution of equivalent phenotypes among these lake systems. However, we also find the genetic signature of human-mediated gene flow and diversity loss within many lakes, highlighting the fragility of recent radiations.

Space trajectories optimization using variable-chromosome-length genetic algorithms

The problem of optimal design of a multi-gravity-assist space trajectories, with free number of deep space maneuvers (MGADSM) poses multi-modal cost functions. In the general form of the problem, the number of design variables is solution dependent. To handle global optimization problems where the number of design variables varies from one solution to another, two novel genetic-based techniques are introduced: hidden genes genetic algorithm (HGGA) and dynamic-size multiple population genetic algorithm (DSMPGA). In HGGA, a fixed length for the design variables is assigned for all solutions. Independent variables of each solution are divided into effective and ineffective (hidden) genes. Hidden genes are excluded in cost function evaluations. Full-length solutions undergo standard genetic operations. In DSMPGA, sub-populations of fixed size design spaces are randomly initialized. Standard genetic operations are carried out for a stage of generations. A new population is then created by reproduction from all members based on their relative fitness. The resulting sub-populations have different sizes from their initial sizes. The process repeats, leading to increasing the size of sub-populations of more fit solutions. Both techniques are applied to several MGADSM problems. They have the capability to determine the number of swing-bys, the planets to swing by, launch and arrival dates, and the number of deep space maneuvers as well as their locations, magnitudes, and directions in an optimal sense. The results show that solutions obtained using the developed tools match known solutions for complex case studies. The HGGA is also used to obtain the asteroids sequence and the mission structure in the global trajectory optimization competition (GTOC) problem. As an application of GA optimization to Earth orbits, the problem of visiting a set of ground sites within a constrained time frame is solved. The J2 perturbation and zonal coverage are considered to design repeated Sun-synchronous orbits. Finally, a new set of orbits, the repeated shadow track orbits (RSTO), is introduced. The orbit parameters are optimized such that the shadow of a spacecraft on the Earth visits the same locations periodically every desired number of days.

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.