Design and development of exome capture sequencing for the domestic pig (Sus scrofa)

Background

The domestic pig (Sus scrofa) is both an important livestock species and a model for biomedical research. Exome sequencing has accelerated identification of protein-coding variants underlying phenotypic traits in human and mouse. We aimed to develop and validate a similar resource for the pig.

Results

We developed probe sets to capture pig exonic sequences based upon the current Ensembl pig gene annotation supplemented with mapped expressed sequence tags (ESTs) and demonstrated proof-of-principle capture and sequencing of the pig exome in 96 pigs, encompassing 24 capture experiments. For most of the samples at least 10x sequence coverage was achieved for more than 90% of the target bases. Bioinformatic analysis of the data revealed over 236,000 high confidence predicted SNPs and over 28,000 predicted indels.

Conclusions

We have achieved coverage statistics similar to those seen with commercially available human and mouse exome kits. Exome capture in pigs provides a tool to identify coding region variation associated with production traits, including loss of function mutations which may explain embryonic and neonatal losses, and to improve genomic assemblies in the vicinity of protein coding genes in the pig.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-550) contains supplementary material, which is available to authorized users.     

Similarities and Differences between Exome Sequences Found in a Variety of Tissues from the Same Individual

DNA is the most stable nucleic acid and most important store of genetic information. DNA sequences are conserved in virtually all the cells of a multicellular organism. To analyze the sequences of various individuals with distinct pathological disorders, DNA is routinely isolated from blood, independently of the tissue that is the target of the disease. This approach has proven useful for the identification of familial diseases where mutations are present in parental germinal cells. With the capacity to compare DNA sequences from distinct tissues or cells, present technology can be used to study whether DNA sequences in tissues are invariant. Here we explored the presence of specific SNVs (Single Nucleotide Variations) in various tissues of the same individual. We tested for the presence of tissue-specific exonic SNVs, taking blood exome as a control. We analyzed the chromosomal location of these SNVs. The number of SNVs per chromosome was found not to depend on chromosome length, but mainly on the number of protein-coding genes per chromosome. Although similar but not identical patterns of chromosomal distribution of tissue-specific SNVs were found, clear differences were detected. This observation supports the notion that each tissue has a specific SNV exome signature.     

Copy number variants in Italian Large White pigs detected using high‐density single nucleotide polymorphisms and their association with back fat thickness

The aim of this study was to identify copy number variants (CNVs) in Italian Large White pigs and test them for association with back fat thickness (BFT). Within a population of 12 000 performance‐tested pigs, two groups of animals with extreme and divergent BFT estimated breeding values (EBVs; 147 with negative and 150 with positive EBVs) were genotyped with the Illumina Porcine SNP60 BeadChip. CNVs were detected with penncnv software. We identified a total of 4146 CNV events in 170 copy number variation regions (CNVRs) located on 15 porcine autosomes. Validation of detected CNVRs was carried out (i) by comparing CNVRs already detected by other studies and (ii) by semiquantitative fluorescent multiplex (SQFM) PCR of a few CNVRs. Most of CNVRs detected in Italian Large White pigs (71.2%) were already reported in other pig breeds/populations, and 82.1% of the CNV events detected by penncnv were confirmed by SQFM PCR. For each CNVR, we compared the occurrence of CNV events between the pigs of the high and low BFT EBV tails. Sixteen regions showed significance at P < 0.10, and seven were significant at P < 0.05 but were not significant after Bonferroni correction (Fisher's exact test). These results indicated that CNVs could explain a limited fraction of the genetic variability of fat deposition in Italian Large White pigs. However, it was interesting to note that one of these CNVRs encompassed the ZPLD1 gene. In humans, a rare CNV event including this gene is associated with obesity. Studies identifying CNVs in pigs could assist in elucidating the genetic mechanisms underlying human obesity.     

Association study highlights the influence of ELOVL fatty acid elongase 6 gene region on backfat fatty acid composition in Large White pig breed

Dietary fatty acid (FA) composition has an impact on human health. There is an increasing request from consumers for healthier food and pork industry must respond to it without worsening performance and the technological properties of pork products. The inclusion of genetic markers for carcass FA composition in pig selection schemes could be a useful tool to reach the right balance between unsaturated and saturated FAs to satisfy market demands. With the aim of finding genomic regions associated with porcine backfat FA composition, a genome-wide association study was performed on 798 Italian Large White pigs genotyped using Illumina PorcineSNP60 k. The strongest associations with backfat contents of palmitic, palmitoleic, oleic, medium-chain and long-chain FAs were found for the Sus scrofa chromosome (SSC) 8 region located at 119 to 122 Mb, where the gene ELOVL FA elongase 6 is mapped. Palmitic, palmitoleic, stearic and oleic acid contents were also found associated with SSC14, in particular with the genomic region at 121 to 124 Mb, where stearoyl-CoA desaturase Δ9 gene lies. On the other hand, the genomic regions associated with backfat contents of arachidic, arachidonic, n-6 and n-3 FAs showed to harbour mainly genes involved in dietary lipids and carbohydrates digestion, absorption and utilisation. To our knowledge, this is the first study performed in Large White pigs identifying markers and genomic regions associated with backfat FA composition. The results validate in Large White some associations previously detected in other pig breeds and indicate the involvement of distinct metabolic pathways in the deposition pattern of essential and non-essential FAs.     

A genome-wide association study of production traits in a commercial population of Large White pigs: evidence of haplotypes affecting meat quality

Background

Numerous quantitative trait loci (QTL) have been detected in pigs over the past 20 years using microsatellite markers. However, due to the low density of these markers, the accuracy of QTL location has generally been poor. Since 2009, the dense genome coverage provided by the Illumina PorcineSNP60 BeadChip has made it possible to more accurately map QTL using genome-wide association studies (GWAS). Our objective was to perform high-density GWAS in order to identify genomic regions and corresponding haplotypes associated with production traits in a French Large White population of pigs.

Methods

Animals (385 Large White pigs from 106 sires) were genotyped using the PorcineSNP60 BeadChip and evaluated for 19 traits related to feed intake, growth, carcass composition and meat quality. Of the 64 432 SNPs on the chip, 44 412 were used for GWAS with an animal mixed model that included a regression coefficient for the tested SNPs and a genomic kinship matrix. SNP haplotype effects in QTL regions were then tested for association with phenotypes following phase reconstruction based on the Sscrofa10.2 pig genome assembly.

Results

Twenty-three QTL regions were identified on autosomes and their effects ranged from 0.25 to 0.75 phenotypic standard deviation units for feed intake and feed efficiency (four QTL), carcass (12 QTL) and meat quality traits (seven QTL). The 10 most significant QTL regions had effects on carcass (chromosomes 7, 10, 16, 17 and 18) and meat quality traits (two regions on chromosome 1 and one region on chromosomes 8, 9 and 13). Thirteen of the 23 QTL regions had not been previously described. A haplotype block of 183 kb on chromosome 1 (six SNPs) was identified and displayed three distinct haplotypes with significant (0.0001 < P < 0.03) associations with all evaluated meat quality traits.

Conclusions

GWAS analyses with the PorcineSNP60 BeadChip enabled the detection of 23 QTL regions that affect feed consumption, carcass and meat quality traits in a LW population, of which 13 were novel QTL. The proportionally larger number of QTL found for meat quality traits suggests a specific opportunity for improving these traits in the pig by genomic selection.     

Development of Genetic Markers in Eucalyptus Species by Target Enrichment and Exome Sequencing

The advent of next-generation sequencing has facilitated large-scale discovery, validation and assessment of genetic markers for high density genotyping. The present study was undertaken to identify markers in genes supposedly related to wood property traits in three Eucalyptus species. Ninety four genes involved in xylogenesis were selected for hybridization probe based nuclear genomic DNA target enrichment and exome sequencing. Genomic DNA was isolated from the leaf tissues and used for on-array probe hybridization followed by Illumina sequencing. The raw sequence reads were trimmed and high-quality reads were mapped to the E. grandis reference sequence and the presence of single nucleotide variants (SNVs) and insertions/ deletions (InDels) were identified across the three species. The average read coverage was 216X and a total of 2294 SNVs and 479 InDels were discovered in E. camaldulensis, 2383 SNVs and 518 InDels in E. tereticornis, and 1228 SNVs and 409 InDels in E. grandis. Additionally, SNV calling and InDel detection were conducted in pair-wise comparisons of E. tereticornis vs. E. grandis, E. camaldulensis vs. E. tereticornis and E. camaldulensis vs. E. grandis. This study presents an efficient and high throughput method on development of genetic markers for family– based QTL and association analysis in Eucalyptus.     

Genome-wide detection of copy number variants in European autochthonous and commercial pig breeds by whole-genome sequencing of DNA pools identified breed-characterising copy number states

In this study, we identified copy number variants (CNVs) in 19 European autochthonous pig breeds and in two commercial breeds (Italian Large White and Italian Duroc) that represent important genetic resources for this species. The genome of 725 pigs was sequenced using a breed-specific DNA pooling approach (30–35 animals per pool) obtaining an average depth per pool of 42×. This approach maximised CNV discovery as well as the related copy number states characterising, on average, the analysed breeds. By mining more than 17.5 billion reads, we identified a total of 9592 CNVs (~683 CNVs per breed) and 3710 CNV regions (CNVRs; 1.15% of the reference pig genome), with an average of 77 CNVRs per breed that were considered as private. A few CNVRs were analysed in more detail, together with other information derived from sequencing data. For example, the CNVR encompassing the KIT gene was associated with coat colour phenotypes in the analysed breeds, confirming the role of the multiple copies in determining breed-specific coat colours. The CNVR covering the MSRB3 gene was associated with ear size in most breeds. The CNVRs affecting the ELOVL6 and ZNF622 genes were private features observed in the Lithuanian Indigenous Wattle and in the Turopolje pig breeds respectively. Overall, the genome variability unravelled here can explain part of the genetic diversity among breeds and might contribute to explain their origin, history and adaptation to a variety of production systems.     

Whole Exome Re-Sequencing Implicates CCDC38 and Cilia Structure and Function in Resistance to Smoking Related Airflow Obstruction

Chronic obstructive pulmonary disease (COPD) is a leading cause of global morbidity and mortality and, whilst smoking remains the single most important risk factor, COPD risk is heritable. Of 26 independent genomic regions showing association with lung function in genome-wide association studies, eleven have been reported to show association with airflow obstruction. Although the main risk factor for COPD is smoking, some individuals are observed to have a high forced expired volume in 1 second (FEV₁) despite many years of heavy smoking. We hypothesised that these “resistant smokers” may harbour variants which protect against lung function decline caused by smoking and provide insight into the genetic determinants of lung health. We undertook whole exome re-sequencing of 100 heavy smokers who had healthy lung function given their age, sex, height and smoking history and applied three complementary approaches to explore the genetic architecture of smoking resistance. Firstly, we identified novel functional variants in the “resistant smokers” and looked for enrichment of these novel variants within biological pathways. Secondly, we undertook association testing of all exonic variants individually with two independent control sets. Thirdly, we undertook gene-based association testing of all exonic variants. Our strongest signal of association with smoking resistance for a non-synonymous SNP was for rs10859974 (P = 2.34×10⁻⁴) in CCDC38, a gene which has previously been reported to show association with FEV₁/FVC, and we demonstrate moderate expression of CCDC38 in bronchial epithelial cells. We identified an enrichment of novel putatively functional variants in genes related to cilia structure and function in resistant smokers. Ciliary function abnormalities are known to be associated with both smoking and reduced mucociliary clearance in patients with COPD. We suggest that genetic influences on the development or function of cilia in the bronchial epithelium may affect growth of cilia or the extent of damage caused by tobacco smoke.     

MuPIT interactive: webserver for mapping variant positions to annotated, interactive 3D structures

Mutation position imaging toolbox (MuPIT) interactive is a browser-based application for single-nucleotide variants (SNVs), which automatically maps the genomic coordinates of SNVs onto the coordinates of available three-dimensional (3D) protein structures. The application is designed for interactive browser-based visualization of the putative functional relevance of SNVs by biologists who are not necessarily experts either in bioinformatics or protein structure. Users may submit batches of several thousand SNVs and review all protein structures that cover the SNVs, including available functional annotations such as binding sites, mutagenesis experiments, and common polymorphisms. Multiple SNVs may be mapped onto each structure, enabling 3D visualization of SNV clusters and their relationship to functionally annotated positions. We illustrate the utility of MuPIT interactive in rationalizing the impact of selected polymorphisms in the PharmGKB database, somatic mutations identified in the Cancer Genome Atlas study of invasive breast carcinomas, and rare variants identified in the exome sequencing project. MuPIT interactive is freely available for non-profit use at http://mupit.icm.jhu.edu.     

A single nucleotide polymorphism in the porcine cathepsin K (<Emphasis Type="Italic">CTSK</Emphasis>) gene is associated with back fat thickness and production traits in Italian Duroc pigs

Cathepsin K (CTSK) was selected as a candidate gene for fat deposition in pigs because recently, in human and mouse, it was shown that this lysosomal proteinase is an obesity marker. A single nucleotide polymorphism (SNP) was identified in intron 4 of the porcine CTSK gene (g.15G>A; FM209043). Allele frequencies of this polymorphism were analysed in seven pig breeds. Radiation hybrid mapping confirmed the localization of CTSK to porcine chromosome 4, close to the FAT1 QTL region. Three populations of pigs (one Italian Large White and two Italian Duroc groups of pigs) were selected for association analysis. In the Italian Large White breed the g.15G>A SNP was not informative. Association analysis including all Italian Duroc pigs showed that the CTSK marker was associated with back fat thickness and lean cuts (P < 0.01), and average daily gain and feed:gain ratio (P < 0.05) estimated breeding values.     

Metabolomics evidences plasma and serum biomarkers differentiating two heavy pig breeds

In pigs, many production traits are known to vary among breeds or lines. These traits can be considered end phenotypes or external traits as they are the final results of complex biological interactions and processes whose fine biological mechanisms are still largely unknown. This study was designed to compare plasma and serum metabolomic profiles between animals of two heavy pig breeds (12 Italian Large White and 12 Italian Duroc), testing indirectly the hypothesis that different genetic backgrounds might be the determining factors of differences observed on the level of metabolites in the analyzed biofluids between breeds. We used a targeted metabolomic approach based on mass spectrometric detection of about 180 metabolites and applied a statistical validation pipeline to identify differences in the metabolomic profiles of the two heavy pig breeds. Blood samples were collected after jugulation at the slaughterhouse and prepared for metabolomics analysis that was carried out using the Biocrates AbsoluteIDQ p180 Kit, covering five different biochemical classes: glycerophospholipids, amino acids, biogenic amines, hexoses and acylcarnitines. A statistical pipeline that included the selection of the most relevant metabolites differentiating the two breeds by sparse Partial Least Squares Discriminant Analysis (sPLS-DA) was coupled with a stability test and significance test determined with leave one out and permutation procedures. sPLS-DA plots clearly separated the pigs of the two investigated breeds. A few metabolites (a total of five metabolites considering the two biofluids) involved in key metabolic pathways largely contributed to these differences between breeds. In particular, a higher level of the sphingomyelins SM (OH) C14:1 (both in plasma and serum), SM (OH) C16:1 (in serum) and SM C16:0 (in serum) were observed in Italian Duroc than in Italian Large White pigs and the inverse was for the biogenic amine kynurenine (in plasma). The level of another biogenic amine (acetylornithine) was higher in Italian Large White than in Italian Duroc pigs in both analysed biofluids. These results provided biomarkers that could be important to understand the biological differences between these two heavy pig breeds. In particular, according to the functional role played by sphingomyelins in obesity-induced inflammatory responses, it could be possible to speculate that a higher level of sphingomyelins in Italian Duroc might be related to the higher interrmuscular fat deposition of this breed compared with the Italian Large White. Additional studies will be needed to evaluate the relevance of these biomarkers for practical applications in pig breeding and nutrition.     

Comparison and evaluation of two exome capture kits and sequencing platforms for variant calling

Background

To promote the clinical application of next-generation sequencing, it is important to obtain accurate and consistent variants of target genomic regions at low cost. Ion Proton, the latest updated semiconductor-based sequencing instrument from Life Technologies, is designed to provide investigators with an inexpensive platform for human whole exome sequencing that achieves a rapid turnaround time. However, few studies have comprehensively compared and evaluated the accuracy of variant calling between Ion Proton and Illumina sequencing platforms such as HiSeq 2000, which is the most popular sequencing platform for the human genome. The Ion Proton sequencer combined with the Ion TargetSeq™ Exome Enrichment Kit together make up TargetSeq-Proton, whereas SureSelect-Hiseq is based on the Agilent SureSelect Human All Exon v4 Kit and the HiSeq 2000 sequencer.

Results

Here, we sequenced exonic DNA from four human blood samples using both TargetSeq-Proton and SureSelect-HiSeq. We then called variants in the exonic regions that overlapped between the two exome capture kits (33.6 Mb). The rates of shared variant loci called by two sequencing platforms were from 68.0 to 75.3 % in four samples, whereas the concordance of co-detected variant loci reached 99 %. Sanger sequencing validation revealed that the validated rate of concordant single nucleotide polymorphisms (SNPs) (91.5 %) was higher than the SNPs specific to TargetSeq-Proton (60.0 %) or specific to SureSelect-HiSeq (88.3 %). With regard to 1-bp small insertions and deletions (InDels), the Sanger sequencing validated rates of concordant variants (100.0 %) and SureSelect-HiSeq-specific (89.6 %) were higher than those of TargetSeq-Proton-specific (15.8 %).

Conclusions

In the sequencing of exonic regions, a combination of using of two sequencing strategies (SureSelect-HiSeq and TargetSeq-Proton) increased the variant calling specificity for concordant variant loci and the sensitivity for variant loci called by any one platform. However, for the sequencing of platform-specific variants, the accuracy of variant calling by HiSeq 2000 was higher than that of Ion Proton, specifically for the InDel detection. Moreover, the variant calling software also influences the detection of SNPs and, specifically, InDels in Ion Proton exome sequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1796-6) contains supplementary material, which is available to authorized users.     

New SNP of the porcine Perilipin 2 (<Emphasis Type="Italic">PLIN2</Emphasis>) gene,association with carcass traits and expression analysis in skeletal muscle

PLIN2 (perilipin 2) is a cytosolic protein that promotes the formation and stabilization of the intracellular lipid droplets, organelles involved in the storage of lipid depots. Porcine PLIN2 gene represents a biological and positional candidate for fat deposition, a polygenic trait that affects carcass and meat quality. The aim of the present study was to screen PLIN2 gene for polymorphisms, to evaluate the association with carcass quality traits, and to investigate the gene expression in skeletal muscle. Six new single nucleotide polymorphisms (SNP) were detected by sequencing 32 samples from five pig breeds (Italian Large White, Italian Duroc, Italian Landrace, Belgian Landrace, Pietrain). Two SNP localized in introns, two in the 3′-untranslated region (UTR), and two missense SNP were found in exons. A 3′-UTR mutation (GU461317:g.98G>A), genotyped in 290 Italian Duroc pigs by High Resolution Melting, resulted significantly associated (P < 0.01) with average daily gain, feed conversion ratio, lean cuts and hams weight estimated breeding values. PLIN2 gene expression analysis in skeletal muscle of Italian Large White and Italian Duroc pigs divergent for backfat thickness and visible intermuscular fat showed a trend of higher expression level in pigs with higher intermuscular fat. These results suggest that PLIN2 can be a marker for carcass quality in pigs. Further investigation at both gene and protein level could elucidate its role on fat deposition.     

Widespread Genomic Signatures of Natural Selection in Hominid Evolution

Selection acting on genomic functional elements can be detected by its indirect effects on population diversity at linked neutral sites. To illuminate the selective forces that shaped hominid evolution, we analyzed the genomic distributions of human polymorphisms and sequence differences among five primate species relative to the locations of conserved sequence features. Neutral sequence diversity in human and ancestral hominid populations is substantially reduced near such features, resulting in a surprisingly large genome average diversity reduction due to selection of 19–26% on the autosomes and 12–40% on the X chromosome. The overall trends are broadly consistent with “background selection” or hitchhiking in ancestral populations acting to remove deleterious variants. Average selection is much stronger on exonic (both protein-coding and untranslated) conserved features than non-exonic features. Long term selection, rather than complex speciation scenarios, explains the large intragenomic variation in human/chimpanzee divergence. Our analyses reveal a dominant role for selection in shaping genomic diversity and divergence patterns, clarify hominid evolution, and provide a baseline for investigating specific selective events.     

The porcine <Emphasis Type="Italic">TBC1D1</Emphasis> gene: mapping,SNP identification,and association study with meat,carcass and production traits in Italian heavy pigs

TBC1D1 [TBC1 (tre-2/USP6, BUB2, cdc16) domain family, member 1] is a Rab-GTPase-activating related protein implicated in regulating the trafficking of glucose transporter 4 (GLUT4 or SLC2A4) storage vesicles to the cell surface in response to insulin and AMPK-activating stimuli in skeletal muscle. Mutations in the human and mouse TBC1D1 genes confer risk of obesity or leanness. We identified five single nucleotide polymorphisms (SNPs) in the porcine TBC1D1 gene. One of them (FN677935:g.219G>A) was genotyped either by high resolution melting and PCR-RFLP analyses to study allele frequencies in a few pig breeds and evaluate association with meat production and carcass traits in five groups of sib-tested pigs of Italian Large White and Italian Duroc breeds. The g.219G>A SNP was associated (P < 0.05) with ham weight, back fat thickness and lean cuts content in Italian Large White and with visible intermuscular fat in Italian Duroc pigs.     

猪ACTA2基因的克隆、表达分析及其与生产性状的关联

为鉴定对猪生产性状有重要影响的新分子标记,文章采用电子克隆结合PCR方法获得了猪肌动蛋白α2(Actin alpha 2, ACTA2)基因编码区序列和部分基因组序列,建立了第2内含子C1554T替换的PCR-HinfⅠ- RFLP基因分型方法,在所检测的7个不同猪群中除大白和梅大群体外,其他群体中均是C等位基因频率高于T等位基因的频率。标记与性状关联分析发现ACTA2基因型与肩部背膘厚、臀部背膘厚、肥肉率、瘦肉率、股二头肌pH和肌内脂肪显著或极显著相关,TT基因型与CC基因型相比具有更高的瘦肉率以及更低的肥肉率和背膘厚。通过Real-time RT-PCR分析发现ACTA2在大白和梅山两个品种猪骨骼肌中的表达量都随着日龄的增加而降低,在各个阶段,梅山猪中的表达量都比大白猪中的表达量高。     

Exonic Splicing Mutations Are More Prevalent than Currently Estimated and Can Be Predicted by Using In Silico Tools

The identification of a causal mutation is essential for molecular diagnosis and clinical management of many genetic disorders. However, even if next-generation exome sequencing has greatly improved the detection of nucleotide changes, the biological interpretation of most exonic variants remains challenging. Moreover, particular attention is typically given to protein-coding changes often neglecting the potential impact of exonic variants on RNA splicing. Here, we used the exon 10 of MLH1, a gene implicated in hereditary cancer, as a model system to assess the prevalence of RNA splicing mutations among all single-nucleotide variants identified in a given exon. We performed comprehensive minigene assays and analyzed patient’s RNA when available. Our study revealed a staggering number of splicing mutations in MLH1 exon 10 (77% of the 22 analyzed variants), including mutations directly affecting splice sites and, particularly, mutations altering potential splicing regulatory elements (ESRs). We then used this thoroughly characterized dataset, together with experimental data derived from previous studies on BRCA1, BRCA2, CFTR and NF1, to evaluate the predictive power of 3 in silico approaches recently described as promising tools for pinpointing ESR-mutations. Our results indicate that ΔtESRseq and ΔHZ_EI-based approaches not only discriminate which variants affect splicing, but also predict the direction and severity of the induced splicing defects. In contrast, the ΔΨ-based approach did not show a compelling predictive power. Our data indicates that exonic splicing mutations are more prevalent than currently appreciated and that they can now be predicted by using bioinformatics methods. These findings have implications for all genetically-caused diseases.     

Artificial selection on introduced Asian haplotypes shaped the genetic architecture in European commercial pigs

Early pig farmers in Europe imported Asian pigs to cross with their local breeds in order to improve traits of commercial interest. Current genomics techniques enabled genome-wide identification of these Asian introgressed haplotypes in modern European pig breeds. We propose that the Asian variants are still present because they affect phenotypes that were important for ancient traditional, as well as recent, commercial pig breeding. Genome-wide introgression levels were only weakly correlated with gene content and recombination frequency. However, regions with an excess or absence of Asian haplotypes (AS) contained genes that were previously identified as phenotypically important such as FASN, ME1, and KIT. Therefore, the Asian alleles are thought to have an effect on phenotypes that were historically under selection. We aimed to estimate the effect of AS in introgressed regions in Large White pigs on the traits of backfat (BF) and litter size. The majority of regions we tested that retained Asian deoxyribonucleic acid (DNA) showed significantly increased BF from the Asian alleles. Our results suggest that the introgression in Large White pigs has been strongly determined by the selective pressure acting upon the introgressed AS. We therefore conclude that human-driven hybridization and selection contributed to the genomic architecture of these commercial pigs.     

Somatic mtDNA Mutation Spectra in the Aging Human Putamen

The accumulation of heteroplasmic mitochondrial DNA (mtDNA) deletions and single nucleotide variants (SNVs) is a well-accepted facet of the biology of aging, yet comprehensive mutation spectra have not been described. To address this, we have used next generation sequencing of mtDNA-enriched libraries (Mito-Seq) to investigate mtDNA mutation spectra of putamen from young and aged donors. Frequencies of the “common” deletion and other “major arc” deletions were significantly increased in the aged cohort with the fold increase in the frequency of the common deletion exceeding that of major arc deletions. SNVs also increased with age with the highest rate of accumulation in the non-coding control region which contains elements necessary for translation and replication. Examination of predicted amino acid changes revealed a skew towards pathogenic SNVs in the coding region driven by mutation bias. Levels of the pathogenic m.3243A>G tRNA mutation were also found to increase with age. Novel multimeric tandem duplications that resemble murine control region multimers and yeast ρ⁻ mtDNAs, were identified in both young and aged specimens. Clonal ∼50 bp deletions in the control region were found at high frequencies in aged specimens. Our results reveal the complex manner in which the mitochondrial genome alters with age and provides a foundation for studies of other tissues and disease states.