Silvia T. Rodríguez-Ramilo

The effect of close relatives on unsupervised Bayesian clustering algorithms in population genetic structure analysis.

The inference of population genetic structures is essential in many research areas in population genetics, conservation biology and evolutionary biology. Recently, unsupervised Bayesian clustering algorithms have been developed to detect a hidden population structure from genotypic data, assuming among others that individuals taken from the population are unrelated. Under this assumption, markers in a sample taken from a subpopulation can be considered to be in Hardy–Weinberg and linkage equilibrium. However, close relatives might be sampled from the same subpopulation, and consequently, might cause Hardy–Weinberg and linkage disequilibrium and thus bias a population genetic structure analysis. In this study, we used simulated and real data to investigate the impact of close relatives in a sample on Bayesian population structure analysis. We also showed that, when close relatives were identified by a pedigree reconstruction approach and removed, the accuracy of a population genetic structure analysis can be greatly improved. The results indicate that unsupervised Bayesian clustering algorithms cannot be used blindly to detect genetic structure in a sample with closely related individuals. Rather, when closely related individuals are suspected to be frequent in a sample, these individuals should be first identified and removed before conducting a population structure analysis.

QTL detection for Aeromonas salmonicida resistance related traits in turbot ( Scophthalmus maximus )

BackgroundInteractions between fish and pathogens, that may be harmless under natural conditions, often result in serious diseases in aquaculture systems. This is especially important due to the fact that the strains used in aquaculture are derived from wild strains that may not have had enough time to adapt to new disease pressures. The turbot is one of the most promising European aquaculture species. Furunculosis, caused by the bacterium Aeromonas salmonicida, produces important losses to turbot industry. An appealing solution is to achieve more robust broodstock, which can prevent or diminish the devastating effects of epizooties. Genomics strategies have been developed in turbot to look for candidate genes for resistance to furunculosis and a genetic map with appropriate density to screen for genomic associations has been also constructed. In the present study, a genome scan for QTL affecting resistance and survival to A. salmonicida in four turbot families was carried out. The objectives were to identify consistent QTL using different statistical approaches (linear regression and maximum likelihood) and to locate the tightest associated markers for their application in genetic breeding strategies.ResultsSignificant QTL for resistance were identified by the linear regression method in three linkage groups (LGs 4, 6 and 9) and for survival in two LGs (6 and 9). The maximum likelihood methodology identified QTL in three LGs (5, 6 and 9) for both traits. Significant association between disease traits and genotypes was detected for several markers, some of them explaining up to 17% of the phenotypic variance. We also identified candidate genes located in the detected QTL using data from previously mapped markers.ConclusionsSeveral regions controlling resistance to A. salmonicida in turbot have been detected. The observed concordance between different statistical methods at particular linkage groups gives consistency to our results. The detected associated markers could be useful for genetic breeding strategies. A finer mapping will be necessary at the detected QTL intervals to narrow associations and around the closely associated markers to look for candidate genes through comparative genomics or positional cloning strategies. The identification of associated variants at specific genes will be essential, together with the QTL associations detected in this study, for future marker assisted selection programs.

Management of genetic diversity of subdivided populations in conservation programmes

Population subdivision must be explicitly considered in the management of conservation programmes, as most populations of wild species at risk of extinction and those kept in captivity are spatially structured. The partition of gene and allelic diversity in within- and between-subpopulation components allows for the integral management of populations. We summarise the main aspects of this partition and some of its applications in terms of priorisation of populations for conservation and establishment of synthetic populations. The procedures for the maintenance of diversity in subdivided populations making use of molecular markers and its implementation by the software METAPOP are illustrated with empirical data.

Identification of Quantitative Trait Loci Associated with Resistance to Viral Haemorrhagic Septicaemia (VHS) in Turbot (Scophthalmus maximus): A Comparison Between Bacterium, Parasite and Virus Diseases

One of the main objectives of genetic breeding programs in turbot industry is to reduce disease-related mortality. In the present study, a genome scan to detect quantitative trait loci (QTL) affecting resistance and survival to viral haemorrhagic septicaemia (VHS) was carried out. Three full-sib families with approximately 90 individuals each were genotyped and evaluated by linear regression and maximum likelihood approaches. In addition, a comparison between QTL detected for resistance and survival time to other important bacterial and parasite diseases affecting turbot (furunculosis and scuticociliatosis) was also carried out. Finally, the relationship between QTL affecting resistance/survival time to the virus and growth-related QTL was also evaluated. Several genomic regions controlling resistance and survival time to VHS were detected. Also significant associations between the evaluated traits and genotypes at particular markers were identified, explaining up to 14xa0% of the phenotypic variance. Several genomic regions controlling general and specific resistance to different diseases in turbot were detected. A preliminary gene mining approach identified candidate genes related to general or specific immunity. This information will be valuable to develop marker-assisted selection programs and to discover candidate genes related to disease resistance to improve turbot production.

Uncovering QTL for resistance and survival time to Philasterides dicentrarchi in turbot (Scophthalmus maximus).

Disease resistance-related traits have received increasing importance in aquaculture breeding programs worldwide. Currently, genomic information offers new possibilities in breeding to address the improvement of this kind of traits. The turbot is one of the most promising European aquaculture species, and Philasterides dicentrarchi is a scuticociliate parasite causing fatal disease in farmed turbot. An appealing approach to fight against disease is to achieve a more robust broodstock, which could prevent or diminish the devastating effects of scuticociliatosis on farmed individuals. In the present study, a genome scan for quantitative trait loci (QTL) affecting resistance and survival time to P.xa0dicentrarchi in four turbot families was carried out. The objectives were to identify QTL using different statistical approaches [linear regression (LR) and maximum likelihood (ML)] and to locate significantly associated markers for their application in genetic breeding strategies. Several genomic regions controlling resistance and survival time to P.xa0dicentrarchi were detected. When analyzing each family separately, significant QTL for resistance were identified by the LR method in two linkage groups (LG1 and LG9) and for survival time in LG1, while the ML methodology identified QTL for resistance in LG9 and LG23 and for survival time in LG6 and LG23. The analysis of the total data set identified an additional significant QTL for resistance and survival time in LG3 with the LR method. Significant association between disease resistance-related traits and genotypes was detected for several markers, a single one explaining up to 22% of the phenotypic variance. Obtained results will be essential to identify candidate genes for resistance and to apply them in marker-assisted selection programs to improve turbot production.

Relaxation of Selection With Equalization of Parental Contributions in Conservation Programs: An Experimental Test With Drosophila melanogaster

Equalization of parental contributions is one of the most simple and widely recognized methods to maintain genetic diversity in conservation programs, as it halves the rate of increase in inbreeding and genetic drift. It has, however, the negative side effect of implying a reduced intensity of natural selection so that deleterious genes are less efficiently removed from the population with possible negative consequences on the reproductive capacity of the individuals. Theoretical results suggest that the lower fitness resulting from equalization of family sizes relative to that for free contribution schemes is expected to be substantial only for relatively large population sizes and after many generations. We present a long-term experiment with Drosophila melanogaster, comparing the fitness performance of lines maintained with equalization of contributions (EC) and others maintained with no management (NM), allowing for free matings and contributions from parents. Two (five) replicates of size N = 100 (20) individuals of each type of line were maintained for 38 generations. As expected, EC lines retained higher gene diversity and allelic richness for four microsatellite markers and a higher heritability for sternopleural bristle number. Measures of life-history traits, such as egg-to-adult viability, mating success, and global fitness declined with generations, but no significant differences were observed between EC and NM lines. Our results, therefore, provide no evidence to suggest that equalization of family sizes entails a disadvantage on the reproductive capacity of conserved populations in comparison with no management procedures, even after long periods of captivity.

A new method for the partition of allelic diversity within and between subpopulations

A method is proposed for the analysis of allelic diversity in the context of subdivided populations. The definition of an allelic distance between subpopulations allows for the partition of total allelic diversity into within- and between-subpopulation components, in a way analogous to the classical partition of gene diversity. A new definition of allelic differentiation, AST, between subpopulations results from this partition, and is contrasted with the concept of allelic richness differentiation. The partition of allelic diversity makes it possible to establish the relative contribution of each subpopulation to within and between-subpopulation components of diversity with implications in priorisation for conservation. A comparison between this partition and that corresponding to allelic richness is illustrated with an example. Computer simulations are used to investigate the behaviour of the new statistic AST in comparison with FST for a finite island model under a range of mutation and migration rates. AST has less dependence on migration rate than FST for large values of migration rate, but the opposite occurs for low migration rates. In addition, the variance in the estimates of AST is higher than that of FST for low mutation rates, but the opposite for high mutation rates.

Mutation-selection balance accounting for genetic variation for viability in Drosophila melanogaster as deduced from an inbreeding and artificial selection experiment.

We carried out an experiment of inbreeding and upward artificial selection for egg‐to‐adult viability in a recently captured population of Drosophila melanogaster, as well as computer simulations of the experimental design, in order to obtain information on the nature of genetic variation for this important fitness component. The inbreeding depression was linear with a rate of 0.70u2003±u20030.11% of the initial mean per 1% increase in inbreeding coefficient, and the realized heritability was 0.06u2003±u20030.07. We compared the empirical observations of inbreeding depression and selection response with computer simulations assuming a balance between the occurrence of partially recessive deleterious mutations and their elimination by selection. Our results suggest that a model assuming mutation‐selection balance with realistic mutational parameters can explain the genetic variation for viability in the natural population studied. Several mutational models are incompatible with some observations and can be discarded. Mutational models assuming a low rate of mutations of large average effect and highly recessive gene action, and others assuming a high rate of mutations of small average effect and close to additive gene action, are compatible with all the observations.

LACK OF NONADDITIVE GENETIC EFFECTS ON EARLY FECUNDITY IN DROSOPHILA MELANOGASTER

Abstract Fecundity is usually considered as a trait closely connected to fitness and is expected to exhibit substantial nonadditive genetic variation and inbreeding depression. However, two independent experiments, using populations of different geographical origin, indicate that early fecundity inDrosophila melanogaster behaves as a typical additive trait of low heritability. The first experiment involved artificial selection in inbred and non‐inbred lines, all of them started from a common base population previously maintained in the laboratory for about 35 generations. The realized heritability estimate was 0.151 ± 0.075 and the inbreeding depression was very small and nonsignificant (0.09 ± 0.09% of the non‐inbred mean per 1% increase in inbreeding coefficient). With inbreeding, the observed decrease in the within‐line additive genetic variance and the corresponding increase of the between‐line variance were very close to their expected values for pure additive gene action. This result is at odds with previous studies showing inbreeding depression and, therefore, directional dominance for the same trait and species. All experiments, however, used laboratory populations, and it is possible that the original genetic architecture of the trait in nature was subsequently altered by the joint action of random drift and adaptation to captivity. Thus, we carried out a second experiment, involving inbreeding without artificial selection in a population recently collected from the wild. In this case we obtained, again, a maximum‐likelihood heritability estimate of 0.210 ± 0.027 and very little nonsignificant inbreeding depression (0.06 ± 0.12%). The results suggest that, for fitness‐component traits, low levels of additive genetic variance are not necessarily associated with large inbreeding depression or high levels of nonadditive genetic variance.

The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged‐Banded (RB) and Smooth‐Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid‐shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome‐wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow.