SERPINB11 frameshift variant associated with novel hoof specific phenotype in Connemara ponies.

Horses belong to the order Perissodactyla and bear the majority of their weight on their third toe; therefore, tremendous force is applied to each hoof. An inherited disease characterized by a phenotype restricted to the dorsal hoof wall was identified in the Connemara pony. Hoof wall separation disease (HWSD) manifests clinically as separation of the dorsal hoof wall along the weight-bearing surface of the hoof during the first year of life. Parents of affected ponies appeared clinically normal, suggesting an autosomal recessive mode of inheritance. A case-control allelic genome wide association analysis was performed (ncases = 15, ncontrols = 24). Population stratification (λ = 1.48) was successfully improved by removing outliers (ncontrols = 7) identified on a multidimensional scaling plot. A genome-wide significant association was detected on chromosome 8 (praw = 1.37x10-10, pgenome = 1.92x10-5). A homozygous region identified in affected ponies spanned from 79,936,024-81,676,900 bp and contained a family of 13 annotated SERPINB genes. Whole genome next-generation sequencing at 6x coverage of two cases and two controls revealed 9,758 SNVs and 1,230 indels within the ~1.7-Mb haplotype, of which 17 and 5, respectively, segregated with the disease and were located within or adjacent to genes. Additional genotyping of these 22 putative functional variants in 369 Connemara ponies (ncases = 23, ncontrols = 346) and 169 horses of other breeds revealed segregation of three putative variants adjacent or within four SERPIN genes. Two of the variants were non-coding and one was an insertion within SERPINB11 that introduced a frameshift resulting in a premature stop codon. Evaluation of mRNA levels at the proximal hoof capsule (ncases = 4, ncontrols = 4) revealed that SERPINB11 expression was significantly reduced in affected ponies (p<0.001). Carrier frequency was estimated at 14.8%. This study describes the first genetic variant associated with a hoof wall specific phenotype and suggests a role of SERPINB11 in maintaining hoof wall structure.

Whole genome sequencing of peach (Prunus persica L.) for SNP identification and selection.

BACKGROUND: The application of next generation sequencing technologies and bioinformatic scripts to identify high frequency SNPs distributed throughout the peach genome is described. Three peach genomes were sequenced using Roche 454 and Illumina/Solexa technologies to obtain long contigs for alignment to the draft 'Lovell' peach sequence as well as sufficient depth of coverage for 'in silico' SNP discovery. DESCRIPTION: The sequences were aligned to the 'Lovell' peach genome released April 01, 2010 by the International Peach Genome Initiative (IPGI). 'Dr. Davis', 'F8, 1-42' and 'Georgia Belle' were sequenced to add SNPs segregating in two breeding populations, Pop DF ('Dr. Davis' × 'F8, 1-42') and Pop DG ('Dr. Davis' × 'Georgia Belle'). Roche 454 sequencing produced 980,000 total reads with 236 Mb sequence for 'Dr. Davis' and 735,000 total reads with 172 Mb sequence for 'F8, 1-42'. 84 bp × 84 bp paired end Illumina/Solexa sequences yielded 25.5, 21.4, 25.5 million sequences for 'Dr. Davis', 'F8, 1-42' and 'Georgia Belle', respectively. BWA/SAMtools were used for alignment of raw reads and SNP detection, with custom PERL scripts for SNP filtering. Velvet's Columbus module was used for sequence assembly. Comparison of aligned and overlapping sequences from both Roche 454 and Illumina/Solexa resulted in the selection of 6654 high quality SNPs for 'Dr. Davis' vs. 'F8, 1-42' and 'Georgia Belle', distributed on eight major peach genome scaffolds as defined from the 'Lovell' assembly. CONCLUSION: The eight scaffolds contained about 215-225 Mb of peach genomic sequences with one SNP/~ 40,000 bases. All sequences from Roche 454 and Illumina/Solexa have been submitted to NCBI for public use in the Short Read Archive database. SNPs have been deposited in the NCBI SNP database.

Genome-Scale Analysis of Programmed DNA Elimination Sites in Tetrahymena thermophila.

Genetically programmed DNA rearrangements can regulate mRNA expression at an individual locus or, for some organisms, on a genome-wide scale. Ciliates rely on a remarkable process of whole-genome remodeling by DNA elimination to differentiate an expressed macronucleus (MAC) from a copy of the germline micronucleus (MIC) in each cycle of sexual reproduction. Here we describe results from the first high-throughput sequencing effort to investigate ciliate genome restructuring, comparing Sanger long-read sequences from a Tetrahymena thermophila MIC genome library to the MAC genome assembly. With almost 25% coverage of the unique-sequence MAC genome by MIC genome sequence reads, we created a resource for positional analysis of MIC-specific DNA removal that pinpoints MAC genome sites of DNA elimination at nucleotide resolution. The widespread distribution of internal eliminated sequences (IES) in promoter regions and introns suggests that MAC genome restructuring is essential not only for what it removes (for example, active transposons) but also for what it creates (for example, splicing-competent introns). Consistent with the heterogeneous boundaries and epigenetically modulated efficiency of individual IES deletions studied to date, we find that IES sites are dramatically under-represented in the ∼25% of the MAC genome encoding exons. As an exception to this general rule, we discovered a previously unknown class of small (<500 bp) IES with precise elimination boundaries that can contribute the 3' exon of an mRNA expressed during genome restructuring, providing a new mechanism for expanding mRNA complexity in a developmentally regulated manner.