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.

Virulence of broad- and narrow-host-range Salmonella enterica serovars in the streptomycin-pretreated mouse model.

Salmonella enterica subspecies I serovars are common bacterial pathogens causing diseases ranging from enterocolitis to systemic infections. Some serovars are adapted to specific hosts, whereas others have a broad host range. The molecular mechanisms defining the virulence characteristics and the host range of a given S. enterica serovar are unknown. Streptomycin pretreated mice provide a surrogate host model for studying molecular aspects of the intestinal inflammation (colitis) caused by serovar Typhimurium (S. Hapfelmeier and W. D. Hardt, Trends Microbiol. 13:497-503, 2005). Here, we studied whether this animal model is also useful for studying other S. enterica subspecies I serovars. All three tested strains of the broad-host-range serovar Enteritidis (125109, 5496/98, and 832/99) caused pronounced colitis and systemic infection in streptomycin pretreated mice. Different levels of virulence were observed among three tested strains of the host-adapted serovar Dublin (SARB13, SD2229, and SD3246). Several strains of host restricted serovars were also studied. Two serovar Pullorum strains (X3543 and 449/87) caused intermediate levels of colitis. No intestinal inflammation was observed upon infection with three different serovar Paratyphi A strains (SARB42, 2804/96, and 5314/98) and one serovar Gallinarum strain (X3796). A second serovar Gallinarum strain (287/91) was highly virulent and caused severe colitis. This strain awaits future analysis. In conclusion, the streptomycin pretreated mouse model can provide an additional tool to study virulence factors (i.e., those involved in enteropathogenesis) of various S. enterica subspecies I serovars. Five of these strains (125109, 2229, 287/91, 449/87, and SARB42) are subject of Salmonella genome sequencing projects. The streptomycin pretreated mouse model may be useful for testing hypotheses derived from this genomic data.

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.

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

Proteomic analysis of Plasmodium in the mosquito: progress and pitfalls

Here we discuss proteomic analyses of whole cell preparations of the mosquito stages of malaria parasite development (i.e. gametocytes, microgamete, ookinete, oocyst and sporozoite) of Plasmodium berghei. We also include critiques of the proteomes of two cell fractions from the purified ookinete, namely the micronemes and cell surface. Whereas we summarise key biological interpretations of the data, we also try to identify key methodological constraints we have met, only some of which we were able to resolve. Recognising the need to translate the potential of current genome sequencing into functional understanding, we report our efforts to develop more powerful combinations of methods for the in silico prediction of protein function and location. We have applied this analysis to the proteome of the male gamete, a cell whose very simple structural organisation facilitated interpretation of data. Some of the in silico predictions made have now been supported by ongoing protein tagging and genetic knockout studies. We hope this discussion may assist future studies.

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.

New insights into the performance of human whole-exome capture platforms.

Whole exome sequencing (WES) is increasingly used in research and diagnostics. WES users expect coverage of the entire coding region of known genes as well as sufficient read depth for the covered regions. It is, however, unknown which recent WES platform is most suitable to meet these expectations. We present insights into the performance of the most recent standard exome enrichment platforms from Agilent, NimbleGen and Illumina applied to six different DNA samples by two sequencing vendors per platform. Our results suggest that both Agilent and NimbleGen overall perform better than Illumina and that the high enrichment performance of Agilent is stable among samples and between vendors, whereas NimbleGen is only able to achieve vendor- and sample-specific best exome coverage. Moreover, the recent Agilent platform overall captures more coding exons with sufficient read depth than NimbleGen and Illumina. Due to considerable gaps in effective exome coverage, however, the three platforms cannot capture all known coding exons alone or in combination, requiring improvement. Our data emphasize the importance of evaluation of updated platform versions and suggest that enrichment-free whole genome sequencing can overcome the limitations of WES in sufficiently covering coding exons, especially GC-rich regions, and in characterizing structural variants.

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.

Genetic redundancy among durum wheat accessions as assessed with SSRs and endosperm proteins

Reducing duplication in ex-situ collections is complicated and requires good quality genetic markers. This study was conducted to assess the value of endosperm proteins and SSRs for validation of potential duplicates and monitoring intra-accession variability. Fifty durum wheat (Triticum turgidum ssp. durum) accessions grouped in 23 potential duplicates, and previously characterised for 30 agro-morphological traits, were analysed for gliadin and high molecular weight glutenin (HMWG) subunit alleles, total protein, and 24 SSRs, covering a wide genome area. Similarity and dissimilarity matrices were generated based on protein and SSRs alleles. For heterogeneous accessions at gliadins the percent pattern homology (PH) between gliadin patterns and the Nei’s coefficient of genetic identity (I) were computed. Eighteen duplicates identical for proteins showed none or less than 3 unshared SSRs alleles. For heterogeneous accessions PH and I values lower than 80 identified clearly off-types with more than 3 SSRs unshared. Only those biotypes differing in no more than one protein-coding locus were confirmed with SSRs. A good concordance among proteins, morphological traits, and SSR were detected. However, the discrepancy in similarity detected in some cases showed that it is advisable to evaluate redundancy through distinct approaches. The analysis in proteins together with SSRs data are very useful to identify duplicates, biotypes, close related genotypes, and contaminations

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.

Identification of a family of zinc transporter genes from Arabidopsis that respond to zinc deficiency

Millions of people worldwide suffer from nutritional imbalances of essential metals like zinc. These same metals, along with pollutants like cadmium and lead, contaminate soils at many sites around the world. In addition to posing a threat to human health, these metals can poison plants, livestock, and wildlife. Deciphering how metals are absorbed, transported, and incorporated as protein cofactors may help solve both of these problems. For example, edible plants could be engineered to serve as better dietary sources of metal nutrients, and other plant species could be tailored to remove metal ions from contaminated soils. We report here the cloning of the first zinc transporter genes from plants, the ZIP1, ZIP2, and ZIP3 genes of Arabidopsis thaliana. Expression in yeast of these closely related genes confers zinc uptake activities. In the plant, ZIP1 and ZIP3 are expressed in roots in response to zinc deficiency, suggesting that they transport zinc from the soil into the plant. Although expression of ZIP2 has not been detected, a fourth related Arabidopsis gene identified by genome sequencing, ZIP4, is induced in both shoots and roots of zinc-limited plants. Thus, ZIP4 may transport zinc intracellularly or between plant tissues. These ZIP proteins define a family of metal ion transporters that are found in plants, protozoa, fungi, invertebrates, and vertebrates, making it now possible to address questions of metal ion accumulation and homeostasis in diverse organisms.