Data compression can discriminate broilers by selection line, detect haplotypes, and estimate genetic potential for complex phenotypes.

Accurately establishing the relationships among individuals lays the foundation for genetic analyses such as genome-wide association studies and identification of selection signatures. Of particular interest to the poultry industry are estimates of genetic merit based on molecular data. These estimates can be commercially exploited in marker-assisted breeding programs to accelerate genetic improvement. Here, we test the utility of a new method we have recently developed to estimate animal relatedness and applied it to genetic parameter estimation in commercial broilers. Our approach is based on the concept of data compression from information theory. Using the real-world compressor gzip to estimate normalized compression distance (NCD) we have built compression-based relationship matrices (CRM) for 988 chickens from 4 commercial broiler lines-2 male and 2 female lines. For all pairs of individuals, we found a strong negative relationship between the commonly used genomic relationship matrix (GRM) and NCD. This reflects the fact that "similarity" is the inverse of "distance." The CRM explained more genetic variation than the corresponding GRM in 2 of 3 phenotypes, with corresponding improvements in accuracy of genomic-enabled predictions of breeding value. A sliding-window version of the analysis highlighted haplotype regions of the genome apparently under selection in a line-specific manner. In the male lines, we retrieved high population-specific scores for IGF-1 and a cognate receptor, INSR. For the female lines, we detected an extreme score for a region containing a reproductive hormone receptor (GNRHR). We conclude that our compression-based method is a valid approach to established relationships and identify regions under selective pressure in commercial lines of broiler chickens.

Prednisolone-induced predisposition to femoral head separation and the accompanying plasma protein changes in chickens.

Femoral head separation (FHS) is an idiopathic bone problem that causes lameness and production losses in commercial poultry. In a model of prednisolone-induced susceptibility to FHS, the changes in plasma proteins and peptides were analyzed to find possible biomarkers. Plasma samples from control and FHS-susceptible birds were depleted of their high abundance proteins by acetonitrile precipitation and were then subjected to cation exchange and reverse-phase (RP) fractionations. Analysis with matrix assisted laser desorption ionization-time-of-flight mass spectrometry (MALDI-TOF-MS) showed several differentially expressed peptides, two of which were isolated by RP-HPLC and identified as the fragments of apolipoprotein A-I. The acetonitrile fractionated plasma proteins were subjected to reduction/alkylation and trypsin digestion followed by liquid chromatography and tandem mass spectrometry, which showed the absence of protocadherin 15, vascular endothelial growth factor-C, and certain transcription and ubiquitin-mediated proteolytic factors in FHS-prone birds. It appears that prednisolone-induced dyslipidemia, vascular, and tissue adhesion problems may be consequential to FHS. Validity of these biomarkers in our model and the natural disease must be verified in future using traditional approaches. BIOMARKER INSIGHTS: Lameness because of femoral head separation (FHS) is a production and welfare problem in the poultry industry. Selection against FHS requires identification of the birds with subclinical disease with biomarkers from a source such as blood. Prednisolone can induce femoral head problems and predisposition to FHS. Using this experimental model, we analyzed the plasma peptides and proteins from normal and FHS-prone chickens by mass spectrometry to identify differentially expressed peptides and proteins. We found two peptides, both derived from apolipoprotein A-I, quantitatively elevated and two proteins, protocadherin 15 and VEGF-C, that were conspicuously absent in FHS-susceptible birds.

Associations of chicken Mx1 polymorphism with antiviral responses in avian influenza virus infected embryos and broilers.

Avian influenza virus (AIV) is a major respiratory disease of poultry that causes catastrophic losses to the poultry industry. The Mx protein has been shown to confer antiviral responses to influenza viruses in mice. One nonsynonymous substitution (S631N) in the chicken Mx protein is reported to be associated with resistance to AIV infection in vitro. The previous studies suggested controversy over whether this substitution in the Mx protein plays an important antiviral role in AIV infection in the chicken. It would be intriguing to investigate if the substitution is associated with resistance to AIV infection both in ovo and in vivo in chickens. In this study, the embryos and young chicks were generated from the cross of Mx1 heterozygous (S631N) parents with an expected segregating ratio of 1:2:1 in the progeny. A PCR length polymorphism was developed to genotype the Mx1 gene from 119 embryos and 48 chickens. The embryonated chicken eggs were inoculated with 10(6) 50% embryo infectious dose (EID(50)) H5N9 AIV on d 13. Hemagglutinating units in allantoic fluid were determined at 48 h postinoculation. For the in vivo study, twenty-four 1-wk-old broilers were inoculated with 10(6) EID(50) H5N3, and virus titers in lungs were evaluated at d 4 postinoculation. This is the first report revealing no significant association between Mx1 genotypes and low pathogenesis AIV infection both in ovo and in vivo in the chicken. Total RNA samples were isolated from chicken lung tissues in the in vivo study, and the Mx1 mRNA expression assay among 3 genotypes also suggested that only heterozygote birds had significantly greater expression with AIV infection than noninfected birds. A recombination breakpoint within Mx1 gene was also first identified, which has laid a solid foundation for further understanding biological function of the Mx1 gene in chickens. The current study provides valuable information on the effect of the Mx1 gene on the genetic resistance to AIV in chickens, and Mx1 will not be applicable for enhancing genetic resistance to AIV infection in chickens.

Effects of high fat diets or prednisolone treatment on femoral head separation in chickens.

1. The effects of high fat diets and prednisolone treatment were studied to understand the etiology of femoral head separation (FHS) in fast growing broiler chickens. Dietary effects on production parameters such as growth, feed conversion ratio (FCR) and blood chemistry were also measured. 2. Three groups of chickens, consisting of 30 birds each, in two replicate pens, were fed isonitrogenous diets containing 40 (control), 60, or 80 g poultry fat supplements per kg feed. The birds were fed a starter diet containing the fat supplements for the first three weeks, then switched to a grower diet containing the same supplements for the rest of the experimental period. Two groups of birds were also raised with the control diets, but were administered either cholesterol or prednisolone intramuscularly at 30 and 32 days of age to evaluate their effects on FHS incidences. 3. The chickens were euthanised and necropsied at 37 d of age. The presence of femoral head weakness was determined by applying mild pressure on the pelvic joint to cause the growth plate to become detached from its articular cartilage in affected cases. 4. High fat diets did not change FHS incidences, but increased 28 d body weights (BW) and FCR. At 37 d of age the BW differences were not significant but the FCR (gain: feed ratio) remained higher in high fat fed groups. Prednisolone treatment, by contrast, resulted in decreased BW, decreased feed efficiency, increased FHS index, and elevated blood lipid levels. 5. The results suggest that high dietary fats do not affect FHS incidence in broilers. Prednisolone treatment causes hyperlipidaemia and increases FHS index, and may therefore provide a suitable experimental model of FHS pathogenesis in growing chickens.

Histopathology and serum clinical chemistry evaluation of broilers with femoral head separation disorder.

Femoral head separation (FHS) and necrosis is a sporadic leg problem of unknown etiology in broiler breeders. To determine the underlying physiology of FHS, the blood chemistry and histopathology of the femoral growth plates of the affected chickens were compared with their age-matched controls and with birds having tibial dyschondroplasia. Femoral problems were categorized on the basis of 1) femoral head separating from articular cartilage without any visible damage to the growth plate (FHS) and 2) FHS with significant tearing and lesions in the growth plate (FHSL). Tibial dyschondroplasia was identified by a widening of the growth plate with an unresorbed plug of cartilage at the proximal end of the tibia. Control birds were without any femoral or tibial problems. The histopathology of FHSL growth plates revealed occasional chondrocyte death, hypocellularity, dysplasia in the prehypertrophic zones, and the absence of inflammatory infiltrates in the lesion areas. Hematoxylin and eosin staining showed brown chromogenic deposits in the metaphyseal bone marrow areas. Blood chemistry of chickens with FHSL showed a modest but significant elevation of cholesterol, triglycerides, and low-density lipoproteins. Only cholesterol and low-density lipoproteins were moderately elevated in FHS-affected chickens. Other blood parameters, such as protein, magnesium, and iron levels, showed differential changes in birds with leg problems, but there were no specific trends. Neither blood ovotransferrin, a marker of chronic inflammation, nor corticosterone, a marker of stress, showed any significant differences from the controls. These results indicate that FHS may be a metabolic problem in poultry, one that is related to fat metabolism disorders, possibly contributing to an unbalanced growth in the articular-epiphyseal complex that leads to its separation under sheer stress.

Polymorphisms in uncoupling protein, melanocortin 3 receptor, melanocortin 4 receptor, and pro-opiomelanocortin genes and association with production traits in a commercial broiler line.

Because avian uncoupling protein (avUCP), melanocortin 3 receptor (MC3R), melanocortin 4 receptor (MC4R), and pro-opiomelanocortin (POMC) genes may be associated with production traits [e.g., BW, weight gain (WG), and feed conversion ratio (FCR)], male and female broilers from an elite broiler line were screened for polymorphisms in these genes. The PCR-restriction fragment length polymorphism (RFLP) tests were developed to type the missense polymorphisms UCPAla118Val, MC4RSer76Leu, MC3R-Met54Leu, and Gly104Ser and POMCPro61Leu. Of 39 single nucleotide polymorphisms identified in all 4 genes, 24/39 were transitions with 11 having a C to T change. Of the 23 polymorphisms in UCP, 17 represented at least 7 haplotypes in this pedigreed broiler line. The UCP Ala-118Val allele was associated with a) high feed efficiency (FE; P = 0.03) and WG (P = 0.053) in selected males, and b) high BW in selected females (P = 0.07) and unselected males (P = 0.015). The UCPVal118Val allele was found in approximately 10% of the birds that were screened. Five silent substitutions, 3 in MC3R and 2 in MC4R, were also identified. Thirteen polymorphisms were identified in the POMC gene representing at least 3 different alleles. A missense Pro61Leu heterozygote was associated with greater BW in females. The heterozygote MC3R Gly104Ser polymorphism was associated with greater FE in selected males (P = 0.03) and greater BW in unselected males (P = 0.007). The MC4R Ser76Leu heterozygote polymorphism was associated with greater BW than the Leu76 homozygote in females (P = 0.05). From these findings, we hypothesize that UCP, MC3R, MC4R and POMC genes may play important roles and could be candidate loci for production traits such as feed conversion and BW in commercial broiler breeding stock.

Genetic variation at the tumour virus B locus in commercial and laboratory chicken populations assessed by a medium-throughput or a high-throughput assay.

The tumour virus B (TVB) locus encodes cellular receptors mediating infection by three subgroups of avian leukosis virus (B, D, and E). Three major alleles, TVB*S1, TVB*S3, and TVB*R, have been described. TVB*S1 encodes a cellular receptor mediating infection of subgroups B, D, and E. TVB*S3 encodes the receptor for two subgroups, B and D, and TVB*R encodes a dysfunctional receptor that does not permit infection by any of the subgroups, B, D, or E. Genetic diversity at the TVB locus of chickens was investigated in both layer and broiler commercial pure lines and laboratory lines. Genotyping assays were developed for both medium-throughput and high-throughput analysis. Of the 36 broiler lines sampled, 14 were fixed for the susceptible allele TVB*S1. Across all broiler lines, 83% of chickens were typed as TVB*S1/*S1, 3% as TVB*R/*R, and 14% as TVB*S1/*R. In the egg-layer lines, five of the 16 tested were fixed for TVB*S1/*S1. About 44% of egg-layers were typed as TVB*S1/*S1, 15% as TVB*R/*R, with the rest segregating for two or three of the alleles. In the laboratory chickens, 60% were fixed for TVB*S1/*S1, 6% for TVB*S3/*S3, 14% for TVB*R/*R, and the rest were heterozygotes (TVB*S1/*S3 or TVB*S1/*R). All commercial pure lines examined in this study carry the TVB*S1 allele that sustains the susceptibility to avian leukosis viruses B, D, and E. More importantly, the TVB*R allele was identified in multiple populations, thus upholding the opportunities for genetic improvement through selection.

Association of mitochondrial function with feed efficiency within a single genetic line of male broilers.

Studies were conducted to determine relationships between feed efficiency and mitochondrial function and biochemistry. After feed efficiency (FE; gain:feed) was determined in broiler breeder males between 6 and 7 wk of age, mitochondria were isolated from breast and leg muscle from birds with high FE (0.83+/-0.01, n = 6) and low FE (0.64+/-0.01, n = 7). Respiratory chain coupling, assessed by the respiratory control ratio (RCR), was greater in high FE breast, and leg mitochondria provided NADH-linked, but not FADH-linked, energy substrates. There were no differences, however, in the adenosine diphosphate to oxygen (ADP:O) ratio (an index of oxidative phosphorylation) when mitochondria were provided either energy substrate. Electron leak, as determined by generation of H202, was greater in the low FE than in high FE breast mitochondria. Electron leak increased following inhibition of electron transport at Complex I (with rotenone) and Complex III (with antimycin A) in low FE but not in high FE breast mitochondria. There were no differences in basal electron leak in leg mitochondria between groups, but H202 generation was elevated (P < 0.07) compared to basal values in low FE leg mitochondria after Complex I inhibition. The activities of Complexes I and II were greater in high FE breast and leg muscle mitochondria compared to those in low FE mitochondria. The results indicate that lower respiratory chain coupling in low FE muscle mitochondria may be due to lower activities of Complexes I and II and defects in electron leak and provide insight into cellular mechanisms associated with the phenotypic expression of feed efficiency in broilers.

Mitochondrial DNA of Hydra attenuata (Cnidaria): a sequence that includes an end of one linear molecule and the genes for l-rRNA, tRNA(f-Met), tRNA(Trp), COII, and ATPase8.

The 3231-nucleotide-pair (ntp) sequence of one end of one of the two linear mitochondrial (mt) DNA molecules of Hydra attenuata (phylum Cnidaria, class Hydrozoa, order Anthomedusae) has been determined. This segment contains complete genes for tRNA(f-Met), l-rRNA, tRNA(Trp), subunit 2 of cytochrome c oxidase (COII), subunit 8 of ATP synthetase (ATPase8), and the 5' 136 ntp of ATPase6. These genes are arranged in the order given and are transcribed from the same strand of the molecule. As in two other cnidarians, the hexacorallian anthozoan Metridium senile and the octocorallian anthozoan Sarcophyton glaucum, the mt-genetic code of H. attenuata is near standard. The only modification appears to be that TGA specifies tryptophan rather than termination. Also as in M. senile and S. glaucum, the encoded H. attenuata mt-tRNA(f-Met) has primary and secondary structural features resembling those of Escherichia coli initiator tRNA(t-Met). As the encoded mt-tRNA(Trp) cannot be folded into a totally orthodox secondary structure, two alternative forms are suggested. The encoded H. attenuata mt-l-rRNA is 1738 nt, which is 451 nt shorter than the M. senile mt-l-rRNA. Comparisons of secondary structure models of these two mt-l-rRNAs indicate that most of the size difference results from loss of nucleotides in the H. attenuata molecule at a minimum of 46 locations, which includes elimination of six distinct helical elements.

Mytilus mitochondrial DNA contains a functional gene for a tRNASer(UCN) with a dihydrouridine arm-replacement loop and a pseudo-tRNASer(UCN) gene.

A 2500-nucleotide pair (ntp) sequence of F-type mitochondrial (mt) DNA of the Pacific Rim mussel Mytilus californianus (class Bivalvia, phylum Mollusca) that contains two complete (ND2 and ND3) and two partial (COI and COIII) protein genes and nine tRNA genes is presented. Seven of the encoded tRNAs (Ala, Arg, His, Met(AUA), Pro, Ser(UCN), and Trp) have the potential to fold into the orthodox four-armed tRNA secondary structure, while two [tRNASer(AGN) and a second tRNASer(UCN)] will fold only into tRNAs with a dihydrouridine (DHU) arm-replacement loop. Comparison of these mt-tRNA gene sequences with previously published, corresponding M. edulis F-type mtDNA indicates that similarity between the four-armed tRNASer(UCN) genes is only 63.8% compared with an average of 92.1% (range 86.2-98. 5%) for the remaining eight tRNA genes. Northern blot analysis indicated that mature tRNAs encoded by the DHU arm-replacement loop-containing tRNASer(UCN), tRNASer(AGN), tRNAMet(AUA), tRNATrp, and tRNAPro genes occur in M. californianus mitochondria, strengthening the view that all of these genes are functional. However, Northern blot and 5' RACE (rapid amplification of cDNA ends) analyses indicated that the four-armed tRNASer(UCN) gene is transcribed into a stable RNA that includes the downstream COI sequence and is not processed into a mature tRNA. On the basis of these observations the M. californianus and M. edulis four-armed tRNASer(UCN) sequences are interpreted as pseudo-tRNASer(UCN) genes.

The mitochondrial genome of the sea anemone Metridium senile (Cnidaria): introns, a paucity of tRNA genes, and a near-standard genetic code.

The circular, 17,443 nucleotide-pair mitochondrial (mt) DNA molecule of the sea anemone, Metridium senile (class Anthozoa, phylum Cnidaria) is presented. This molecule contains genes for 13 energy pathway proteins and two ribosomal (r) RNAs but, relative to other metazoan mtDNAs, has two unique features: only two transfer RNAs (tRNA(f-Met) and tRNA(Trp)) are encoded, and the cytochrome c oxidase subunit I (COI) and NADH dehydrogenase subunit 5 (ND5) genes each include a group I intron. The COI intron encodes a putative homing endonuclease, and the ND5 intron contains the molecule's ND1 and ND3 genes. Most of the unusual characteristics of other metazoan mtDNAs are not found in M. senile mtDNA: unorthodox translation initiation codons and partial translation termination codons are absent, the use of TGA to specify tryptophan is the only genetic code modification, and both encoded tRNAs have primary and secondary structures closely resembling those of standard tRNAs. Also, with regard to size and secondary structure potential, the mt-s-rRNA and mt-1-rRNA have the least deviation from Escherichia coli 16S and 23S rRNAs of all known metazoan mt-rRNAs. These observations indicate that most of the genetic variations previously reported in metazoan mtDNAs developed after Cnidaria diverged from the common ancestral line of all other Metazoa.

DNA marker technology: a revolution in animal genetics.

The development of DNA-based markers has had a revolutionary impact on gene mapping and, more generally, on all of animal and plant genetics. With DNA-based markers, it is theoretically possible to exploit the entire diversity in DNA sequence that exists in any cross. For this reason, high resolution genetic maps are being developed at an unprecedented speed. The most commonly used DNA-based markers include those based on a cloned and (usually) sequenced DNA fragment and other, more random, assays for genetic polymorphism that can be grouped under the heading of fingerprint markers. The advantages and disadvantages of the various marker types are discussed, along with their application to the reference chicken genetic linkage maps and to the search for quantitative trait loci (QTL). The prospects for the use of DNA-based markers in marker-assisted selection are considered, along with likely future trends in poultry gene mapping. Further development of both physical and linkage genome maps of the chicken will allow animal scientists to more efficiently detect and characterize QTL and will provide them access to the wealth of genetic information that is being generated about the human genome and the genomes of model species, such as the mouse and Drosophila.

Characterization of CR1 repeat random PCR markers for mapping the chicken genome.

Polymerase chain reaction (PCR) primers complementary to portions of the chicken repetitive element CR1 have been used previously to generate useful markers on the chicken genome linkage map. To understand better the genetic basis for this technique and to convert CR1-PCR loci to markers useful in physical genome mapping, five polymorphic CR1-PCR-generated DNAs were cloned and partially sequenced. Inverse PCR was then employed to clone the corresponding region of the genomes of both the Jungle Fowl (JF) and White Leghorn (WL) parental DNA templates. Our results demonstrate that some of the CR1-PCR-generated DNAs arise from priming at an endogenous CR1 element, whereas others are due to chance complementarity between the CR1-PCR primer in use and random annealing sites in the genome, unrelated to a demonstrable CR1 element. In all five instances, it was possible to identify the sequence difference between the JF and WL parental DNAs that gave rise to the initial polymorphism and design allele-specific PCR primer sets that uniquely detect that polymorphism. In four of the five instances, the polymorphism was a one or two basepair sequence difference within the primer annealing site, but in the fifth case the responsible difference was outside, but very close to, the annealing site. In all instances the allele-specific PCR for the sequence polymorphism mapped identically with the corresponding CR1-PCR amplification polymorphism. We conclude that CR1-PCR provides an efficient and reliable mechanism for genome mapping in avians that can correlate linkage and physical mapping approaches.

Genomic organization and evolution of the soybean SB92 satellite sequence.

Repetitive DNA sequences comprise a large percentage of plant genomes, and their characterization provides information about both species and genome evolution. We have isolated a recombinant clone containing a highly repeated DNA element (SB92) that is homologous to ca. 0.9% of the soybean genome or about 10(5) copies. This repeated sequence is tandemly arranged and is found in four or five major genomic locations. FISH analysis of metaphase chromosomes suggests that two of these locations are centromeric. We have determined the sequence of two cloned repeats and performed genomic sequencing to obtain a consensus sequence. The consensus repeat size was 92 bp and exhibited an average of 10% nucleotide substitution relative to the two cloned repeats. This high level of sequence diversity suggests an ancient origin but is inconsistent with the limited phylogenetic distribution of SB92, which is found at high copy number only in the annual soybeans. It therefore seems likely that this sequence is undergoing very rapid evolution.

Mapping quantitative trait loci controlling milk production in dairy cattle by exploiting progeny testing.

We have exploited "progeny testing" to map quantitative trait loci (QTL) underlying the genetic variation of milk production in a selected dairy cattle population. A total of 1,518 sires, with progeny tests based on the milking performances of > 150,000 daughters jointly, was genotyped for 159 autosomal microsatellites bracketing 1645 centimorgan or approximately two thirds of the bovine genome. Using a maximum likelihood multilocus linkage analysis accounting for variance heterogeneity of the phenotypes, we identified five chromosomes giving very strong evidence (LOD score > or = 3) for the presence of a QTL controlling milk production: chromosomes 1, 6, 9, 10 and 20. These findings demonstrate that loci with considerable effects on milk production are still segregating in highly selected populations and pave the way toward marker-assisted selection in dairy cattle breeding.

The mitochondrial ribosomal RNA genes of the nematodes Caenorhabditis elegans and Ascaris suum: consensus secondary-structure models and conserved nucleotide sets for phylogenetic analysis.

The small- and large-subunit mitochondrial ribosomal RNA genes (mt-s-rRNA and mt-l-rRNA) of the nematode worms Caenorhabditis elegans and Ascaris suum encode the smallest rRNAs so far reported for metazoa. These size reductions correlate with the previously described, smaller, structurally anomalous mt-tRNAs of C. elegans and A. suum. Using primer extension analysis, the 5' end nucleotides of the mt-s-rRNA and mt-l-rRNA genes were determined to be adjacent to the 3' end nucleotides of the tRNA(Glu) and tRNA(His) genes, respectively. Detailed, consensus secondary-structure models were constructed for the mt-s-rRNA genes and the 3' 64% of mt-l-rRNA genes of the two nematodes. The mt-s-rRNA secondary-structure model bears a remarkable resemblance to the previously defined universal core structure of E. coli 16S rRNA: most of the nucleotides that have been classified as variable or semiconserved in the E. coli model appear to have been eliminated from the C. elegans and A. suum sequences. Also, the secondary structure model constructed for the 3' 64% of the mt-l-rRNA is similar to the corresponding portion of the previously defined E. coli 23S rRNA core secondary structure. The proposed C. elegans/A. suum mt-s-rRNA and mt-l-rRNA models include all of the secondary-structure element-forming sequences that in E. coli rRNAs contain nucleotides important for A-site and P-site (but not E-site) interactions with tRNAs. Sets of apparently homologous sequences within the mt-s-rRNA and mt-l-rRNA core structures, derived by alignment of the C. elegans and A. suum mt-rRNAs to the corresponding mt-rRNAs of other eukaryotes, and E. coli rRNAs were used in maximum-likelihood analyses. The patterns of divergence of metazoan phyla obtained show considerable agreement with the most prevalent metazoan divergence patterns derived from more classical, morphological, and developmental data.