A de novo mutation in CACNA1A is associated with autosomal dominant bovine familial convulsions and ataxia in Angus cattle.

Bovine familial convulsions and ataxia (BFCA) is considered an autosomal dominant syndrome with incomplete penetrance. Nine Angus calves from the same herd were diagnosed with BFCA within days of birth. Necropsy revealed cerebellar and spinal cord lesions associated with the condition. Parentage testing confirmed that all affected calves had a common sire. The sire was then bred to 36 cows across two herds using artificial insemination, producing an additional 14 affected calves. The objective of this investigation was to identify hypothesized dominant genetic variation underlying the condition. Whole-genome sequencing was performed on the sire, six affected and seven unaffected paternal half-sibling calves and combined with data from 135 unrelated controls. The sire and five of the six affected calves were heterozygous for a nonsense variant (Chr7 g.12367906C>T, c.5073C>T, p.Arg1681*) in CACNA1A. The other affected calves (N = 8) were heterozygous for the variant but it was absent in the other unaffected calves (N = 7) and parents of the sire. This variant was also absent in sequence data from over 6500 other cattle obtained via public repositories and collaborator projects. The variant in CACNA1A is expressed in the cerebellum of the ataxic calves as detected in the transcriptome and was not differentially expressed compared with controls. The CACNA1A protein is part of a highly expressed cerebellar calcium voltage gated channel. The nonsense variant is proposed to cause haploinsufficiency, preventing proper transmission of neuronal signals through the channel and resulting in BFCA.

A KRT71 Loss-of-Function Variant Results in Inner Root Sheath Dysplasia and Recessive Congenital Hypotrichosis of Hereford Cattle.

Genodermatoses, such as heritable skin disorders, mostly represent Mendelian conditions. Congenital hypotrichosis (HY) characterize a condition of being born with less hair than normal. The purpose of this study was to characterize the clinicopathological phenotype of a breed-specific non-syndromic form of HY in Hereford cattle and to identify the causative genetic variant for this recessive disorder. Affected calves showed a very short, fine, wooly, kinky and curly coat over all parts of the body, with a major expression in the ears, the inner part of the limbs, and in the thoracic-abdominal region. Histopathology showed a severely altered morphology of the inner root sheath (IRS) of the hair follicle with abnormal Huxley and Henle's layers and severely dysplastic hair shafts. A genome-wide association study revealed an association signal on chromosome 5. Homozygosity mapping in a subset of cases refined the HY locus to a 690 kb critical interval encompassing a cluster of type II keratin encoding genes. Protein-coding exons of six positional candidate genes with known hair or hair follicle function were re-sequenced. This revealed a protein-changing variant in the KRT71 gene that encodes a type II keratin specifically expressed in the IRS of the hair follicle (c.281delTGTGCCCA; p.Met94AsnfsX14). Besides obvious phenocopies, a perfect concordance between the presence of this most likely pathogenic loss-of-function variant located in the head domain of KRT71 and the HY phenotype was found. This recessive KRT71-related form of hypotrichosis provides a novel large animal model for similar human conditions. The results have been incorporated in the Online Mendelian Inheritance in Animals (OMIA) database (OMIA 002114-9913).

A Nonsense Variant in Hephaestin Like 1 (HEPHL1) Is Responsible for Congenital Hypotrichosis in Belted Galloway Cattle.

Genodermatosis such as hair disorders mostly follow a monogenic mode of inheritance. Congenital hypotrichosis (HY) belong to this group of disorders and is characterized by abnormally reduced hair since birth. The purpose of this study was to characterize the clinical phenotype of a breed-specific non-syndromic form of HY in Belted Galloway cattle and to identify the causative genetic variant for this recessive disorder. An affected calf born in Switzerland presented with multiple small to large areas of alopecia on the limbs and on the dorsal part of the head, neck, and back. A genome-wide association study using Swiss and US Belted Galloway cattle encompassing 12 cases and 61 controls revealed an association signal on chromosome 29. Homozygosity mapping in a subset of cases refined the HY locus to a 1.5 Mb critical interval and subsequent Sanger sequencing of protein-coding exons of positional candidate genes revealed a stop gain variant in the HEPHL1 gene that encodes a multi-copper ferroxidase protein so-called hephaestin like 1 (c.1684A>T; p.Lys562*). A perfect concordance between the homozygous presence of this most likely pathogenic loss-of-function variant and the HY phenotype was found. Genotyping of more than 700 purebred Swiss and US Belted Galloway cattle showed the global spread of the mutation. This study provides a molecular test that will permit the avoidance of risk matings by systematic genotyping of relevant breeding animals. This rare recessive HEPHL1-related form of hypotrichosis provides a novel large animal model for similar human conditions. The results have been incorporated in the Online Mendelian Inheritance in Animals (OMIA) database (OMIA 002230-9913).

Genome-wide association study for feed efficiency and growth traits in U.S. beef cattle.

BACKGROUND: Single nucleotide polymorphism (SNP) arrays for domestic cattle have catalyzed the identification of genetic markers associated with complex traits for inclusion in modern breeding and selection programs. Using actual and imputed Illumina 778K genotypes for 3887 U.S. beef cattle from 3 populations (Angus, Hereford, SimAngus), we performed genome-wide association analyses for feed efficiency and growth traits including average daily gain (ADG), dry matter intake (DMI), mid-test metabolic weight (MMWT), and residual feed intake (RFI), with marker-based heritability estimates produced for all traits and populations. RESULTS: Moderate and/or large-effect QTL were detected for all traits in all populations, as jointly defined by the estimated proportion of variance explained (PVE) by marker effects (PVE ≥ 1.0%) and a nominal P-value threshold (P ≤ 5e-05). Lead SNPs with PVE ≥ 2.0% were considered putative evidence of large-effect QTL (n = 52), whereas those with PVE ≥ 1.0% but < 2.0% were considered putative evidence for moderate-effect QTL (n = 35). Identical or proximal lead SNPs associated with ADG, DMI, MMWT, and RFI collectively supported the potential for either pleiotropic QTL, or independent but proximal causal mutations for multiple traits within and between the analyzed populations. Marker-based heritability estimates for all investigated traits ranged from 0.18 to 0.60 using 778K genotypes, or from 0.17 to 0.57 using 50K genotypes (reduced from Illumina 778K HD to Illumina Bovine SNP50). An investigation to determine if QTL detected by 778K analysis could also be detected using 50K genotypes produced variable results, suggesting that 50K analyses were generally insufficient for QTL detection in these populations, and that relevant breeding or selection programs should be based on higher density analyses (imputed or directly ascertained). CONCLUSIONS: Fourteen moderate to large-effect QTL regions which ranged from being physically proximal (lead SNPs ≤ 3Mb) to fully overlapping for RFI, DMI, ADG, and MMWT were detected within and between populations, and included evidence for pleiotropy, proximal but independent causal mutations, and multi-breed QTL. Bovine positional candidate genes for these traits were functionally conserved across vertebrate species.

Precise gene editing paves the way for derivation of Mannheimia haemolytica leukotoxin-resistant cattle.

Signal peptides of membrane proteins are cleaved by signal peptidase once the nascent proteins reach the endoplasmic reticulum. Previously, we reported that, contrary to the paradigm, the signal peptide of ruminant CD18, the ß subunit of ß2 integrins, is not cleaved and hence remains intact on mature CD18 molecules expressed on the surface of ruminant leukocytes. Leukotoxin secreted by Mannheimia (Pasteurella) haemolytica binds to the intact signal peptide and causes cytolysis of ruminant leukocytes, resulting in acute inflammation and lung tissue damage. We also demonstrated that site-directed mutagenesis leading to substitution of cleavage-inhibiting glutamine (Q), at amino acid position 5 upstream of the signal peptide cleavage site, with cleavage-inducing glycine (G) results in the cleavage of the signal peptide and abrogation of leukotoxin-induced cytolysis of target cells. In this proof-of-principle study, we used precise gene editing to induce Q(‒5)G substitution in both alleles of CD18 in bovine fetal fibroblast cells. The gene-edited fibroblasts were used for somatic nuclear transfer and cloning to produce a bovine fetus homozygous for the Q(‒5)G substitution. The leukocyte population of this engineered ruminant expressed CD18 without the signal peptide. More importantly, these leukocytes were absolutely resistant to leukotoxin-induced cytolysis. This report demonstrates the feasibility of developing lines of cattle genetically resistant to M. haemolytica-caused pneumonia, which inflicts an economic loss of over $1 billion to the US cattle industry alone.

QTLs associated with dry matter intake, metabolic mid-test weight, growth and feed efficiency have little overlap across 4 beef cattle studies.

BACKGROUND: The identification of genetic markers associated with complex traits that are expensive to record such as feed intake or feed efficiency would allow these traits to be included in selection programs. To identify large-effect QTL, we performed a series of genome-wide association studies and functional analyses using 50 K and 770 K SNP genotypes scored in 5,133 animals from 4 independent beef cattle populations (Cycle VII, Angus, Hereford and Simmental×Angus) with phenotypes for average daily gain, dry matter intake, metabolic mid-test body weight and residual feed intake. RESULTS: A total of 5, 6, 11 and 10 significant QTL (defined as 1-Mb genome windows with Bonferroni-corrected P-value<0.05) were identified for average daily gain, dry matter intake, metabolic mid-test body weight and residual feed intake, respectively. The identified QTL were population-specific and had little overlap across the 4 populations. The pleiotropic or closely linked QTL on BTA 7 at 23 Mb identified in the Angus population harbours a promising candidate gene ACSL6 (acyl-CoA synthetase long-chain family member 6), and was the largest effect QTL associated with dry matter intake and mid-test body weight explaining 10.39% and 14.25% of the additive genetic variance, respectively. Pleiotropic or closely linked QTL associated with average daily gain and mid-test body weight were detected on BTA 6 at 38 Mb and BTA 7 at 93 Mb confirming previous reports. No QTL for residual feed intake explained more than 2.5% of the additive genetic variance in any population. Marker-based estimates of heritability ranged from 0.21 to 0.49 for residual feed intake across the 4 populations. CONCLUSIONS: This GWAS study, which is the largest performed for feed efficiency and its component traits in beef cattle to date, identified several large-effect QTL that cumulatively explained a significant percentage of additive genetic variance within each population. Differences in the QTL identified among the different populations may be due to differences in power to detect QTL, environmental variation, or differences in the genetic architecture of trait variation among breeds. These results enhance our understanding of the biology of growth, feed intake and utilisation in beef cattle.

Bivariate genome-wide association analysis of the growth and intake components of feed efficiency.

Single nucleotide polymorphisms (SNPs) associated with average daily gain (ADG) and dry matter intake (DMI), two major components of feed efficiency in cattle, were identified in a genome-wide association study (GWAS). Uni- and multi-SNP models were used to describe feed efficiency in a training data set and the results were confirmed in a validation data set. Results from the univariate and bivariate analyses of ADG and DMI, adjusted by the feedlot beef steer maintenance requirements, were compared. The bivariate uni-SNP analysis identified (P-value <0.0001) 11 SNPs, meanwhile the univariate analyses of ADG and DMI identified 8 and 9 SNPs, respectively. Among the six SNPs confirmed in the validation data set, five SNPs were mapped to KDELC2, PHOX2A, and TMEM40. Findings from the uni-SNP models were used to develop highly accurate predictive multi-SNP models in the training data set. Despite the substantially smaller size of the validation data set, the training multi-SNP models had slightly lower predictive ability when applied to the validation data set. Six Gene Ontology molecular functions related to ion transport activity were enriched (P-value <0.001) among the genes associated with the detected SNPs. The findings from this study demonstrate the complementary value of the uni- and multi-SNP models, and univariate and bivariate GWAS analyses. The identified SNPs can be used for genome-enabled improvement of feed efficiency in feedlot beef cattle, and can aid in the design of empirical studies to further confirm the associations.

Single nucleotide polymorphisms and haplotypes associated with feed efficiency in beef cattle.

BACKGROUND: General, breed- and diet-dependent associations between feed efficiency in beef cattle and single nucleotide polymorphisms (SNPs) or haplotypes were identified on a population of 1321 steers using a 50 K SNP panel. Genomic associations with traditional two-step indicators of feed efficiency - residual feed intake (RFI), residual average daily gain (RADG), and residual intake gain (RIG) - were compared to associations with two complementary one-step indicators of feed efficiency: efficiency of intake (EI) and efficiency of gain (EG). Associations uncovered in a training data set were evaluated on independent validation data set. A multi-SNP model was developed to predict feed efficiency. Functional analysis of genes harboring SNPs significantly associated with feed efficiency and network visualization aided in the interpretation of the results. RESULTS: For the five feed efficiency indicators, the numbers of general, breed-dependent, and diet-dependent associations with SNPs (P-value < 0.0001) were 31, 40, and 25, and with haplotypes were six, ten, and nine, respectively. Of these, 20 SNP and six haplotype associations overlapped between RFI and EI, and five SNP and one haplotype associations overlapped between RADG and EG. This result confirms the complementary value of the one and two-step indicators. The multi-SNP models included 89 SNPs and offered a precise prediction of the five feed efficiency indicators. The associations of 17 SNPs and 7 haplotypes with feed efficiency were confirmed on the validation data set. Nine clusters of Gene Ontology and KEGG pathway categories (mean P-value < 0.001) including, 9nucleotide binding; ion transport, phosphorous metabolic process, and the MAPK signaling pathway were overrepresented among the genes harboring the SNPs associated with feed efficiency. CONCLUSIONS: The general SNP associations suggest that a single panel of genomic variants can be used regardless of breed and diet. The breed- and diet-dependent associations between SNPs and feed efficiency suggest that further refinement of variant panels require the consideration of the breed and management practices. The unique genomic variants associated with the one- and two-step indicators suggest that both types of indicators offer complementary description of feed efficiency that can be exploited for genome-enabled selection purposes.

Association of the porcine transforming growth factor beta type I receptor (TGFBR1) gene with growth and carcass traits.

BACKGROUND: Growth and carcass traits are of great economic importance in livestock production. A large number of quantitative trait loci (QTL) have been identified for growth and carcass traits on porcine chromosome one (SSC1). A key positional candidate for this chromosomal region is TGFBR1 (transforming growth factor beta type I receptor). This gene plays a key role in inherited disorders at cardiovascular, craniofacial, neurocognitive, and skeletal development in mammals. RESULTS: In this study, 27 polymorphic SNPs in the porcine TGFBR1 gene were identified on the University of Illinois Yorkshire × Meishan resource population. Three SNPs (SNP3, SNP43, SNP64) representing major polymorphic patterns of the 27 SNPs in F1 and F0 individuals of the Illinois population were selected for analyses of QTL association and genetic diversity. An association analysis for growth and carcass traits was completed using these three representative SNPs in the Illinois population with 298 F2 individuals and a large commercial population of 1008 animals. The results indicate that the TGFBR1 gene polymorphism (SNP64) is significantly associated (p < 0.05) with growth rates including average daily gains between birth and 56 kg (p = 0.049), between 5.5 and 56 kg (p = 0.024), between 35 and 56 kg (p = 0.021). Significant associations (p < 0.05) were also identified between TGFBR1 gene polymorphisms (SNP3/SNP43) and carcass traits including loin-eye-area (p = 0.022) in the Illinois population, and back-fat thickness (p = 0.0009), lean percentage (p = 0.0023) and muscle color (p = 0.021) in the commercial population. These three SNPs were also used to genotype a diverse panel of 130 animals representing 11 pig breeds. Alleles SNP3_T and SNP43_G were fixed in Pietrain and Sinclair pig breeds. SNP64_G allele was uniquely identified in Chinese Meishan pigs. Strong evidence of association (p < 0.01) between both SNP3 and SNP64 alleles and reproductive traits including gestation length and number of corpora lutea were also observed in the Illinois population. CONCLUSION: This study gives the first evidence of association between the porcine TGFBR1 gene and traits of economic importance and provides support for using TGFBR1 markers for pig breeding and selection programs. The genetic diversities in different pig breeds would be helpful to understand the genetic background and migration of the porcine TGFBR1 gene.

Integration of physical and genetic maps in apple confirms whole-genome and segmental duplications in the apple genome.

A total of 355 simple sequence repeat (SSR) markers were developed, based on expressed sequence tag (EST) and bacterial artificial chromosome (BAC)-end sequence databases, and successfully used to construct an SSR-based genetic linkage map of the apple. The consensus linkage map spanned 1143 cM, with an average density of 2.5 cM per marker. Newly developed SSR markers along with 279 SSR markers previously published by the HiDRAS project were further used to integrate physical and genetic maps of the apple using a PCR-based BAC library screening approach. A total of 470 contigs were unambiguously anchored onto all 17 linkage groups of the apple genome, and 158 contigs contained two or more molecular markers. The genetically mapped contigs spanned ∼421 Mb in cumulative physical length, representing 60.0% of the genome. The sizes of anchored contigs ranged from 97 kb to 4.0 Mb, with an average of 995 kb. The average physical length of anchored contigs on each linkage group was ∼24.8 Mb, ranging from 17.0 Mb to 37.73 Mb. Using BAC DNA as templates, PCR screening of the BAC library amplified fragments of highly homologous sequences from homoeologous chromosomes. Upon integrating physical and genetic maps of the apple, the presence of not only homoeologous chromosome pairs, but also of multiple locus markers mapped to adjacent sites on the same chromosome was detected. These findings demonstrated the presence of both genome-wide and segmental duplications in the apple genome and provided further insights into the complex polyploid ancestral origin of the apple.

Cell cycle and aging, morphogenesis, and response to stimuli genes are individualized biomarkers of glioblastoma progression and survival.

BACKGROUND: Glioblastoma is a complex multifactorial disorder that has swift and devastating consequences. Few genes have been consistently identified as prognostic biomarkers of glioblastoma survival. The goal of this study was to identify general and clinical-dependent biomarker genes and biological processes of three complementary events: lifetime, overall and progression-free glioblastoma survival. METHODS: A novel analytical strategy was developed to identify general associations between the biomarkers and glioblastoma, and associations that depend on cohort groups, such as race, gender, and therapy. Gene network inference, cross-validation and functional analyses further supported the identified biomarkers. RESULTS: A total of 61, 47 and 60 gene expression profiles were significantly associated with lifetime, overall, and progression-free survival, respectively. The vast majority of these genes have been previously reported to be associated with glioblastoma (35, 24, and 35 genes, respectively) or with other cancers (10, 19, and 15 genes, respectively) and the rest (16, 4, and 10 genes, respectively) are novel associations. Pik3r1, E2f3, Akr1c3, Csf1, Jag2, Plcg1, Rpl37a, Sod2, Topors, Hras, Mdm2, Camk2g, Fstl1, Il13ra1, Mtap and Tp53 were associated with multiple survival events.Most genes (from 90 to 96%) were associated with survival in a general or cohort-independent manner and thus the same trend is observed across all clinical levels studied. The most extreme associations between profiles and survival were observed for Syne1, Pdcd4, Ighg1, Tgfa, Pla2g7, and Paics. Several genes were found to have a cohort-dependent association with survival and these associations are the basis for individualized prognostic and gene-based therapies. C2, Egfr, Prkcb, Igf2bp3, and Gdf10 had gender-dependent associations; Sox10, Rps20, Rab31, and Vav3 had race-dependent associations; Chi3l1, Prkcb, Polr2d, and Apool had therapy-dependent associations. Biological processes associated glioblastoma survival included morphogenesis, cell cycle, aging, response to stimuli, and programmed cell death. CONCLUSIONS: Known biomarkers of glioblastoma survival were confirmed, and new general and clinical-dependent gene profiles were uncovered. The comparison of biomarkers across glioblastoma phases and functional analyses offered insights into the role of genes. These findings support the development of more accurate and personalized prognostic tools and gene-based therapies that improve the survival and quality of life of individuals afflicted by glioblastoma multiforme.

A deletion mutation in bovine SLC4A2 is associated with osteopetrosis in Red Angus cattle.

BACKGROUND: Osteopetrosis is a skeletal disorder of humans and animals characterized by the formation of overly dense bones, resulting from a deficiency in the number and/or function of bone-resorbing osteoclast cells. In cattle, osteopetrosis can either be induced during gestation by viral infection of the dam, or inherited as a recessive defect. Genetically affected calves are typically aborted late in gestation, display skull deformities and exhibit a marked reduction of osteoclasts. Although mutations in several genes are associated with osteopetrosis in humans and mice, the genetic basis of the cattle disorder was previously unknown. RESULTS: We have conducted a whole-genome association analysis to identify the mutation responsible for inherited osteopetrosis in Red Angus cattle. Analysis of >54,000 SNP genotypes for each of seven affected calves and nine control animals localized the defective gene to the telomeric end of bovine chromosome 4 (BTA4). Homozygosity analysis refined the interval to a 3.4-Mb region containing the SLC4A2 gene, encoding an anion exchanger protein necessary for proper osteoclast function. Examination of SLC4A2 from normal and affected animals revealed a approximately 2.8-kb deletion mutation in affected calves that encompasses exon 2 and nearly half of exon 3, predicted to prevent normal protein function. Analysis of RNA from a proven heterozygous individual confirmed the presence of transcripts lacking exons 2 and 3, in addition to normal transcripts. Genotyping of additional animals demonstrated complete concordance of the homozygous deletion genotype with the osteopetrosis phenotype. Histological examination of affected tissues revealed scarce, morphologically abnormal osteoclasts displaying evidence of apoptosis. CONCLUSIONS: These results indicate that a deletion mutation within bovine SLC4A2 is associated with osteopetrosis in Red Angus cattle. Loss of SLC4A2 function appears to induce premature cell death, and likely results in cytoplasmic alkalinization of osteoclasts which, in turn, may disrupt acidification of resorption lacunae.

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.

An integrated RH map of porcine chromosome 10.

BACKGROUND: Whole genome radiation hybrid (WG-RH) maps serve as "scaffolds" to significantly improve the orientation of small bacterial artificial chromosome (BAC) contigs, order genes within the contigs and assist assembly of a sequence-ready map for virtually any species. Here, we report the construction of a porcine: human comparative map for pig (Sus scrofa) chromosome 10 (SSC10) using the IMNpRH2(12,000-rad) porcine WG-RH panel, integrated with the IMpRH(7000-rad) WG-RH, genetic and BAC fingerprinted contig (FPC) maps. RESULTS: Map vectors from the IMNpRH2(12,000-rad) and IMpRH(7,000-rad) panels were merged to construct parallel framework (FW) maps, within which FW markers common to both panels have an identical order. This strategy reduced map discrepancies between the two panels and significantly improved map accuracy. A total of 216 markers, including 50 microsatellites (MSs), 97 genes and ESTs, and 69 BAC end sequences (BESs), were ordered within two linkage groups at two point (2 pt) LOD score of 8. One linkage group covers SSC10p with accumulated map distances of 738.2 cR(7,000) and 1814.5 cR(12,000), respectively. The second group covers SSC10q at map distances of 1336.9 cR(7,000) and 3353.6 cR(12,000), yielding an overall average map resolution of 16.4 kb/cR(12,000) or 393.5 kb per marker on SSC10. This represents an approximately 2.5-fold increase in map resolution over the IMpRH(7,000-rad) panel. Based on 127 porcine markers that have homologous sequences in the human genome, a detailed comparative map between SSC10 and human (Homo sapiens) chromosome (HSA) 1, 9 and 10 was built. CONCLUSION: This initial comparative RH map of SSC10 refines the syntenic regions between SSC10 and HSA1, 9 and 10. It integrates the IMNpRH2(12,000-rad) and IMpRH(7,000-rad), genetic and BAC FPC maps and provides a scaffold to close potential gaps between contigs prior to genome sequencing and assembly. This map is also useful in fine mapping of QTLs on SSC10.

Diversification of porcine MHC class II genes: evidence for selective advantage.

The major histocompatibility complex (MHC) is an immunological gene-dense region of high diversity in mammalian species. Sus scrofa was domesticated by at least six independent events over Eurasia during the Holocene period. It has been hypothesized that the level and distribution of MHC variation in pig populations reflect genetic selection and environmental influences. In an effort to define the complexity of MHC polymorphisms and the role of selection in the generation of class II gene diversity (DQB, DRB1, and pseudogene PsiDRB3), DNA from globally distributed unrelated domestic pigs of European and Asian origins and a Suidae out-group was analyzed. The number of pseudogene alleles identified (PsiDRB3 33) was greater than those found in the expressed genes (DQB 20 and DRB1 23) but the level of observed heterozygosity (PsiDRB3 0.452, DQB 0.732, and DRB1 0.767) and sequence diversity (PsiDRB3 0.029, DQB 0.062, and DRB1 0.074) were significantly lower in the pseudogene, respectively. The substitution ratios reflected an excess of d (N) (DQB 1.476, DRB1 1.724, and PsiDRB3 0.508) and the persistence of expressed gene alleles suggesting the influence of balancing selection, while the pseudogene was undergoing purifying selection. The lack of a clear MHC phylogeographic tree, coupled with close genetic distances observed between the European and Asian populations (DQB 0.047 and DRB1 0.063) suggested that unlike observations using mtDNA, the MHC diversity lacks phylogeographic structure and appears to be globally uniform. Taken together, these results suggest that, despite regional differences in selective breeding and environments, no skewing of MHC diversity has occurred.

BAC-end sequence-based SNPs and Bin mapping for rapid integration of physical and genetic maps in apple.

A genome-wide BAC physical map of the apple, Malus x domestica Borkh., has been recently developed. Here, we report on integrating the physical and genetic maps of the apple using a SNP-based approach in conjunction with bin mapping. Briefly, BAC clones located at ends of BAC contigs were selected, and sequenced at both ends. The BAC end sequences (BESs) were used to identify candidate SNPs. Subsequently, these candidate SNPs were genetically mapped using a bin mapping strategy for the purpose of mapping the physical onto the genetic map. Using this approach, 52 (23%) out of 228 BESs tested were successfully exploited to develop SNPs. These SNPs anchored 51 contigs, spanning approximately 37 Mb in cumulative physical length, onto 14 linkage groups. The reliability of the integration of the physical and genetic maps using this SNP-based strategy is described, and the results confirm the feasibility of this approach to construct an integrated physical and genetic maps for apple.

Piggy-BACing the human genome I: constructing a porcine BAC physical map through comparative genomics.

Availability of the human genome sequence and high similarity between humans and pigs at the molecular level provides an opportunity to use a comparative mapping approach to piggy-BAC the human genome. In order to advance the pig genome sequencing initiative, sequence similarity between large-scale porcine BAC-end sequences (BESs) and human genome sequence was used to construct a comparatively-anchored porcine physical map that is a first step towards sequencing the pig genome. A total of 50,300 porcine BAC clones were end-sequenced, yielding 76,906 BESs after trimming with an average read length of 538 bp. To anchor the porcine BACs on the human genome, these BESs were subjected to BLAST analysis using the human draft sequence, revealing 31.5% significant hits (E < e(-5)). Both genic and non-genic regions of homology contributed to the alignments between the human and porcine genomes. Porcine BESs with unique homology matches within the human genome provided a source of markers spaced approximately 70 to 300 kb along each human chromosome. In order to evaluate the utility of piggy-BACing human genome sequences, and confirm predictions of orthology, 193 evenly spaced BESs with similarity to HSA3 and HSA21 were selected and then utilized for developing a high-resolution (1.22 Mb) comparative radiation hybrid map of SSC13 that represents a fusion of HSA3 and HSA21. Resulting RH mapping of SSC13 covers 99% and 97% of HSA3 and HSA21, respectively. Seven evolutionary conserved blocks were identified including six on HSA3 and a single syntenic block corresponding to HSA21. The strategy of piggy-BACing the human genome described in this study demonstrates that through a directed, targeted comparative genomics approach construction of a high-resolution anchored physical map of the pig genome can be achieved. This map supports the selection of BACs to construct a minimal tiling path for genome sequencing and targeted gap filling. Moreover, this approach is highly relevant to other genome sequencing projects.

Fine-mapping of a QTL influencing pork tenderness on porcine chromosome 2.

BACKGROUND: In a previous study, a quantitative trait locus (QTL) exhibiting large effects on both Instron shear force and taste panel tenderness was detected within the Illinois Meat Quality Pedigree (IMQP). This QTL mapped to the q arm of porcine chromosome 2 (SSC2q). Comparative analysis of SSC2q indicates that it is orthologous to a segment of human chromosome 5 (HSA5) containing a strong positional candidate gene, calpastatin (CAST). CAST polymorphisms have recently been shown to be associated with meat quality characteristics; however, the possible involvement of other genes and/or molecular variation in this region cannot be excluded, thus requiring fine-mapping of the QTL. RESULTS: Recent advances in porcine genome resources, including high-resolution radiation hybrid and bacterial artificial chromosome (BAC) physical maps, were utilized for development of novel informative markers. Marker density in the ~30-Mb region surrounding the most likely QTL position was increased by addition of eighteen new microsatellite markers, including nine publicly-available and nine novel markers. Two newly-developed markers were derived from a porcine BAC clone containing the CAST gene. Refinement of the QTL position was achieved through linkage and haplotype analyses. Within-family linkage analyses revealed at least two families segregating for a highly-significant QTL in strong positional agreement with CAST markers. A combined analysis of these two families yielded QTL intervals of 36 cM and 7 cM for Instron shear force and taste panel tenderness, respectively, while haplotype analyses suggested further refinement to a 1.8 cM interval containing CAST markers. The presence of additional tenderness QTL on SSC2q was also suggested. CONCLUSION: These results reinforce CAST as a strong positional candidate. Further analysis of CAST molecular variation within the IMQP F1 boars should enhance understanding of the molecular basis of pork tenderness, and thus allow for genetic improvement of pork products. Furthermore, additional resources have been generated for the targeted investigation of other putative QTL on SSC2q, which may lead to further advancements in pork quality.

Characterization of the porcine ATM gene: towards the generation of a novel non-murine animal model for Ataxia-Telangiectasia.

Ataxia-Telangiectasia (A-T) is a genetic disorder causing cerebellar degeneration, immune deficiency, cancer predisposition, chromosomal instability and radiation sensitivity. Among the mutations responsible for A-T, 85% represent truncating mutations that result in the production of shorter, highly unstable forms of ATM (AT-mutated) protein leading to a null ATM phenotype. Several ATM-deficient mice have been created however none reflects the extent of neurological degeneration observed in humans. In an attempt to identify an alternative animal model, we have characterized the porcine ortholog of ATM (pATM). When compared to the human ATM (hATM), the pATM showed a high level of homology in the coding region, particularly in the regions coding for functional domains, and had extensive alternative splicing of the 5'UTR, characteristic for the human ATM mRNA. Six different 5'UTRs resulting from alternative splicing of the first three exons were identified. The porcine 5'UTRs varied in size, had multiple ATG codons and different secondary structures, supporting the possibility of complex transcriptional regulation. Three of the six transcripts demonstrated alternative splicing of exon 3, the first putative coding exon, altering the translation start and giving rise to a putative protein lacking the N-terminus substrate binding domain (82-89 aa) involved in activation of human p53 and BRCA1 pathways. Real time-PCR analysis revealed variable expression levels of total ATM transcripts in individual tissues. Although each splice variant was ubiquitously expressed among the tissues, differences in the relative abundances of specific 5'UTRs were observed. The extensive alternative splicing of the pATM gene resembles the complex splicing observed in the hATM and could provide insights for differences observed between mice and humans with regards to the onset of A-T. Thus, the pig may provide a more relevant clinical model of A-T.

Transcriptome profiling of the small intestinal epithelium in germfree versus conventional piglets.

BACKGROUND: To gain insight into host-microbe interactions in a piglet model, a functional genomics approach was used to address the working hypothesis that transcriptionally regulated genes associated with promoting epithelial barrier function are activated as a defensive response to the intestinal microbiota. Cesarean-derived germfree (GF) newborn piglets were colonized with adult swine feces, and villus and crypt epithelial cell transcriptomes from colonized and GF neonatal piglets were compared using laser-capture microdissection and high-density porcine oligonucleotide microarray technology. RESULTS: Consistent with our hypothesis, resident microbiota induced the expression of genes contributing to intestinal epithelial cell turnover, mucus biosynthesis, and priming of the immune system. Furthermore, differential expression of genes associated with antigen presentation (pan SLA class I, B2M, TAP1 and TAPBP) demonstrated that microbiota induced immune responses using a distinct regulatory mechanism common for these genes. Specifically, gene network analysis revealed that microbial colonization activated both type I (IFNAR) and type II (IFNGR) interferon receptor mediated signaling cascades leading to enhanced expression of signal transducer and activator of transcription 1 (STAT1), STAT2 and IFN regulatory factor 7 (IRF7) transcription factors and the induction of IFN-inducible genes as a reflection of intestinal epithelial inflammation. In addition, activated RNA expression of NF-kappa-B inhibitor alpha (NFkappaBIA; a.k.a I-kappa-B-alpha, IKBalpha) and toll interacting protein (TOLLIP), both inhibitors of inflammation, along with downregulated expression of the immunoregulatory transcription factor GATA binding protein-1 (GATA1) is consistent with the maintenance of intestinal homeostasis. CONCLUSION: This study supports the concept that the intestinal epithelium has evolved to maintain a physiological state of inflammation with respect to continuous microbial exposure, which serves to sustain a tight intestinal barrier while preventing overt inflammatory responses that would compromise barrier function.