Spontaneous gliomas in dogs occur at a frequency similar to that in humans and may provide a translational model for therapeutic development and comparative biological investigations. Copy number alterations in 38 canine gliomas, including diffuse astrocytomas, glioblastomas, oligodendrogliomas, and mixed oligoastrocytomas, were defined using an Illumina 170K single nucleotide polymorphism array. Highly recurrent alterations were seen in up to 85% of some tumor types, most notably involving chromosomes 13, 22, and 38, and gliomas clustered into 2 major groups consisting of high-grade IV astrocytomas, or oligodendrogliomas and other tumors. Tumor types were characterized by specific broad and focal chromosomal events including focal loss of the INK4A/B locus in glioblastoma and loss of the RB1 gene and amplification of the PDGFRA gene in oligodendrogliomas. Genes associated with the 3 critical pathways in human high-grade gliomas (TP53, RB1, and RTK/RAS/PI3K) were frequently associated with canine aberrations. Analysis of oligodendrogliomas revealed regions of chromosomal losses syntenic to human 1p involving tumor suppressor genes, such as CDKN2C, as well as genes associated with apoptosis, autophagy, and response to chemotherapy and radiation. Analysis of high frequency chromosomal aberrations with respect to human orthologues may provide insight into both novel and common pathways in gliomagenesis and response to therapy.

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.

Abstract 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.

Dilated cardiomyopathy (DCM) is characterized by decreased systolic function and dilation of one or both ventricles, often leading to heart failure or sudden death. Two 10-month-old sibling Nova Scotia Duck Tolling Retrievers (NSDTR) died acutely with evidence of dilated cardiomyopathy with myocardial fibrosis. Association analysis using two cases and 35 controls identified three candidate regions homozygous in the two cases. Whole genome sequencing identified a frameshift deletion in the LMNA gene (NC_049228.1:g.41688530del, NP_001274080:p.(Asp576ThrfsTer124)). Three retrospectively identified NSDTRs with sudden death before 2 years of age and severe myocardial fibrosis were also homozygous for the deletion. One 5 year old with sudden death and myocardial fibrosis was heterozygous for the deletion. This variant was not identified in 722 dogs of other breeds, nor was it identified to be homozygous in 784 NSDTR. LMNA codes for lamin A/C proteins, which are type V intermediate filaments that provide structural support to the nuclear membrane. In humans, LMNA variants can cause DCM with sudden death as well as diseases of striated muscles, lipodystrophy, neuropathies, and accelerated aging disorders. This frameshift deletion is predicted to affect processing of prelamin A into lamin A. Pedigree analysis in the NSDTR and functional evaluation of heterozygotes is consistent with a predominantly recessive mode of inheritance and possibly low penetrance in heterozygotes in contrast to people, where most pathogenic LMNA variants are dominantly inherited.

Female parenting is obligate in mammals, but fathering behavior among mammals is rare. Only 3-5% of mammalian species exhibit biparental care, including humans, and mechanisms of fathering behavior remain sparsely studied. However, in species where it does exist, paternal care is often crucial to the survivorship of offspring. The present study is the first to identify new gene targets linked to the experience of fathering behavior in a biparental species using RNA sequencing. In order to determine the pattern of gene expression within the medial preoptic area that is specifically associated with fathering behavior, we identified genes in male prairie voles (Microtus ochrogaster) that experienced one of three social conditions: virgin males, pair bonded males, and males with fathering experience. A list of genes exhibiting different expression patterns in each comparison (i.e. Virgin vs Paired, Virgin vs Fathers, and Paired vs Fathers) was evaluated using the gene ontology enrichment analysis, and Kyoto Encyclopedia of Genes and Genomes pathways analysis to reveal metabolic pathways associated with specific genes. Using these tools, we generated a filtered list of genes that exhibited altered patterns of expression in voles with different amounts of social experience. Finally, we used NanoString to quantify differences in the expression of these selected genes. These genes are involved in a variety of processes, with enrichment in genes associated with immune function, metabolism, synaptic plasticity, and the remodeling of dendritic spines. The identification of these genes and processes will lead to novel insights into the biological basis of fathering behavior.

Background & aims

The pathogenesis of Wilson disease (WD) involves hepatic and brain copper accumulation resulting from pathogenic variants affecting the ATP7B gene and downstream epigenetic and metabolic mechanisms. Prior methylome investigations in human WD liver and blood and in the Jackson Laboratory (Bar Harbor, ME) C3He-Atp7b^tx-j/J (tx-j) WD mouse model revealed an epigenetic signature of WD, including changes in histone deacetylase (HDAC) 5. We tested the hypothesis that histone acetylation is altered with respect to copper overload and aberrant DNA methylation in WD.

Methods

We investigated class IIa HDAC4 and HDAC5 and H3K9/H3K27 histone acetylation in tx-j mouse livers compared with C3HeB/FeJ (C3H) control in response to 3 treatments: 60% kcal fat diet, D-penicillamine (copper chelator), and choline (methyl group donor). Experiments with copper-loaded hepatoma G2 cells were conducted to validate in vivo studies.

Results

In 9-week tx-j mice, HDAC5 levels increased significantly after 8 days of a 60% kcal fat diet compared with chow. In 24-week tx-j mice, HDAC4/5 levels were reduced 5- to 10-fold compared with C3H, likely through mechanisms involving HDAC phosphorylation. HDAC4/5 levels were affected by disease progression and accompanied by increased acetylation. D-penicillamine and choline partially restored HDAC4/5 and H3K9ac/H3K27ac to C3H levels. Integrated RNA and chromatin immunoprecipitation sequencing analyses revealed genes regulating energy metabolism and cellular stress/development, which, in turn, were regulated by histone acetylation in tx-j mice compared with C3H mice, with Pparα and Pparγ among the most relevant targets.

Conclusions

These results suggest dietary modulation of class IIa HDAC4/5, and subsequent H3K9/H3K27 acetylation/deacetylation can regulate gene expression in key metabolic pathways in the pathogenesis of WD.

Abstract Female parenting is obligate in mammals, but fathering behavior among mammals is rare. Only 3–5% of mammalian species exhibit biparental care, including humans, and mechanisms of fathering behavior remain sparsely studied. However, in species where it does exist, paternal care is often crucial to the survivorship of offspring. The present study is the first to identify new gene targets linked to the experience of fathering behavior in a biparental species using RNA sequencing. In order to determine the pattern of gene expression within the medial preoptic area that is specifically associated with fathering behavior, we identified differentially expressed genes in male prairie voles ( Microtus ochrogaster) that experienced one of three social conditions: virgin males, pair bonded males, and males with fathering experience. Differentially expressed genes from each comparison (i.e., Virgin vs Paired, Virgin vs Fathers, and Paired vs Fathers) were evaluated using the Gene Ontology enrichment analysis, and Kegg pathways analysis to reveal metabolic pathways associated with specific differentially expressed genes. Using these tools, we identified a group of genes that are differentially expressed in voles with different amounts of social experience. These genes are involved in a variety of processes, with particular enrichment in genes associated with immune function, metabolism, synaptic plasticity, and the remodeling of dendritic spines. The identification of these genes and processes will lead to novel insights into the biological basis of fathering behavior.

Squamous cell carcinoma (SCC) is the most common periocular cancer in horses and the second most common tumor of the horse overall. A missense mutation in damage-specific DNA-binding protein 2 (DDB2, c.1012 C>T, p.Thr338Met) was previously found to be strongly associated with ocular SCC in Haflinger and Belgian horses, explaining 76% of cases across both breeds. To determine if this same variant in DDB2 contributes to risk for ocular SCC in the Arabian, Appaloosa, and Percheron breeds and to determine if the variant contributes to risk for oral or urogenital SCC, histologically confirmed SCC cases were genotyped for the DDB2 variant and associations were investigated. Horses with urogenital SCC that were heterozygous for the DDB2 risk allele were identified in the Appaloosa breed, but a significant association between the DDB2 variant and SCC occurring at any location in this breed was not detected. The risk allele was not identified in Arabians, and no Percherons were homozygous for the risk allele. High-throughput sequencing data from six Haflingers were analyzed to ascertain if any other variant from the previously associated 483 kb locus on ECA12 was more concordant with the SCC phenotype than the DDB2 variant. Sixty polymorphisms were prioritized for evaluation, and no other variant from this locus explained the genetic risk better than the DDB2 allele (P = 3.39 × 10-17, n = 118). These data provide further support of the DDB2 variant contributing to risk for ocular SCC, specifically in the Haflinger and Belgian breeds.

Infiltration of regulatory T (Treg) cells into many tumor types correlates with poor patient prognoses. However, mechanisms of intratumoral Treg cell function remain to be elucidated. We investigated Treg cell function in a genetically engineered mouse model of lung adenocarcinoma and found that Treg cells suppressed anti-tumor responses in tumor-associated tertiary lymphoid structures (TA-TLSs). TA-TLSs have been described in human lung cancers, but their function remains to be determined. TLSs in this model were spatially associated with >90% of tumors and facilitated interactions between T cells and tumor-antigen-presenting dendritic cells (DCs). Costimulatory ligand expression by DCs and T cell proliferation rates increased in TA-TLSs upon Treg cell depletion, leading to tumor destruction. Thus, we propose that Treg cells in TA-TLSs can inhibit endogenous immune responses against tumors, and targeting these cells might provide therapeutic benefit for cancer patients.

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.