Population structure, habitat use and conservation of short-finned pilot whales Globicephala macrorhynchus in the Archipelago of Madeira

This thesis provides information on the grouping structure, survival, abundance, dive characteristics and habitat preferences of short-finned pilot whales occurring in the oceanic archipelago of Madeira (Portugal, NE Atlantic), based on data collected between 2001-2011, and contributes for its conservation. Photo-identification methods and genetic analyses demonstrated that there is a large degree of variability in site fidelity, including resident, regular visitor and transient whales, and that they may not be genetically isolated. It is proposed that the pilot whales encountered in Madeira belong to a single population encompassing several clans, possibly three clans of island-associated (i.e. resident and regular visitor) whales and others of transients, each containing two to three matrilineal pods. Mark-recapture methods estimated that the island-associated community is composed of less than 150 individuals and that their survival rate is within the range of other long-lived cetacean species, and that around 300 whales of different residency patterns uses the southern area of the island of Madeira from mid-summer to mid-autumn. No significant trend was observed between years. Time-depth recorders deployed in adult whales during daytime revealed that they spend over ¾ of their time at the surface, that they have a low diving rate, and that transient whales also forage during their passage. The analyses of visual data collected from nautical and aerial line-transect surveys indicate a core/preferred habitat area in the south-east of the island of Madeira. That area is used for resting, socializing, foraging, breeding, calving and birthing. Thus, that area should be considered as an important habitat for this species, at least seasonally (during autumn) when the species is more abundant, and included in conservation plans. No direct threat needing urgent measures was identified, although the impact of some activities like whale-watching or marine traffic should be assessed.

Population structure and effects of inbreeding on milk yield and quality of Murrah buffaloes

To provide data for conservation, selection, and expansion programs of buffalo herds, this study evaluated the history of a population of Murrah buffaloes based on population structure and the effect of inbreeding on accumulated 305-d milk yield (MY), fat yield (FY), protein yield (PY), mozzarella production (MProd), and somatic cell score (SCS). The usefulness of including the individual inbreeding coefficient (F) or individual increase in inbreeding coefficient (Delta F) in the model to describe inbreeding depression was evaluated. Pedigree information from 8,054 animals born between 1976 and 2008 and 4,497 lactation records obtained from 12 herds were used. The realized effective population size was 40.10 +/- 1.27, and the mean F of the entire population was 2.14%. The ratio between the number of founders and ancestors demonstrated the existence of a bottleneck in the pedigree of this population, which may contribute to a reduction of genetic diversity. The effect of F on MY, FY, PY, MProd, and SCS was -1.005 kg, -0.299 kg, -0.246 kg, -1.201 kg, and -0.002 units, and the effect of Delta F transformed to equivalent F (%) for a mean of 2.57 equivalent generations was -4.287 kg, -0.581 kg, -0.383 kg, -2.001 kg, and -0.007 units, respectively. The inbreeding depression observed may have important economic repercussions for production systems. The Delta F can be considered the better of the two indicators of inbreeding depression due to its properties that prevent underestimation of this effect. A designed mating system to avoid inbreeding may be applied to this population to maintain genetic diversity.

Influence of data structure on the estimation of the additive genetic direct and maternal covariance for early growth traits in Nellore cattle

The objective of the present study was to investigate the effect of data structure on estimated genetic parameters and predicted breeding values of direct and maternal genetic effects for weaning weight (WW) and weight gain from birth to weaning (BWG), including or not the genetic covariance between direct and maternal effects. Records of 97,490 Nellore animals born between 1993 and 2006, from the Jacarezinho cattle raising farm, were used. Two different data sets were analyzed: DI_all, which included all available progenies of dams without their own performance; DII_all, which included DI_all + 20% of recorded progenies with maternal phenotypes. Two subsets were obtained from each data set (DI_all and DII_all): DI_1 and DII_1, which included only dams with three or fewer progenies; DI_5 and DII_5, which included only dams with five or more progenies. (Co)variance components and heritabilities were estimated by Bayesian inference through Gibbs sampling using univariate animal models. In general, for the population and traits studied, the proportion of dams with known phenotypic information and the number of progenies per dam influenced direct and maternal heritabilities, as well as the contribution of maternal permanent environmental variance to phenotypic variance. Only small differences were observed in the genetic and environmental parameters when the genetic covariance between direct and maternal effects was set to zero in the data sets studied. Thus, the inclusion or not of the genetic covariance between direct and maternal effects had little effect on the ranking of animals according to their breeding values for WW and BWG. Accurate estimation of genetic correlations between direct and maternal genetic effects depends on the data structure. Thus, this covariance should be set to zero in Nellore data sets in which the proportion of dams with phenotypic information is low, the number of progenies per dam is small, and pedigree relationships are poorly known. (c) 2012 Elsevier B.V. All rights reserved.

Estimates of genetic parameters, and cluster and principal components analyses of breeding values related to egg production traits in a White Leghorn population

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Atlantic forest fragmentation and genetic diversity of an isolated population of the Blue-manakin, Chiroxiphia caudata (Pipridae), assessed by microsatellite analyses

Habitat fragmentation is predicted to restrict gene flow, which can result in the loss of genetic variation and inbreeding depression. The Brazilian Atlantic forest has experienced extensive loss of habitats since European settlement five centuries ago, and many bird populations and species are vanishing. Genetic variability analysis in fragmented populations could be important in determining their long-term viability and for guiding management plans. Here we analyzed genetic diversity of a small understory bird, the Blue-manakins Chiroxiphia caudata (Pipridae), from an Atlantic forest fragment (112 ha) isolated 73 years ago, and from a 10,000 ha continuous forest tract (control), using orthologous microsatellite loci. Three of the nine loci tested were polymorphic. No statistically significant heterozygote loss was detected for the fragment population. Although genetic diversity, which was estimated by expected heterozygosity and allelic richness, has been lower in the fragment population in relation to the control, it was not statistically significant, suggesting that this 112 ha fragment can be sufficient to maintain a blue-manakin population large enough to avoid stochastic effects, such as inbreeding and/or genetic drift. Alternatively, it is possible that 73 years of isolation did not accumulate sufficient generations for these effects to be detected. However, some alleles have been likely lost, specially the rare ones, what is expected from genetic drift for such a small and isolated population. A high genetic differentiation was detected between populations by comparing both allelic and genotypic distributions. Only future studies in continuous areas are likely to answer if such a structure was caused by the isolation resulted from the forest fragmentation or by natural population structure.

Evaluating genetic algorithms with different population structures on a lot sizing and scheduling problem

This paper studies the use of different population structures in a Genetic Algorithm (GA) applied to lot sizing and scheduling problems. The population approaches are divided into two types: single-population and multi-population. The first type has a non-structured single population. The multi-population type presents non-structured and structured populations organized in binary and ternary trees. Each population approach is tested on lot sizing and scheduling problems found in soft drink companies. These problems have two interdependent levels with decisions concerning raw material storage and soft drink bottling. The challenge is to simultaneously determine the lot sizing and scheduling of raw materials in tanks and products in lines. Computational results are reported allowing determining the better population structure for the set of problem instances evaluated. Copyright 2008 ACM.

Population structure of Lutjanus purpureus (Lutjanidae - Perciformes) on the Brazilian coast: further existence evidence of a single species of red snapper in the western Atlantic

ABSTRACT: The present study focus on the mitochondrial control region to investigate phylogeographic patterns and population structure in Lutjanus purpureus, and to evaluate the genetic similarity between L. purpureus and L. campechanus. For the initial analysis, 810 base pairs sequence from control region were obtained from 239 specimens of L. purpureus collected from four localities off the Brazilian coast. The results revealed the presence of a single panmictic population characterized by high values of genetic diversity. The 299 base pairs hypervariable portion were used for the combined analysis of L. purpureus and L. campechanus, being 275 haplotypes identified in the 414 specimens. Phylogenetic tree and haplotype network did not indicate phylogeographic substructuring between the two species, but rather an intense intermingling of individuals. Considering their marked morphological similarity, the molecular data presented here indicatethat only one species of red snapper exists in the western Atlantic.

Global optimization using a genetic algorithm with hierarchically structured population

This paper applies a genetic algorithm with hierarchically structured population to solve unconstrained optimization problems. The population has individuals distributed in several overlapping clusters, each one with a leader and a variable number of support individuals. The hierarchy establishes that leaders must be fitter than its supporters with the topological organization of the clusters following a tree. Computational tests evaluate different population structures, population sizes and crossover operators for better algorithm performance. A set of known benchmark test problems is solved and the results found are compared with those obtained from other methods described in the literature, namely, two genetic algorithms, a simulated annealing, a differential evolution and a particle swarm optimization. The results indicate that the method employed is capable of achieving better performance than the previous approaches in regard as the two criteria usually employed for comparisons: the number of function evaluations and rate of success. The method also has a superior performance if the number of problems solved is taken into account. (C) 2013 Elsevier B.V. All rights reserved.

History, structure, and genetic diversity of Brazilian Gir cattle

The Gir cattle breed (Bos indicus) is an important genetic resource for milk production throughout the tropics. The small number of Gir animals introduced in Brazil, rapid dissemination of this breed recently, and the intensification of selection practices could contribute to increase of inbreeding level and reduce genetic diversity in this population. The population was analyzed in terms of pedigree completeness level, inbreeding coefficient, coancestry, generation interval, effective population size, effective number of founders and ancestors, among others. Despite the low mean inbreeding (around 2%), minor problems were identified in the population structure of the Brazilian Gir cattle, e.g., trend of narrower bottlenecks in the pedigree in recent years. The effective population sizes based on inbreeding (94) or coancestry (165.9) as well as the effective number of ancestors (76) and founders (143) were relativity high. The major subdivision of the Gir breed was observed between 1993 and 2002 (dairy and dual-purpose herds, wide use of within-herd matings). In this period the level of inbreeding remained at a higher level, there was a small increase in coancestry and the number of equivalent subpopulations was approximately 6. After 2002, there was genetic exchange between subpopulations, reduction in the average inbreeding, pronounced increase in the average coancestry, and the number of equivalent subpopulations was about 2. Furthermore, it was found that the mean generation interval of the population tended to increase in recent years (around 9 years). About 23% of genetic diversity has been lost since the first generation of founders. Based on the effective population size, number of equivalent subpopulations, inbreeding, coancestry, and loss of genetic diversity, the Gir population is still highly structured, but there is ample room for artificial selection. The results regarding the effective number of founders and ancestors in the present population demonstrate the existence of bottlenecks in the pedigree and indicate the need for population structure monitoring. Nevertheless, the Brazilian Gir breed can perfectly face a breeding program with high selection intensity. (C) 2014 Elsevier B.V. All rights reserved.

Mycobacterium tuberculosis population structure shift in a 5-year molecular epidemiology surveillance follow-up study in a low endemic agro-industrial setting in São Paulo, Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Conservation genetics of European wildcat (Felis silvestris silvestris): a wide and integrating analysis protocol for admixture inferences and population structure

Introgression of domestic cat genes into European wildcat (Felis silvestris silvestris) populations and reduction of wildcats’ range in Europe, leaded by habitat loss and fragmentation, are considered two of the main conservation problems for this endangered feline. This thesis addressed the questions related with the artificial hybridization and populations’ fragmentation, using a conservation genetics perspective. We combined the use of highly polymorphic loci, Bayesian statistical inferences and landscape analyses tools to investigate the origin of the geographic-genetic substructure of European wildcats (Felis silvestris silvestris) in Italy and Europe. The genetic variability of microsatellites evidenced that European wildcat populations currently distributed in Italy differentiated in, and expanded from two distinct glacial refuges during the Last Glacial Maximum. The genetic and geographic substructure detected between the eastern and western sides of the Apennine ridge, resulted by adaptation to specific ecological conditions of the Mediterranean habitats. European wildcat populations in Europe are strongly structured into 5 geographic-genetic macro clusters corresponding to: the Italian peninsular & Sicily; Balkans & north-eastern Italy; Germany eastern; central Europe; and Iberian Peninsula. Central European population might have differentiated in the extra-Mediterranean Würm ice age refuge areas (Northern Alps, Carpathians, and the Bulgarian mountain systems), while the divergence among and within the southern European populations might have resulted by the Pleistocene bio geographical framework of Europe, with three southern refugia localized in the Balkans, Italian Peninsula and Iberia Peninsula. We further combined the use of most informative autosomal SNPs with uniparental markers (mtDNA and Y-linked) for accurately detecting parental genotypes and levels of introgressive hybridization between European wild and domestic cats. A total of 11 hybrids were identified. The presence of domestic mitochondrial haplotypes shared with some wild individuals led us to hypnotize the possibility that ancient introgressive events might have occurred and that further investigation should be recommended.

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.

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.

DOUBLE-STRAND BREAK REPAIR PATHWAYS IN DNA STRUCTURE-INDUCED GENETIC INSTABILITY

Genetic instability in mammalian cells can occur by many different mechanisms. In the absence of exogenous sources of DNA damage, the DNA structure itself has been implicated in genetic instability. When the canonical B-DNA helix is naturally altered to form a non-canonical DNA structure such as a Z-DNA or H-DNA, this can lead to genetic instability in the form of DNA double-strand breaks (DSBs) (1, 2). Our laboratory found that the stability of these non-B DNA structures was different in mammals versus Escherichia coli (E.coli) bacteria (1, 2). One explanation for the difference between these species may be a result of how DSBs are repaired within each species. Non-homologous end-joining (NHEJ) is primed to repair DSBs in mammalian cells, while bacteria that lack NHEJ (such as E.coli), utilize homologous recombination (HR) to repair DSBs. To investigate the role of the error-prone NHEJ repair pathway in DNA structure-induced genetic instability, E.coli cells were modified to express genes to allow for a functional NHEJ system under different HR backgrounds. The Mycobacterium tuberculosis NHEJ sufficient system is composed of Ku and Ligase D (LigD) (3). These inducible NHEJ components were expressed individually and together in E.coli cells, with or without functional HR (RecA/RecB), and the Z-DNA and H-DNA-induced mutations were characterized. The Z-DNA structure gave rise to higher mutation frequencies compared to the controls, regardless of the DSB repair pathway(s) available; however, the type of mutants produced after repair was greatly dictated on the available DSB repair system, indicated by the shift from 2% large-scale deletions in the total mutant population to 24% large-scale deletions when NHEJ was present (4). This suggests that NHEJ has a role in the large deletions induced by Z-DNA-forming sequences. H-DNA structure, however, did not exhibit an increase in mutagenesis in the newly engineered E.coli environment, suggesting the involvement of other factors in regulating H-DNA formation/stability in bacterial cells. Accurate repair by established DNA DSB repair pathways is essential to maintain the stability of eukaryotic and prokaryotic genomes and our results suggest that an error-prone NHEJ pathway was involved in non-B DNA structure-induced mutagenesis in both prokaryotes and eukaryotes.

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.