Landscape genomics reveals regions associated with adaptive phenotypic and genetic variation in Ethiopian indigenous chickens

Climate change is a threat to sustainable livestock production and livelihoods in the tropics. It has adverse impacts on feed and water availability, disease prevalence, production, environmental temperature, and biodiversity. Unravelling the drivers of local adaptation and understanding the underlying genetic variation in random mating indigenous livestock populations informs the design of genetic improvement programmes that aim to increase productivity and resilience. In the present study, we combined environmental, genomic, and phenotypic information of Ethiopian indigenous chickens to investigate their environmental adaptability. Through a hybrid sampling strategy, we captured wide biological and ecological variabilities across the country. Our environmental dataset comprised mean values of 34 climatic, vegetation and soil variables collected over a thirty-year period for 260 geolocations. Our biological dataset included whole genome sequences and quantitative measurements (on eight traits) from 513 individuals, representing 26 chicken populations spread along 4 elevational gradients (6-7 populations per gradient). We performed signatures of selection analyses ([Formula: see text] and XP-EHH) to detect footprints of natural selection, and redundancy analyses (RDA) to determine genotype-environment and genotype-phenotype-associations. RDA identified 1909 outlier SNPs linked with six environmental predictors, which have the highest contributions as ecological drivers of adaptive phenotypic variation. The same method detected 2430 outlier SNPs that are associated with five traits. A large overlap has been observed between signatures of selection identified by[Formula: see text]and XP-EHH showing that both methods target similar selective sweep regions. Average genetic differences measured by [Formula: see text] are low between gradients, but XP-EHH signals are the strongest between agroecologies. Genes in the calcium signalling pathway, those associated with the hypoxia-inducible factor (HIF) transcription factors, and sports performance (GALNTL6) are under selection in high-altitude populations. Our study underscores the relevance of landscape genomics as a powerful interdisciplinary approach to dissect adaptive phenotypic and genetic variation in random mating indigenous livestock populations.

Genotype by Environment Interactions (G*E) of Chickens Tested in Ethiopia Using Body Weight as a Performance Trait

Ethiopia is a developing nation that could highly benefit from securing food via improved smallholder poultry farming. To support farmer and breeding decisions regarding which chicken strain to use in which Ethiopian environment, G*E analyses for body weight (BW) of growing male and female chickens were conducted. Research questions were (1) if a G*E is present for BW and (2) which strain performs best in which environment in terms of predicted BW. Analyses were performed using predicted BW at four different ages (90, 120, 150, and 180 days) of five strains (Horro, Koekoek, Kuroiler, Sasso-Rhode Island Red (S-RIR), and Sasso) tested in five Ethiopian regions (Addis Ababa, Amhara, Oromia, South Region, and Tigray) that are part of three Agro-Ecological Zones (AEZ) (cool humid, cool sub-humid, and warm semi-arid). The indigenous Horro strain was used as a control group to compare four other introduced tropically adapted strains. The dataset consisted of 999 female and 989 male farm-average BW measurements. G*E was strongly present (p < 0.001) for all combinations of strain and region analyzed. In line with previous research, Sasso was shown to have the highest predicted BW, especially at an early age, followed by Kuroiler. Horro had the lowest predicted BW at most ages and in most regions, potentially due to its young breeding program. The highest predicted BW were observed in Tigray, Oromia, and Amhara regions, which are in the main part of the cool sub-humid AEZ.

The induction of oocyte maturation and ovulation in the European eel (Anguilla anguilla): in vitro and in vivo comparison of progesterone with 17α,20β-dihydroxy-4-pregnen-3-one

Ovulation in European eel is induced by injection of 17α,20β-dihydroxy-4-pregnen-3-one (DHP) as the maturation-inducing hormone (MIH). Female eels need to ovulate within 18 h after injection to release good quality eggs. Progesterone (P), as an upstream precursor of DHP, may promote endogenous DHP production and improve egg quality. The purpose of this study was therefore to compare treatment of P with DHP on batch level, in vitro, to determine dose-response effects, and in vivo, at a single dose. For the in vitro experiment, ovarian tissue was extracted and placed in culture plates containing hormone-free medium and media supplemented with the treatment: DHP at 1, 10 and 100 ng mL-1, or P at 10, 100 and 1,000 ng mL-1. At the start of incubation, the folliculated oocytes were sampled for histology, microscopy and qPCR. After incubation for 12 and 18 h, the oocytes were sampled for microscopy and qPCR analysis. For the in vivo experiment, females were either injected with DHP or P at a dose of 2 mg kg-1 to assess their effects on ovulation and reproductive success. At the moment of release, eggs were sampled for RNA sequencing to compare effects of DHP and P on the expression of genes involved in egg quality aspects. Remaining eggs were fertilized and larval viability was recorded. Both DHP and P were able to induce GVBD (DHP at 10 and 100 ng mL-1, P at 100 and 1,000 ng mL-1) in vitro. Expression of genes involved in oocyte maturation and ovulation was similar in vitro for both DHP and P treatments. Regarding the in vivo results, RNAseq results reflected similar DHP and P effects on the expression of genes involved in egg quality aspects. Females injected with either DHP or P ovulated, released eggs, and were equally able to produce larvae without any differences in reproductive success. Our results support the conclusion that DHP and P work equally well in vitro and in vivo. P is more attractive to apply as the price is 3,000 times lower than the price of DHP.

Species and Phenotypic Distribution Models Reveal Population Differentiation in Ethiopian Indigenous Chickens

Smallholder poultry production dominated by indigenous chickens is an important source of livelihoods for most rural households in Ethiopia. The long history of domestication and the presence of diverse agroecologies in Ethiopia create unique opportunities to study the effect of environmental selective pressures. Species distribution models (SDMs) and Phenotypic distribution models (PDMs) can be applied to investigate the relationship between environmental variation and phenotypic differentiation in wild animals and domestic populations. In the present study we used SDMs and PDMs to detect environmental variables related with habitat suitability and phenotypic differentiation among nondescript Ethiopian indigenous chicken populations. 34 environmental variables (climatic, soil, and vegetation) and 19 quantitative traits were analyzed for 513 adult chickens from 26 populations. To have high variation in the dataset for phenotypic and ecological parameters, animals were sampled from four spatial gradients (each represented by six to seven populations), located in different climatic zones and geographies. Three different ecotypes are proposed based on correlation test between habitat suitability maps and phenotypic clustering of sample populations. These specific ecotypes show phenotypic differentiation, likely in response to environmental selective pressures. Nine environmental variables with the highest contribution to habitat suitability are identified. The relationship between quantitative traits and a few of the environmental variables associated with habitat suitability is non-linear. Our results highlight the benefits of integrating species and phenotypic distribution modeling approaches in characterization of livestock populations, delineation of suitable habitats for specific breeds, and understanding of the relationship between ecological variables and quantitative traits, and underlying evolutionary processes.

Genome-wide association analysis of adaptation to oxygen stress in Nile tilapia (Oreochromis niloticus)

Background

Tilapia is one of the most abundant species in aquaculture. Hypoxia is known to depress growth rate, but the genetic mechanism by which this occurs is unknown. In this study, two groups consisting of 3140 fish that were raised in either aerated (normoxia) or non-aerated pond (nocturnal hypoxia). During grow out, fish were sampled five times to determine individual body weight (BW) gains. We applied a genome-wide association study to identify SNPs and genes associated with the hypoxic and normoxic environments in the 16th generation of a Genetically Improved Farmed Tilapia population.

Results

In the hypoxic environment, 36 SNPs associated with at least one of the five body weight measurements (BW1 till BW5), of which six, located between 19.48 Mb and 21.04 Mb on Linkage group (LG) 8, were significant for body weight in the early growth stage (BW1 to BW2). Further significant associations were found for BW in the later growth stage (BW3 to BW5), located on LG1 and LG8. Analysis of genes within the candidate genomic region suggested that MAPK and VEGF signalling were significantly involved in the later growth stage under the hypoxic environment. Well-known hypoxia-regulated genes such as igf1rb, rora, efna3 and aurk were also associated with growth in the later stage in the hypoxic environment. Conversely, 13 linkage groups containing 29 unique significant and suggestive SNPs were found across the whole growth period under the normoxic environment. A meta-analysis showed that 33 SNPs were significantly associated with BW across the two environments, indicating a shared effect independent of hypoxic or normoxic environment. Functional pathways were involved in nervous system development and organ growth in the early stage, and oocyte maturation in the later stage.

Conclusions

There are clear genotype-growth associations in both normoxic and hypoxic environments, although genome architecture involved changed over the growing period, indicating a transition in metabolism along the way. The involvement of pathways important in hypoxia especially at the later growth stage indicates a genotype-by-environment interaction, in which MAPK and VEGF signalling are important components.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07486-5.

What Goes Wrong during Early Development of Artificially Reproduced European Eel Anguilla anguilla? Clues from the Larval Transcriptome and Gene Expression Patterns

Simple Summary

Closing the life cycle of the European eel in captivity is urgently needed to gain perspective for the commercial production of juvenile glass eels. Larvae are produced weekly at our facilities, but large variations in larval mortality are observed during the first week after hatching. Although much effort has been devoted to investigating ways to prevent early larval mortality, it remains unclear what the causes are. The aim of this study was to perform a transcriptomic study on European eel larvae in order to identify genes and physiological pathways that are differentially regulated in the comparison of larvae from batches that did not survive for longer than three days vs. larvae from batches that survived for at least a week up to 22 days after hatching (non-viable vs. viable larvae). In contrast to earlier published studies on European eel, we conclude that larvae exhibit immune competency. Non-viable larvae initiated an inflammatory and host protection immune response and tried to maintain osmoregulatory homeostasis. As a perspective, microbial control and salinity reduction might benefit eel larvae in terms of lower mortality and improved development by lowering the costs of immune functioning and osmoregulation.

Abstract

In eels, large variations in larval mortality exist, which would impede the viable production of juvenile glass eels in captivity. The transcriptome of European eel larvae was investigated to identify physiological pathways and genes that show differential regulation between non-viable vs. viable larvae. Expression of genes involved in inflammation and host protection was higher, suggesting that non-viable larvae suffered from microbial infection. Expression of genes involved in osmoregulation was also higher, implying that non-viable larvae tried to maintain homeostasis by strong osmoregulatory adaptation. Expression of genes involved in myogenesis, neural, and sensory development was reduced in the non-viable larvae. Expression of the major histocompatibility complex class-I (mhc1) gene, M-protein (myom2), the dopamine 2B receptor (d2br), the melatonin receptor (mtr1), and heat-shock protein beta-1 (hspb1) showed strong differential regulation and was therefore studied in 1, 8, and 15 days post-hatch (dph) larvae by RT-PCR to comprehend the roles of these genes during ontogeny. Expression patterning of these genes indicated the start of active swimming (8 dph) and feed searching behavior (15 dph) and confirmed immunocompetence immediately after hatching. This study revealed useful insights for improving larval survival by microbial control and salinity reduction.

Swimming Performance and Oxygen Consumption as Non-lethal Indicators of Production Traits in Atlantic Salmon and Gilthead Seabream

The aim of this study was to investigate swimming performance and oxygen consumption as non-lethal indicator traits of production parameters in Atlantic salmon Salmo salar L. and Gilthead seabream Sparus aurata L. A total of 34 individual fish of each species were subjected to a series of experiments: (1) a critical swimming speed (Ucrit) test in a swim-gutter, followed by (2) two starvation-refeeding periods of 42 days, and (3) swimming performance experiments coupled to respirometry in swim-tunnels. Ucrit was assessed first to test it as a predictor trait. Starvation-refeeding traits included body weight; feed conversion ratio based on dry matter; residual feed intake; average daily weight gain and loss. Swim-tunnel respirometry provided oxygen consumption in rest and while swimming at the different speeds, optimal swim speed and minimal cost of transport (COT). After experiments, fish were dissected and measured for tissue weights and body composition in terms of dry matter, ash, fat, protein and moist, and energy content. The Ucrit test design was able to provide individual Ucrit values in high throughput manner. The residual Ucrit (RUcrit) should be considered in order to remove the size dependency of swimming performance. Most importantly, RUcrit predicted filet yield in both species. The minimal COT, the oxygen consumption when swimming at Uopt, added predictive value to the seabream model for feed intake.

Optimizing design to estimate genetic correlations between environments with common environmental effects

Breeding programs for different species aim to improve performance by testing members of full-sib (FS) and half-sib (HS) families in different environments. When genotypes respond differently to changes in the environment, this is defined as genotype by environment (G × E) interaction. The presence of common environmental effects within families generates covariance between siblings, and these effects should be taken into account when estimating a genetic correlation. Therefore, an optimal design should be established to accurately estimate the genetic correlation between environments in the presence of common environmental effects. We used stochastic simulation to find the optimal population structure using a combination of FS and HS groups with different levels of common environmental effects. Results show that in a population with a constant population size of 2,000 individuals per environment, ignoring common environmental effects when they are present in the population will lead to an upward bias in the estimated genetic correlation of on average 0.3 when the true genetic correlation is 0.5. When no common environmental effects are present in the population, the lowest standard error (SE) of the estimated genetic correlation was observed with a mating ratio of one dam per sire, and 10 offspring per sire per environment. When common environmental effects are present in the population and are included in the model, the lowest SE is obtained with mating ratios of at least 5 dams per sire and with a minimum number of 10 offspring per sire per environment. We recommend that studies that aim to estimate the magnitude of G × E in pigs, chicken, and fish should acknowledge the potential presence of common environmental effects and adjust the mating ratio accordingly.

The genetic correlation between feed conversion ratio and growth rate affects the design of a breeding program for more sustainable fish production

BACKGROUND

Most fish breeding programs aim at improving growth rate and include feed conversion ratio (FCR) neither in the breeding goal nor in the selection index, although decreasing FCR is known to increase farm profit and decrease environmental impacts. This is because FCR is difficult to measure in fish that live in groups and FCR is assumed to have a favourable (negative) genetic correlation with growth, although the magnitude of this correlation is unknown. We investigated the effect of the genetic correlation between growth and FCR on the economic and environmental responses of a two-trait breeding goal (growth and FCR), compared to a single-trait breeding goal (growth only). Next, we evaluated the weights to assign to growth and FCR in a two-trait breeding goal to maximize sustainability of fish production.

METHODS

We used pseudo-best linear unbiased prediction (BLUP) index calculations to simulate a breeding program for sea bass. For the single-trait breeding goal, the trait in the breeding goal and in the index was thermal growth coefficient (TGC) and for the two-trait breeding goal, the traits in the breeding goal were TGC and FCR and the traits in the index were TGC and percentage of fat in the dorsal muscle (an indirect measure of FCR). We simulated responses to selection for genetic and phenotypic correlations between TGC and FCR ranging from 0 to - 0.8. Then, in the two-trait breeding goal, we calculated the economic return and the change in eutrophication when using economic values (EV) or environmental values (ENV).

RESULTS

When the genetic correlation between TGC and FCR was lower than - 0.45, we found major differences in economic returns and in eutrophication between single and two-trait breeding programs. At a correlation of - 0.25, the two-trait breeding goal based on EV increased economic return by 25% compared to the single-trait breeding goal, while using ENV decreased eutrophication by 1.34% per ton of fish produced after one generation of selection.

CONCLUSIONS

The genetic correlation between TGC and FCR affects the magnitude of economic losses due to omitting FCR in the breeding program. In addition, the genetic correlation affects the importance of choosing EV or ENV to reduce eutrophication and increase profit.

Chronic simultaneous exposure of common carp (Cyprinus carpio) from embryonic to juvenile stage to drospirenone and gestodene at low ng/L level caused intersex

Synthetic progestins are emerging contaminants of the aquatic environment with endocrine disrupting potential. The main aim of the present study was to investigate the effects of the synthetic progestins gestodene, and drospirenone on sex differentiation in common carp (Cyprinus carpio) by histological analysis. To gain insights into the mechanisms behind the observations from the in vivo experiment on sex differentiation, we analyzed expression of genes involved in hypothalamus-pituitary-gonad (HPG) and hypothalamus-pituitary-thyroid (HPT) axes, histology of hepatopancreas, and in vitro bioassays. Carp were continuously exposed to concentrations of 2 ng/L of single progestins (gestodene or drospirenone) or to their mixture at concentration 2 ng/L of each. The exposure started 24 h after fertilization of eggs and concluded 160 days post-hatching. Our results showed that exposure of common carp to a binary mixture of drospirenone and gestodene caused increased incidence of intersex (32%) when compared to clean water and solvent control groups (both 3%). Intersex most probably was induced by a combination of multiple modes of action of the studied substances, namely anti-gonadotropic activity, interference with androgen receptor, and potentially also with HPT axis or estrogen receptor.

Optimized grouping to increase accuracy of prediction of breeding values based on group records in genomic selection breeding programs

BACKGROUND

Phenotypic records of group means or group sums are a good alternative to individual records for some difficult to measure, but economically important traits such as feed efficiency or egg production. Accuracy of predicted breeding values based on group records increases with increasing relationships between group members. The classical way to form groups with more closely-related animals is based on pedigree information. When genotyping information is available before phenotyping, its use to form groups may further increase the accuracy of prediction from group records. This study analyzed two grouping methods based on genomic information: (1) unsupervised clustering implemented in the STRUCTURE software and (2) supervised clustering that models genomic relationships.

RESULTS

Using genomic best linear unbiased prediction (GBLUP) models, estimates of the genetic variance based on group records were consistent with those based on individual records. When genomic information was available to constitute the groups, genomic relationship coefficients between group members were higher than when random grouping of paternal half-sibs and of full-sibs was applied. Grouping methods that are based on genomic information resulted in higher accuracy of genomic estimated breeding values (GEBV) prediction compared to random grouping. The increase was ~ 1.5% for full-sibs and ~ 11.5% for paternal half-sibs. In addition, grouping methods that are based on genomic information led to lower coancestry coefficients between the top animals ranked by GEBV. Of the two proposed methods, supervised clustering was superior in terms of accuracy, computation requirements and applicability. By adding surplus genotyped offspring (more genotyped offspring than required to fill the groups), the advantage of supervised clustering increased by up to 4.5% compared to random grouping of full-sibs, and by 14.7% compared to random grouping of paternal half-sibs. This advantage also increased with increasing family sizes or decreasing genome sizes.

CONCLUSIONS

The use of genotyping information for grouping animals increases the accuracy of selection when phenotypic group records are used in genomic selection breeding programs.

Genetic correlations between growth performance and carcass traits of purebred and crossbred pigs raised in tropical and temperate climates1

In pig breeding, selection commonly takes place in purebred (PB) pigs raised mainly in temperate climates (TEMP) under optimal environmental conditions in nucleus farms. However, pork production typically makes use of crossbred (CB) animals raised in nonstandardized commercial farms, which are located not only in TEMP regions but also in tropical and subtropical regions (TROP). Besides the differences in the genetic background of PB and CB, differences in climate conditions, and differences between nucleus and commercial farms can lower the genetic correlation between the performance of PB in the TEMP (PBTEMP) and CB in the TROP (CBTROP). Genetic correlations (rg) between the performance of PB and CB growing-finishing pigs in TROP and TEMP environments have not been reported yet, due to the scarcity of data in both CB and TROP. Therefore, the present study aimed 1) to verify the presence of genotype × environment interaction (G × E) and 2) to estimate the rg for carcass and growth performance traits when PB and 3-way CB pigs are raised in 2 different climatic environments (TROP and TEMP). Phenotypic records of 217,332 PB and 195,978 CB, representing 2 climatic environments: TROP (Brazil) and TEMP (Canada, France, and the Netherlands) were available for this study. The PB population consisted of 2 sire lines, and the CB population consisted of terminal 3-way cross progeny generated by crossing sires from one of the PB sire lines with commercially available 2-way maternal sow crosses. G × E appears to be present for average daily gain, protein deposition, and muscle depth given the rg estimates between PB in both environments (0.64 to 0.79). With the presence of G × E, phenotypes should be collected in TROP when the objective is to improve the performance of CB in the TROP. Also, based on the rg estimates between PBTEMP and CBTROP (0.22 to 0.25), and on the expected responses to selection, selecting based only on the performance of PBTEMP would give limited genetic progress in the CBTROP. The rg estimates between PBTROP and CBTROP are high (0.80 to 0.99), suggesting that combined crossbred-purebred selection schemes would probably not be necessary to increase genetic progress in CBTROP. However, the calculated responses to selection show that when the objective is the improvement of CBTROP, direct selection based on the performance of CBTROP has the potential to lead to the higher genetic progress compared with indirect selection on the performance of PBTROP.

Economic optimization of full-sib test group size and genotyping effort in a breeding program for Atlantic salmon

BACKGROUND

Breeding companies may want to maximize the rate of genetic gain from their breeding program within a limited budget. In salmon breeding programs, full-sibs of selection candidates are subjected to performance tests for traits that cannot be recorded on selection candidates. While marginal gains in the aggregate genotype from phenotyping and genotyping more full-sibs per candidate decrease, costs increase linearly, which suggests that there is an optimum in the allocation of the budget among these activities. Here, we studied how allocation of the fixed budget to numbers of phenotyped and genotyped test individuals in performance tests can be optimized.

METHODS

Gain in the aggregate genotype was a function of the numbers of full-sibs of selection candidates that were (1) phenotyped in a challenge test for sea lice resistance (2) phenotyped in a slaughter test (3) genotyped in the challenge test, and (4) genotyped in the slaughter test. Each of these activities was subject to budget constraints. Using a grid search, we optimized allocation of the budget among activities to maximize gain in the aggregate genotype. We performed sensitivity analyses on the maximum gain in the aggregate genotype and on the relative allocation of the budget among activities at the optimum.

RESULTS

Maximum gain in the aggregate genotype was €386/ton per generation. The response surface for gain in the aggregate genotype was rather flat around the optimum, but it curved strongly near the extremes. Maximum gain was sensitive to the size of the budget and the relative emphasis on breeding goal traits, but less sensitive to the accuracy of genomic prediction and costs of phenotyping and genotyping. The relative allocation of budget among activities at the optimum was sensitive to costs of phenotyping and genotyping and the relative emphasis on breeding goal traits, but was less sensitive to the accuracy of genomic prediction and the size of the budget.

CONCLUSIONS

There is an optimum allocation of budget to the numbers of full-sibs of selection candidates that are phenotyped and genotyped in performance tests that maximizes gain in the aggregate genotype. Although potential gains from optimizing group sizes and genotyping effort may be small, they come at no extra cost.

Derivation of the economic value of R0 for macroparasitic diseases and application to sea lice in salmon

BACKGROUND

Macroparasites, such as ticks, lice, and helminths, are a concern in livestock and aquaculture production, and can be controlled by genetic improvement of the host population. Genetic improvement should aim at reducing the rate at which parasites spread across the farmed population. This rate is determined by the basic reproduction ratio, i.e. [Formula: see text], which is the appropriate breeding goal trait. This study aims at providing a method to derive the economic value of [Formula: see text].

METHODS

Costs of a disease are the sum of production losses and expenditures on disease control. Genetic improvement of [Formula: see text] lowers the loss-expenditure frontier. Its economic effect depends on whether the management strategy is optimized or not. The economic value may be derived either from the reduction in losses with constant expenditures or from the reduction in expenditures with constant losses.

RESULTS

When [Formula: see text] ≤ 1, the economic value of a further reduction is zero because there is no risk of a major epidemic. When [Formula: see text] > 1 and management is optimized, the economic value increases with decreasing values of [Formula: see text], because both the mean number of parasites per host and frequency of treatments decrease at an increasing rate when [Formula: see text] decreases. When [Formula: see text] > 1 and management is not optimized, the economic value depends on whether genetic improvement is used for reducing expenditures or losses. For sea lice in salmon, the economic value depends on a reduction in expenditures with constant losses, and is estimated to be 0.065€/unit [Formula: see text]/kg production.

DISCUSSION

Response to selection for measures of disease prevalence cannot be predicted from quantitative genetic theory alone. Moreover, many studies fail to address the issue of whether genetic improvement results in reduced losses or expenditures. Using [Formula: see text] as the breeding goal trait, weighed by its appropriate economic value, avoids these issues.

CONCLUSION

When management is optimized, the economic value increases with decreasing values of [Formula: see text] (until the threshold of [Formula: see text], where it drops to zero). When management is not optimized, the economic value depends on whether genetic improvement is used for reduced expenditures or production losses. For sea lice in salmon, the economic value is estimated to be 0.065 €/unit [Formula: see text]/kg production.

Cost-benefit analysis of aquaculture breeding programs

Background

Profitability of breeding programs is a key determinant in the adoption of selective breeding, and can be evaluated using cost-benefit analysis. There are many options to design breeding programs, with or without a multiplier tier. Our objectives were to evaluate different breeding program designs for aquaculture and to optimize the number of selection candidates for these programs.

Methods

The baseline was based on an existing breeding program for gilthead seabream, where improvement of the nucleus had priority over improvement of the multiplier tier, which was partly replaced once every 3 years. Alternative breeding programs considered were annual multiplier tier replacement, annual multiplier tier replacement with priority on improvement of the multiplier tier, and a program without a multiplier tier. Cost-benefit analyses were performed to compare breeding programs. The outcomes were used to describe relationships between profitability and the number of selection candidates, length of the time horizon, and production output, and to estimate the optimum numbers of selection candidates.

Results

The baseline breeding program was profitable after 5 years and reached a net present value of 2.9 million euro in year 10. All alternative programs were more profitable up to year 17. The program without a multiplier tier was the most profitable one up to year 22, followed by the program with annual multiplier tier replacement and nucleus priority. The optimum number of selection candidates increased with the length of the time horizon and production output.

Conclusions

The baseline breeding program was profitable after 5 years. For a short time horizon, putting priority on improvement of the multiplier tier over the nucleus is more profitable than putting priority on nucleus improvement, and vice versa for a long time horizon. Use of a multiplier tier increases the delay between costs made for selection and resulting benefits. Thus, avoiding the use of a multiplier tier will increase the profitability of the breeding program in the short term. The optimum number of selection candidates increases with the length of the time horizon and production output. Using too many selection candidates relative to the optimum leads to less reduction in profitability than using too few selection candidates.

Multiple criteria decision-making process to derive consensus desired genetic gains for a dairy cattle breeding objective for diverse production systems

Dairy cattle industries contribute to food and nutrition security and are a source of income for numerous households in many developing countries. Selective breeding can enhance efficiency in these industries. Developing dairy industries are characterized by diverse production and marketing systems. In this paper, we use weighted goal aggregating procedure to derive consensus trait preferences for different producer categories and processors. We based the study on the dairy industry in Kenya. The analytic hierarchy process was used to derive individual preferences for milk yield (MY), calving interval (CIN), production lifetime (PLT), mature body weight (MBW), and fat yield (FY). Results show that classical classification of production systems into large-scale and smallholder systems does not capture all differences in trait preferences. These differences became apparent when classification was based on productivity at the individual animal level, with high and low intensity producers and processors as the most important groups. High intensity producers had highest preferences for PLT and MY, whereas low intensity producers had highest preference for CIN and PLT; processors preferred MY and FY the most. The highest disagreements between the groups were observed for FY, PLT, and MY. Individual and group preferences were aggregated into consensus preferences using weighted goal programming. Desired gains were obtained as a product of consensus preferences and percentage genetic gains (G%). These were 2.42, 0.22, 2.51, 0.15, and 0.87 for MY, CIN, PLT, MBW, and FY, respectively. Consensus preferences can be used to derive a single compromise breeding objective for situations where the same genetic resources are used in diverse production and marketing circumstances.

Economic evaluation of progeny-testing and genomic selection schemes for small-sized nucleus dairy cattle breeding programs in developing countries

In developing countries minimal and erratic performance and pedigree recording impede implementation of large-sized breeding programs. Small-sized nucleus programs offer an alternative but rely on their economic performance for their viability. We investigated the economic performance of 2 alternative small-sized dairy nucleus programs [i.e., progeny testing (PT) and genomic selection (GS)] over a 20-yr investment period. The nucleus was made up of 453 male and 360 female animals distributed in 8 non-overlapping age classes. Each year 10 active sires and 100 elite dams were selected. Populations of commercial recorded cows (CRC) of sizes 12,592 and 25,184 were used to produce test daughters in PT or to create a reference population in GS, respectively. Economic performance was defined as gross margins, calculated as discounted revenues minus discounted costs following a single generation of selection. Revenues were calculated as cumulative discounted expressions (CDE, kg) × 0.32 (€/kg of milk) × 100,000 (size commercial population). Genetic superiorities, deterministically simulated using pseudo-BLUP index and CDE, were determined using gene flow. Costs were for one generation of selection. Results show that GS schemes had higher cumulated genetic gain in the commercial cow population and higher gross margins compared with PT schemes. Gross margins were between 3.2- and 5.2-fold higher for GS, depending on size of the CRC population. The increase in gross margin was mostly due to a decreased generation interval and lower running costs in GS schemes. In PT schemes many bulls are culled before selection. We therefore also compared 2 schemes in which semen was stored instead of keeping live bulls. As expected, semen storage resulted in an increase in gross margins in PT schemes, but gross margins remained lower than those of GS schemes. We conclude that implementation of small-sized GS breeding schemes can be economically viable for developing countries.

Derivation of economic values for production traits in aquaculture species

BACKGROUND

In breeding programs for aquaculture species, breeding goal traits are often weighted based on the desired gains but economic gain would be higher if economic values were used instead. The objectives of this study were: (1) to develop a bio-economic model to derive economic values for aquaculture species, (2) to apply the model to determine the economic importance and economic values of traits in a case-study on gilthead seabream, and (3) to validate the model by comparison with a profit equation for a simplified production system.

METHODS

A bio-economic model was developed to simulate a grow-out farm for gilthead seabream, and then used to simulate gross margin at the current levels of the traits and after one genetic standard deviation change in each trait with the other traits remaining unchanged. Economic values were derived for the traits included in the breeding goal: thermal growth coefficient (TGC), thermal feed intake coefficient (TFC), mortality rate (M), and standard deviation of harvest weight ([Formula: see text]). For a simplified production system, improvement in TGC was assumed to affect harvest weight instead of growing period. Using the bio-economic model and a profit equation, economic values were derived for harvest weight, cumulative feed intake at harvest, and overall survival.

RESULTS

Changes in gross margin showed that the order of economic importance of the traits was: TGC, TFC, M, and [Formula: see text]. Economic values in € (kg production)^-1 (trait unit)^-1 were: 0.40 for TGC, -0.45 for TFC, -7.7 for M, and -0.0011 to -0.0010 for [Formula: see text]. For the simplified production system, similar economic values were obtained with the bio-economic model and the profit equation. The advantage of the profit equation is its simplicity, while that of the bio-economic model is that it can be applied to any aquaculture species, because it can include any limiting factor and/or environmental condition that affects production.

CONCLUSIONS

We confirmed the validity of the bio-economic model. TGC is the most important trait to improve, followed by TFC and M, and the effect of [Formula: see text] on gross margin is small.

Economic values of growth and feed efficiency for fish farming in recirculating aquaculture system with density and nitrogen output limitations: a case study with African catfish (Clarias gariepinus)

In fish farming, economic values (EV) of breeding goal traits are lacking, even though they are key parameters when defining selection objectives. The aim of this study was to develop a bioeconomic model to estimate EV of 2 traits representing production performances in fish farming: the thermal growth coefficient (TGC) and the feed conversion ratio (FCR). This approach was applied to a farm producing African catfish (Clarias gariepinus) in a recirculating aquaculture system (RAS). In the RAS, 2 factors could limit production level: the nitrogen treatment capacity of the biofilter or the fish density in rearing tanks at harvest. Profit calculation includes revenue from fish sales, cost of juveniles, cost of feed, cost of waste water treatment, and fixed costs. In the reference scenario, profit was modeled to zero. EV were calculated as the difference in profit per kilogram of fish between the current population mean for both traits (µt) and the next generation of selective breeding (µt+Δt) for either TGC or FCR. EV of TGC and FCR were calculated for three generations of hypothetical selection on either TGC or FCR (respectively 6.8% and 7.6% improvement per generation). The results show that changes in TGC and FCR can affect both the number of fish that can be stocked (number of batches per year and number of fish per batch) and the factor limiting production. The EV of TGC and FCR vary and depend on the limiting factors. When dissolved NH3-N is the limiting factor for both µt and µt+Δt, increasing TGC decreases the number of fish that can be stocked but increases the number of batches that can be grown. As a result, profit remains constant and EVTGC is zero. Increasing FCR, however, increases the number of fish stocked and the ratio of fish produced per kilogram of feed consumed ("economic efficiency"). The EVFCR is 0.14 €/kg of fish, and profit per kilogram of fish increases by about 10%. When density is the limiting factor for both µt and µt+Δt, the number of fish stocked per batch is fixed; therefore, extra profit is obtained by increasing either TGC, which increases the annual number of batches, or by decreasing FCR, which decreases annual feed consumption. EVTGC is 0.03 €/kg of fish and EVFCR is 0.05-0.06 €/kg of fish. These results emphasize the importance of calculating economic values in the right context to develop efficient future breeding programs in aquaculture.

Rates of inbreeding and genetic adaptation for populations managed as herds in zoos with a rotational mating system or with optimized contribution of parents

This study compares two genetic management scenarios for species kept in herds, such as deer. The simulations were designed so that their results can be extended to a wide range of zoo populations. In the first scenario, the simulated populations of size 3 × 20, 6 × 40 or 20 × 60 (herds × animals in herd) were managed with a rotational mating (RM) scheme in which 10%, 20% or 50% of males were selected for breeding and moved between herds in a circular fashion. The second scenario was based on optimal contribution theory (OC). OC requires an accurate pedigree to calculate kinship; males were selected and assigned numbers of offspring to minimize kinship in the next generation. RM was efficient in restriction of inbreeding and produced results comparable with OC. However, RM can result in genetic adaptation of the population to the zoo environment, in particular when 20% or less males are selected for rotation and selection of animals is not random. Lowest rates of inbreeding were obtained by combining OC with rotation of males as in the RM scheme. RM is easy to implement in practice and does not require pedigree data. When full pedigree is available, OC management is preferable.

A genetic linkage map of sole (Solea solea): a tool for evolutionary and comparative analyses of exploited (flat)fishes

Linkage maps based on markers derived from genes are essential evolutionary tools for commercial marine fish to help identify genomic regions associated with complex traits and subject to selective forces at play during exploitation or selective breeding. Additionally, they allow the use of genomic information from other related species for which more detailed information is available. Sole (solea solea L.) is a commercially important flatfish species in the North Sea, subject to overexploitation and showing evidence of fisheries-induced evolutionary changes in growth- and maturation-related traits. Sole would definitely benefit from a linkage map to better understand how evolution has shaped its genome structure. This study presents a linkage map of sole based on 423 single nucleotide polymorphisms derived from expressed sequence tags and 8 neutral microsatellite markers. The total map length is 1233.8 cM and consists of 38 linkage groups with a size varying between 0 to 92.1 cM. Being derived from expressed sequence tags allowed us to align the map with the genome of four model fish species, namely medaka (Oryzias latipes), Nile tilapia (Oreochromis niloticus), three-spined stickleback (Gasterosteus aculeatus) and green spotted pufferfish (Tetraodon nigroviridis). This comparison revealed multiple conserved syntenic regions with all four species, and suggested that the linkage groups represent 21 putative sole chromosomes. The map was also compared to the linkage map of turbot (Scophthalmus maximus), another commercially important flatfish species and closely related to sole. For all putative sole chromosomes (except one) a turbot homolog was detected, confirming the even higher degree of synteny between these two flatfish species.

Identifying environmental variables explaining genotype-by-environment interaction for body weight of rainbow trout (Onchorynchus mykiss): reaction norm and factor analytic models

BACKGROUND

Identifying the relevant environmental variables that cause GxE interaction is often difficult when they cannot be experimentally manipulated. Two statistical approaches can be applied to address this question. When data on candidate environmental variables are available, GxE interaction can be quantified as a function of specific environmental variables using a reaction norm model. Alternatively, a factor analytic model can be used to identify the latent common factor that explains GxE interaction. This factor can be correlated with known environmental variables to identify those that are relevant. Previously, we reported a significant GxE interaction for body weight at harvest in rainbow trout reared on three continents. Here we explore their possible causes.

METHODS

Reaction norm and factor analytic models were used to identify which environmental variables (age at harvest, water temperature, oxygen, and photoperiod) may have caused the observed GxE interaction. Data on body weight at harvest was recorded on 8976 offspring reared in various locations: (1) a breeding environment in the USA (nucleus), (2) a recirculating aquaculture system in the Freshwater Institute in West Virginia, USA, (3) a high-altitude farm in Peru, and (4) a low-water temperature farm in Germany. Akaike and Bayesian information criteria were used to compare models.

RESULTS

The combination of days to harvest multiplied with daily temperature (Day*Degree) and photoperiod were identified by the reaction norm model as the environmental variables responsible for the GxE interaction. The latent common factor that was identified by the factor analytic model showed the highest correlation with Day*Degree. Day*Degree and photoperiod were the environmental variables that differed most between Peru and other environments. Akaike and Bayesian information criteria indicated that the factor analytical model was more parsimonious than the reaction norm model.

CONCLUSIONS

Day*Degree and photoperiod were identified as environmental variables responsible for the strong GxE interaction for body weight at harvest in rainbow trout across four environments. Both the reaction norm and the factor analytic models can help identify the environmental variables responsible for GxE interaction. A factor analytic model is preferred over a reaction norm model when limited information on differences in environmental variables between farms is available.

Genotype-by-environment interaction of growth traits in rainbow trout (Oncorhynchus mykiss): a continental scale study

Rainbow trout is a globally important fish species for aquaculture. However, fish for most farms worldwide are produced by only a few breeding companies. Selection based solely on fish performance recorded at a nucleus may lead to lower-than-expected genetic gains in other production environments when genotype-by-environment (G × E) interaction exists. The aim was to quantify the magnitude of G × E interaction of growth traits (tagging weight; BWT, harvest weight; BWH, and growth rate; TGC) measured across 4 environments, located in 3 different continents, by estimating genetic correlations between environments. A total of 100 families, of at least 25 in size, were produced from the mating 58 sires and 100 dams. In total, 13,806 offspring were reared at the nucleus (selection environment) in Washington State (NUC) and in 3 other environments: a recirculating aquaculture system in Freshwater Institute (FI), West Virginia; a high-altitude farm in Peru (PE), and a cold-water farm in Germany (GER). To account for selection bias due to selective mortality, a multitrait multienvironment animal mixed model was applied to analyze the performance data in different environments as different traits. Genetic correlation (rg) of a trait measured in different environments and rg of different traits measured in different environments were estimated. The results show that heterogeneity of additive genetic variances was mainly found for BWH measured in FI and PE. Additive genetic coefficient of variation for BWH in NUC, FI, PE, and GER were 7.63, 8.36, 8.64, and 9.75, respectively. Genetic correlations between the same trait in different environments were low, indicating strong reranking (BWT: rg = 0.15 to 0.37, BWH: rg = 0.19 to 0.48, TGC: rg = 0.31 to 0.36) across environments. The rg between BWT in NUC and BWH in both FI (0.31) and GER (0.36) were positive, which was also found between BWT in NUC and TGC in both FI (0.10) and GER (0.20). However, rg were negative between BWT in NUC and both BWH (-0.06) and TGC (-0.20) in PE. Correction for selection bias resulted in higher additive genetic variances. In conclusion, strong G × E interaction was found for BWT, BWH, and TGC. Accounting for G × E interaction in the breeding program, either by using sib information from testing stations or environment-specific breeding programs, would increase genetic gains for environments that differ significantly from NUC.

Defining desired genetic gains for rainbow trout breeding objective using analytic hierarchy process

Distributing animals from a single breeding program to a global market may not satisfy all producers, as they may differ in market objectives and farming environments. Analytic hierarchy process (AHP) is used to estimate preferences, which can be aggregated to consensus preference values using weighted goal programming (WGP). The aim of this study was to use an AHP-WGP based approach to derive desired genetic gains for rainbow trout breeding and to study whether breeding trait preferences vary depending on commercial products and farming environments. Two questionnaires were sent out. Questionnaire-A (Q-A) was distributed to 178 farmers from 5 continents and used to collect information on commercial products and farming environments. In this questionnaire, farmers were asked to rank the 6 most important traits for genetic improvement from a list of 13 traits. Questionnaire B (Q-B) was sent to all farmers who responded to Q-A (53 in total). For Q-B, preferences of the 6 traits were obtained using pairwise comparison. Preference intensity was given to quantify (in % of a trait mean; G%) the degree to which 1 trait is preferred over the other. Individual preferences, social preferences, and consensus preferences (Con-P) were estimated using AHP and WGP. Desired gains were constructed by multiplying Con-P by G%. The analysis revealed that the 6 most important traits were thermal growth coefficient (TGC), survival (Surv), feed conversion ratio (FCR), condition factor (CF), fillet percentage (FIL%), and late maturation (LMat). Ranking of traits based on average Con-P values were Surv (0.271), FCR (0.246), TGC (0.246), LMat (0.090), FIL% (0.081), and CF (0.067). Corresponding desired genetic gains (in % of trait mean) were 1.63, 1.87, 1.67, 1.29, 0.06, and 0.33%, respectively. The results from Con-P values show that trait preferences may vary for different types of commercial production or farming environments. This study demonstrated that combination of AHP and WGP can be used to derive desired gains for a breeding program and to quantify differences due to variations market demand or production environment.

Feeding behaviour, swimming activity and boldness explain variation in feed intake and growth of sole (Solea solea) reared in captivity

The major economic constraint for culturing sole (Solea solea) is its slow and variable growth. The objective was to study the relationship between feed intake/efficiency, growth, and (non-) feeding behaviour of sole. Sixteen juveniles with an average (SD) growth of 2.7 (1.9) g/kg(0.8)/d were selected on their growth during a 4-week period in which they were housed communally with 84 other fish. Selected fish were housed individually during a second 4-week period to measure individual feed intake, growth, and behaviour. Fish were hand-fed three times a day during the dark phase of the day until apparent satiation. During six different days, behaviour was recorded twice daily during 3 minutes by direct observations. Total swimming activity, frequency of burying and of escapes were recorded. At the beginning and end of the growth period, two sequential behavioural tests were performed: "Novel Environment" and "Light Avoidance". Fish housed individually still exhibited pronounced variation in feed intake (CV = 23%), growth (CV = 25%) and behavior (CV = 100%). Differences in feed intake account for 79% of the observed individual differences in growth of sole. Fish with higher variation in feed intake between days and between meals within days had significantly a lower total feed intake (r = -0.65 and r = -0.77) and growth. Active fish showed significantly higher feed intake (r = 0.66) and growth (r = 0.58). Boldness during both challenge tests was related to fast growth: (1) fish which reacted with a lower latency time to swim in a novel environment had significantly higher feed intake (r = -0.55) and growth (r = -0.66); (2) fish escaping during the light avoidance test tended to show higher feed intake (P<0.1) and had higher growth (P<0.05). In conclusion, feeding consistency, swimming activity in the tank, and boldness during behavioral tests are related to feed intake and growth of sole in captivity.

Conservation priorities for Ethiopian sheep breeds combining threat status, breed merits and contributions to genetic diversity

Prioritizing livestock breeds for conservation needs to incorporate both genetic and non-genetic aspects important for the survival of the breeds. Here, we apply a maximum-utility-strategy to prioritize 14 traditional Ethiopian sheep breeds based on their threat status, contributions to farmer livelihoods (current breed merits) and contributions to genetic diversity. Contributions of the breeds to genetic diversity were quantified using Eding's marker-estimated kinship approaches. Non-genetic aspects included threats (e.g. low population size, low preferences by farmers) and current merits (economic, ecological and cultural merits). Threat analysis identified eight of the 14 breeds as threatened. Analysis of current merits showed that sub-alpine and arid-lowland breeds contribute most to farmer livelihoods in comparison to other breeds. The highest contribution to the genetic diversity conserved was from the Simien breed. Simien showed high between-breed (low between-breed kinship = 0.04) as well as high within-breed diversity (low within-breed kinship = 0.09 and high H_E = 0.73 and allelic richness = 6.83). We combined the results on threat status, current breed merits and contributions to genetic diversity to produce a ranking of the 14 breeds for conservation purposes. Our results balance the trade-offs between conserving breeds as insurance against future uncertainties and current sustainable utilization. The ranking of breeds provides a basis for conservation strategies for Ethiopian sheep and contributes to a regional or global conservation plan.

Population structure, genetic variation and morphological diversity in indigenous sheep of Ethiopia

We investigated genetic and morphological diversity and population structure of 14 traditional sheep populations originating from four ecological zones in Ethiopia (sub-alpine, wet highland, sub-humid lowland and arid lowland). All animals (n = 672) were genotyped for 17 microsatellite markers and scored for 12 morphological characters. The sheep were initially classified as fat-tailed (11 populations), thin-tailed (one population) and fat-rumped sheep (two populations). These classifications are thought to correspond to three consecutive introduction events of sheep from the Near-East into East Africa. For the 14 populations, allelic richness ranged from 5.87 to 7.51 and expected heterozygosity (H(E)) from 0.66 to 0.75. Genetic differentiations (F(ST) values) between all pairs of populations, except between sub-alpine populations, were significantly different from zero (P < 0.001). Cluster analysis of morphological characters and a dendrogram constructed from genetic distances were broadly consistent with the classification into fat-tailed, thin-tailed and fat-rumped sheep. Bayesian cluster analysis using microsatellite markers indicated that there has been further genetic differentiation after the initial introduction of sheep into Ethiopia. Investigation of factors associated with genetic variation showed that an isolation-by-distance model, independently of other factors, explained most of the observed genetic variation. We also obtained a strong indication of adaptive divergence in morphological characters, patterns of morphological variation being highly associated with ecology even when the effect of neutral genetic divergence (F(ST)) was parcelled out in partial Mantel tests. Using a combination of F(ST) values, Bayesian clustering analysis and morphological divergence, we propose a classification of Ethiopian sheep into six breed groups and nine breeds.

Comparative genomics in cyprinids: common carp ESTs help the annotation of the zebrafish genome

Background

Automatic annotation of sequenced eukaryotic genomes integrates a combination of methodologies such as ab-initio methods and alignment of homologous genes and/or proteins. For example, annotation of the zebrafish genome within Ensembl relies heavily on available cDNA and protein sequences from two distantly related fish species and other vertebrates that have diverged several hundred million years ago. The scarcity of genomic information from other cyprinids provides the impetus to leverage EST collections to understand gene structures in this diverse teleost group.

Results

We have generated 6,050 ESTs from the differentiating testis of common carp (Cyprinus carpio) and clustered them with 9,303 non-gonadal ESTs from CarpBase as well as 1,317 ESTs and 652 common carp mRNAs from GenBank. Over 28% of the resulting 8,663 unique transcripts are exclusively testis-derived ESTs. Moreover, 974 of these transcripts did not match any sequence in the zebrafish or fathead minnow EST collection.

A total of 1,843 unique common carp sequences could be stringently mapped to the zebrafish genome (version 5), of which 1,752 matched coding sequences of zebrafish genes with or without potential splice variants. We show that 91 common carp transcripts map to intergenic and intronic regions on the zebrafish genome assembly and regions annotated with non-teleost sequences. Interestingly, an additional 42 common carp transcripts indicate the potential presence of new splicing variants not found in zebrafish databases so far. The fact that common carp transcripts help the identification or confirmation of these coding regions in zebrafish exemplifies the usefulness of sequences from closely related species for the annotation of model genomes.

We also demonstrate that 5' UTR sequences of common carp and zebrafish orthologs share a significant level of similarity based on preservation of motif arrangements for as many as 10 ab-initio motifs.

Conclusion

Our data show that there is sufficient homology between the transcribed sequences of common carp and zebrafish to warrant an even deeper cyprinid transcriptome comparison. On the other hand, the comparative analysis illustrates the value in utilizing partially sequenced transcriptomes to understand gene structure in this diverse teleost group. We highlight the need for integrated resources to leverage the wealth of fragmented genomic data.

Genetic characterization of four strains of Nile tilapia (Oreochromis niloticus L.) using microsatellite markers

Four domesticated strains of Nile tilapia (Oreochromis niloticus L.) were genetically characterized using 14 microsatellite markers and 64 animals per strain. Two strains, Chitralada (AIT) and International Development Research Centers (IDRC) were obtained from the AIT institute, Bangkok, Thailand. The GIFT strain (5th generation) came from NAGRI, Thailand, and the GOTT strain was supplied by the University of Göttingen, Germany. The average numbers of alleles per marker were 5.0 (GOTT), 5.4 (AIT), 5.6 (IDRC) and 7.5 (GIFT). Private alleles were found at all markers with the exception of two. No fixation of alleles was found at any marker. Population differentiation, FST, was 0.178 (great genetic differentiation) and confirmed grouping of the animals in strains. The expected level of heterozygosity ranged from 0.624 to 0.711, but the observed level of heterozygosity significantly deviated from the expected level in three strains. This was probably because of small population size. Moderate to great genetic differentiation was found between strains. A phylogenetic tree reflected the strains known histories. Application of the Weitzman approach showed that all strains have added value for the total genetic diversity and thus should be retained.

Corticosteroids affect the testicular androgen production in male common carp (Cyprinus carpio L.)

Our previous experiments to study the effect of stress adaptation on pubertal development in carp showed that repeated temperature stress and prolonged feeding with cortisol-containing food pellets, which mimics the endocrine stress effects, retarded the first waves of spermatogenesis and decreased 11-ketotestosterone (11KT) plasma levels. The objective of the present study was to investigate whether the decrease in plasma 11KT is caused by a direct effect of cortisol on the steroid-producing capacity of the testis or by an indirect effect, such as a decrease in plasma LH. Pubertal and adolescent isogenic male common carp (Cyprinus carpio L.) were fed with either cortisol-containing food pellets or control food pellets over a prolonged period. Our results indicate that cortisol has a direct inhibitory effect on the testicular androgen secretion independent of the LH secretion. Furthermore, the pubertal period is critical to the influence of cortisol regarding testicular androgen secretion, because the effect is no longer observed at adolescence.

In vitro vitellogenin production by carp (Cyprinus carpio) hepatocytes as a screening method for determining (anti)estrogenic activity of xenobiotics

The yolk protein precursor vitellogenin (Vtg) is secreted by the liver of female as well as male fish, in response to estrogenic compounds. In this study, an in vitro assay was developed for measuring Vtg induction, using cultured primary hepatocytes from genetically uniform strains of carp (Cyprinus carpio). Vtg production was measured by indirect competitive ELISA, using a polyclonal antiserum against goldfish Vtg that cross-reacts with carp Vtg. Vtg was dose-dependently induced by 17beta-estradiol (E2) in hepatocytes of both sexes. E2 had a lowest observed effect concentration (LOEC) for Vtg induction of 2 nM, an EC50 between 50 and 150 nM, and a maximum response at 2 microM. The plasticizer and xenoestrogen bisphenol-A induced Vtg secretion by hepatocytes of both sexes at 50 and 100 microM. This carp hepatocyte (CARP-HEP) assay can also be used to detect antiestrogenic activity, which was measured as the reduction of E2-stimulated Vtg synthesis. Two well-known antiestrogenic compounds, tamoxifen and 2,3,7, 8-tetrachlorodibenzo-p-dioxin (TCDD), were tested. TCDD caused a reduction in Vtg synthesis in female hepatocytes at concentrations <0.1 nM, making it approximately 10,000-fold more potent than tamoxifen. Carp hepatocytes were also sensitive to induction of cytochrome P4501A (CYP1A) activity, measured as ethoxyresorufin O-deethylase (EROD). Depending on the exposure time, 18 or 96 h, EROD EC50 values for TCDD were 27 or 6 pM, respectively. The CARP-HEP assay, using the 96-well plate format, offers good possibilities to screen large numbers of compounds for (anti)estrogenic properties. In addition, it can simultaneously determine aryl hydrocarbon receptor agonist properties, measured as CYP1A induction.

Precocious onset of spermatogenesis in juvenile carp (Cyprinus carpio L., Teleostei) following treatment with low doses of L-thyroxine

Two strains of carp (Cyprinus carpio L.), a heterogenic and an all male isogeneic, were used as experimental models to assess the effects of low doses of L-thyroxine (T4) on testicular development in juvenile fish. In four experiments (two with each strain) precocious onset of spermatogenesis was induced following exposure to 15 micrograms T4 l-1, but not to 8 micrograms or 30 micrograms T4 l-1. No significant increase in body weight and testis diameter (week 16 excepted) or numbers of primary spermatogonia was observed. Total T4, measured in the isogeneic strain, increased significantly in the serum of experimental specimens within two weeks after the onset of exposure, evidencing a readily absorption of the exogeneous T4. The contents of free T4 and total triiodothyronine (T3) in the serum were measured too, showing different levels in the two experiments, as also was observed with total T4. The content of free T3, which is formed in the cells of peripheral tissues by deiodination of T4 and which supposedly exerted the observed effect on spermatogenesis, could not be determined by the used standard Amerlite assays. It is hypothesized that T3 enhances in males the activity of gonadotropic hormones, which in turn will stimulate steroid production, similar to findings in females.