Assessing Pathogenicity for Novel Mutation/Sequence Variants: The Value of Healthy Older Individuals

Improvement in DNA technology is increasingly revealing unexpected/unknown mutations in healthy persons and generating anxiety due to their still unknown health consequences. We report a 44-year-old healthy father of a 10-year-old daughter with bilateral coloboma and hearing loss, but without muscle weakness, in whom a whole-genome CGH revealed a deletion of exons 38-44 in the dystrophin gene. This mutation was inherited from her asymptomatic father, who was further clinically and molecularly evaluated for prognosis and genetic counseling (GC). This deletion was never identified by us in 982 Duchenne/Becker patients. To assess whether the present case represents a rare case of non-penetrance, and aiming to obtain more information for prognosis and GC, we suggested that healthy older relatives submit their DNA for analysis, to which several complied. Mutation analysis revealed that his mother, brother, and 56-year-old maternal uncle also carry the 38-44 deletion, suggesting it an unlikely cause of muscle weakness. Genome sequencing will disclose mutations and variants whose health impact are still unknown, raising important problems in interpreting results, defining prognosis, and discussing GC. We suggest that, in addition to family history, keeping the DNA of older relatives could be very informative, in particular for those interested in having their genome sequenced.

Xylella fastidiosa: An in vivo system to study possible survival strategies within citrus xylem vessels based on global gene expression analysis

Xylella fastidiosa inhabits the plant xylem, a nutrient-poor environment, so that mechanisms to sense and respond to adverse environmental conditions are extremely important for bacterial survival in the plant host. Although the complete genome sequences of different Xylella strains have been determined, little is known about stress responses and gene regulation in these organisms. In this work, a DNA microarray was constructed containing 2,600 ORFs identified in the genome sequencing project of Xylella fastidiosa 9a5c strain, and used to check global gene expression differences in the bacteria when it is infecting a symptomatic and a tolerant citrus tree. Different patterns of expression were found in each variety, suggesting that bacteria are responding differentially according to each plant xylem environment. The global gene expression profile was determined and several genes related to bacterial survival in stressed conditions were found to be differentially expressed between varieties, suggesting the involvement of different strategies for adaptation to the environment. The expression pattern of some genes related to the heat shock response, toxin and detoxification processes, adaptation to atypical conditions, repair systems as well as some regulatory genes are discussed in this paper. DNA microarray proved to be a powerful technique for global transcriptome analyses. This is one of the first studies of Xylella fastidiosa gene expression in vivo which helped to increase insight into stress responses and possible bacterial survival mechanisms in the nutrient-poor environment of xylem vessels.

Trypanosoma cruzi mitochondrial maxicircles display species- and strain-specific variation and a conserved element in the non-coding region

Abstract Background The mitochondrial DNA of kinetoplastid flagellates is distinctive in the eukaryotic world due to its massive size, complex form and large sequence content. Comprised of catenated maxicircles that contain rRNA and protein-coding genes and thousands of heterogeneous minicircles encoding small guide RNAs, the kinetoplast network has evolved along with an extreme form of mRNA processing in the form of uridine insertion and deletion RNA editing. Many maxicircle-encoded mRNAs cannot be translated without this post-transcriptional sequence modification. Results We present the complete sequence and annotation of the Trypanosoma cruzi maxicircles for the CL Brener and Esmeraldo strains. Gene order is syntenic with Trypanosoma brucei and Leishmania tarentolae maxicircles. The non-coding components have strain-specific repetitive regions and a variable region that is unique for each strain with the exception of a conserved sequence element that may serve as an origin of replication, but shows no sequence identity with L. tarentolae or T. brucei. Alternative assemblies of the variable region demonstrate intra-strain heterogeneity of the maxicircle population. The extent of mRNA editing required for particular genes approximates that seen in T. brucei. Extensively edited genes were more divergent among the genera than non-edited and rRNA genes. Esmeraldo contains a unique 236-bp deletion that removes the 5'-ends of ND4 and CR4 and the intergenic region. Esmeraldo shows additional insertions and deletions outside of areas edited in other species in ND5, MURF1, and MURF2, while CL Brener has a distinct insertion in MURF2. Conclusion The CL Brener and Esmeraldo maxicircles represent two of three previously defined maxicircle clades and promise utility as taxonomic markers. Restoration of the disrupted reading frames might be accomplished by strain-specific RNA editing. Elements in the non-coding region may be important for replication, transcription, and anchoring of the maxicircle within the kinetoplast network.

New approaches to open problems in gene expression microarray data

In the past decade, the advent of efficient genome sequencing tools and high-throughput experimental biotechnology has lead to enormous progress in the life science. Among the most important innovations is the microarray tecnology. It allows to quantify the expression for thousands of genes simultaneously by measurin the hybridization from a tissue of interest to probes on a small glass or plastic slide. The characteristics of these data include a fair amount of random noise, a predictor dimension in the thousand, and a sample noise in the dozens. One of the most exciting areas to which microarray technology has been applied is the challenge of deciphering complex disease such as cancer. In these studies, samples are taken from two or more groups of individuals with heterogeneous phenotypes, pathologies, or clinical outcomes. these samples are hybridized to microarrays in an effort to find a small number of genes which are strongly correlated with the group of individuals. Eventhough today methods to analyse the data are welle developed and close to reach a standard organization (through the effort of preposed International project like Microarray Gene Expression Data -MGED- Society [1]) it is not unfrequant to stumble in a clinician's question that do not have a compelling statistical method that could permit to answer it.The contribution of this dissertation in deciphering disease regards the development of new approaches aiming at handle open problems posed by clinicians in handle specific experimental designs. In Chapter 1 starting from a biological necessary introduction, we revise the microarray tecnologies and all the important steps that involve an experiment from the production of the array, to the quality controls ending with preprocessing steps that will be used into the data analysis in the rest of the dissertation. While in Chapter 2 a critical review of standard analysis methods are provided stressing most of problems that In Chapter 3 is introduced a method to adress the issue of unbalanced design of miacroarray experiments. In microarray experiments, experimental design is a crucial starting-point for obtaining reasonable results. In a two-class problem, an equal or similar number of samples it should be collected between the two classes. However in some cases, e.g. rare pathologies, the approach to be taken is less evident. We propose to address this issue by applying a modified version of SAM [2]. MultiSAM consists in a reiterated application of a SAM analysis, comparing the less populated class (LPC) with 1,000 random samplings of the same size from the more populated class (MPC) A list of the differentially expressed genes is generated for each SAM application. After 1,000 reiterations, each single probe given a "score" ranging from 0 to 1,000 based on its recurrence in the 1,000 lists as differentially expressed. The performance of MultiSAM was compared to the performance of SAM and LIMMA [3] over two simulated data sets via beta and exponential distribution. The results of all three algorithms over low- noise data sets seems acceptable However, on a real unbalanced two-channel data set reagardin Chronic Lymphocitic Leukemia, LIMMA finds no significant probe, SAM finds 23 significantly changed probes but cannot separate the two classes, while MultiSAM finds 122 probes with score >300 and separates the data into two clusters by hierarchical clustering. We also report extra-assay validation in terms of differentially expressed genes Although standard algorithms perform well over low-noise simulated data sets, multi-SAM seems to be the only one able to reveal subtle differences in gene expression profiles on real unbalanced data. In Chapter 4 a method to adress similarities evaluation in a three-class prblem by means of Relevance Vector Machine [4] is described. In fact, looking at microarray data in a prognostic and diagnostic clinical framework, not only differences could have a crucial role. In some cases similarities can give useful and, sometimes even more, important information. The goal, given three classes, could be to establish, with a certain level of confidence, if the third one is similar to the first or the second one. In this work we show that Relevance Vector Machine (RVM) [2] could be a possible solutions to the limitation of standard supervised classification. In fact, RVM offers many advantages compared, for example, with his well-known precursor (Support Vector Machine - SVM [3]). Among these advantages, the estimate of posterior probability of class membership represents a key feature to address the similarity issue. This is a highly important, but often overlooked, option of any practical pattern recognition system. We focused on Tumor-Grade-three-class problem, so we have 67 samples of grade I (G1), 54 samples of grade 3 (G3) and 100 samples of grade 2 (G2). The goal is to find a model able to separate G1 from G3, then evaluate the third class G2 as test-set to obtain the probability for samples of G2 to be member of class G1 or class G3. The analysis showed that breast cancer samples of grade II have a molecular profile more similar to breast cancer samples of grade I. Looking at the literature this result have been guessed, but no measure of significance was gived before.

Novel isolates of Cydia pomonella granulovirus (CpGV): deciphering the molecular mechanism for overcoming CpGV resistance in codling moth (Cydia pomonella)

During the last twenty years, Cydia pomonolla granulovirus (CpGV, Baculoviridae) has become the most important biological control agent for the codling moth (CM) in organic and integrated apple production. All registered products in Europe are based on the isolate CpGV-M, which was discovered 1964 in Mexico. A serious threat to future application of CpGV is the occurrence of CM field populations resistant to CpGV. Since 2003, populations with up to 10,000-fold reduced susceptibility were reported from orchards in Germany, France, Italy, Switzerland, Austria and the Netherlands. A putative alternative to CpGV-M are novel CpGV isolates which are able to overcome CM resistance. This thesis focuses on the identification and characterisation of resistance overcoming CpGV isolates and the analysis of their molecular difference to CpGV-M.rnSixteen CpGV isolates were tested against CM lab strains in bioassays. Hereby, five isolates were identified which were able to completely overcome resistance. The genomes of these isolates were compared to CpGV-M by restriction fragment length polymorphism (RFLP) analysis. To identify the molecular factor responsible for improved virulence of some CpGV isolates, major genomic differences were sequenced and analysed. A 0.7 kb insertion was found in CpGV-I01, -I12 and -E2, but not in other resistance overcoming isolates. Analysis of the insertions sequence revealed that it might be due to a transposition event, but not involved in overcoming resistance. rnFor unequivocal identification of CpGV isolates, a new method based on molecular analysis was established. Partial sequencing of the conserved polyhedrin/granulin (polh/gran), late expression factor-8 (lef-8) and late expression factor-9 (lef-9) genes revealed single nucleotide polymorphisms (SNPs). SNP analysis correlated with the grouping obtained by RFLP analysis. A phylogenetic classification due to different genome types A-E is proposed. Phylogenetic analysis suggested that CpGV-M was the phylogenetically youngest of the tested CpGV isolates.rnWhole genome sequencing of two resistance overcoming isolates CpGV-I12 (type D genome) and -S (type E genome) and CpGV-M (type A genome) was performed. Comparison of the three genomes revealed a high sequence identity. Several insertions and deletions ranging from 1-700 nucleotides (nt) were found. Comparison on open reading frame (ORF) level revealed that CpGV-I12 and -S shared only one protein alteration when compared to CpGV-M: a stretch of 24 nt present in ORF cp24 was not found in any of the resistance overcoming isolates. Cp24 codes for the early gene pe38. Combined with the results of phylogenetic analysis, it is proposed that these 24 nt are a recent insertion into the CpGV-M genome. The role of pe38 in overcoming resistance was investigated by knocking out pe38 of a CpGV-M based bacmid and swapping of CpGV-I12 pe38 of into the k.o. bacmid. When pe38 of CpGV-I12 was inserted into the k.o. bacmid, the infectivity could not be rescued, suggesting that the genomic portion of pe38 might play a role in its function.rnIt can be concluded that the recently observed CpGV resistance in CM is only related to type A genomes. RFLP and SNP analysis provide tools for identifying and characterising different CpGV isolates reliably, a pre-condition for a future registration of CpGV products based on novel CpGV isolates.rnrnrn

Holocentric chromosomes: convergent evolution, meiotic adaptations, and genomic analysis

In most eukaryotes, the kinetochore protein complex assembles at a single locus termed the centromere to attach chromosomes to spindle microtubules. Holocentric chromosomes have the unusual property of attaching to spindle microtubules along their entire length. Our mechanistic understanding of holocentric chromosome function is derived largely from studies in the nematode Caenorhabditis elegans, but holocentric chromosomes are found over a broad range of animal and plant species. In this review, we describe how holocentricity may be identified through cytological and molecular methods. By surveying the diversity of organisms with holocentric chromosomes, we estimate that the trait has arisen at least 13 independent times (four times in plants and at least nine times in animals). Holocentric chromosomes have inherent problems in meiosis because bivalents can attach to spindles in a random fashion. Interestingly, there are several solutions that have evolved to allow accurate meiotic segregation of holocentric chromosomes. Lastly, we describe how extensive genome sequencing and experiments in nonmodel organisms may allow holocentric chromosomes to shed light on general principles of chromosome segregation.

Bioinformatics tools and databases for the study of human growth hormone

Advances in novel molecular biological diagnostic methods are changing the way of diagnosis and study of metabolic disorders like growth hormone deficiency. Faster sequencing and genotyping methods require strong bioinformatics tools to make sense of the vast amount of data generated by modern laboratories. Advances in genome sequencing and computational power to analyze the whole genome sequences will guide the diagnostics of future. In this chapter, an overview of some basic bioinformatics resources that are needed to study metabolic disorders are reviewed and some examples of bioinformatics analysis of human growth hormone gene, protein and structure are provided.

A Staphylococcus xylosus isolate with a new mecC allotype

Recently, a novel variant of mecA known as mecC (mecA(LGA251)) was identified in Staphylococcus aureus isolates from both humans and animals. In this study, we identified a Staphylococcus xylosus isolate that harbors a new allotype of the mecC gene, mecC1. Whole-genome sequencing revealed that mecC1 forms part of a class E mec complex (mecI-mecR1-mecC1-blaZ) located at the orfX locus as part of a likely staphylococcal cassette chromosome mec element (SCCmec) remnant, which also contains a number of other genes present on the type XI SCCmec.

Novel pseudo-staphylococcal cassette chromosome mec element (ψSCCmec57395) in methicillin-resistant Staphylococcus pseudintermedius CC45

Genetic characterization of methicillin-resistant Staphylococcus pseudintermedius (MRSP) from Thailand and Israel revealed the presence of a predominant atypical clonal lineage which was not typeable by SmaI-PFGE and SCCmec typing. All the atypical isolates (n = 34) belonged to CC45 (30 ST45 and 2 ST179 isolates, 1 ST57 isolate, and 1 ST85 isolate). The isolates originated from healthy and diseased dogs and cats, as well as from the environment of one clinic. Cfr9I-pulsed-field gel electrophoresis (Cfr9I-PFGE) and dru typing permitted the further distinction of CC45 isolates from the two different countries. Microarray analysis identified genes that confer resistance to β-lactams (mecA; blaZ), aminoglycosides [aac(6')-Ie-aph(2')-Ia; aph(3')-III; ant(6)-Ia], macrolides and lincosamides [erm(B)], tetracyclines [tet(M)], trimethoprim [dfr(G)], streptothricin (sat4), and chloramphenicol (catpC221). Fluoroquinolone resistance was attributed to specific amino acid substitutions, i.e., Ser84Leu in GyrA and Ser80Ile and Asp84Asn in GrlA. A novel pseudo-staphylococcal cassette chromosome (ΨSCCmec57395) element was identified in MRSP strain 57395 (sequence type ST45) by whole-genome sequencing. The 12,282-bp ΨSCCmec57395 element contained a class C1 mec gene complex but no ccr genes. In addition to the methicillin resistance gene mecA, ΨSCCmec57395 also carried determinants of resistance to heavy metals, such as arsenic, cadmium, and copper. Bsu36I restriction analysis of the ΨSCCmec57395 element amplified by long-range PCR revealed the presence of ΨSCCmec57395 in the 33 additional isolates of MRSP CC45. The ΨSCCmec57395 element represents a new class of SCCmec and has been identified in MRSP of CC45, which is a predominant clonal lineage in Israel and Thailand.

A point mutation in cpsE renders Streptococcus pneumoniae nonencapsulated and enhances its growth, adherence and competence

BackgroundThe polysaccharide capsule is a major virulence factor of the important human pathogen Streptococcus pneumoniae. However, S. pneumoniae strains lacking capsule do occur.ResultsHere, we report a nasopharyngeal isolate of Streptococcus pneumoniae composed of a mixture of two phenotypes; one encapsulated (serotype 18C) and the other nonencapsulated, determined by serotyping, electron microscopy and fluorescence isothiocyanate dextran exclusion assay.By whole genome sequencing, we demonstrated that the phenotypes differ by a single nucleotide base pair in capsular gene cpsE (C to G change at gene position 1135) predicted to result in amino acid change from arginine to glycine at position 379, located in the cytoplasmic, enzymatically active, region of this transmembrane protein. This SNP is responsible for loss of capsule production as the phenotype is transferred with the capsule operon. The nonencapsulated variant is superior in growth in vitro and is also 117-fold more adherent to and more invasive into Detroit 562 human epithelial cells than the encapsulated variant.Expression of six competence pathway genes and one competence-associated gene was 11 to 34-fold higher in the nonencapsulated variant than the encapsulated and transformation frequency was 3.7-fold greater.ConclusionsWe identified a new single point mutation in capsule gene cpsE of a clinical S. pneumoniae serotype 18C isolate sufficient to cause loss of capsule expression resulting in the co-existence of the encapsulated and nonencapsulated phenotype. The mutation caused phenotypic changes in growth, adherence to epithelial cells and transformability. Mutation in capsule gene cpsE may be a way for S. pneumoniae to lose its capsule and increase its colonization potential.

Global phylogenomic analysis of nonencapsulated Streptococcus pneumoniae reveals a deep-branching classic lineage that is distinct from multiple sporadic lineages.

The surrounding capsule of Streptococcus pneumoniae has been identified as a major virulence factor and is targeted by pneumococcal conjugate vaccines (PCV). However, nonencapsulated Streptococcus pneumoniae (Non-Ec-Sp) have also been isolated globally, mainly in carriage studies. It is unknown if Non-Ec-Sp evolve sporadically, if they have high antibiotic non-susceptiblity rates and a unique, specific gene content. Here, whole genome sequencing of 131 Non-Ec-Sp isolates sourced from 17 different locations around the world was performed. Results revealed a deep-branching classic lineage that is distinct from multiple sporadic lineages. The sporadic lineages clustered with a previously sequenced, global collection of encapsulated S. pneumoniae (Ec-Sp) isolates while the classic lineage is comprised mainly of the frequently identified multi-locus sequences types ST344 (n=39) and ST448 (n=40). All ST344 and nine ST448 isolates had high non-susceptiblity rates to β-lactams and other antimicrobials. Analysis of the accessory genome reveals that the classic Non-Ec-Sp contained an increased number of mobile elements, than Ec-Sp and sporadic Non-Ec-Sp. Performing adherence assays to human epithelial cells for selected classic and sporadic Non-Ec-Sp revealed that the presence of a integrative conjugative element (ICE) results in increased adherence to human epithelial cells (P=0.005). In contrast, sporadic Non-Ec-Sp lacking the ICE had greater growth in vitro possibly resulting in improved fitness. In conclusion, Non-Ec-Sp isolates from the classic lineage have evolved separately. They have spread globally, are well adapted to nasopharyngeal carriage and are able to coexist with Ec-Sp. Due to continued use of pneumococcal conjugate vaccines, Non-Ec-Sp may become more prevalent.

Independent polled mutations leading to complex gene expression differences in cattle

The molecular regulation of horn growth in ruminants is still poorly understood. To investigate this process, we collected 1019 hornless (polled) animals from different cattle breeds. High-density SNP genotyping confirmed the presence of two different polled associated haplotypes in Simmental and Holstein cattle co-localized on BTA 1. We refined the critical region of the Simmental polled mutation to 212 kb and identified an overlapping region of 932 kb containing the Holstein polled mutation. Subsequently, whole genome sequencing of polled Simmental and Holstein cows was used to determine polled associated genomic variants. By genotyping larger cohorts of animals with known horn status we found a single perfectly associated insertion/deletion variant in Simmental and other beef cattle confirming the recently published possible Celtic polled mutation. We identified a total of 182 sequence variants as candidate mutations for polledness in Holstein cattle, including an 80 kb genomic duplication and three SNPs reported before. For the first time we showed that hornless cattle with scurs are obligate heterozygous for one of the polled mutations. This is in contrast to published complex inheritance models for the bovine scurs phenotype. Studying differential expression of the annotated genes and loci within the mapped region on BTA 1 revealed a locus (LOC100848215), known in cow and buffalo only, which is higher expressed in fetal tissue of wildtype horn buds compared to tissue of polled fetuses. This implicates that the presence of this long noncoding RNA is a prerequisite for horn bud formation. In addition, both transcripts associated with polledness in goat and sheep (FOXL2 and RXFP2), show an overexpression in horn buds confirming their importance during horn development in cattle.

The novel macrolide-Lincosamide-Streptogramin B resistance gene erm(44) is associated with a prophage in Staphylococcus xylosus.

A novel erythromycin ribosome methylase gene, erm(44), that confers resistance to macrolide, lincosamide, and streptogramin B (MLSB) antibiotics was identified by whole-genome sequencing of the chromosome of Staphylococcus xylosus isolated from bovine mastitis milk. The erm(44) gene is preceded by a regulatory sequence that encodes two leader peptides responsible for the inducible expression of the methylase gene, as demonstrated by cloning in Staphylococcus aureus. The erm(44) gene is located on a 53-kb putative prophage designated ΦJW4341-pro. The 56 predicted open reading frames of ΦJW4341-pro are structurally organized into the five functional modules found in members of the family Siphoviridae. ΦJW4341-pro is site-specifically integrated into the S. xylosus chromosome, where it is flanked by two perfect 19-bp direct repeats, and exhibits the ability to circularize. The presence of erm(44) in three additional S. xylosus strains suggests that this putative prophage has the potential to disseminate MLSB resistance.

Bacillus Calmette-Guérin strain differences have an impact on clinical outcome in bladder cancer immunotherapy

BACKGROUND Whether the commonly used bacillus Calmette-Guérin (BCG) strains Connaught and Tice confer different treatment responses in non-muscle-invasive bladder cancer (NMIBC) is unknown. OBJECTIVES To compare clinical efficacy, immunogenicity, and genetics of BCG Connaught and Tice. DESIGN, SETTING, AND PARTICIPANTS A prospective randomized single-institution trial with treatment of 142 high-risk NMIBC patients with BCG Connaught or Tice. INTERVENTION Patients were randomized to receive six instillations of BCG Connaught or Tice. For experimental studies, BCG strains were compared in C57Bl/6 mice. Bladders and lymphoid tissues were analyzed by cytometry and the latter cultivated to detect live BCG. BCG genomic DNA was sequenced and compared with reference genomes. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS Recurrence-free survival was the primary end point of the clinical study. The Kaplan-Meier estimator was used for estimating survival and time-to-event end points. Nonparametric tests served for the analysis of the in vivo results. RESULTS AND LIMITATIONS Treatment with BCG Connaught conferred significantly greater 5-yr recurrence-free survival compared with treatment with BCG Tice (p=0.0108). Comparable numbers of patients experienced BCG therapy-related side effects in each treatment group (p=0.09). In mice, BCG Connaught induced stronger T-helper cell 1-biased responses, greater priming of BCG-specific CD8(+) T cells, and more robust T-cell recruitment to the bladder than BCG Tice. Genome sequencing of the BCG strains revealed candidate genes potentially involved in the differential clinical responses. CONCLUSIONS BCG strain may have an impact on treatment outcome in NMIBC immunotherapy. PATIENT SUMMARY We compared the efficacy of two commonly used bacillus Calmette-Guérin (BCG) strains for the treatment of NMIBC and found that treatment with BCG Connaught prevented recurrences more efficiently than BCG Tice. Comparison of the immunogenicity of the two strains in mice indicated superior immunogenicity of BCG Connaught. We also identified genetic differences that may explain the differential efficacy of the Connaught and Tice BCG strains. TRIAL REGISTRATION NCT00003779.

Epidermolysis bullosa in Danish Hereford calves is caused by a deletion in LAMC2 gene.

BACKGROUND Heritable forms of epidermolysis bullosa (EB) constitute a heterogeneous group of skin disorders of genetic aetiology that are characterised by skin and mucous membrane blistering and ulceration in response to even minor trauma. Here we report the occurrence of EB in three Danish Hereford cattle from one herd. RESULTS Two of the animals were necropsied and showed oral mucosal blistering, skin ulcerations and partly loss of horn on the claws. Lesions were histologically characterized by subepidermal blisters and ulcers. Analysis of the family tree indicated that inbreeding and the transmission of a single recessive mutation from a common ancestor could be causative. We performed whole genome sequencing of one affected calf and searched all coding DNA variants. Thereby, we detected a homozygous 2.4 kb deletion encompassing the first exon of the LAMC2 gene, encoding for laminin gamma 2 protein. This loss of function mutation completely removes the start codon of this gene and is therefore predicted to be completely disruptive. The deletion co-segregates with the EB phenotype in the family and absent in normal cattle of various breeds. Verifying the homozygous private variants present in candidate genes allowed us to quickly identify the causative mutation and contribute to the final diagnosis of junctional EB in Hereford cattle. CONCLUSIONS Our investigation confirms the known role of laminin gamma 2 in EB aetiology and shows the importance of whole genome sequencing in the analysis of rare diseases in livestock.