Pathogenesis and shedding of Usutu virus in juvenile chickens

Usutu virus (USUV; family: Flaviviridae, genus: Flavivirus), is an emerging zoonotic arbovirus that causes severe neuroinvasive disease in humans and has been implicated in the loss of breeding bird populations in Europe. USUV is maintained in an enzootic cycle between ornithophilic mosquitos and wild birds. As a member of the Japanese encephalitis serocomplex, USUV is closely related to West Nile virus (WNV) and St. Louis encephalitis virus (SLEV), both neuroinvasive arboviruses endemic in wild bird populations in the United States. An avian model for USUV is essential to understanding zoonotic transmission. Here we describe the first avian models of USUV infection with the development of viremia. Juvenile commercial ISA Brown chickens were susceptible to infection by multiple USUV strains with evidence of cardiac lesions. Juvenile chickens from two chicken lines selected for high (HAS) or low (LAS) antibody production against sheep red blood cells showed markedly different responses to USUV infection. Morbidity and mortality were observed in the LAS chickens, but not HAS chickens. LAS chickens had significantly higher viral titers in blood and other tissues, as well as oral secretions, and significantly lower development of neutralizing antibody responses compared to HAS chickens. Mathematical modelling of virus-host interactions showed that the viral clearance rate is a stronger mitigating factor for USUV viremia than neutralizing antibody response in this avian model. These chicken models provide a tool for further understanding USUV pathogenesis in birds and evaluating transmission dynamics between avian hosts and mosquito vectors.This work was supported by National Science Foundation [grant number: 1813011]; Virginia Polytechnic Institute and State University [grant number: Data and Decisions Proposal, VMCVM Internal Research Competition].S

Mitogenomic analysis of a 50-generation chicken pedigree reveals a rapid rate of mitochondrial evolution and evidence for paternal mtDNA inheritance

Mitochondrial genomes represent a valuable source of data for evolutionary research, but studies of their short-term evolution have typically been limited to invertebrates, humans and laboratory organisms. Here we present a detailed study of 12 mitochondrial genomes that span a total of 385 transmissions in a well-documented 50-generation pedigree in which two lineages of chickens were selected for low and high juvenile body weight. These data allowed us to test the hypothesis of time-dependent evolutionary rates and the assumption of strict maternal mitochondrial transmission, and to investigate the role of mitochondrial mutations in determining phenotype. The identification of a non-synonymous mutation in ND4L and a synonymous mutation in CYTB, both novel mutations in Gallus, allowed us to estimate a molecular rate of 3.13 × 10(-7) mutations/site/year (95% confidence interval 3.75 × 10(-8)-1.12 × 10(-6)). This is substantially higher than avian rate estimates based upon fossil calibrations. Ascertaining which of the two novel mutations was present in an additional 49 individuals also revealed an instance of paternal inheritance of mtDNA. Lastly, an association analysis demonstrated that neither of the point mutations was strongly associated with the phenotypic differences between the two selection lines. Together, these observations reveal the highly dynamic nature of mitochondrial evolution over short time periods

Gut Microbiota Co-microevolution with Selection for Host Humoral Immunity

To explore coevolution between the gut microbiota and the humoral immune system of the host, we used chickens as the model organism. The host populations were two lines (HAS and LAS) developed from a common founder that had undergone 40 generations of divergent selection for antibody titers to sheep red blood cells (SRBC) and two relaxed sublines (HAR and LAR). Analysis revealed that microevolution of host humoral immunity contributed to the composition of gut microbiota at the taxa level. Relaxing selection enriched some microorganisms whose functions were opposite to host immunity. Particularly, Ruminococcaceae and Oscillospira enriched in high antibody relaxed (HAR) and contributed to reduction in antibody response, while Lactobacillus increased in low antibody relaxed (LAR) and elevated the antibody response. Microbial functional analysis showed that alterations were involved in pathways relating to the immune system and infectious diseases. Our findings demonstrated co-microevolution relationships of host-microbiota and that gut microorganisms influenced host immunity

Bidirectional Selection for Body Weight on Standing Genetic Variation in a Chicken Model

Experimental populations of model organisms provide valuable opportunities to unravel the genomic impact of selection in a controlled system. The Virginia body weight chicken lines represent a unique resource to investigate signatures of selection in a system where long-term, single-trait, bidirectional selection has been carried out for more than 60 generations. At 55 generations of divergent selection, earlier analyses of pooled genome resequencing data from these lines revealed that 14.2% of the genome showed extreme differentiation between the selected lines, contained within 395 genomic regions. Here, we report more detailed analyses of these data exploring the regions displaying within- and between-line genomic signatures of the bidirectional selection applied in these lines. Despite the strict selection regime for opposite extremes in body weight, this did not result in opposite genomic signatures between the lines. The lines often displayed a duality of the sweep signatures, where an extended region of homozygosity in one line, in contrast to mosaic pattern of heterozygosity in the other line. These haplotype mosaics consisted of short, distinct haploblocks of variable between-line divergence, likely the results of a complex demographic history involving bottlenecks, introgressions and moderate inbreeding. We demonstrate this using the example of complex haplotype mosaicism in the growth1 QTL. These mosaics represent the standing genetic variation available at the onset of selection in the founder population. Selection on standing genetic variation can thus result in different signatures depending on the intensity and direction of selection

Genomic signatures of 60 years of bidirectional selection for 8-week body weight in chickens

Sixty years, constituting 60 generations, have passed since the founding of the Virginia body weight lines, an experimental population of White Plymouth Rock chickens. Using a stringent breeding scheme for divergent 8-week body weight, the lines, which originated from a common founder population, have responded to bidirectional selection with an approximate 15-fold difference in the selected trait. They provide a model system to study the genetics of complex traits in general and the influences of artificial selection on quantitative genetic architectures in particular. As we reflect on the 60th anniversary of the initiation of the Virginia body weight lines, there is opportunity to discuss the findings obtained using different analytical and experimental genetic and genomic strategies and integrate them with a recent pooled genome resequencing dataset. Hundreds of regions across the genome show differentiation between the 2 lines, reinforcing previous findings that response to selection relied on standing variation across many genes and giving insights into the haplotype complexity underlying regions associated with body weight

Haplotype Purging after Relaxation of Selection in Lines of Chickens That Had Undergone Long-Term Selection for High and Low Body Weight

Bi-directional selection for increased and decreased 56-day body weights (BW56) has been applied to two lines of White Plymouth Rock chickens-the Virginia high (HWS) and low (LWS) body weight lines. Correlated responses have been observed, including negative effects on traits related to fitness. Here, we use high and low body weight as proxies for fitness. On a genome-wide level, relaxed lines (HWR, LWR) bred from HWS and LWS purged some genetic variants in the selected lines. Whole-genome re-sequencing was here used to identify individual loci where alleles that accumulated during directional selection were purged when selection was relaxed. In total, 11 loci with significant purging signals were identified, five in the low (LW) and six in the high (HW) body weight lineages. Associations between purged haplotypes in these loci and BW56 were tested in an advanced intercross line (AIL). Two loci with purging signals and haplotype associations to BW56 are particularly interesting for further functional characterization, one locus on chromosome 6 in the LW covering the sour-taste receptor genePKD2L1, a functional candidate gene for the decreased appetite observed in the LWS and a locus on chromosome 20 in the HW containing a skeletal muscle hypertrophy gene,DNTTIP1

Imputation-Based Fine-Mapping Suggests That Most QTL in an Outbred Chicken Advanced Intercross Body Weight Line Are Due to Multiple, Linked Loci

The Virginia chicken lines have been divergently selected for juvenile body weight for more than 50 generations. Today, the high- and low-weight lines show a >12-fold difference for the selected trait, 56-d body weight. These lines provide unique opportunities to study the genetic architecture of long-term, single-trait selection. Previously, several quantitative trait loci (QTL) contributing to weight differences between the lines were mapped in an F2-cross between them, and these were later replicated and fine-mapped in a nine-generation advanced intercross of them. Here, we explore the possibility to further increase the fine-mapping resolution of these QTL via a pedigree-based imputation strategy that aims to better capture the genetic diversity in the divergently selected, but outbred, founder lines. The founders of the intercross were high-density genotyped, and then pedigree-based imputation was used to assign genotypes throughout the pedigree. Imputation increased the marker density 20-fold in the selected QTL, providing 6911 markers for the subsequent analysis. Both single-marker association and multi-marker backward-elimination analyses were used to explore regions associated with 56-d body weight. The approach revealed several statistically and population structure independent associations and increased the mapping resolution. Further, most QTL were also found to contain multiple independent associations to markers that were not fixed in the founder populations, implying a complex underlying architecture due to the combined effects of multiple, linked loci perhaps located on independent haplotypes that still segregate in the selected lines

Genotyping by low-coverage whole-genome sequencing in intercross pedigrees from outbred founders : a cost-efficient approach

Background Experimental intercrosses between outbred founder populations are powerful resources for mapping loci that contribute to complex traits i.e. quantitative trait loci (QTL). Here, we present an approach and its accompanying software for high-resolution reconstruction of founder mosaic genotypes in the intercross offspring from such populations using whole-genome high-coverage sequence data on founder individuals (similar to 30x) and very low-coverage sequence data on intercross individuals (< 0.5x). Sets of founder-line informative markers were selected for each full-sib family and used to infer the founder mosaic genotypes of the intercross individuals. The application of this approach and the quality of the estimated genome-wide genotypes are illustrated in a large F-2 pedigree between two divergently selected lines of chickens. Results We describe how we obtained whole-genome genotype data for hundreds of individuals in a cost- and time-efficient manner by using a Tn5-based library preparation protocol and an imputation algorithm that was optimized for this application. In total, 7.6 million markers segregated in this pedigree and, within each full-sib family, between 10.0 and 13.7% of these were fully informative, i.e. fixed for alternative alleles in the founders from the divergent lines, and were used for reconstruction of the offspring mosaic genotypes. The genotypes that were estimated based on the low-coverage sequence data were highly consistent (> 95% agreement) with those obtained using individual single nucleotide polymorphism (SNP) genotyping. The estimated resolution of the inferred recombination breakpoints was relatively high, with 50% of them being defined on regions shorter than 10 kb. Conclusions A method and software for inferring founder mosaic genotypes in intercross offspring from low-coverage whole-genome sequencing in pedigrees from heterozygous founders are described. They provide high-quality, high-resolution genotypes in a time- and cost-efficient manner. The software is freely available at https://github.com/CarlborgGenomics/Stripes.

Glucose Tolerance and Plasma Non-Esterified Fatty Acid Levels in Chickens Selected for Low Body Weight, Red Junglefowl, and their Reciprocal Cross

Responses of an individual to food deprivation, such as a 16-h fast, are complex, and are influenced by environmental and genetic factors. Domestication is an ongoing process during which adaptations to changing environments occur over generations. Food deprivation by their caretakers is less for domestic chickens than for their junglefowl ancestors. Unlike domestic chicken, the junglefowl adapted over generations to periods of food deprivation, which may be reflected in differences in metabolic responses to brief periods without food. Here, we compared the blood glucose and plasma levels of non-esterified fatty acids (NEFA) among four populations when deprived of feed for 16 h. The four populations included a domestic White Rock experimental line (LWS) maintained for generations under ad libitum feeding, adult red junglefowl (RJF), and a reciprocal cross of the lines. Although there were significant differences in adult (31-week) body weight between the RJF (683 g) and LWS (1282 g), with the weight of F-1 crosses being intermediate, the amount of abdominal fat relative to body weight was similar for all populations. Patterns for blood glucose responses to a glucose bolus after a 16-h fast were similar for the initial and final points in the parental and cross populations. However, RJF reached their peak faster than LWS, with the reciprocal cross intermediate to the parental populations. Plasma NEFA concentrations were higher after the 16-h fast than in fed states, with no population differences for the fasting state. However, in the fed state, NEFA levels were lesser for LWS than for others, which was reflected further in percentage change from fed to fasted. This larger change in LWS suggests differences in mobilization of energy substrates and implies that during domestication or development of the LWS line, thresholds for responses to acute stressors may have increased