Breeding value prediction for production traits in layer chickens using pedigree or genomic relationships in a reduced animal model

Background

Genomic selection involves breeding value estimation of selection candidates based on high-density SNP genotypes. To quantify the potential benefit of genomic selection, accuracies of estimated breeding values (EBV) obtained with different methods using pedigree or high-density SNP genotypes were evaluated and compared in a commercial layer chicken breeding line.

Methods

The following traits were analyzed: egg production, egg weight, egg color, shell strength, age at sexual maturity, body weight, albumen height, and yolk weight. Predictions appropriate for early or late selection were compared. A total of 2,708 birds were genotyped for 23,356 segregating SNP, including 1,563 females with records. Phenotypes on relatives without genotypes were incorporated in the analysis (in total 13,049 production records).The data were analyzed with a Reduced Animal Model using a relationship matrix based on pedigree data or on marker genotypes and with a Bayesian method using model averaging. Using a validation set that consisted of individuals from the generation following training, these methods were compared by correlating EBV with phenotypes corrected for fixed effects, selecting the top 30 individuals based on EBV and evaluating their mean phenotype, and by regressing phenotypes on EBV.

Results

Using high-density SNP genotypes increased accuracies of EBV up to two-fold for selection at an early age and by up to 88% for selection at a later age. Accuracy increases at an early age can be mostly attributed to improved estimates of parental EBV for shell quality and egg production, while for other egg quality traits it is mostly due to improved estimates of Mendelian sampling effects. A relatively small number of markers was sufficient to explain most of the genetic variation for egg weight and body weight.     

Enrichment of DNRA bacteria in a continuous culture

Denitrification and dissimilatory nitrate reduction to ammonium (DNRA) are competing microbial nitrate-reduction processes. The occurrence of DNRA has been shown to be effected qualitatively by various parameters in the environment. A more quantitative understanding can be obtained using enrichment cultures in a laboratory reactor, yet no successful DNRA enrichment culture has been described. We showed that a stable DNRA-dominated enrichment culture can be obtained in a chemostat system. The enrichment was based on the hypothesis that nitrate limitation is the dominant factor in selecting for DNRA. First, a conventional denitrifying culture was enriched from activated sludge, with acetate and nitrate as substrates. Next, the acetate concentration in the medium was increased to obtain nitrate-limiting conditions. As a result, conversions shifted from denitrification to DNRA. In this selection of a DNRA culture, two important factors were the nitrate limitation and a relatively low dilution rate (0.026 h⁻¹). The culture was a highly enriched population of Deltaproteobacteria most closely related to Geobacter lovleyi, based on 16S rRNA gene sequencing (97% similarity). We established a stable and reproducible cultivation method for the enrichment of DNRA bacteria in a continuously operated reactor system. This enrichment method allows to further investigate the DNRA process and address the factors for competition between DNRA and denitrification, or other N-conversion pathways.     

Genomic selection of purebred animals for crossbred performance in the presence of dominant gene action

Background

Genomic selection is an appealing method to select purebreds for crossbred performance. In the case of crossbred records, single nucleotide polymorphism (SNP) effects can be estimated using an additive model or a breed-specific allele model. In most studies, additive gene action is assumed. However, dominance is the likely genetic basis of heterosis. Advantages of incorporating dominance in genomic selection were investigated in a two-way crossbreeding program for a trait with different magnitudes of dominance. Training was carried out only once in the simulation.

Results

When the dominance variance and heterosis were large and overdominance was present, a dominance model including both additive and dominance SNP effects gave substantially greater cumulative response to selection than the additive model. Extra response was the result of an increase in heterosis but at a cost of reduced purebred performance. When the dominance variance and heterosis were realistic but with overdominance, the advantage of the dominance model decreased but was still significant. When overdominance was absent, the dominance model was slightly favored over the additive model, but the difference in response between the models increased as the number of quantitative trait loci increased. This reveals the importance of exploiting dominance even in the absence of overdominance. When there was no dominance, response to selection for the dominance model was as high as for the additive model, indicating robustness of the dominance model. The breed-specific allele model was inferior to the dominance model in all cases and to the additive model except when the dominance variance and heterosis were large and with overdominance. However, the advantage of the dominance model over the breed-specific allele model may decrease as differences in linkage disequilibrium between the breeds increase. Retraining is expected to reduce the advantage of the dominance model over the alternatives, because in general, the advantage becomes important only after five or six generations post-training.

Conclusion

Under dominance and without retraining, genomic selection based on the dominance model is superior to the additive model and the breed-specific allele model to maximize crossbred performance through purebred selection.     

Metastasis of bronchogenic carcinoma to the thumb

Bone metastasis in the hand is rare. The etiology is quite different from that of metastasis to other bones; bronchogenic carcinoma is by far the most frequent case. Distal phalanges are mainly involved with irregular osteolysis and cortical destruction. Differential diagnosis of phalangeal metastasis includes osteomyelitis, rheumatoid arthritis and gout. The prognosis is always that of metastatic bronchial cancer with an average survival of three months. Treatment may involve distal digital amputation or antalgic radiotherapy. A case of bronchogenic carcinoma with metastasis to the thumb is presented. The metastasis was located in the distal phalanx of the left thumb. The primary tumor was located in the lung. Treatment consisted of amputation. The overall survival was five months.     

Can Nep1-like proteins form oligomers?

Nep1-like proteins (NLPs) are a novel family of microbial elicitors of plant necrosis that induce a hypersensitive-like response in dicot plants. The spatial structure and role of these proteins are yet unknown. In a paper published in BMC Plant Biology (2008; 8:50) we have proposed that the core region of Nep1-like proteins (NLPs) belong to the Cupin superfamily. Based on what is known about the Cupin superfamily, in this addendum to the paper we discuss how NLPs could form oligomers.Key words: quaternary structure, necrosis and ethylene inducing proteins, NLPs, MpNEP1, MpNEP2, NPP1, Moniliophthora perniciosa, Phytophthora parasiticaCupins may be organized as monomers, dimers, hexamers and octamers of β-barrel domains.¹ To the best of our knowledge trimers have not been detected yet. The interaction of two monomers building up a dimeric structure is basically performed by three types of interactions: hydrophobic interactions between β-strands in different subunits, salt bridges and hydrogen bonds between β-strands. In cupin dimers, the hydrophobic interactions occur between two βI strands in different subunits (Fig. 1A and B). This strand represents the central axis of rotation of the dimer as one residue in βI interacts with the corresponding residue in the other subunit (Fig. 1B). Therefore, all residues in βI must be hydrophobic, as one residue interacts with the other subunit and the next one in the sequence interacts with the interior of the protein. Charged residues in βI would disrupt such interactions. Most cupin dimers have strong hydrophobic residues such as tryptophan (W), phenylalanine (F) and methionine (M) pointing towards the own subunit (↓), while small hydrophobic residues such as leucine (L), isoleucine (I), and valine (V) point to the other subunit (↑). A particular case is leucine that interacts with other subunits, for instance, βI = liaW (positions 217–220 in Fig. 1B) and βI = LVsw of type I and II NLP consensuses, respectively. Therefore, the pattern of hydropathicity suggests that the side chain orientation is βI = l217 ↑ i218 ↓ a219 ↑ W220 ↓ d221 ↑. However we observe that just after βI there is a charged residue (aspartate D221) which would point outwards disrupting the dimer or at least making it less stable. It is interesting to observe that the requirement for a negatively charged residue at this last position is very high: 96% of all type I NLPs contains an aspartate (D) or glutamate (E) indicating an important role for it, maybe in avoiding dimerization of the NLPs. A second interesting hypothesis is as follows: several cupins are oxygenases, decarboxylases, etc. and use a negatively charged residue, such as aspartate or glutamate as proton donor.¹ Now, if the alternate pattern of side chains of the residues is βI = l217 ↓ i218 ↑ a219 ↓ W220 ↑ d221 ↓, instead of the previous one, then the aspartate or glutamate residue would point to the hydrophobic pocket and would be positioned to interact with the metal ion, as in cupins with enzymatic activity. However, there are no experimental evidences that the NLPs have enzymatic activity.Open in a separate windowFigure 1(A) Three-dimensional structure prediction for type I NLP consensus, (B) Interface between two βI strands in type I NLP consensus. From the left to the right: EF-coil with the conserved residue H162, βC and βH strands (superposed) with the conserved histidines H133 and H135 in βC, H193 and leucine L195 in βH, W220 in βI and W118 in βB. The strands in the right subunit follow the same pattern but rotated.The second type of interaction is salt bridges between charged residues in different subunits. Analyzing all interacting side chains in the 1VJ2 protein (dimer), we verify that the charged side chains of N35 and E57 (numbers in original structure) are only 2.72 Å apart. In the NLPs, this corresponds to N108^36% (Q108^60%) at the border of βB and E138^98%. The negatively charged residue D125 helps to correct the orientation of the subunits in relation to each other avoiding any disorientation. The high conservation level of these residues suggests that NLPs are dimeric structures. However, as we will see next, only hydrophobic and charged interactions are not enough to build a dimer.Garcia et al. (2007)² have used small angle X-ray scattering (SAXS) to show that, in solution, at low concentrations (<2 mg/ml) the two copies of the NLPs of Moniliophthora perniciosa, MpNEP1 and MpNEP2, exist as dimers and monomers, respectively. The same technique showed that at higher concentrations, >5 mg/ml, both proteins exist as dimers, as is the case for PpNPP1.² They also reported, based on electrophoresis analysis, that PpNPP1 and MpNEP1 exist as oligomers and MpNEP2 as monomers.² However, experiments with the PpNPP1 in size exclusion chromatography using myoglobin as size standard suggest that PpNPP1 is a monomer.³ Figure 2 compares MpNEP1, MpNEP2 and PpNPP1, where the most relevant differences in sequence are marked with asterisks (*) and are possibly related to the differences in oligomeric properties between MpNEP1 and PpNPP1 with MpNEP2. These positions are methionine M27 and leucine L35, which occur only in MpNEP2, glycine G250, which occurs only in MpNEP2 and NEP1 (Fusarium oxysporum) and lysine K31, which occurs only MpNEP2, {"type":"entrez-protein","attrs":{"text":"BAB04114","term_id":"10173008"}}BAB04114 (Bacillus halodurans) and {"type":"entrez-protein","attrs":{"text":"AAU23136","term_id":"52003194"}}AAU23136 (Bacillus licheniformis). The other residues are aspartate D28, which occurs 9 times and alanine A37 which occurs 7 times of all investigated NLPs. Thus, the sequence mdHDkiakl at the start of the NLPs seems to explain the monomeric state of MpNEP2, although at higher concentrations they form dimers. Besides the weak hydrophobic interactions, dimeric cupins and bicupins (two β barrels in the same sequence building up a dimeric-like 4d-structure) are stable structures (see Fig. 1 in ref. ⁴). By aggregating the first β-strand in the start domain of one β-barrel to the ABIDG β-sheet of the other β-barrel, composing a big ABIDGY β-sheet (Y is the first β-strand). For instance, using the bicupin 1L3J (oxalate decarboxylase) as template, the low confidence level β-strand at position 26–33 (v in H29D30 avv) in type I NLPs corresponds to the first β-strand. Since this proceeds from both barrels they can build a stable structure (see Fig. 1 in ref. ⁴). The quaternary structure is related to the presence of interaction residues in the BID β-sheet of the cupin structure. These are present in the NLPs and would enable them to form dimers.Open in a separate windowFigure 2Alignment of type I NLP consensus, PpNPP1, MpNEP1 and MpNEP2. Solid line boxes are β-strands, double line boxes are α-helices. The sequence positions marked with asterisks (*) are possibly related to the differences in oligomeric properties between MpNEP1 and PpNPP1 with MpNEP2.     

The multiple sex chromosomes of platypus and echidna are not completely identical and several share homology with the avian Z

Background

Sex-determining systems have evolved independently in vertebrates. Placental mammals and marsupials have an XY system, birds have a ZW system. Reptiles and amphibians have different systems, including temperature-dependent sex determination, and XY and ZW systems that differ in origin from birds and placental mammals. Monotremes diverged early in mammalian evolution, just after the mammalian clade diverged from the sauropsid clade. Our previous studies showed that male platypus has five X and five Y chromosomes, no SRY, and DMRT1 on an X chromosome. In order to investigate monotreme sex chromosome evolution, we performed a comparative study of platypus and echidna by chromosome painting and comparative gene mapping.

Results

Chromosome painting reveals a meiotic chain of nine sex chromosomes in the male echidna and establishes their order in the chain. Two of those differ from those in the platypus, three of the platypus sex chromosomes differ from those of the echidna and the order of several chromosomes is rearranged. Comparative gene mapping shows that, in addition to bird autosome regions, regions of bird Z chromosomes are homologous to regions in four platypus X chromosomes, that is, X₁, X₂, X₃, X₅, and in chromosome Y₁.

Conclusion

Monotreme sex chromosomes are easiest to explain on the hypothesis that autosomes were added sequentially to the translocation chain, with the final additions after platypus and echidna divergence. Genome sequencing and contig anchoring show no homology yet between platypus and therian Xs; thus, monotremes have a unique XY sex chromosome system that shares some homology with the avian Z.     

Accuracies of genomic breeding values in American Angus beef cattle using K-means clustering for cross-validation

Background

Genomic selection is a recently developed technology that is beginning to revolutionize animal breeding. The objective of this study was to estimate marker effects to derive prediction equations for direct genomic values for 16 routinely recorded traits of American Angus beef cattle and quantify corresponding accuracies of prediction.

Methods

Deregressed estimated breeding values were used as observations in a weighted analysis to derive direct genomic values for 3570 sires genotyped using the Illumina BovineSNP50 BeadChip. These bulls were clustered into five groups using K-means clustering on pedigree estimates of additive genetic relationships between animals, with the aim of increasing within-group and decreasing between-group relationships. All five combinations of four groups were used for model training, with cross-validation performed in the group not used in training. Bivariate animal models were used for each trait to estimate the genetic correlation between deregressed estimated breeding values and direct genomic values.

Results

Accuracies of direct genomic values ranged from 0.22 to 0.69 for the studied traits, with an average of 0.44. Predictions were more accurate when animals within the validation group were more closely related to animals in the training set. When training and validation sets were formed by random allocation, the accuracies of direct genomic values ranged from 0.38 to 0.85, with an average of 0.65, reflecting the greater relationship between animals in training and validation. The accuracies of direct genomic values obtained from training on older animals and validating in younger animals were intermediate to the accuracies obtained from K-means clustering and random clustering for most traits. The genetic correlation between deregressed estimated breeding values and direct genomic values ranged from 0.15 to 0.80 for the traits studied.

Conclusions

These results suggest that genomic estimates of genetic merit can be produced in beef cattle at a young age but the recurrent inclusion of genotyped sires in retraining analyses will be necessary to routinely produce for the industry the direct genomic values with the highest accuracy.     

Prospects for estimating nucleotide divergence with RAPDs

The technique of random amplification of polymorphic DNA (RAPD), which is simply polymerase chain reaction (PCR) amplification of genomic DNA by a single short oligonucleotide primer, produces complex patterns of anonymous polymorphic DNA fragments. The information provided by these banding patterns has proved to be of great utility for mapping and for verification of identity of bacterial strains. Here we consider whether the degree of similarity of the banding patterns can be used to estimate nucleotide diversity and nucleotide divergence. With haploid data, fragments generated by RAPD-PCR can be treated in a fashion very similar to that for restriction-fragment data. Amplification of diploid samples, on the other hand, requires consideration of the fact that presence of a band is dominant to absence of the band. After describing a method for estimating nucleotide divergence on the basis of diploid samples, we summarize the restrictions and criteria that must be met when RAPD data are used for estimating population genetic parameters.