Design of reference populations for genomic selection in crossbreeding programs

Background

In crossbreeding programs, genomic selection offers the opportunity to make efficient use of information on crossbred (CB) individuals in the selection of purebred (PB) candidates. In such programs, reference populations often contain genotyped PB animals, although the breeding objective is usually more focused on CB performance. The question is what would be the benefit of including a larger proportion of CB individuals in the reference population.

Methods

In a deterministic simulation study, we evaluated the benefit of including various proportions of CB animals in a reference population for genomic selection of PB animals in a crossbreeding program. We used a pig breeding scheme with selection for a moderately heritable trait and a size of 6000 for the reference population.

Results

Applying genomic selection to improve the performance of CB individuals, with a genetic correlation between PB and CB performance (r_PC) of 0.7, selection accuracy of PB candidates increased from 0.49 to 0.52 if the reference population consisted of PB individuals, it increased to 0.55 if the reference population consisted of the same number of CB individuals, and to 0.60 if the size of the CB reference population was twice that of the reference population for each PB line. The advantage of using CB rather than PB individuals increased linearly with the proportion of CB individuals in the reference population. This advantage disappeared quickly if r_PC was higher or if the breeding objective put some emphasis on PB performance. The benefit of adding CB individuals to an existing PB reference population was limited for high r_PC.

Conclusions

Using CB rather than PB individuals in a reference population for genomic selection can provide substantial advantages, but only when correlations between PB and CB performances are not high and PB performance is not part of the breeding objective.     

Changes in Uric Acid Levels following Bariatric Surgery Are Not Associated with SLC2A9 Variants in the Swedish Obese Subjects Study

Context and Objective

Obesity and SLC2A9 genotype are strong determinants of uric acid levels. However, data on SLC2A9 variants and weight loss induced changes in uric acid levels are missing. We examined whether the changes in uric acid levels two- and ten-years after weight loss induced by bariatric surgery were associated with SLC2A9 single nucleotide polymorphisms (SNPs) in the Swedish Obese Subjects study.

Methods

SNPs (N = 14) identified by genome-wide association studies and exonic SNPs in the SLC2A9 gene locus were genotyped. Cross-sectional associations were tested before (N = 1806), two (N = 1664) and ten years (N = 1201) after bariatric surgery. Changes in uric acid were compared between baseline and Year 2 (N = 1660) and years 2 and 10 (N = 1172). A multiple testing corrected threshold of P = 0.007 was used for statistical significance.

Results

Overall, 11 of the 14 tested SLC2A9 SNPs were significantly associated with cross-sectional uric acid levels at all three time points, with rs13113918 showing the strongest association at each time point (R² = 3.7−5.2%, 3.9×10⁻²²≤p≤7.7×10⁻¹¹). One SNP (rs737267) showed a significant association (R² = 0.60%, P = 0.002) with change in uric acid levels from baseline to Year 2, as common allele homozygotes (C/C, N = 957) showed a larger decrease in uric acid (−61.4 µmol/L) compared to minor allele carriers (A/X: −51.7 µmol/L, N = 702). No SNPs were associated with changes in uric acid from years 2 to 10.

Conclusions

SNPs in the SLC2A9 locus contribute significantly to uric acid levels in obese individuals, and the associations persist even after considerable weight loss due to bariatric surgery. However, we found little evidence for an interaction between genotype and weight change on the response of uric acid to bariatric surgery over ten years. Thus, the fluctuations in uric acid levels among the surgery group appear to be driven by the weight losses and gains, independent of SLC2A9 genotypes.     

Anti-citrullinated fibronectin antibodies in rheumatoid arthritis are associated with human leukocyte antigen-DRB1 shared epitope alleles

Introduction

Fibronectin is one of the most abundant proteins present in the inflamed joint. Here, we characterized the citrullination of fibronectin in the joints of rheumatoid arthritis (RA) patients and studied the prevalence, epitope specificity and human leukocyte antigen (HLA) association of autoantibodies against citrullinated fibronectin in RA.

Methods

Citrullinated residues in fibronectin isolated from RA patient synovial fluid were identified by mass spectrometry. The corresponding citrullinated and non-citrullinated peptides were synthesized and used to analyze the presence of autoantibodies to these peptides in RA sera and sera from other diseases and healthy controls by ELISA. The data were compared with risk factors like shared epitope HLA alleles and smoking, and with clinical features.

Results

Five citrullinated residues were identified in fibronectin from RA synovial fluid. RA sera reacted in a citrulline-dependent manner with two out of four citrullinated fibronectin peptides, one of which contains two adjacent citrulline residues, in contrast to non-RA sera, which were not reactive. The most frequently recognized peptide (FN-Cit_1035,1036, LTVGLTXXGQPRQY, in which × represents citrulline) was primarily targeted by anti-CCP (cyclic citrullinated peptide) 2-positive RA patients. Anti-FN-Cit_1035,1036 autoantibodies were detected in 50% of established anti-CCP2-positive RA patients and in 45% of such patients from a early arthritis clinic. These antibodies appeared to be predominantly of the immunoglobulin G (IgG) isotype and to be associated with HLA shared epitope alleles (odds ratio = 2.11).

Conclusions

Fibronectin in the inflamed synovia of RA patients can be citrullinated at least at five positions. Together with the flanking amino acids, three of these citrullinated residues comprise two epitopes recognized by RA autoantibodies. Anti-citrullinated fibronectin peptide antibodies are associated with HLA shared epitope alleles.     

Automatically extracting functionally equivalent proteins from SwissProt

Background  

There is a frequent need to obtain sets of functionally equivalent homologous proteins (FEPs) from different species. While it is usually the case that orthology implies functional equivalence, this is not always true; therefore datasets of orthologous proteins are not appropriate. The information relevant to extracting FEPs is contained in databanks such as UniProtKB/Swiss-Prot and a manual analysis of these data allow FEPs to be extracted on a one-off basis. However there has been no resource allowing the easy, automatic extraction of groups of FEPs – for example, all instances of protein C.     

Improved genome assembly and evidence-based global gene model set for the chordate Ciona intestinalis: new insight into intron and operon populations

Background

The draft genome sequence of the ascidian Ciona intestinalis, along with associated gene models, has been a valuable research resource. However, recently accumulated expressed sequence tag (EST)/cDNA data have revealed numerous inconsistencies with the gene models due in part to intrinsic limitations in gene prediction programs and in part to the fragmented nature of the assembly.

Results

We have prepared a less-fragmented assembly on the basis of scaffold-joining guided by paired-end EST and bacterial artificial chromosome (BAC) sequences, and BAC chromosomal in situ hybridization data. The new assembly (115.2 Mb) is similar in length to the initial assembly (116.7 Mb) but contains 1,272 (approximately 50%) fewer scaffolds. The largest scaffold in the new assembly incorporates 95 initial-assembly scaffolds. In conjunction with the new assembly, we have prepared a greatly improved global gene model set strictly correlated with the extensive currently available EST data. The total gene number (15,254) is similar to that of the initial set (15,582), but the new set includes 3,330 models at genomic sites where none were present in the initial set, and 1,779 models that represent fusions of multiple previously incomplete models. In approximately half, 5'-ends were precisely mapped using 5'-full-length ESTs, an important refinement even in otherwise unchanged models.

Conclusion

Using these new resources, we identify a population of non-canonical (non-GT-AG) introns and also find that approximately 20% of Ciona genes reside in operons and that operons contain a high proportion of single-exon genes. Thus, the present dataset provides an opportunity to analyze the Ciona genome much more precisely than ever.     

Monitoring of residues of ochratoxin A in blood and kidney of chicken using HPLC-FLD

Detoxification of ochratoxin A can be achieved by chemical or enzymatic hydrolyzation, the products of such reactions are ochratoxin α and phenylalanine. Ochratoxin α like ochratoxin A, is a fluorescing molecule, therefore sensitive analysis is possible at very low concentration levels. Methods have been established that make it possible to look for residues of ochratoxin A and its main metabolite ochratoxin α in blood and tissues at very low concentration levels. Plasma is extracted by the use of small amounts of chloroform; the extract is cleaned with water and afterwards evaporated to dryness]. The residue is re-dissolved and analysed by HPLC-FLD. Using this method a limit of detection of 0.5μg/l for both ochratoxin A and ochratoxin α can be reached.