A comparison of alternative methods to compute conditional genotype probabilities for genetic evaluation with finite locus models

An increased availability of genotypes at marker loci has prompted the development of models that include the effect of individual genes. Selection based on these models is known as marker-assisted selection (MAS). MAS is known to be efficient especially for traits that have low heritability and non-additive gene action. BLUP methodology under non-additive gene action is not feasible for large inbred or crossbred pedigrees. It is easy to incorporate non-additive gene action in a finite locus model. Under such a model, the unobservable genotypic values can be predicted using the conditional mean of the genotypic values given the data. To compute this conditional mean, conditional genotype probabilities must be computed. In this study these probabilities were computed using iterative peeling, and three Markov chain Monte Carlo (MCMC) methods – scalar Gibbs, blocking Gibbs, and a sampler that combines the Elston Stewart algorithm with iterative peeling (ESIP). The performance of these four methods was assessed using simulated data. For pedigrees with loops, iterative peeling fails to provide accurate genotype probability estimates for some pedigree members. Also, computing time is exponentially related to the number of loci in the model. For MCMC methods, a linear relationship can be maintained by sampling genotypes one locus at a time. Out of the three MCMC methods considered, ESIP, performed the best while scalar Gibbs performed the worst.     

Haplotype inference in crossbred populations without pedigree information

Background

Current methods for haplotype inference without pedigree information assume random mating populations. In animal and plant breeding, however, mating is often not random. A particular form of nonrandom mating occurs when parental individuals of opposite sex originate from distinct populations. In animal breeding this is called crossbreeding and hybridization in plant breeding. In these situations, association between marker and putative gene alleles might differ between the founding populations and origin of alleles should be accounted for in studies which estimate breeding values with marker data. The sequence of alleles from one parent constitutes one haplotype of an individual. Haplotypes thus reveal allele origin in data of crossbred individuals.

Results

We introduce a new method for haplotype inference without pedigree that allows nonrandom mating and that can use genotype data of the parental populations and of a crossbred population. The aim of the method is to estimate line origin of alleles. The method has a Bayesian set up with a Dirichlet Process as prior for the haplotypes in the two parental populations. The basic idea is that only a subset of the complete set of possible haplotypes is present in the population.

Conclusion

Line origin of approximately 95% of the alleles at heterozygous sites was assessed correctly in both simulated and real data. Comparing accuracy of haplotype frequencies inferred with the new algorithm to the accuracy of haplotype frequencies inferred with PHASE, an existing algorithm for haplotype inference, showed that the DP algorithm outperformed PHASE in situations of crossbreeding and that PHASE performed better in situations of random mating.     

You can't always get what you want: size assortative mating by mutual mate choice as a resolution of sexual conflict

Background  

Assortative mating patterns for mate quality traits like body size are often observed in nature. However, the underlying mechanisms that cause assortative mating patterns are less well known. Sexual selection is one important explanation for assortment, suggesting that i) one (usually the female) or both sexes could show preferences for mates of similar size or ii) mutual mate choice could resolve sexual conflict over quality traits into assortment. We tested these hypotheses experimentally in the socially monogamous cichlid fish Pelvicachromis taeniatus, in which mate choice is mutual.     

Genome-wide association and genomic prediction for host response to porcine reproductive and respiratory syndrome virus infection

Background

Host genetics has been shown to play a role in porcine reproductive and respiratory syndrome (PRRS), which is the most economically important disease in the swine industry. A region on Sus scrofa chromosome (SSC) 4 has been previously reported to have a strong association with serum viremia and weight gain in pigs experimentally infected with the PRRS virus (PRRSV). The objective here was to identify haplotypes associated with the favorable phenotype, investigate additional genomic regions associated with host response to PRRSV, and to determine the predictive ability of genomic estimated breeding values (GEBV) based on the SSC4 region and based on the rest of the genome. Phenotypic data and 60 K SNP genotypes from eight trials of ~200 pigs from different commercial crosses were used to address these objectives.

Results

Across the eight trials, heritability estimates were 0.44 and 0.29 for viral load (VL, area under the curve of log-transformed serum viremia from 0 to 21 days post infection) and weight gain to 42 days post infection (WG), respectively. Genomic regions associated with VL were identified on chromosomes 4, X, and 1. Genomic regions associated with WG were identified on chromosomes 4, 5, and 7. Apart from the SSC4 region, the regions associated with these two traits each explained less than 3% of the genetic variance. Due to the strong linkage disequilibrium in the SSC4 region, only 19 unique haplotypes were identified across all populations, of which four were associated with the favorable phenotype. Through cross-validation, accuracies of EBV based on the SSC4 region were high (0.55), while the rest of the genome had little predictive ability across populations (0.09).

Conclusions

Traits associated with response to PRRSV infection in growing pigs are largely controlled by genomic regions with relatively small effects, with the exception of SSC4. Accuracies of EBV based on the SSC4 region were high compared to the rest of the genome. These results show that selection for the SSC4 region could potentially reduce the effects of PRRS in growing pigs, ultimately reducing the economic impact of this disease.     

Targeted overactivity of beta cell K(ATP) channels induces profound neonatal diabetes

A paradigm for control of insulin secretion is that glucose metabolism elevates cytoplasmic [ATP]/[ADP] in beta cells, closing K(ATP) channels and causing depolarization, Ca2+ entry, and insulin release. Decreased responsiveness of K(ATP) channels to elevated [ATP]/[ADP] should therefore lead to decreased insulin secretion and diabetes. To test this critical prediction, we generated transgenic mice expressing beta cell K(ATP) channels with reduced ATP sensitivity. Animals develop severe hyperglycemia, hypoinsulinemia, and ketoacidosis within 2 days and typically die within 5. Nevertheless, islet morphology, insulin localization, and alpha and beta cell distributions were normal (before day 3), pointing to reduced insulin secretion as causal. The data indicate that normal K(ATP) channel activity is critical for maintenance of euglycemia and that overactivity can cause diabetes by inhibiting insulin secretion.     

Apoptosis,cytokine and chemokine induction by non-structural 1 (NS1) proteins encoded by different influenza subtypes

Background

Influenza pandemic remains a serious threat to human health. Viruses of avian origin, H5N1, H7N7 and H9N2, have repeatedly crossed the species barrier to infect humans. Recently, a novel strain originated from swine has evolved to a pandemic. This study aims at improving our understanding on the pathogenic mechanism of influenza viruses, in particular the role of non-structural (NS1) protein in inducing pro-inflammatory and apoptotic responses.

Methods

Human lung epithelial cells (NCI-H292) was used as an in-vitro model to study cytokine/chemokine production and apoptosis induced by transfection of NS1 mRNA encoded by seven infleunza subtypes (seasonal and pandemic H1, H2, H3, H5, H7, and H9), respectively.

Results

The results showed that CXCL-10/IP10 was most prominently induced (> 1000 folds) and IL-6 was slightly induced (< 10 folds) by all subtypes. A subtype-dependent pattern was observed for CCL-2/MCP-1, CCL3/MIP-1α, CCL-5/RANTES and CXCL-9/MIG; where induction by H5N1 was much higher than all other subtypes examined. All subtypes induced a similar temporal profile of apoptosis following transfection. The level of apoptosis induced by H5N1 was remarkably higher than all others. The cytokine/chemokine and apoptosis inducing ability of the 2009 pandemic H1N1 was similar to previous seasonal strains.

Conclusions

In conclusion, the NS1 protein encoded by H5N1 carries a remarkably different property as compared to other avian and human subtypes, and is one of the keys to its high pathogenicity. NCI-H292 cells system proves to be a good in-vitro model to delineate the property of NS1 proteins.

     

Expression of the cDNA for NAD(+)-dependent 15-hydroxyprostaglandin dehydrogenase as a catalytically active enzyme in Escherichia coli.

NAD(+)-dependent 15-hydroxyprostaglandin dehydrogenase (15-PGDH) is a key enzyme involved in the catabolism of the prostaglandins. The cDNA for human placental 15-PGDH has been expressed in Escherichia coli as a catalytically active protein. The polymerase chain reaction was used to introduce restriction endonuclease sites at each end of the 15-PGDH coding sequence. The 15-PGDH DNA was then inserted into the bacterial expression plasmids pUC-18 and pUC-19 which contain the isopropyl-l-thio-beta-D-galactopyranoside (IPTG) inducible lacZ promoter. Extracts from E. coli containing these expression plasmids exhibited 15-PGDH activity which was inducible with (IPTG). Crude extracts from E. coli expressing 15-PGDH activity were found to contain proteins of the predicted sizes in stained SDS-polyacrylamide gels and in Western blots using human placental 15-PGDH antiserum. The specific activity in E. coli extracts was several hundred-fold higher than that seen in extracts from human placenta.     

Neurotoxic,cytotoxic, apoptotic and antiproliferative effects of some marine algae extracts on the NA2B cell line

Oxidative stress contributes to cancer pathologies and to apoptosis. Marine algae exhibit cytotoxic, antiproliferative and apoptotic effects; their metabolites have been used to treat many types of cancer. We investigated in culture extracts of Petalonia fascia, Jania longifurca and Halimeda tuna to determine their effects on mouse neuroblastoma cell line, NA2B. NA2B cells were treated with algae extracts, and the survival and proliferation of NA2B cells were assessed using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. The effects of algae extracts on oxidative stress in NA2B cells also were investigated using nitric oxide synthase (NOS) immunocytochemistry and apoptosis was assessed using terminal deoxynucleotidyl transferase dUTP nick end labeling. We observed significant neurite inhibition with moderate damage by the neurotoxicity-screening test (NST) at IC₅₀ dilutions of the extracts. MTT demonstrated that J. longifurca extracts were more toxic than P. fascia and H. tuna extracts. We found an increase of endothelial and inducible NOS immunostaining for oxidative stress and TUNEL analysis revealed increased apoptosis after application of extract. Our findings suggest that the algae we tested may have potential use for treatment of cancer.     

Regulators of complement activity mediate inhibitory mechanisms through a common C3b‐binding mode

Regulators of complement activation (RCA) inhibit complement‐induced immune responses on healthy host tissues. We present crystal structures of human RCA (MCP, DAF, and CR1) and a smallpox virus homolog (SPICE) bound to complement component C3b. Our structural data reveal that up to four consecutive homologous CCP domains (i–iv), responsible for inhibition, bind in the same orientation and extended arrangement at a shared binding platform on C3b. Large sequence variations in CCP domains explain the diverse C3b‐binding patterns, with limited or no contribution of some individual domains, while all regulators show extensive contacts with C3b for the domains at the third site. A variation of ~100° rotation around the longitudinal axis is observed for domains binding at the fourth site on C3b, without affecting the overall binding mode. The data suggest a common evolutionary origin for both inhibitory mechanisms, called decay acceleration and cofactor activity, with variable C3b binding through domains at sites ii, iii, and iv, and provide a framework for understanding RCA disease‐related mutations and immune evasion.