Temporal variation of genetic diversity and differentiation of albacore, Thunnus alalunga, in the Mediterranean Sea

This study is on albacore (Thunnus alalunga, Bonnaterre 1788), an epi- and mesopelagic oceanic tuna species cosmopolitan in the tropical and temperate waters of all oceans including the Mediterranean Sea, extending in a broad band between 40°N and 40°S. What it’s known about albacore population structure is based on different studies that used fisheries data, RFLP, mtDNA control region and nuDNA markers, blood lectins analysis, individual tags and microsatellite. At the moment, for T. alalunga six management units are recognized: the North Pacific, South Pacific, Indian, North Atlantic, South Atlantic and Mediterranean stocks. In this study I have done a temporal and spatial comparison of genetic variability between different Mediterranean populations of Thunnus alalunga matching an historical dataset ca. from 1920s composed of 43 individuals divided in 3 populations (NADR, SPAIN and CMED) with a modern dataset composed of 254 individuals and 7 populations (BAL, CYP, LIG, TYR, TUR, ADR, ALB). The investigation was possible using a panel of 94 nuclear SNPs, built specifically for the target species at the University of Basque Country UPV/EHU. First analysis done was the Hardy-Weinberg, then the number of clusters (K) was determined using STRUCTURE and to assess the genetic variability, allele frequencies, the average number of alleles per locus, expected (He) and observed (Ho) heterozygosis, and the index of polymorphism (P) was used the software Genetix. Historical and modern samples gives different results, showing a clear loss of genetic diversity over time leading to a single cluster in modern albacore instead of the two found in historical samples. What this study reveals is very important for conservation concerns, and additional research endeavours are needed.

Are there genetic breaks between Atlantic and Pacific Yellowfin tuna (Thunnus albacares) populations? A preliminary study based on microsatellites gene variation

Yellowfin tuna (Thunnus albacares, YFT, Bonnaterre 1788) is one of the most important market tuna species in the world. The high mortality of juveniles is in part caused by their bycatch. Indeed, if unregulated, it could permanently destabilize stocks health. For this reason investigating and better knowing the stock boundaries represent a crucial concern. Aim of this thesis was to preliminary investigate the YFT population structure within and between Atlantic and Pacific Oceans through the analysis of genetic variation at eight microsatellite loci and assess the occurrence of barriers to the gene flow between Oceans. For this propouse we collected 4 geographical samples coming from Atlantic and Pacific Ocean and selected a panel of 8 microsatellites loci developped by Antoni et al., (2014). Samples 71-2-Y and 77-2-Y, came from rispectively west central pacific ocean (WCPO) and east central pacific ocean (ECPO), instead samples 41-1-Y and 34-2-Y derive from west central atlantic ocean (WCAO) and east central atlantic ocean (ECAO). Total 160 specimens were analyzed (40 per sample) and were carried out several genetic information as allele frequencies, allele number, allelic richness, HWE (using He and Ho) and pairwise Fst genetic distance. Results obtained, may support the panmictic theory of this species, only one of pairwise Fst obtained is statistically significant (Fst= 0.00927; pV= 0.00218) between 41-1-Y and 71-2-Y samples. Results suggest low genetic differentiation and consequent high level of gene flow between Atlantic and Pacific populations. Furthermore, we performed an analysis of molecular taxonomy through the use of ATCO (the flaking region between ATPse6 and cytochrome oxidase subunit III genes mt DNA, to discriminate within the gener Thunnus two of the related species (Yellofin and bigeye tuna) according with their difficult recognition at certain size (<40 cm). ATCO analysis in this thesis, has provided strong discriminate evidence between the target species proving to be one of the most reliable genetic tools capable to indagate within the genus Thunnus. Thus, our study has provided useful information for possible use of this protocol for conservation plans and management of this fish stocks.

Inference of Raja miraletus population structure using mitochondrial and nuclear DNA: comparing the resolution power of molecular markers in exploring species and population boundaries

This study focused on the role of oceanographic discontinuities and the presence of transitional areas in shaping the population structure and the phylogeography of the Raja miraletus species complex, coupled with the test of the effective occurrence of past speciation events. The comparisons between the Atlantic African and the North-Eastern Atlantic-Mediterranean geographic populations were unravelled using both Cytochrome Oxidase I and eight microsatellite loci. This approach guaranteed a robust dataset for the identification of a speciation event between the Atlantic African clade, corresponding to the ex Raja ocellifera nominal species, and the NE Atlantic-Mediterranean R. miraletus clade. As a matter of fact, the origin of the Atlantic Africa and the NE Atlantic-Mediterranean deep split dated about 11.74MYA and was likely due to the synergic influence currents and two upwelling areas crossing the Western African Waters. Within the Mediterranean Sea, particular attention was also paid to the transitional area represented by Adventura and Maltese Bank, that might have contributed in sustaining the connectivity of the Western and the Eastern Mediterranean geographical populations. Furthermore, the geology of the easternmost part of Sicily and the geo-morphological depression of the Calabrian Arc could have driven the differentiation of the Eastern Mediterranean Sea. Although bathymetric and oceanographic discontinuity could represent barriers to dispersal and migration between Eastern and Western Mediterranean samples, a clear and complete genetic separation among them was not detected. Results produced by this work identified a speciation event defining Raja ocellifera and R. miraletus as two different species, and describing the R. miraletus species complex as the most ancient cryptic speciation event in the family Rajidae, representing another example of how strictly connected the environment, the behavioural habits and the evolutionary and ecologic drivers are.

Evidence of neutral and adaptive genetic divergence between European trout populations sampled along altitudinal gradients

Species with a wide geographical distribution are often composed of distinct subgroups which may be adapted to their local environment. European trout (Salmo trutta species complex) provide an example of such a complex consisting of several genetically and ecologically distinct forms. However, trout populations are strongly influenced by human activities, and it is unclear to what extent neutral and adaptive genetic differences have persisted. We sampled 30 Swiss trout populations from heterogeneous environments along replicated altitudinal gradients in three major European drainages. More than 850 individuals were genotyped at 18 microsatellite loci which included loci diagnostic for evolutionary lineages and candidate markers associated with temperature tolerance, reproductive timing and immune defence. We find that the phylogeographic structure of Swiss trout populations has not been completely erased by stocking. Distinct genetic clusters corresponding to the different drainages could be identified, although nonindigenous alleles were clearly present, especially in the two Mediterranean drainages. We also still detected neutral genetic differentiation within rivers which was often associated with the geographical distance between populations. Five loci showed evidence of divergent selection between populations with several drainage-specific patterns. Lineage-diagnostic markers, a marker linked to a quantitative trait locus for upper temperature tolerance in other salmonids and a marker linked to the major histocompatibility class I gene were implicated in local adaptation and some patterns were associated with altitude. In contrast, tentative evidence suggests a signal of balancing selection at a second immune relevant gene (TAP2). Our results confirm the persistence of both neutral and potentially adaptive genetic differences between trout populations in the face of massive human-mediated dispersal.

Mapping genetic variants associated with beta-adrenergic responses in inbred mice

β-blockers and β-agonists are primarily used to treat cardiovascular diseases. Inter-individual variability in response to both drug classes is well recognized, yet the identity and relative contribution of the genetic players involved are poorly understood. This work is the first genome-wide association study (GWAS) addressing the values and susceptibility of cardiovascular-related traits to a selective β(1)-blocker, Atenolol (ate), and a β-agonist, Isoproterenol (iso). The phenotypic dataset consisted of 27 highly heritable traits, each measured across 22 inbred mouse strains and four pharmacological conditions. The genotypic panel comprised 79922 informative SNPs of the mouse HapMap resource. Associations were mapped by Efficient Mixed Model Association (EMMA), a method that corrects for the population structure and genetic relatedness of the various strains. A total of 205 separate genome-wide scans were analyzed. The most significant hits include three candidate loci related to cardiac and body weight, three loci for electrocardiographic (ECG) values, two loci for the susceptibility of atrial weight index to iso, four loci for the susceptibility of systolic blood pressure (SBP) to perturbations of the β-adrenergic system, and one locus for the responsiveness of QTc (p<10(-8)). An additional 60 loci were suggestive for one or the other of the 27 traits, while 46 others were suggestive for one or the other drug effects (p<10(-6)). Most hits tagged unexpected regions, yet at least two loci for the susceptibility of SBP to β-adrenergic drugs pointed at members of the hypothalamic-pituitary-thyroid axis. Loci for cardiac-related traits were preferentially enriched in genes expressed in the heart, while 23% of the testable loci were replicated with datasets of the Mouse Phenome Database (MPD). Altogether these data and validation tests indicate that the mapped loci are relevant to the traits and responses studied.

Analysis of pedigree and conformation data to explain genetic variability of the horse breed Franches-Montagnes

Franches-Montagnes is the only native horse breed in Switzerland, therefore special efforts should be made for ensuring its survival. The objectives of this study were to characterize the structure of this population as well as genetic variability with pedigree data, conformation traits and molecular markers. Studies were focused to clarify if this population is composed of a heavy- and a light-type subpopulation. Extended pedigree records of 3-year-old stallions (n = 68) and mares (n = 108) were available. Evaluations of body conformation traits as well as pedigree data and molecular markers did not support the two-subpopulation hypothesis. The generation interval ranged from 7.8 to 9.3 years. The complete generation equivalent was high (>12). The number of effective ancestors varied between 18.9 and 20.1, whereof 50% of the genetic variability was attributed to seven of them. Genetic contribution of Warmblood horses ranged from 36% to 42% and that of Coldblood horses from 4% to 6%. The average inbreeding coefficient reached 6%. Inbreeding effective population size was 114.5 when the average increase of the inbreeding coefficient per year since 1910 was taken. Our results suggest that bottleneck situations occurred because of selection of a small number of sire lines. Promotion of planned matings between parents that are less related is recommended in order to avoid a reduction of the genetic diversity.

Genetic diversity among horse populations with a special focus on the Franches-Montagnes breed

Genetic characterization helps to assure breed integrity and to assign individuals to defined populations. The objective of this study was to characterize genetic diversity in six horse breeds and to analyse the population structure of the Franches-Montagnes breed, especially with regard to the degree of introgression with Warmblood. A total of 402 alleles from 50 microsatellite loci were used. The average number of alleles per locus was significantly lower in Thoroughbreds and Arabians. Average heterozygosities between breeds ranged from 0.61 to 0.72. The overall average of the coefficient of gene differentiation because of breed differences was 0.100, with a range of 0.036-0.263. No significant correlation was found between this parameter and the number of alleles per locus. An increase in the number of homozygous loci with increasing inbreeding could not be shown for the Franches-Montagnes horses. The proportion of shared alleles, combined with the neighbour-joining method, defined clusters for Icelandic Horse, Comtois, Arabians and Franches-Montagnes. A more disparate clustering could be seen for European Warmbloods and Thoroughbreds, presumably from frequent grading-up of Warmbloods with Thoroughbreds. Grading-up effects were also observed when Bayesian and Monte Carlo resampling approaches were used for individual assignment to a given population. Individual breed assignments to defined reference populations will be very difficult when introgression has occurred. The Bayesian approach within the Franches-Montagnes breed differentiated individuals with varied proportions of Warmblood.

Genetic association studies under the population stratification, family pedigree and application to genome-wide association studies

This dissertation has three separate parts: the first part deals with the general pedigree association testing incorporating continuous covariates; the second part deals with the association tests under population stratification using the conditional likelihood tests; the third part deals with the genome-wide association studies based on the real rheumatoid arthritis (RA) disease data sets from Genetic Analysis Workshop 16 (GAW16) problem 1. Many statistical tests are developed to test the linkage and association using either case-control status or phenotype covariates for family data structure, separately. Those univariate analyses might not use all the information coming from the family members in practical studies. On the other hand, the human complex disease do not have a clear inheritance pattern, there might exist the gene interactions or act independently. In part I, the new proposed approach MPDT is focused on how to use both the case control information as well as the phenotype covariates. This approach can be applied to detect multiple marker effects. Based on the two existing popular statistics in family studies for case-control and quantitative traits respectively, the new approach could be used in the simple family structure data set as well as general pedigree structure. The combined statistics are calculated using the two statistics; A permutation procedure is applied for assessing the p-value with adjustment from the Bonferroni for the multiple markers. We use simulation studies to evaluate the type I error rates and the powers of the proposed approach. Our results show that the combined test using both case-control information and phenotype covariates not only has the correct type I error rates but also is more powerful than the other existing methods. For multiple marker interactions, our proposed method is also very powerful. Selective genotyping is an economical strategy in detecting and mapping quantitative trait loci in the genetic dissection of complex disease. When the samples arise from different ethnic groups or an admixture population, all the existing selective genotyping methods may result in spurious association due to different ancestry distributions. The problem can be more serious when the sample size is large, a general requirement to obtain sufficient power to detect modest genetic effects for most complex traits. In part II, I describe a useful strategy in selective genotyping while population stratification is present. Our procedure used a principal component based approach to eliminate any effect of population stratification. The paper evaluates the performance of our procedure using both simulated data from an early study data sets and also the HapMap data sets in a variety of population admixture models generated from empirical data. There are one binary trait and two continuous traits in the rheumatoid arthritis dataset of Problem 1 in the Genetic Analysis Workshop 16 (GAW16): RA status, AntiCCP and IgM. To allow multiple traits, we suggest a set of SNP-level F statistics by the concept of multiple-correlation to measure the genetic association between multiple trait values and SNP-specific genotypic scores and obtain their null distributions. Hereby, we perform 6 genome-wide association analyses using the novel one- and two-stage approaches which are based on single, double and triple traits. Incorporating all these 6 analyses, we successfully validate the SNPs which have been identified to be responsible for rheumatoid arthritis in the literature and detect more disease susceptibility SNPs for follow-up studies in the future. Except for chromosome 13 and 18, each of the others is found to harbour susceptible genetic regions for rheumatoid arthritis or related diseases, i.e., lupus erythematosus. This topic is discussed in part III.

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.

Genetic variation and demographic history of the Haplochromis laparogramma group of Lake Victoria-An analysis based on SINEs and mitochondrial DNA

More than 500 endemic haplochromine cichlid species inhabit Lake Victoria. This striking species diversity is a classical example of recent explosive adaptive radiation thought to have happened within the last similar to 15,000 years. In this study, we examined the population structure and historical demography of 3 pelagic haplochromine cichlid species that resemble in morphology and have similar niche, Haplochromis (Yssichromis) laparogramma, Haplochromis (Y.) pyrrhocephalus, and Haplochromis (Y.) sp. "glaucocephalus". We investigated the sequences of the mitochondrial DNA control region and the insertion patterns of short interspersed elements (SINEs) of 759 individuals. We show that sympatric forms are genetically differentiated in 4 of 6 cases, but we also found apparent weakening of the genetic differentiation in areas with turbid water. We estimated the timings of population expansion and species divergence to coincide with the refilling of the lake at the Pleistocene/Holocene boundary. We also found that estimates can be altered significantly by the choice of the shape of the molecular clock. If we employ the nonlinear clock model of evolutionary rates in which the rates are higher towards the recent, the population expansion was dated at around the event of desiccation of the lake ca. 17,000 YBP. Thus, we succeeded in clarifying the species and population structure of closely related Lake Victoria cichlids and in showing the importance of applying appropriate clock calibrations in elucidating recent evolutionary events. (C) 2009 Elsevier B.V. All rights reserved.

Genetic diversity and ecological success of Staphylococcus aureus strains colonizing humans

The genetic determinants and phenotypic traits which make a Staphylococcus aureus strain a successful colonizer are largely unknown. The genetic diversity and population structure of 133 S. aureus isolates from healthy, generally risk-free adult carriers were investigated using four different typing methods: multilocus sequence typing (MLST), amplified fragment length polymorphism analysis (AFLP), double-locus sequence typing (DLST), and spa typing were compared. Carriage isolates displayed great genetic diversity which could only be revealed fully by DLST. Results of AFLP and MLST were highly concordant in the delineation of genotypic clusters of closely related isolates, roughly equivalent to clonal complexes. spa typing and DLST provided considerably less phylogenetic information. The resolution of spa typing was similar to that of AFLP and inferior to that of DLST. AFLP proved to be the most universal method, combining a phylogeny-building capacity similar to that of MLST with a much higher resolution. However, it had a lower reproducibility than sequencing-based MLST, DLST, and spa typing. We found two cases of methicillin-resistant S. aureus colonization, both of which were most likely associated with employment at a health service. Of 21 genotypic clusters detected, 2 were most prevalent: cluster 45 and cluster 30 each colonized 24% of the carrier population. The number of bacteria found in nasal samples varied significantly among the clusters, but the most prevalent clusters were not particularly numerous in the nasal samples. We did not find much evidence that genotypic clusters were associated with different carrier characteristics, such as age, sex, medical conditions, or antibiotic use. This may provide empirical support for the idea that genetic clusters in bacteria are maintained in the absence of adaptation to different niches. Alternatively, carrier characteristics other than those evaluated here or factors other than human hosts may exert selective pressure maintaining genotypic clusters.

Genetic diversity of Swiss sheep breeds in the focus of conservation research

There is constant pressure to improve evaluation of animal genetic resources in order to prevent their erosion. Maintaining the integrity of livestock species as well as their genetic diversity is of paramount interest for long-term agricultural policies. One major use of DNA techniques in conservation is to reveal genetic diversity within and between populations. Forty-one microsatellites were analysed to assess genetic diversity in nine Swiss sheep breeds and to measure the loss of the overall diversity when one breed would become extinct. The expected heterozygosities varied from 0.65 to 0.74 and 10.8% of the total genetic diversity can be explained by the variation among breeds. Based on the proportion of shared alleles, each of the nine breeds were clearly defined in their own cluster in the neighbour-joining tree describing the relationships among the breeds. Bayesian clustering methods assign individuals to groups based on their genetic similarity and infer the number of populations. In STRUCTURE, this approach pooled the Valais Blacknose and the Valais Red. With BAPS method the two Valais sheep breeds could be separated. Caballero & Toro approach (2002) was used to calculate the loss or gain of genetic diversity when each of the breeds would be removed from the set. The changes in diversity based on between-breed variation ranged from -12.2% (Valais Blacknose) to 0% (Swiss Black Brown Mountain and Mirror Sheep); based on within-breed diversity the removal of a breed could also produce an increase in diversity (-0.6% to + 0.6%). Allelic richness ranged from 4.9 (Valais Red) to 6.7 (Brown Headed Meat sheep and Red Engadine Sheep). Breed conservation decisions cannot be limited to genetic diversity alone. In Switzerland, conservation goals are embedded in the desire to carry the cultural legacy over to future generations.

A genome-to-genome analysis of associations between human genetic variation, HIV-1 sequence diversity, and viral control

HIV-1 sequence diversity is affected by selection pressures arising from host genomic factors. Using paired human and viral data from 1071 individuals, we ran >3000 genome-wide scans, testing for associations between host DNA polymorphisms, HIV-1 sequence variation and plasma viral load (VL), while considering human and viral population structure. We observed significant human SNP associations to a total of 48 HIV-1 amino acid variants (p<2.4 × 10−12). All associated SNPs mapped to the HLA class I region. Clinical relevance of host and pathogen variation was assessed using VL results. We identified two critical advantages to the use of viral variation for identifying host factors: (1) association signals are much stronger for HIV-1 sequence variants than VL, reflecting the ‘intermediate phenotype’ nature of viral variation; (2) association testing can be run without any clinical data. The proposed genome-to-genome approach highlights sites of genomic conflict and is a strategy generally applicable to studies of host–pathogen interaction.

Population structure, inbreeding and local adaptation within an endangered riverine specialist: the nase (Chondrostoma nasus)

Chondrostoma nasus is a cyprinid fish with highly specialized, ecologically and geographically distinct, ontogenetic trophic niches. Nase population numbers across their Swiss range have shown massive declines and many localized extinctions. Here we integrate data from different genetic markers with phenotypic and demographic data to survey patterns of neutral and adaptive genetic diversity in all extant (and one extinct) Swiss nase populations, with the aim to delineate intraspecific conservation units (CUs) and to inform future population management strategies. We discovered two major genetically and geographically distinct population groupings. The first population grouping comprises nase inhabiting rivers flowing into Lake Constance; the second comprises nase populations from Rhine drainages below Lake Constance. Within these clusters there is generally limited genetic differentiation among populations. Genomic outlier scans based on 256–377 polymorphic AFLP loci revealed little evidence of local adaptation both within and among population clusters, with the exception of one candidate locus identified in scans involving the inbred Schanzengraben population. However, significant phenotypic differentiation in body shape between certain populations suggests a need for more intensive future studies of local adaptation. Our data strongly suggests that the two major population groups should be treated as distinct CUs, with any supplemental stocking and reintroductions sourced only from within the range of the CU concerned.

Genetic and physical analysis of the bacteriophage P22 tail protein

A plasmid based genetic system was developed for the tail protein of the Salmonella typhimurium bacteriophage P22 and used to isolate and characterize tail protein mutants. The tail protein is a trimeric structural protein of the phage and an endorhamnosidase whose activity is essential for infection. The gene for the tail protein has previously been cloned into a plasmid expression vector and sequenced. A plate complementation assay for tail protein produced from the cloned gene was developed and used to isolate 27 tail protein mutants following mutagenesis of the cloned gene. These mutations were mapped into 12 deletion intervals using deletions which were made on plasmids in vitro and crossed onto P22. The base substitutions were determined by DNA sequencing. The majority of mutants had missense or nonsense mutations in the protein coding portion of the gene; however four of the mutants were in the putative transcription terminator. The oligomeric state of tail protein from the 15 missense mutants was investigated using SDS and nondenaturing polyacrylamide gel electrophoresis of cell lysates. Wild-type tail protein retains its trimeric structure in SDS gels at room temperature. Two of the mutant proteins also migrated as trimers in SDS gels, yet one of these had a considerably faster mobility than wild-type trimer. Its migration was the same as wild-type in a nondenaturing gel, so it is thought to be a trimer which is partially denatured by SDS. Four of the mutants produced proteins which migrate at the position of a monomer in an SDS gel but cannot be seen on a nondenaturing gel. These proteins are thought to be either monomers or soluble aggregates which cannot enter the nondenaturing gel. The remainder of mutants produce protein which is degraded. The mutant tail protein which had normal trimeric mobility on SDS and nondenaturing gels was purified. This protein has essentially wild-type ability to attach to phage capsids, but its endorhamnosidase activity is only 4% of wild-type. ^