Imputation of sequence level genotypes in the Franches-Montagnes horse breed

BACKGROUND A cost-effective strategy to increase the density of available markers within a population is to sequence a small proportion of the population and impute whole-genome sequence data for the remaining population. Increased densities of typed markers are advantageous for genome-wide association studies (GWAS) and genomic predictions. METHODS We obtained genotypes for 54 602 SNPs (single nucleotide polymorphisms) in 1077 Franches-Montagnes (FM) horses and Illumina paired-end whole-genome sequencing data for 30 FM horses and 14 Warmblood horses. After variant calling, the sequence-derived SNP genotypes (~13 million SNPs) were used for genotype imputation with the software programs Beagle, Impute2 and FImpute. RESULTS The mean imputation accuracy of FM horses using Impute2 was 92.0%. Imputation accuracy using Beagle and FImpute was 74.3% and 77.2%, respectively. In addition, for Impute2 we determined the imputation accuracy of all individual horses in the validation population, which ranged from 85.7% to 99.8%. The subsequent inclusion of Warmblood sequence data further increased the correlation between true and imputed genotypes for most horses, especially for horses with a high level of admixture. The final imputation accuracy of the horses ranged from 91.2% to 99.5%. CONCLUSIONS Using Impute2, the imputation accuracy was higher than 91% for all horses in the validation population, which indicates that direct imputation of 50k SNP-chip data to sequence level genotypes is feasible in the FM population. The individual imputation accuracy depended mainly on the applied software and the level of admixture.

Hairless Streaks in Cattle Implicate TSR2 in Early Hair Follicle Formation

Four related cows showed hairless streaks on various parts of the body with no correlation to the pigmentation pattern. The stripes occurred in a consistent pattern resembling the lines of Blaschko. The non-syndromic hairlessness phenotype observed occurred across three generations of a single family and was compatible with an X-linked mode of inheritance. Linkage analysis and subsequent whole genome sequencing of one affected female identified two perfectly associated non-synonymous sequence variants in the critical interval on bovine chromosome X. Both variants occurred in complete linkage disequilibrium and were absent in more than 3900 controls. An ERCC6L missense mutation was predicted to cause an amino acid substitution of a non-conserved residue. Analysis in mice showed no specific Ercc6l expression pattern related to hair follicle development and therefore ERCC6L was not considered as causative gene. A point mutation at the 5'-splice junction of exon 5 of the TSR2, 20S rRNA accumulation, homolog (S. cerevisiae), gene led to the production of two mutant transcripts, both of which contain a frameshift and generate a premature stop codon predicted to truncate approximately 25% of the protein. Interestingly, in addition to the presence of both physiological TSR2 transcripts, the two mutant transcripts were predominantly detected in the hairless skin of the affected cows. Immunohistochemistry, using an antibody against the N-terminal part of the bovine protein demonstrated the specific expression of the TSR2 protein in the skin and the hair of the affected and the control cows as well as in bovine fetal skin and hair. The RNA hybridization in situ showed that Tsr2 was expressed in pre- and post-natal phases of hair follicle development in mice. Mammalian TSR2 proteins are highly conserved and are known to be broadly expressed, but their precise in vivo functions are poorly understood. Thus, by dissecting a naturally occurring mutation in a domestic animal species, we identified TSR2 as a regulator of hair follicle development.

Reclassification of Actinobacillus muris as Muribacter muris gen. nov. comb. nov.

To reinvestigate the taxonomy of [Actinobacillus] muris, 474 strains mainly from mice and rats were characterized by phenotype and 130 strains selected for genotypic characterization by 16S rRNA and partial rpoB gene sequencing. The type strain was further investigated by whole genome sequencing. Phylogenetic analysis of the DNA sequences showed one monophyletic group with intra group similarities of 96.7 % and 97.2 % for 16S rRNA and rpoB genes, respectively. The lowest 16S rRNA similarity to the closest related valid named taxon outside the group was 95.9 % to the type strain of [Pasteurella] pneumotropica. The closest related taxon based on rpoB sequence comparison was 'Haemophilus influenzae-murium' with 88.4 %. A new genus, Muribacter is proposed based on a distinct phylogenetic position based on 16S rRNA and rpoB gene sequence comparisons with major divergence to the existing genera of Pasteurellaceae. The new genus includes the characteristics of [Actinobacillus] muris with the emendation that acid formation from (-)-D-mannitol is variable as well the hydrolysis of esculin while the α-glucosidase test is positive. There is no requirement for exogenously supplied nicotinamide adenine dinucleotide (V factor) for the majority of strains investigated, however, one strain was found positive. The major fatty acids of the type strain of Muribacter muris were C 14:0, C 14:0 3OH/C 16:1 ISOI, C 16:1 ω7c and C 16:0 which is in line with most genera of Pasteurellaceae. The type strain of Muribacter muris is CCUG 16938T ( = NCTC 12432T = ATCC 49577T).

Bone-Conditioned Medium Changes Gene Expression in Bone-Derived Fibroblasts.

PURPOSE Autologous bone is used for augmentation in the course of oral implant placement. Bone grafts release paracrine signals that can modulate mesenchymal cell differentiation in vitro. The detailed genetic response of the bone-derived fibroblasts to these paracrine signals has remained elusive. Paracrine signals accumulate in bone-conditioned medium (BCM) prepared from porcine cortical bone chips. MATERIALS AND METHODS In this study, bone-derived fibroblasts were exposed to BCM followed by a whole genome expression profiling and downstream quantitative reverse transciptase polymerase chain reaction of the most strongly regulated genes. RESULTS The data show that ADM, IL11, IL33, NOX4, PRG4, and PTX3 were differentially expressed in response to BCM in bone-derived fibroblasts. The transforming growth factor beta (TGF-β) receptor 1 antagonist SB431542 blocked the effect of BCM on the expression of the gene panel, except for IL33. CONCLUSION These in vitro results extend existing evidence that cortical bone chips release paracrine signals that provoke a robust genetic response in mesenchymal cells that is not exclusively mediated via the TGF-β receptor. The present data provide further insights into the process of graft consolidation.

TECPR2 Associated Neuroaxonal Dystrophy in Spanish Water Dogs.

Clinical, pathological and genetic examination revealed an as yet uncharacterized juvenile-onset neuroaxonal dystrophy (NAD) in Spanish water dogs. Affected dogs presented with various neurological deficits including gait abnormalities and behavioral deficits. Histopathology demonstrated spheroid formation accentuated in the grey matter of the cerebral hemispheres, the cerebellum, the brain stem and in the sensory pathways of the spinal cord. Iron accumulation was absent. Ultrastructurally spheroids contained predominantly closely packed vesicles with a double-layered membrane, which were characterized as autophagosomes using immunohistochemistry. The family history of the four affected dogs suggested an autosomal recessive inheritance. SNP genotyping showed a single genomic region of extended homozygosity of 4.5 Mb in the four cases on CFA 8. Linkage analysis revealed a maximal parametric LOD score of 2.5 at this region. By whole genome re-sequencing of one affected dog, a perfectly associated, single, non-synonymous coding variant in the canine tectonin beta-propeller repeat-containing protein 2 (TECPR2) gene affecting a highly conserved region was detected (c.4009C>T or p.R1337W). This canine NAD form displays etiologic parallels to an inherited TECPR2 associated type of human hereditary spastic paraparesis (HSP). In contrast to the canine NAD, the spinal cord lesions in most types of human HSP involve the sensory and the motor pathways. Furthermore, the canine NAD form reveals similarities to cases of human NAD defined by widespread spheroid formation without iron accumulation in the basal ganglia. Thus TECPR2 should also be considered as candidate gene for human NAD. Immunohistochemistry and the ultrastructural findings further support the assumption, that TECPR2 regulates autophagosome accumulation in the autophagic pathways. Consequently, this report provides the first genetic characterization of juvenile canine NAD, describes the histopathological features associated with the TECPR2 mutation and provides evidence to emphasize the association between failure of autophagy and neurodegeneration.

Spatial genomic heterogeneity within localized, multifocal prostate cancer

Herein we provide a detailed molecular analysis of the spatial heterogeneity of clinically localized, multifocal prostate cancer to delineate new oncogenes or tumor suppressors. We initially determined the copy number aberration (CNA) profiles of 74 patients with index tumors of Gleason score 7. Of these, 5 patients were subjected to whole-genome sequencing using DNA quantities achievable in diagnostic biopsies, with detailed spatial sampling of 23 distinct tumor regions to assess intraprostatic heterogeneity in focal genomics. Multifocal tumors are highly heterogeneous for single-nucleotide variants (SNVs), CNAs and genomic rearrangements. We identified and validated a new recurrent amplification of MYCL, which is associated with TP53 deletion and unique profiles of DNA damage and transcriptional dysregulation. Moreover, we demonstrate divergent tumor evolution in multifocal cancer and, in some cases, tumors of independent clonal origin. These data represent the first systematic relation of intraprostatic genomic heterogeneity to predicted clinical outcome and inform the development of novel biomarkers that reflect individual prognosis.

A transposable element insertion in APOB causes cholesterol deficiency in Holstein cattle

Cholesterol deficiency, a new autosomal recessive inherited genetic defect in Holstein cattle, has been recently reported to have an influence on the rearing success of calves. The affected animals show unresponsive diarrhea accompanied by hypocholesterolemia and usually die within the first weeks or months of life. Here, we show that whole genome sequencing combined with the knowledge about the pedigree and inbreeding status of a livestock population facilitates the identification of the causative mutation. We resequenced the entire genomes of an affected calf and a healthy partially inbred male carrying one copy of the critical 2.24-Mb chromosome 11 segment in its ancestral state and one copy of the same segment with the cholesterol deficiency mutation. We detected a single structural variant, homozygous in the affected case and heterozygous in the non-affected carrier male. The genetic makeup of this key animal provides extremely strong support for the causality of this mutation. The mutation represents a 1.3kb insertion of a transposable LTR element (ERV2-1) in the coding sequence of the APOB gene, which leads to truncated transcripts and aberrant splicing. This finding was further supported by RNA sequencing of the liver transcriptome of an affected calf. The encoded apolipoprotein B is an essential apolipoprotein on chylomicrons and low-density lipoproteins, and therefore, the mutation represents a loss of function mutation similar to autosomal recessive inherited familial hypobetalipoproteinemia-1 (FHBL1) in humans. Our findings provide a direct gene test to improve selection against this deleterious mutation in Holstein cattle.

Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage.

Mycobacterium tuberculosis strains of the Beijing lineage are globally distributed and are associated with the massive spread of multidrug-resistant (MDR) tuberculosis in Eurasia. Here we reconstructed the biogeographical structure and evolutionary history of this lineage by genetic analysis of 4,987 isolates from 99 countries and whole-genome sequencing of 110 representative isolates. We show that this lineage initially originated in the Far East, from where it radiated worldwide in several waves. We detected successive increases in population size for this pathogen over the last 200 years, practically coinciding with the Industrial Revolution, the First World War and HIV epidemics. Two MDR clones of this lineage started to spread throughout central Asia and Russia concomitantly with the collapse of the public health system in the former Soviet Union. Mutations identified in genes putatively under positive selection and associated with virulence might have favored the expansion of the most successful branches of the lineage.

Taxonomy, phylogeny and molecular epidemiology of Echinococcus multilocularis: From fundamental knowledge to health ecology.

Alveolar echinococcosis, caused by the tapeworm Echinococcus multilocularis, is one of the most severe parasitic diseases in humans and represents one of the 17 neglected diseases prioritised by the World Health Organisation (WHO) in 2012. Considering the major medical and veterinary importance of this parasite, the phylogeny of the genus Echinococcus is of considerable importance; yet, despite numerous efforts with both mitochondrial and nuclear data, it has remained unresolved. The genus is clearly complex, and this is one of the reasons for the incomplete understanding of its taxonomy. Although taxonomic studies have recognised E. multilocularis as a separate entity from the Echinococcus granulosus complex and other members of the genus, it would be premature to draw firm conclusions about the taxonomy of the genus before the phylogeny of the whole genus is fully resolved. The recent sequencing of E. multilocularis and E. granulosus genomes opens new possibilities for performing in-depth phylogenetic analyses. In addition, whole genome data provide the possibility of inferring phylogenies based on a large number of functional genes, i.e. genes that trace the evolutionary history of adaptation in E. multilocularis and other members of the genus. Moreover, genomic data open new avenues for studying the molecular epidemiology of E. multilocularis: genotyping studies with larger panels of genetic markers allow the genetic diversity and spatial dynamics of parasites to be evaluated with greater precision. There is an urgent need for international coordination of genotyping of E. multilocularis isolates from animals and human patients. This could be fundamental for a better understanding of the transmission of alveolar echinococcosis and for designing efficient healthcare strategies.

Standard genotyping overestimates transmission of Mycobacterium tuberculosis among immigrants in a low incidence country.

Immigrants from high tuberculosis (TB) incidence regions are a risk group for TB in low-incidence countries such as Switzerland. In a previous analysis of a nationwide collection of 520 Mycobacterium tuberculosis isolates from 2000-2008, we identified 35 clusters comprising 90 patients based on standard genotyping (24-loci MIRU-VNTR and spoligotyping). Here, we used whole genome sequencing (WGS) to revisit these transmission clusters. Genome-based transmission clusters were defined as isolate pairs separated by ≤12 single nucleotide polymorphisms (SNPs). WGS confirmed 17/35 (49%) MIRU-VNTR clusters; the other 18 clusters contained pairs separated by >12 SNPs. Most transmission clusters (3/4) of Swiss-born patients were confirmed by WGS, as opposed to 25% (4/16) of clusters involving only foreign-born patients. The overall clustering proportion using standard genotyping was 17% (90 patients, 95% confidence interval [CI]: 14-21%), but only 8% (43 patients, 95% CI: 6-11%) using WGS. The clustering proportion was 17% (67/401, 95% CI: 13-21%) using standard genotyping and 7% (26/401, 95% CI: 4-9%) using WGS among foreign-born patients, and 19% (23/119, 95% CI: 13-28%) and 14% (17/119, 95% CI: 9-22%), respectively, among Swiss-born patients. Using weighted logistic regression, we found weak evidence for an association between birth origin and transmission (aOR 2.2, 95% CI: 0.9-5.5, comparing Swiss-born patients to others). In conclusion, standard genotyping overestimated recent TB transmission in Switzerland when compared to WGS, particularly among immigrants from high TB incidence regions, where genetically closely related strains often predominate. We recommend the use of WGS to identify transmission clusters in low TB incidence settings.

Initial characterization of stiff skin-like syndrome in West Highland white terriers

BACKGROUND Stiff skin syndrome and systemic or localized scleroderma are cutaneous disorders characterized by dermal fibrosis and present clinically with induration of the skin, with or without joint, internal organ or vascular involvement. OBJECTIVES To provide clinical, histological and preliminary genetic analysis of two West Highland white terrier siblings presenting with indurated skin resembling stiff skin syndrome in humans. ANIMALS Two client owned full sibling West Highland white terriers from two different litters. METHODS Clinical examination, histopathological examination and whole genome sequencing analysis of affected and unaffected West Highland white terriers. RESULTS Affected dogs exhibited markedly indurated skin that was attached firmly to the underlying tissue and incomplete closure of the mouth and eyes. No abnormalities were found by neurological or orthopaedic examination, radiographs of the head or whole body computed tomography. Histologically, the dermis and pannicular septa were thickened by a marked increase in coarse collagen fibres and a mild to moderate increase in collagen fibre diameter. The syndrome most likely follows an autosomal recessive mode of inheritance. The sequence analysis did not reveal any obvious causative variant in the investigated candidate genes ADAMTSL2 and FBN1. CONCLUSION AND CLINICAL IMPORTANCE The clinical phenotype and histopathological features of two West Highland white terrier siblings resembled stiff skin syndrome in humans. Unlike in humans, or previously described beagles with stiff skin, there was no restriction of joint mobility. Genetic analysis did not detect a candidate causative variant and warrants further research.

Insight into the evolution of the Solanaceae from the parental genomes of Petunia hybrida

Petunia hybrida is a popular bedding plant that has a long history as a genetic model system. We report the whole-genome sequencing and assembly of inbred derivatives of its two wild parents, P. axillaris N and P. inflata S6. The assemblies include 91.3% and 90.2% coverage of their diploid genomes (1.4 Gb; 2n = 14) containing 32,928 and 36,697 protein-coding genes, respectively. The genomes reveal that the Petunia lineage has experienced at least two rounds of hexaploidization: the older gamma event, which is shared with most Eudicots, and a more recent Solanaceae event that is shared with tomato and other solanaceous species. Transcription factors involved in the shift from bee to moth pollination reside in particularly dynamic regions of the genome, which may have been key to the remarkable diversity of floral colour patterns and pollination systems. The high-quality genome sequences will enhance the value of Petunia as a model system for research on unique biological phenomena such as small RNAs, symbiosis, self-incompatibility and circadian rhythms.

Lethal chondrodysplasia in a family of Holstein cattle is associated with a de novo splice site variant of COL2A1

Background Lethal chondrodysplasia (bulldog syndrome) is a well-known congenital syndrome in cattle and occurs sporadically in many breeds. In 2015, it was noticed that about 12 % of the offspring of the phenotypically normal Danish Holstein sire VH Cadiz Captivo showed chondrodysplasia resembling previously reported bulldog calves. Pedigree analysis of affected calves did not display obvious inbreeding to a common ancestor, suggesting the causative allele was not a rare recessive. The normal phenotype of the sire suggested a dominant inheritance with incomplete penetrance or a mosaic mutation. Results Three malformed calves were examined by necropsy, histopathology, radiology, and computed tomography scanning. These calves were morphologically similar and displayed severe disproportionate dwarfism and reduced body weight. The syndrome was characterized by shortening and compression of the body due to reduced length of the spine and the long bones of the limbs. The vicerocranium had severe dysplasia and palatoschisis. The bones had small irregular diaphyses and enlarged epiphyses consisting only of chondroid tissue. The sire and a total of four affected half-sib offspring and their dams were genotyped with the BovineHD SNP array to map the defect in the genome. Significant genetic linkage was obtained for several regions of the bovine genome including chromosome 5 where whole genome sequencing of an affected calf revealed a COL2A1 point mutation (g.32473300 G > A). This private sequence variant was predicted to affect splicing as it altered the conserved splice donor sequence GT at the 5’-end of COL2A1 intron 36, which was changed to AT. All five available cases carried the mutant allele in heterozygous state and all five dams were homozygous wild type. The sire VH Cadiz Captivo was shown to be a gonadal and somatic mosaic as assessed by the presence of the mutant allele at levels of about 5 % in peripheral blood and 15 % in semen. Conclusions The phenotypic and genetic findings are comparable to a previously reported COL2A1 missense mutation underlying lethal chondrodysplasia in the offspring of a mosaic French Holstein sire (Igale Masc). The identified independent spontaneous splice site variant in COL2A1 most likely caused chondrodysplasia and must have occurred during the early foetal development of the sire. This study provides a first example of a dominant COL2A1 splice site variant as candidate causal mutation of a severe lethal chondrodysplasia phenotype. Germline mosaicism is a relatively frequent mechanism in the origin of genetic disorders and explains the prevalence of a certain fraction of affected offspring. Paternal dominant de novo mutations are a risk in cattle breeding, especially because the ratio of defective offspring may be very high and be associated with significant animal welfare problems.

Genomic amplification of the caprine EDNRA locus might lead to a dose dependent loss of pigmentation

The South African Boer goat displays a characteristic white spotting phenotype, in which the pigment is limited to the head. Exploiting the existing phenotype variation within the breed, we mapped the locus causing this white spotting phenotype to chromosome 17 by genome wide association. Subsequent whole genome sequencing identified a 1 Mb copy number variant (CNV) harboring 5 genes including EDNRA. The analysis of 358 Boer goats revealed 3 alleles with one, two, and three copies of this CNV. The copy number is correlated with the degree of white spotting in goats. We propose a hypothesis that ectopic overexpression of a mutant EDNRA scavenges EDN3 required for EDNRB signaling and normal melanocyte development and thus likely lead to an absence of melanocytes in the non-pigmented body areas of Boer goats. Our findings demonstrate the value of domestic animals as reservoir of unique mutants and for identifying a precisely defined functional CNV.

Salivary Pellicle vs Whole Saliva: The Response of Oral Fibroblasts Based on a Genome-Wide Microarray.

PURPOSE Whole saliva comprises components of the salivary pellicle that spontaneously forms on surfaces of implants and teeth. However, there are no studies that functionally link the salivary pellicle with a possible change in gene expression. MATERIALS AND METHODS This study examined the genetic response of oral fibroblasts exposed to the salivary pellicle and whole saliva. Oral fibroblasts were seeded onto a salivary pellicle and the respective untreated surface. Oral fibroblasts were also exposed to freshly harvested sterile-filtered whole saliva. A genome-wide microarray of oral fibroblasts was performed, followed by gene ontology screening with DAVID functional annotation clustering, KEGG pathway analysis, and the STRING functional protein association network. RESULTS Exposure of oral fibroblasts to saliva caused 61 genes to be differentially expressed (P < .05). Gene ontology screening assigned the respective genes into 262 biologic processes, 3 cellular components, 13 molecular functions, and 7 pathways. Most remarkable was the enrichment in the inflammatory response. None of the genes regulated by whole saliva was significantly changed when cells were placed onto a salivary pellicle. CONCLUSION The salivary pellicle per se does not provoke a significant inflammatory response of oral fibroblasts in vitro, whereas sterile-filtered whole saliva does produce a strong inflammatory response.