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.     

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.     

Endogenous chloride channels of insect sf9 cells. Evidence for coordinated activity of small elementary channel units

The endogenous Cl- conductance of Spodoptera frugiperda (Sf9) cells was studied 20-35 h after plating out of either uninfected cells or cells infected by a baculovirus vector carrying the cloned beta-galactosidase gene (beta-Gal cells). With the cation Tris+ in the pipette and Na+ in the bath, the reversal potential of whole-cell currents was governed by the prevailing Cl- equilibrium potential and could be fitted by the Goldman-Hodgkin-Katz equation with similar permeabilities for uninfected and beta-Gal cells. In the frequency range 0.12 < f < 300 Hz, the power density spectrum of whole-cell Cl- currents could be fitted by three Lorentzians. Independent of membrane potential, >50% of the total variance of whole-cell current fluctuations was accounted for by the low frequency Lorentzian (fc = 0.40 +/- 0.03 Hz, n = 6). Single-Cl- channels showed complex gating kinetics with long lasting (seconds) openings interrupted by similar long closures. In the open state, channels exhibited fast burst-like closures. Since the patches normally contained more than a single channel, it was not possible to measure open and closed dwell-time distributions for comparing single-Cl- channel activity with the kinetic features of whole-cell currents. However, the power density spectrum of Cl- currents of cell-attached and excised outside-out patches contained both high and low frequency Lorentzian components, with the corner frequency of the slow component (fc = 0.40 +/- 0.02 Hz, n = 4) similar to that of whole-cell current fluctuations. Chloride channels exhibited multiple conductance states with similar Goldman-Hodgkin-Katz-type rectification. Single-channel permeabilities covered the range from approximately 0.6.10(-14) cm5/s to approximately 6.10(-14) cm3/s, corresponding to a limiting conductance (gamma 150/150) of approximately 3.5 pS and approximately 35 pS, respectively. All states reversed near the same membrane potential, and they exhibited similar halide ion selectivity, P1 > PCl approximately PBr. Accordingly, Cl- current amplitudes larger than current flow through the smallest channel unit resolved seem to result from simultaneous open/shut events of two or more channel units.     

A change in twist of actin provides the force for the extension of the acrosomal process in limulus sperm: the false-discharge reaction

One of the most spectacular motions is the generation of the acrosomal process in the limulus sperm. On contact with the egg, the sperm generates a 60-mum-long process that literally drills its way through the jelly surrounding the egg. This irresversible reaction takes only a few seconds. We suggested earlier that this motion is driven by a change in twist of the actin filaments comprising the acrosomal process. In this paper we analyze the so-called false discharge, a reversible reaction, in which the acrosomal filament bundle extends laterally from the base of the sperm and not anteriorly from the apex. Unlike the true discharge, which is straight, the false discharge is helical. Before extension, the filament bundle is coiled about the base of the sperm. In the coil, the bundle is not smoothly bent but consists of arms (straight segments) and elbows (corners) so that the coil looks like a 14-sided polygon. The extension of the false discharge works as follows: starting at the base of the bundle, the filaments change their twist which concomitantly changes the orientations of the elbows relative to each other; that is, in the coil, the elbows all like in a common plane, but after the change in twist, the plane of each elbow is rotated to be perpendicular to that of its neighbors. This change transforms the bundle from a compact coil into an extended left- handed helix. Because the basal end of the bundle is unconstrained, the extension is lateral. The true discharge works the same way but starts at the apical end of the bundle. The apical end, however, is constrained by its passage through the nuclear canal, which directs the extention anteriorly. Unlike the false discharge, during the true discharge the elbows are melted out, making the reaction irreversible. This study shows that rapid movement can be regenerated by actin without myosin and gives us insight into the molecular mechanism.