D. Chatziplis

Quantitative Trait Loci Involved in Sex Determination and Body Growth in the Gilthead Sea Bream (Sparus aurata L.) through Targeted Genome Scan

Among vertebrates, teleost fish exhibit a considerably wide range of sex determination patterns that may be influenced by extrinsic parameters. However even for model fish species like the zebrafish Danio rerio the precise mechanisms involved in primary sex determination have not been studied extensively. The zebrafish, a gonochoristic species, is lacking discernible sex chromosomes and the sex of juvenile fish is difficult to determine. Sequential protandrous hermaphrodite species provide distinct determination of the gender and allow studying the sex determination process by looking at the mechanism of sex reversal. This is the first attempt to understand the genetic basis of phenotypic variation for sex determination and body weight in a sequential protandrous hermaphrodite species, the gilthead sea bream (Sparus aurata). This work demonstrates a fast and efficient strategy for Quantitative Trait Loci (QTL) detection in the gilthead sea bream, a non-model but target hermaphrodite fish species. Therefore a comparative mapping approach was performed to query syntenies against two other Perciformes, the European sea bass (Dicentrarchus labrax), a gonochoristic species and the Asian sea bass (Lates calcarifer) a protandrous hermaphrodite. In this manner two significant QTLs, one QTL affecting both body weight and sex and one QTL affecting sex, were detected on the same linkage group. The co-segregation of the two QTLs provides a genomic base to the observed genetic correlation between these two traits in sea bream as well as in other teleosts. The identification of QTLs linked to sex reversal and growth, will contribute significantly to a better understanding of the complex nature of sex determination in S. aurata where most individuals reverse to the female sex at the age of two years through development and maturation of the ovarian portion of the gonad and regression of the testicular area. [Genomic sequences reported in this manuscript have been submitted to GenBank under accession numbers HQ021443–HQ021749.]

Quantitative trait loci for body growth and sex determination in the hermaphrodite teleost fish Sparus aurata L.

Gilthead sea bream (Sparus aurata L.) is an important marine fish in Mediterranean aquaculture. Sex determination by age and/or body weight is a critical life-history trait, the genetic basis for which is largely unknown in this sequential hermaphrodite species. Herein, we performed a partial genome scan to map quantitative trait loci (QTL) affecting body weight and sex using 74 informative microsatellite markers from 10 paternal half-sib families to construct nine linkage groups (LG). In total, four growth-related QTL (two chromosome-wide and two genome-wide) and six QTL related to sex determination (three pairs in three different LGs) were detected (two chromosome-wide and one genome-wide). The proportion of phenotypic variation explained by the body-weight QTL ranged from 9.3% to 17.2%, showing their potential for use in marker-assisted selection. The results obtained offer solid ground to investigate the structure and function of the genomic regions involved in the mechanisms of sex reversal.

Segregation analysis of blood oxygen saturation in broilers suggests a major gene influence on ascites.

1. Blood oxygen saturation (SaO) is a potential indicator trait for resistance to ascites in chickens. 2. The objective of the study was to investigate the genetic architecture of SaO in a meat-type chicken line reared in commercial conditions. 3. Data were collected over 15 generations of selection and were divided into two data sets on the basis of a change in recording age from 6 to 5 weeks of age, approximately halfway through the period. The resulting pedigrees comprised in excess of 90 000 birds each and, on average, 12% of these birds had SaO records. 4. Segregation analyses of SaO were carried out assuming a mixed inheritance model that included a major locus segregating in a polygenic background. 5. The analyses suggest that a major gene is involved in the genetic control of SaO in this line. The putative gene acts in a dominant fashion and has an additive effect of around 0·90 σ p , equivalent to a predicted difference in SaO between the two homozygous classes of more than 10%. The frequency of the allele that increases SaO changed from 0·53 to 0·65 from the first to the second set of data, consistent with selection on SaO scores. 6. Using estimated genotype probabilities at the putative major locus, we inferred that it acts in an overdominant fashion on body weight and fleshing score. If the low SaO allele leads to susceptibility to ascites, its combined effects are consistent with it being maintained in the population by a balance of natural selection on fitness and artificial selection on growth and carcase traits. 7. Even with selection on both SaO and growth traits, the combined genotypic effects would make it difficult to remove the unfavourable low-SaO allele by means of traditional selection without the use of genetic markers.

Exploring a Nonmodel Teleost Genome Through RAD Sequencing—Linkage Mapping in Common Pandora, Pagellus erythrinus and Comparative Genomic Analysis

Common pandora (Pagellus erythrinus) is a benthopelagic marine fish belonging to the teleost family Sparidae, and a newly recruited species in Mediterranean aquaculture. The paucity of genetic information relating to sparids, despite their growing economic value for aquaculture, provides the impetus for exploring the genomics of this fish group. Genomic tool development, such as genetic linkage maps provision, lays the groundwork for linking genotype to phenotype, allowing fine-mapping of loci responsible for beneficial traits. In this study, we applied ddRAD methodology to identify polymorphic markers in a full-sib family of common pandora. Employing the Illumina MiSeq platform, we sampled and sequenced a size-selected genomic fraction of 99 individuals, which led to the identification of 920 polymorphic loci. Downstream mapping analysis resulted in the construction of 24 robust linkage groups, corresponding to the karyotype of the species. The common pandora linkage map showed varying degrees of conserved synteny with four other teleost genomes, namely the European seabass (Dicentrarchus labrax), Nile tilapia (Oreochromis niloticus), stickleback (Gasterosteus aculeatus), and medaka (Oryzias latipes), suggesting a conserved genomic evolution in Sparidae. Our work exploits the possibilities of genotyping by sequencing to gain novel insights into genome structure and evolution. Such information will boost the study of cultured species and will set the foundation for a deeper understanding of the complex evolutionary history of teleosts.

European Sea Bass

Aquaculture of European sea bass (Dicentrarchus labrax L.) has taken off in the coastal regions of the Mediterranean Sea and southeastern Atlantic Ocean over the past 25 years and increased to 71,649 metric tons in 2004. Genetic support for this industry was initially limited to cytogenetics and population genetics, but with time it has been complemented with selective breeding, as well as functional and comparative genomics. The haploid genome of sea bass consists of 24 chromosomes, weighing 0.78 pg and containing approximately 1,525 Mb. A number of different types of genetic markers are available. A first-generation linkage map based on 174 microsatellite markers covers 25 linkage groups (815 cM). A draft of an updated linkage map, including 369 microsatellite and AFLP markers, is now available. EST resources based on at least 17 cDNA tissue libraries and surpassing 30,000 sequence traces have been generated. A large insert BAC library has a 13× genomic coverage. Breeding goals have been established and heritability values of various traits measured. Functional genomic analysis in relation to the reproductive biology and stress physiology are in progress. A pilot analysis has detected a QTL for body length on the terminal end of linkage group 1. All these resources bring European sea bass into the group of the top ten genome resource-rich fish species. Additional genomic resources such as EST sequences, macro- and micro-arrays, a second-generation linkage map, and physical maps based on BAC fingerprints and radiation hybrids will become available in the near future. Selective breeding of this species is expected to direct it progressively toward complete domestication.

Second generation genetic linkage map for the gilthead sea bream Sparus aurata L.

An updated second linkage map was constructed for the gilthead sea bream, Sparus aurata L., a fish species of great economic importance for the Mediterranean aquaculture industry. In contrast to the first linkage map which mainly consisted of genomic microsatellites (SSRs), the new linkage map is highly enriched with SSRs found in Expressed Sequence Tags (EST-SSRs), which greatly facilitates comparative mapping with other teleosts. The new map consists of 321 genetic markers in 27 linkage groups (LGs): 232 genomic microsatellites, 85 EST-SSRs and 4 SNPs; of those, 13 markers were linked to LGs but were not ordered. Eleven markers (5 SSRs, 5 EST-SSRs and 1 SNP) are not assigned to any LG. The total length of the sex-averaged map is 1769.7cM, 42% longer than the previously published one, and the number of markers in each LG ranges from 2 to 30. The inter-marker distance varies from 0 to 75.6cM, with an average of 5.75cM. The male and female maps have a length of 1349.2 and 2172.1cM, respectively, and the average distance between markers is 4.38 and 7.05cM, respectively. Comparative mapping with the three-spined stickleback (Gasterosteus acuulatus) chromosomes and scaffolds showed conserved synteny with 132 S. aurata markers (42.9% of those mapped) having a hit on the stickleback genome.

Genetic parameters for blood oxygen saturation, body weight and breast conformation in 4 meat-type chicken lines

1. The objective of the study was to explore the genetic architecture of blood oxygen saturation (SaO) (an indicator trait, negatively correlated with ascites susceptibility), body weight (Weight) and fleshing score (Flesh, a measure of breast conformation) for 4 meat-type chicken lines reared in commercial conditions. 2. Genetic components, including heritabilities and genetic correlations, were estimated by Restricted Maximum likelihood for these traits measured at 6 weeks of age. 3. Data were collected over eight generations of selection and pedigrees comprised in excess of 130 000 birds. 4. Univariate analyses were performed to allow model definition and to obtain starting values for trivariate analyses. The basic model included a random animal effect and, in further models explored, a maternal environmental effect or a genetic maternal effect or both were fitted. Models were compared using likelihood ratio tests. 5. Estimated heritabilities for SaO ranged from 0·1 to 0·2, and there was no evidence of genetic maternal effects for SaO. The environmental maternal component was significant for one of the populations only. Estimated heritabilities for both Weight and Flesh were between 0·2 and 0·4, and there was evidence of environmental and genetic maternal effects for these traits in all populations. 6. Genetic correlations between SaO and Weight and between SaO and Flesh were low and negative. This suggests that, in principle, genetic selection to simultaneously increase SaO, and therefore decrease ascites susceptibility, and Weight and Flesh could be performed using traditional (marker-free) selection methods. We discuss how a putative interaction between ascites and production traits could jeopardise the success of such methods.

Selection and subsequent analysis of sib pair data for QTL detection

Haseman and Elston (1972) developed a robust regression method for the detection of linkage between a marker and a quantitative trait locus (QTL) using sib pair data. The principle underlying this method is that the difference in phenotypes between pairs of sibs becomes larger as they share a decreasing number of alleles at a particular QTL identical by descent (IBD) from their parents. In this case, phenotypically very different sibs will also on average share a proportion of alleles IBD at any marker linked to the QTL that is lower than the expected value of 0.5. Thus, the deviation of the proportion of marker alleles IBD from the expected value in pairs of sibs selected to be phenotypically different (i.e. discordant) can provide a test for the presence of a QTL. A simple regression method for QTL detection in sib pairs selected for high phenotypic differences is presented here. The power of the analytical method was found to be greater than the power obtained using the standard analysis when samples of sib pairs with high phenotypic differences were used. However, the use of discordant sib pairs was found to be less powerful for QTL detection than alternative selective genotyping schemes based on the phenotypic values of the sibs except with intense selection, when its advantage was only marginal. The most effective selection scheme overall was the use of sib pairs from entire families selected on the basis of high within-family variance for the trait in question. There is little effect of selection on QTL position estimates, which are in good agreement with the simulated values. However, QTL variance estimates are biased to a greater or lesser degree, depending on the selection method.

Selective genotyping for QTL detection using sib pair analysis in outbred populations with hierarchical structures.

A simulation study illustrates the effects of the inclusion of half-sib pairs as well as the effects of selective genotyping on the power of detection and the parameter estimates in a sib pair analysis of data from an outbred population. The power of QTL detection obtained from samples of sib pairs selected according to their within family variance or according to the mean within family variance within half sib family was compared and contrasted with the power obtained when only full sib pair analysis was used. There was an increase in power (4–16%) and decrease in the bias of parameter estimates with the use of half-sib information. These improvements in power and parameter estimates depended on the number of the half sib pairs (half sib family size). Almost the same power as that obtained using all the available sib pairs could be achieved by selecting only 50–60% the animals. The most effective method was to select both full and half sib pairs on the basis of high within full sib family variance for the trait in question. The QTL position estimates were in general slightly biased towards the center of the chromosome and the QTL variance estimates were biased upwards, there being quite large differences in bias depending on the selection method.

A practical approach to detect ancestral haplotypes in livestock populations

BackgroundThe effects of different evolutionary forces are expected to lead to the conservation, over many generations, of particular genomic regions (haplotypes) due to the development of linkage disequilibrium (LD). The detection and identification of early (ancestral) haplotypes can be used to clarify the evolutionary dynamics of different populations as well as identify selection signatures and genomic regions of interest to be used both in conservation and breeding programs. The aims of this study were to develop a simple procedure to identify ancestral haplotypes segregating across several generations both within and between populations with genetic links based on whole-genome scanning. This procedure was tested with simulated and then applied to real data from different genotyped populations of Spanish, Fleckvieh, Simmental and Brown-Swiss cattle.ResultsThe identification of ancestral haplotypes has shown coincident patterns of selection across different breeds, allowing the detection of common regions of interest on different bovine chromosomes and mirroring the evolutionary dynamics of the studied populations. These regions, mainly located on chromosomes BTA5, BTA6, BTA7 and BTA21 are related with certain animal traits such as coat colour and milk protein and fat content.ConclusionIn agreement with previous studies, the detection of ancestral haplotypes provides useful information for the development and comparison of breeding and conservation programs both through the identification of selection signatures and other regions of interest, and as indicator of the general genetic status of the populations.