DNA-based diagnosis of rare diseases in veterinary medicine: a 4.4 kb deletion of ITGB4 is associated with epidermolysis bullosa in Charolais cattle.

BACKGROUND Rare diseases in livestock animals are traditionally poorly diagnosed. Other than clinical description and pathological examination, the underlying causes have, for the most part, remained unknown. A single case of congenital skin fragility in cattle was observed, necropsy, histological and ultrastructural examinations were carried out and whole genome sequencing was utilized to identify the causative mutation. RESULTS A single purebred female Charolais calf with severe skin lesions was delivered full-term and died spontaneously after birth. The clinical and pathological findings exactly matched the gross description given by previous reports on epitheliogenesis imperfecta and epidermolysis bullosa (EB) in cattle. Histological and ultrastructural changes were consistent with EB junctionalis (EBJ). Genetic analysis revealed a previously unpublished ITGB4 loss-of-function mutation; the affected calf was homozygous for a 4.4 kb deletion involving exons 17 to 22, and the dam carried a single copy of the deletion indicating recessive inheritance. The homozygous mutant genotype did not occur in healthy controls of various breeds but some heterozygous carriers were found among Charolais cattle belonging to the affected herd. The mutant allele was absent in a representative sample of unrelated sires of the German Charolais population. CONCLUSION This is the first time in which a recessively inherited ITGB4 associated EBJ has been reported in cattle. The identification of heterozygous carriers is of importance in avoiding the transmission of this defect in future. Current DNA sequencing methods offer a powerful tool for understanding the genetic background of rare diseases in domestic animals having a reference genome sequence available.

Mutations upstream of fabI in triclosan resistant Staphylococcus aureus strains are associated with elevated fabI gene expression

Background The enoyl-acyl carrier protein (ACP) reductase enzyme (FabI) is the target for a series of antimicrobial agents including novel compounds in clinical trial and the biocide triclosan. Mutations in fabI and heterodiploidy for fabI have been shown to confer resistance in S. aureus strains in a previous study. Here we further determined the fabI upstream sequence of a selection of these strains and the gene expression levels in strains with promoter region mutations. Results Mutations in the fabI promoter were found in 18% of triclosan resistant clinical isolates, regardless the previously identified molecular mechanism conferring resistance. Although not significant, a higher rate of promoter mutations were found in strains without previously described mechanisms of resistance. Some of the mutations identified in the clinical isolates were also detected in a series of laboratory mutants. Microarray analysis of selected laboratory mutants with fabI promoter region mutations, grown in the absence of triclosan, revealed increased fabI expression in three out of four tested strains. In two of these strains, only few genes other than fabI were upregulated. Consistently with these data, whole genome sequencing of in vitro selected mutants identified only few mutations except the upstream and coding regions of fabI, with the promoter mutation as the most probable cause of fabI overexpression. Importantly the gene expression profiling of clinical isolates containing similar mutations in the fabI promoter also showed, when compared to unrelated non-mutated isolates, a significant up-regulation of fabI. Conclusions In conclusion, we have demonstrated the presence of C34T, T109G, and A101C mutations in the fabI promoter region of strains with fabI up-regulation, both in clinical isolates and/or laboratory mutants. These data provide further observations linking mutations upstream fabI with up-regulated expression of the fabI gene.

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

A RAB3GAP1 SINE Insertion in Alaskan Huskies with Polyneuropathy, Ocular Abnormalities and Neuronal Vacuolation (POANV) Resembling Human Warburg Micro Syndrome 1 (WARBM1)

We observed a hereditary phenotype in Alaskan Huskies, which was characterized by polyneuropathy with ocular abnormalities and neuronal vacuolation (POANV). The affected dogs developed a progressive severe ataxia, which led to euthanasia between 8 and 16 months of age. The pedigrees were consistent with a monogenic autosomal recessive inheritance. We localized the causative genetic defect to a 4 Mb interval on chromosome 19 by a combined linkage and homozygosity mapping approach. Whole genome sequencing of one affected dog, an obligate carrier and an unrelated control revealed a 218 bp SINE insertion into exon 7 of the RAB3GAP1 gene. The SINE insertion was perfectly associated with the disease phenotype in a cohort of 43 Alaskan Huskies and it was absent from 541 control dogs of diverse other breeds. The SINE insertion induced aberrant splicing and led to a transcript with a greatly altered exon 7. RAB3GAP1 loss-of-function variants in humans cause Warburg Micro Syndrome 1 (WARBM1), which is characterized by additional developmental defects compared to canine POANV, whereas Rab3gap1 deficient mice have a much milder phenotype than either humans or dogs. Thus the RAB3GAP1 mutant Alaskan Huskies provide an interesting intermediate phenotype that may help to better understand the function of RAB3GAP1 in development. Furthermore, the identification of the presumed causative genetic variant will enable genetic testing to avoid the non-intentional breeding of affected dogs.

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.

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.

EuroPineDB: a high-coverage web database for maritime pine transcriptome

Pinus pinaster is an economically and ecologically important species that is becoming a woody gymnosperm model. Its enormous genome size makes whole-genome sequencing approaches are hard to apply. Therefore, the expressed portion of the genome has to be characterised and the results and annotations have to be stored in dedicated databases.

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.