Predictable and Divergent Change in the Multivariate P Matrix during Parallel Adaptation.

AbstractAdaptation to replicated environmental conditions can be remarkably predictable, suggesting that parallel evolution may be a common feature of adaptive radiation. An open question, however, is how phenotypic variation itself evolves during repeated adaptation. Here, we use a dataset of morphological measurements from 35 populations of threespine stickleback, consisting of 16 parapatric lake-stream pairs and three marine populations, to understand how phenotypic variation has evolved during transitions from marine to freshwater environments and during subsequent diversification across the lake-stream boundary. We find statistical support for divergent phenotypic covariance (P) across populations, with most diversification of P occurring among freshwater populations. Despite a close correspondence between within-population phenotypic variation and among-population divergence, we find that variation in P is unrelated to total variation in population means across the set of populations. For lake-stream pairs, we find that theoretical predictions for microevolutionary change can explain more than 30% of divergence in P matrices across the habitat boundary. Together, our results indicate that divergence in variance structure occurs primarily in dimensions of trait space with low phenotypic integration, correlated with disparate lake and stream environments. Our findings illustrate how conserved and divergent features of multivariate variation can underlie adaptive radiation.

Genomic Diversity of the Retinta Breed Derived from Two Ancestral Bovine Lineages.

The Retinta breed, an autochthonous type of Spanish beef cattle, is highly adapted to breeding in its natural environment, which is characterized by a Mediterranean climate. The origins of this breed can be traced to two ancestral bovine stocks, which gave rise to distinct morphotypes differentiated primarily by coat color, alongside other significant traits such as growth, morphological conformation and temperament. Specifically, one morphotype comprises blond animals (Rubia Andaluza), genetically resembling the ancestral Bos taurus Aquitanicus, while the other encompasses brown- and red-colored animals (Retinta Extremeña) originating from Bos taurus Turdenatus stock. Over decades, these populations have undergone hybridization, leading to a unified population, albeit with the original subpopulations largely maintaining their genetic integrity. The objective of this study was to undertake genealogical and genomic characterization of these genetic lines, including a particular subpopulation within the blond animals (Tamarona cow). To achieve this, the genealogical records of 22,004 active animals were analyzed, and over 63,000 SNPs from a total of 1030 animals were examined for genomic characterization. Genealogical analysis revealed pedigree completeness and a high level of effective population size (Ne) across the entire population, yet relatively low Ne values within each pure line (ranging from 28.38 to 34.47). These findings underscore the ongoing efforts of the National Association of Retinta Breeders (ACRE) over the past decades to mitigate the loss of variability in this breed. The genomic characterization highlights the persistent differences within the original population and the predominant influence of the Retinto line within the current breed, as evidenced by principal component analysis (PCA) and admixture analysis. Furthermore, the identification of the Tamarona subpopulation within the blond lineage underlines its unique genetic composition, warranting its recognition as an official genetic line within the current Retinta breed. Given the small population size of these lines, particularly the Tamarona subpopulation, protective measures are imperative to preserve this distinct gene pool. Such measures would enhance the genetic diversity of the Retinta breed, which is essential for sustainable breeding practices in its natural habitats.

Laboratory-adapted and wild-type black soldier flies express differential plasticity in bioconversion and nutrition when reared on urban food waste streams.

BACKGROUND: The black soldier fly (BSF) offers a potential solution to address shortages of feed and food sources; however, selecting effective rearing substrates remains a major hurdle in BSF farming. In an urban area like Singapore, current practice is based on rearing BSF on homogeneous waste streams (e.g., spent brewery grains or okara) because heterogeneous food wastes (e.g., mixed kitchen/canteen waste or surplus cooked food) present several operational challenges with respect to the standardization of development, nutritional content, and harvesting. RESULTS: In this study, we compared two genetic strains of BSF larvae (wild-type and laboratory-adapted line) in a bioconversion experiment with diverse types of food waste (homogeneous/heterogeneous; plant/meat) and we quantified the phenotypic plasticity. Our results demonstrate different plasticity in bioconversion performance, larval growth and larval nutrition between the two BSF lines. This difference may be attributed to the selective breeding the laboratory-adapted line has experienced. Notably, larval lipid content displayed little to no genetic variation for plasticity compared with larval protein and carbohydrate content. Despite variation in larval development, heterogeneous food wastes can produce better performance in bioconversion, larval growth, and larval nutrient content than homogeneous food waste. All-meat diets result in high larvae mortality but larval survival could be rescued by mixing meat with plant-based food wastes. CONCLUSION: Overall, we suggest using mixed meals for BSF larvae feeding. Targeted breeding may be a promising strategy for the BSF industry but it is important to consider the selection effects on plasticity in larval nutrition carefully. © 2023 The Authors. Journal of The Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Classifying Genetic Lines in Pork Production by Ileal Crude Protein and Amino Acid Digestibility in Growing Pigs.

The first aim of the study was to evaluate the effect of different dietary lysine (LYS) to energy (DE) ratios on the apparent ileal digestibility (AID) of crude protein (CP) and selected amino acids (AA) in growing pigs (40-60 kg) of different genotypes. The second aim was to classify genotypes into groups based on the AID of CP and AAs. The trials were conducted on a total of 90 cross-bred barrows (30 animals/genotype) in two replicates. Before the trial series, the experimental animals (average initial body weight (BW) = 40.9 ± 8.5 kg) were surgically fitted with post valve T-cannula (PVTC). The diets were formulated with six different total LYS/DE ratios. Titanium dioxide (TiO2) was added to the diets (5 g/kg) as an indigestible marker. Based on our results, it can be concluded that the LYS/DE ratio of the diets affected the AID of the CP and AA in different ways by each genotype (p < 0.05). It can also be concluded that pigs of different genetic potential can be classified with a high accuracy (91.7%) in respect of their CP and AA digestive capacity. Our results indicate the development of genetic-profile-based swine nutrition technologies as a future direction.

Morphostructural Differences between the Historical Genetic Lines of the Spanish Merino Sheep.

The Merino breed, which originates from Spain, is the most emblematic livestock breed in the world, since it is the first with a worldwide extension and has had an important impact on the genetic origin of several of the main current sheep populations. For this reason, it is of vital importance to typify the historical genetic lines of the original Spanish Merino breed and thereby ensure the conservation of its variability. In the present study, we used 337 purebred animals (males and females) registered in the Genealogical Book of the Native Merino Breed. All the animals were descendants of herds from six ancestral genetic lines (Maesso, Egea, Granda, López-Montenegro, Hidalgo, and Donoso). Significant differences were found in all the morphometric traits and indexes between the different genetic lines. Using discriminant analysis, 84% of the animals were classified correctly into their historical genetic lines. Furthermore, the distances between the lines, calculated by a cluster test, showed that Hidalgo, Maesso, and Donoso had the most clearly defined lines, while the Granda, López-Montenegro, and Egea lines were more similar to each other. All this demonstrates the rich genetic variability existing in the genuine gene pool of the Merino sheep breed.

Dynamics of SARS-CoV-2 Major Genetic Lineages in Moscow in the Context of Vaccine Prophylaxis.

Findings collected over two and a half years of the COVID-19 pandemic demonstrated that the level immunity resulting from vaccination and infection is insufficient to stop the circulation of new genetic variants. The short-term decline in morbidity was followed by a steady increase. The early identification of new genetic lineages that will require vaccine adaptation in the future is an important research target. In this study, we summarised data on the variability of genetic line composition throughout the COVID-19 pandemic in Moscow, Russia, and evaluated the virological and epidemiological features of dominant variants in the context of selected vaccine prophylaxes. The prevalence of the Omicron variant highlighted the low effectiveness of the existing immune layer in preventing infection, which points to the necessity of optimising the antigens used in vaccines in Moscow. Logistic growth curves showing the rate at which the new variant displaces the previously dominant variants may serve as early indicators for selecting candidates for updated vaccines, along with estimates of efficacy, reduced viral neutralising activity against the new strains, and viral load in previously vaccinated patients.

Prevalence of epitheliocystis in freshwater Atlantic salmon reared in flow-through and recirculation aquaculture systems.

Epitheliocystis, an intracellular bacterial infection in the gills and skin epithelium, has been frequently reported in Atlantic salmon (Salmo salar) during freshwater production in a number of countries. This study describes the prevalence and intensity of a natural epitheliocystis infection present in the gills of two strains of Atlantic salmon reared in either a flow-through (FT) or a recirculation aquaculture system (RAS) in Ireland. Repeated sampling of gills prior to and throughout seawater transfer, histology and quantitative real-time PCR were used to determine infection prevalence and intensity. Despite no clinical gill disease, and minor histopathological changes, epitheliocystis lesions were identified in histology at all time points. Specific PCR confirmed the presence of Candidatus Clavichlamydia salmonicola in both strains and its number of copies was correlated with intensity of epitheliocystis lesions. A significant interaction between hatchery system and fish strain on the prevalence and intensity of gill epitheliocystis was found both using histological and molecular methods. Specifically, fish from FT had higher prevalence and intensity than RAS reared fish and within FT, the Irish cohort were more affected than Icelandic.

Genomic Population Structure of the Main Historical Genetic Lines of Spanish Merino Sheep.

According to historiographical documentation, the Romans first began to select Merino sheep in the Iberian Peninsula during the first century, with the aim of obtaining a breed appreciated for the quality of its wool. This process continued locally during the Middle Ages, when Spanish sheep were protected, and their export to foreign countries was banned. It was during the 16th century when individual Merino sheep were allowed to spread around the world to be used to improve the wool quality of local breeds. However, the wool crisis of the 1960s shifted the selection criteria of the Merino breed towards meat production at the expenses of wool. Consequently, individuals that display the genetic and phenotypic characteristics of those sheep originally bred in the kingdom of Spain in the Middle Ages are extremely difficult to find in commercial herds. In this study, we characterized the genetic basis of 403 individuals from the main historical Spanish Merino genetic lines (Granda, Hidalgo, Lopez-Montenegro, Maeso, Donoso and Egea), which were bred in isolation over the last 200 years, using a genomic approach based on genotyping data from the Axiom™ Ovine 50K SNP Genotyping Array. Our analysis included measuring population structure, genomic differentiation indexes, runs of homozygosity (ROH) patterns, and an analysis of molecular variance (AMOVA). The results showed large genetic differences between the historical lines, even though they belong to the same breed. In addition, ROH analysis showed differences due to increased inbreeding among the ancient generations compared with the modern Merino lines, confirming the breed's ancestral and closed origin. However, our results also showed a high variability and richness within the Spanish historical Merino lines from a genetic viewpoint. This fact, together with their great ability to produce high-quality wool, suggests that ancestral Merino lines from Spain should be considered a valuable genetic population to be maintained as a resource for the improvement of wool-producing sheep breeds all around the world.

Key parameters of head-heart electrical stunning need to be adapted to improve stunning effectiveness and meat quality in pigs of different genetic lines.

Adjustment of key parameters of electrical stunning of pigs is necessary to better respect animal welfare and avoid impairments of meat quality. In the present study, we analysed the effect of genetic line, stunning program and genetic line × stunning program interactions on stunning effectiveness, carcass and meat quality traits in 10,338 pigs from five different genetic lines in a German slaughterhouse using electrical stunning (two-cycle method head/heart current). Genetic line, stunning program and its interaction significantly influenced stunning effectiveness and pH in the semimembranosus muscle. Backup stunning rates ranged from 1.59 to 5.01% in the five genetic lines, and were significantly decreased (2.85%) in program 1 using a current of 1.3 to 1.35 A, 800 Hz and 1 s for the head-cycle compared to program 2 (4.73%) with 1.8 A, 100 Hz and 4.5 s. We suggest to consider genetic line × stunning program interactions when adjusting officially defined stunning parameters for head-heart electrical stunning.

Performance and processing yield comparisons of Large White male turkeys by genetic lines, sources, and seasonal rearing.

Large White male turkey genetic lines (GL) comparison in performance and processing yields under the same conditions are rare in the literature. Two rearing experiments (EXP) were conducted to accomplish 2 objectives. The first objective was to test the effects of poult source and genetic lines on performance and processing yields. The second objective was to extract season and growth patterns when comparing both EXP common treatments. In EXP 1, male poults from 5 different sources were randomly assigned to 48 concrete: litter-covered floor pens. In EXP 2, male poults from 7 different genetic lines were randomly assigned to 48 concrete: litter-covered floor pens. For both EXP, the experimental design was a completely randomized block design with a one-factor arrangement. Both EXP were placed in the same house with the same management and nutrition in two separate seasons of the same year. Bird performance and carcass processing yield were analyzed in SAS 9.4 or JMP 15.1 in a mixed model. In EXP 1 no significant difference in BW or processing yield was observed. However, a similar GL from a commercial hatchery had an improved feed conversion ratio (FCR) over the same GL sourced directly from the genetic company hatchery. In EXP 2, statistical differences were observed in performance and breast meat yield depending on the GL. A season effect was observed when comparing the two EXP. Birds raised in the fall season had a 2 kg BW increase, on average, over their spring counterparts. This difference in BW can also be observed in a statistically higher breast meat yield by the birds raised in the fall over the ones raised in the spring. In conclusion, a comparison between GL resulted in effects due to genetic line, poult source, and rearing season on bird performance and carcass yield.

Effects of common full-sib families on accuracy of genomic prediction for tagging weight in striped catfish Pangasianodon hypophthalmus.

Common full-sib families (c 2 ) make up a substantial proportion of total phenotypic variation in traits of commercial importance in aquaculture species and omission or inclusion of the c 2 resulted in possible changes in genetic parameter estimates and re-ranking of estimated breeding values. However, the impacts of common full-sib families on accuracy of genomic prediction for commercial traits of economic importance are not well known in many species, including aquatic animals. This research explored the impacts of common full-sib families on accuracy of genomic prediction for tagging weight in a population of striped catfish comprising 11,918 fish traced back to the base population (four generations), in which 560 individuals had genotype records of 14,154 SNPs. Our single step genomic best linear unbiased prediction (ssGLBUP) showed that the accuracy of genomic prediction for tagging weight was reduced by 96.5%-130.3% when the common full-sib families were included in statistical models. The reduction in the prediction accuracy was to a smaller extent in multivariate analysis than in univariate models. Imputation of missing genotypes somewhat reduced the upward biases in the prediction accuracy for tagging weight. It is therefore suggested that genomic evaluation models for traits recorded during the early phase of growth development should account for the common full-sib families to minimise possible biases in the accuracy of genomic prediction and hence, selection response.

Differences in susceptibility to Mycobacterium chelonae in zebrafish (Danio rerio) lines commonly used in scientific research.

Mycobacteriosis is one of the most common diseases encountered in laboratory zebrafish. These infections can present a problem to researchers using zebrafish because they may introduce unknown experimental variables. Whilst differences in severity of infections between species of Mycobacterium infecting zebrafish have been well documented, little is known about differences in susceptibility between zebrafish lines. Previous surveys have found higher prevalence in the TU zebrafish line relative to other lines, suggesting that there may be underlying genetic differences in susceptibility. This study investigates Mycobacterium chelonae H1E2-GFP infections in four different zebrafish lines commonly used in research (AB, 5D, casper and TU). Fish were exposed to a labelled (green-fluorescent protein (GFP)) strain of M. chelonae by intraperitoneal injection, and infection status was evaluated after 10 weeks. Visualization of GFP in euthanized fish and histology were used as endpoints. In GFP images, severity was assessed by image analysis, and in histological sections, counts of granulomas containing acid-fast bacteria were used. Results indicated differences in severity of infections between lines, but no significant differences in prevalence.

Claw Characteristics of Culled Sows from Three Farrow-to-Finish Greek Farms. Part 1: Claw Length Measurements, Lesion Scores and Their Association.

The aim of the study was to investigate variations in lengths and lesions in claws of culled sows and to evaluate their association. All four feet of 185 sows from three Greek farrow-to-finish farms (Farm A: 57 sows; Farm B: 64 sows; Farm C: 64 sows) were examined for lesions and their lengths were measured. All claw lengths were lower in sows of farm C compared to those from sows of B and A. Claw lengths in sows of farm B were lower compared to those from A for all lateral toes of front feet and for all medial and four out of three lateral toes of rear feet. Sum of length measurements of the main toes of the front feet (SLF) associated with lesions on sole, white line and heel of front feet, while sum of length measurements of the main toes of the rear feet (SLR) associated with all lesions of the rear feet. The lengths of the main toes were correlated with the length of dew claws on front and rear feet. Overall, sows' claw lesion severity and claw lengths may differ between farms and frequency of lesions is higher in longer claws.

Genetic Parameters of Effort and Recovery in Sport Horses Assessed with Infrared Thermography.

The way a horse activates (effort phase-EP) and recovers (recovery phase-RP) during a sport event can affect its sport performance. The aim of this manuscript was to test horses' adaptation to sport performance and its genetic basis, using eye temperature assessed with infrared thermography. EP and RP were measured in 495 Spanish Sport Horses, during a performance test, considering sex (2) and genetic lines (5) as fixed effects. The ranking position obtained on an official sport competition was also collected. Differences in variables due to genetic line and sex effects were found, showing that, regardless of the genetic line, stallions tended to recover better than mares after the sport test developed. High positive intra-class correlations (p < 0.001) were found between EP and RP for both fixed effects, so that the higher the EP, the higher the RP. However, for the ranking position, a low negative correlation (p < 0.01) was found, so that the higher the eye temperature increase, the better the position. Heritabilities showed medium-high values with a medium positive genetic correlation between them. Thus, breed origins and sex influence horses' effort and recovery during sport performance, showing a genetic basis adequate for selection.

Gregarious nesting in relation to floor eggs in broiler breeders.

Gregarious nesting has often been observed in laying hens, where hens prefer to visit a nest already occupied by other hens over empty nests. This may result in overcrowding of the nests which is considered a welfare issue and, moreover, can increase the economic issue of floor eggs. This study aimed to describe gregarious nesting and spatial behavior in broiler breeders and how this relates to genetic background, fearfulness and mating behavior. Five commercially available genetic lines of broiler breeders were housed in 21 pens of 550 females and 50 males (six pens for lines 1 and 2, five pens for line 3 and two pens for lines 4 and 5) during the ages 20-60 weeks. Every 10 weeks, the plumage condition and wounds were assessed of 50 random hens per pen. Avoidance distance and novel object tests were performed to assess fearfulness at four time points. Distribution of eggs over nests was observed for 6 weeks at the onset of egg production at 26 weeks of age, and use of space was recorded at four time points, while (floor) egg production was noted daily per pen. We found differences between genetic lines over time in plumage condition and prevalence of wounds. Fear of humans was highest at the earliest age tested and did not correlate with general fearfulness as assessed by the novel object test. The distribution of eggs over nests was related to genetic background and was more uneven at the earliest age compared to later ages, and a more uneven distribution was correlated with an increased percentage of floor eggs. Distribution of birds over the litter area differed between the genetic lines, and less use of the litter area was correlated with an increased fear of humans and presence of wounds, suggesting an association with aggressive mating behavior. This difference in distribution of the birds could also explain the correlation between increased presence of wounds and decreased percentage of floor eggs. It is concluded that broiler breeders do show gregarious nesting, which is affected by genetic background. Both increased gregarious nesting and wounds are related to increased floor egg percentage, which should be studied further in broiler breeder research. Genetic selection for even use of the available nests and of the litter and slatted area would therefore support both broiler breeder welfare and performance.

Differences in Tibia Shape in Organically Reared Chicken Lines Measured by Means of Geometric Morphometrics.

In the present study, the conformation of the tibia of seven genetic lines of broilers was analyzed by Geometric Morphometrics and correlated to carcass weight and walking ability. The used chicken genetic lines were classified as fast, medium, or slow growing and ranked for their walking ability. Six chicken types were reared in an organic farm and slaughtered at 81 days of age while one slow-growing and highly walking line (Naked Neck) was reared in a commercial farm and used as external reference for moving activity and growth speed. A mixed landmarks and semi-landmarks model was applied to the study of tibia shape. Results of this study showed that: (i) body weight gain was positively correlated to the curvature of the antero-posterior axis of the tibia; (ii) the shape of the tibia and the active walking behavior were significantly correlated; (iii) walking and not-walking genetic lines could be discriminated in relation to the overall shape of the tibia; (iv) a prevalence of static behavior was correlated to a more pronounced curvature of the antero-posterior axis of the tibia. Results of this study revealed that the walking genetic types have a more functional and natural tibia conformation. This easy morphologic method for evaluating tibia shape could help to characterize the adaptability of genotypes to organic and outdoor rearing.

Contamination of eggs by Salmonella Enteritidis in experimentally infected laying hens of four commercial genetic lines in conventional cages and enriched colony housing.

Human illness caused by the consumption of eggs contaminated with Salmonella Enteritidis is a continuing international public health concern. This pathogen is deposited inside the edible contents of eggs as a consequence of its ability to colonize reproductive tissues in infected hens. Conditions in the housing environment can influence the persistence and transmission of avian Salmonella infections, but the food safety ramifications of different poultry management systems are not entirely clear. The present study assessed the deposition of S. Enteritidis inside eggs laid by groups of experimentally infected laying hens of 4 commercial genetic lines (designated as white egg lines W1 and W2 and brown egg lines B1 and B2). Groups of hens from each line were housed at 555 cm2 of floor space per bird in both conventional cages and colony units enriched with access to perches and nesting areas. All hens were orally inoculated with 5.75 × 107 cfu of a 2-strain S. Enteritidis mixture, and the internal contents of eggs laid 5 to 24 D post-inoculation were cultured to detect the pathogen. No significant differences in egg contamination frequencies were found between the 2 housing systems for any of the hen lines. Contaminated eggs were laid between 7 and 21 D post-inoculation at an overall frequency of 2.47%, ranging from 0.25 to 4.38% for the 4 hen lines. The frequency of S. Enteritidis recovery from egg samples was significantly (P < 0.05) lower for line B2 than for any of the other lines, and the egg contamination frequency for line W1 was significantly greater than for line W2. The overall incidence of contamination among white eggs (3.38%) was significantly higher than among brown eggs (1.56%). These results demonstrate that S. Enteritidis deposition inside eggs can vary between genetic lines of infected laying hens, but housing these hens in 2 different systems did not affect the production of contaminated eggs.

Colonization of internal organs by Salmonella Enteritidis in experimentally infected laying hens of four commercial genetic lines in conventional cages and enriched colony housing.

The prevalence of Salmonella Enteritidis in commercial egg-laying flocks is a prominent public health concern because contaminated eggs cause human illness. Deposition of this pathogen inside eggs results from bacterial colonization of reproductive tissues in infected hens. Environmental conditions can influence avian Salmonella infections, but the food safety consequences of different poultry housing systems remain uncertain. The present study assessed the invasion of internal organs by Salmonella Enteritidis in groups of experimentally infected laying hens of four commercial genetic lines (designated as white egg lines W1 and W2 and brown egg lines B1 and B2). Groups of hens from each line were housed at 555 cm2 of floor space per bird in both conventional cages and colony units enriched with access to perches and nesting areas. All hens were orally inoculated with 5.75 × 107 colony-forming units of a two-strain Salmonella Enteritidis mixture. At 6 to 7 d post-inoculation, hens were euthanized, and samples of liver, spleen, ovary, oviduct, and intestinal tract were removed for bacteriologic culturing. The frequency of Salmonella Enteritidis recovery from intestinal samples was significantly (P < 0.05) greater for the two white egg lines combined than for the two brown egg lines combined in both conventional cage (72.2% vs. 50.0%) and enriched colony housing systems (66.7% vs. 37.5%). The frequency of intestinal Salmonella Enteritidis isolation from line B1 was significantly higher from hens in conventional cages (47.2%) than in enriched colonies (22.2%), but no differences were observed for other hen lines. Line W1 yielded more positive intestinal samples than either brown egg line in conventional cages, and line B2 had fewer positive intestinal samples than all other lines in enriched colonies. There were no significant differences between hen lines or housing systems in Salmonella Enteritidis isolation from other internal organs. These results demonstrate that Salmonella Enteritidis colonization of the intestinal tract can vary between genetic lines of egg-laying hens and that some lines are subject to housing system influences on Salmonella susceptibility.

Comparative Response of the Hepatic Transcriptomes of Domesticated and Wild Turkey to Aflatoxin B1.

The food-borne mycotoxin aflatoxin B1 (AFB1) poses a significant risk to poultry, which are highly susceptible to its hepatotoxic effects. Domesticated turkeys (Meleagris gallopavo) are especially sensitive, whereas wild turkeys (M. g. silvestris) are more resistant. AFB1 toxicity entails bioactivation by hepatic cytochrome P450s to the electrophilic exo-AFB1-8,9-epoxide (AFBO). Domesticated turkeys lack functional hepatic GST-mediated detoxification of AFBO, and this is largely responsible for the differences in resistance between turkey types. This study was designed to characterize transcriptional changes induced in turkey livers by AFB1, and to contrast the response of domesticated (susceptible) and wild (more resistant) birds. Gene expression responses to AFB1 were examined using RNA-sequencing. Statistically significant differences in gene expression were observed among treatment groups and between turkey types. Expression analysis identified 4621 genes with significant differential expression (DE) in AFB1-treated birds compared to controls. Characterization of DE transcripts revealed genes dis-regulated in response to toxic insult with significant association of Phase I and Phase II genes and others important in cellular regulation, modulation of apoptosis, and inflammatory responses. Constitutive expression of GSTA3 was significantly higher in wild birds and was significantly higher in AFB1-treated birds when compared to controls for both genetic groups. This pattern was also observed by qRT-PCR in other wild and domesticated turkey strains. Results of this study emphasize the differential response of these genetically distinct birds, and identify genes and pathways that are differentially altered in aflatoxicosis.