Design of a High Density SNP Genotyping Assay in the Pig Using SNPs Identified and Characterized by Next Generation Sequencing Technology

Background

The dissection of complex traits of economic importance to the pig industry requires the availability of a significant number of genetic markers, such as single nucleotide polymorphisms (SNPs). This study was conducted to discover several hundreds of thousands of porcine SNPs using next generation sequencing technologies and use these SNPs, as well as others from different public sources, to design a high-density SNP genotyping assay.

Methodology/Principal Findings

A total of 19 reduced representation libraries derived from four swine breeds (Duroc, Landrace, Large White, Pietrain) and a Wild Boar population and three restriction enzymes (AluI, HaeIII and MspI) were sequenced using Illumina''s Genome Analyzer (GA). The SNP discovery effort resulted in the de novo identification of over 372K SNPs. More than 549K SNPs were used to design the Illumina Porcine 60K+SNP iSelect Beadchip, now commercially available as the PorcineSNP60. A total of 64,232 SNPs were included on the Beadchip. Results from genotyping the 158 individuals used for sequencing showed a high overall SNP call rate (97.5%). Of the 62,621 loci that could be reliably scored, 58,994 were polymorphic yielding a SNP conversion success rate of 94%. The average minor allele frequency (MAF) for all scorable SNPs was 0.274.

Conclusions/Significance

Overall, the results of this study indicate the utility of using next generation sequencing technologies to identify large numbers of reliable SNPs. In addition, the validation of the PorcineSNP60 Beadchip demonstrated that the assay is an excellent tool that will likely be used in a variety of future studies in pigs.     

Stable isotope tagging of epitopes: a highly selective strategy for the identification of major histocompatibility complex class I-associated peptides induced upon viral infection

Identification of peptides presented in major histocompatibility complex (MHC) class I molecules after viral infection is of strategic importance for vaccine development. Until recently, mass spectrometric identification of virus-induced peptides was based on comparative analysis of peptide pools isolated from uninfected and virus-infected cells. Here we report on a powerful strategy aiming at the rapid, unambiguous identification of naturally processed MHC class I-associated peptides, which are induced by viral infection. The methodology, stable isotope tagging of epitopes (SITE), is based on metabolic labeling of endogenously synthesized proteins during infection. This is accomplished by culturing virus-infected cells with stable isotope-labeled amino acids that are expected to be anchor residues (i.e. residues of the peptide that have amino acid side chains that bind into pockets lining the peptide-binding groove of the MHC class I molecule) for the human leukocyte antigen allele of interest. Subsequently these cells are mixed with an equal number of non-infected cells, which are cultured in normal medium. Finally peptides are acid-eluted from immunoprecipitated MHC molecules and subjected to two-dimensional nanoscale LC-MS analysis. Virus-induced peptides are identified through computer-assisted detection of characteristic, binomially distributed ratios of labeled and unlabeled molecules. Using this approach we identified novel measles virus and respiratory syncytial virus epitopes as well as infection-induced self-peptides in several cell types, showing that SITE is a unique and versatile method for unequivocal identification of disease-related MHC class I epitopes.     

Heat-induced denaturation and aggregation of ovalbumin at neutral pH described by irreversible first-order kinetics

The heat-induced denaturation kinetics of two different sources of ovalbumin at pH 7 was studied by chromatography and differential scanning calorimetry. The kinetics was found to be independent of protein concentration and salt concentration, but was strongly dependent on temperature. For highly pure ovalbumin, the decrease in nondenatured native protein showed first-order dependence. The activation energy obtained with different techniques varied between 430 and 490 kJ*mole(-1). First-order behavior was studied in detail using differential scanning calorimetry. The calorimetric traces were irreversible and highly scan rate-dependent. The shape of the thermograms as well as the scan rate dependence can be explained by assuming that the thermal denaturation takes place according to a simplified kinetic process where N is the native state, D is denatured (or another final state) and k a first-order kinetic constant that changes with temperature, according to the Arrhenius equation. A kinetic model for the temperature-induced denaturation and aggregation of ovalbumin is presented. Commercially obtained ovalbumin was found to contain an intermediate-stable fraction (IS) of about 20% that was unable to form aggregates. The denaturation of this fraction did not satisfy first-order kinetics.     

Typing single-nucleotide polymorphisms using a gel-based sequencer: a new data analysis tool and suggestions for improved efficiency

Single-nucleotide polymorphisms (SNPs) are increasingly used as genetic markers. Although a high number of SNP-genotyping techniques have been described, most techniques still have low throughput or require major investments. For laboratories that have access to an automated sequencer, a single-base extension (SBE) assay can be implemented using the ABI SNaPshot trade mark kit. Here we present a modified protocol comprising multiplex template generation, multiplex SBE reaction, and multiplex sample analysis on a gel-based sequencer such as the ABI 377. These sequencers run on a Macintosh platform, but on this platform the software available for analysis of data from the ABI 377 has limitations. First, analysis of the size standard included with the kit is not facilitated. Therefore a new size standard was designed. Second, using Genotyper (ABI), the analysis of the data is very tedious and time consuming. To enable automated batch analysis of 96 samples, with 10 SNPs each, we developed SNPtyper. This is a spreadsheet-based tool that uses the data from Genotyper and offers the user a convenient interface to set parameters required for correct allele calling. In conclusion, the method described will enable any lab having access to an ABI sequencer to genotype up to 1000 SNPs per day for a single experimenter, without investing in new equipment.     

Principles of quantitation of viral loads using nucleic acid sequence-based amplification in combination with homogeneous detection using molecular beacons

For quantitative NASBA-based viral load assays using homogeneous detection with molecular beacons, such as the NucliSens EasyQ HIV-1 assay, a quantitation algorithm is required. During the amplification process there is a constant growth in the concentration of amplicons to which the beacon can bind while generating a fluorescence signal. The overall fluorescence curve contains kinetic information on both amplicon formation and beacon binding, but only the former is relevant for quantitation. In the current paper, mathematical modeling of the relevant processes is used to develop an equation describing the fluorescence curve as a function of the amplification time and the relevant kinetic parameters. This equation allows reconstruction of RNA formation, which is characterized by an exponential increase in concentrations as long as the primer concentrations are not rate limiting and by linear growth over time after the primer pool is depleted. During the linear growth phase, the actual quantitation is based on assessing the amplicon formation rate from the viral RNA relative to that from a fixed amount of calibrator RNA. The quantitation procedure has been successfully applied in the NucliSens EasyQ HIV-1 assay.     

Differential expression of inflammatory chemokines by Th1- and Th2-cell promoting dendritic cells: a role for different mature dendritic cell populations in attracting appropriate effector cells to peripheral sites of inflammation

Protective immunity to pathogens depends on efficient immune responses adapted to the type of pathogen and the infected tissue. Dendritic cells (DC) play a pivotal role in directing the effector T cell response to either a protective T helper type 1 (Th1) or type 2 (Th2) phenotype. Human monocyte-derived DC can be differentiated into Th1-, Th2- or Th1/Th2-promoting DC in vitro upon activation with microbial compounds or cytokines. Host defence is highly dependent on mobile leucocytes and cell trafficking is largely mediated by the interactions of chemokines with their specific receptors expressed on the surface of leucocytes. The production of chemokines by mature effector DC remains elusive. Here we assess the differential production of both inflammatory and homeostatic chemokines by monocyte-derived mature Th1/Th2-, Th1- or Th2-promoting DC and its regulation in response to CD40 ligation, thereby mimicking local engagement with activated T cells. We show that mature Th1- and Th1/Th2-, but not Th2-promoting DC, selectively express elevated levels of the inflammatory chemokines CCL2/MCP-1, CCL3/MIP-1alpha, CCL4/MIP-1beta and CCL5/RANTES, as well as the homeostatic chemokine CCL19/MIP-3beta. CCL21/6Ckine is preferentially expressed by Th2-promoting DC. Production of the Th1-attracting chemokines, CXCL9/Mig, CXCL10/IP-10 and CXCL11/I-TAC, is restricted to Th1-promoting DC. In contrast, expression of Th2-associated chemokines does not strictly correlate with the Th2-promoting DC phenotype, except for CCL22/MDC, which is preferentially expressed by Th2-promoting DC. Because inflammatory chemokines and Th1-associated chemokines are constitutively expressed by mature Th1-promoting DC and CCL22/MDC is constitutively expressed by mature Th2-promoting DC, we propose a novel role for mature DC present in inflamed peripheral tissues in orchestrating the immune response by recruiting appropriate leucocyte populations to the site of pathogen entry.     

The chicken leukocyte receptor complex: a highly diverse multigene family encoding at least six structurally distinct receptor types

The chicken Ig-like receptors (CHIR) have been described as two Ig domain molecules with long cytoplasmic tails containing inhibitory motifs. In this study, we demonstrate that CHIR form a large family, with multiple members showing great sequence variability among members as well as a great diversity in domain organization and properties of the transmembrane and cytoplasmic segments. We characterize various novel receptor types with motifs indicative of inhibitory, activating, or both functions. In addition to the inhibitory receptors with two ITIM, receptors with a single immunoreceptor tyrosine-based switch motif or receptors lacking a cytoplasmic domain were isolated. Activating receptors with a short cytoplasmic domain and a transmembrane arginine assembled with the newly identified chicken FcepsilonRIgamma chain. Three bifunctional receptor types were characterized composed of one or two C2-type Ig-like domains, a transmembrane region with a positively charged residue and combinations of cytoplasmic motifs such as ITIM, immunoreceptor tyrosine-based switch motif, and YXXM. RT-PCR revealed distinct expression patterns of individual CHIR. All receptor types shared a conserved genomic architecture, and in single Ig domain receptors a pseudoexon replaced the second Ig exon. Southern blot analyses with probes specific for the Ig1 domain were indicative of a large multigene family. Of 103 sequences from the Ig1 domain of a single animal, 41 unique sequences were obtained that displayed extensive variability within restricted Ig regions. Fluorescence in situ hybridization localized the CHIR gene cluster to microchromosome 31 and identified this region as orthologous to the human leukocyte receptor complex.     

Localization to chicken Chromosome 5 of a novel locus determining salmonellosis resistance

Clear genetic differences in the susceptibility of chickens to visceral infection by Salmonella have been observed and it has been possible to identify resistant and susceptible lines of inbred chickens. We report here the results of experiments to map directly the gene(s) controlling this trait in chickens by examining crosses between highly susceptible and highly resistant lines. In the mapping panel, a region on chicken Chromosome (Chr) 5 was found to have a large effect on resistance, and this effect was observed in three separate resource populations. Mapping of additional marker loci in the region of the resistance gene further localized it to a region of approximately 2 cM, close to the genes for creatine kinase (CKB) and dynein (DNCH1). This region shows conserved synteny with telomeric regions of human Chr 14 and mouse Chr 12. On the basis of this conserved synteny, this resistance gene seems unlikely to correspond to the previously identified salmonellosis resistance genes Lps (located on mouse Chr 4) or Nos(2) (located on mouse Chr 11). There was no association between Nramp1 and resistance in these crosses, although this gene was shown to contribute to resistance in other crosses. The homologous human and mouse regions at present contain no likely candidate genes for this trait. Thus this appears to be a novel resistance gene, which we designate SAL1.     

Mapping the Naked Neck (NA) and Polydactyly (PO) mutants of the chicken with microsatellite molecular markers

The bulked segregant analysis methodology has been used to map, with microsatellite markers, two morphological mutations in the chicken: polydactyly (PO) and naked neck (NA). These autosomal mutations show partial dominance for NA, and dominance with incomplete penetrance for PO. They were mapped previously to different linkage groups of the classical map, PO to the linkage group IV and NA being linked to the erythrocyte antigen CPPP. An informative family of 70 offspring was produced by mating a sire, heterozygous for each of the mutations, to 7 dams homozygous recessive for each locus. Three DNA pools were prepared, pool PO included 20 chicks exhibiting at least one extra-toe, pool NA included 20 non-polydactyly chicks showing the typical phenotype associated with heterozygosity for the naked neck mutation, and pool NP included 20 chicks exhibiting neither of the mutant phenotypes. Typings were done on an ABI-373 automatic sequencer with 147 microsatellite markers covering most of the genome. An unbalanced distribution of sire marker alleles were detected between pool PO, and pools NA and NP, for two markers of chromosome 2p, MCW0082 and MCW0247. A linkage analysis taking into account the incomplete penetrance of polydactyly (80%) was performed with additional markers of this region and showed that the closest marker to the PO locus was MCW0071 (5 cM, lod score = 9). MCW0071 lies within the engrailed gene EN2 in the chicken. In the mouse, the homologous gene maps on chromosome 5, close to the hemimelic extra-toes mutation Hx. In the case of the NA locus, markers of chromosome 3 were selected because CPPP was mapped on this chromosome. Analysis of individual typings showed a linkage of 5.7 cM (lod score = 13) between the NA locus and ADL0237 in the distal region of chromosome 3q. These results contribute to connecting the former classical map to the molecular genetic map of the chicken, and open the way to the identification of the molecular nature of two developmental mutations of the chicken that are known to occur in many breeds of chickens.     

Parallel Genomic Changes Drive Repeated Evolution of Placentas in Live-Bearing Fish

The evolutionary origin of complex organs challenges empirical study because most organs evolved hundreds of millions of years ago. The placenta of live-bearing fish in the family Poeciliidae represents a unique opportunity to study the evolutionary origin of complex organs, because in this family a placenta evolved at least nine times independently. It is currently unknown whether this repeated evolution is accompanied by similar, repeated, genomic changes in placental species. Here, we compare whole genomes of 26 poeciliid species representing six out of nine independent origins of placentation. Evolutionary rate analysis revealed that the evolution of the placenta coincides with convergent shifts in the evolutionary rate of 78 protein-coding genes, mainly observed in transporter- and vesicle-located genes. Furthermore, differences in sequence conservation showed that placental evolution coincided with similar changes in 76 noncoding regulatory elements, occurring primarily around genes that regulate development. The unexpected high occurrence of GATA simple repeats in the regulatory elements suggests an important function for GATA repeats in developmental gene regulation. The distinction in molecular evolution observed, with protein-coding parallel changes more often found in metabolic and structural pathways, compared with regulatory change more frequently found in developmental pathways, offers a compelling model for complex trait evolution in general: changing the regulation of otherwise highly conserved developmental genes may allow for the evolution of complex traits.