Simultaneous testing of rule- and model-based approaches for runs of homozygosity detection opens up a window into genomic footprints of selection in pigs

Background: Past selection events left footprints in the genome of domestic animals, which can be traced back by stretches of homozygous genotypes, designated as runs of homozygosity (ROHs). The analysis of common ROH regions within groups or populations displaying potential signatures of selection requires high-quality SNP data as well as carefully adjusted ROH-defining parameters. In this study, we used a simultaneous testing of rule- and model-based approaches to perform strategic ROH calling in genomic data from different pig populations to detect genomic regions under selection for specific phenotypes. Results: Our ROH analysis using a rule-based approach offered by PLINK, as well as a model-based approach run by RZooRoH demonstrated a high efficiency of both methods. It underlined the importance of providing a high- quality SNP set as input as well as adjusting parameters based on dataset and population for ROH calling. Particularly, ROHs ≤ 20 kb were called in a high frequency by both tools, but to some extent covered different gene sets in sub- sequent analysis of ROH regions common for investigated pig groups. Phenotype associated ROH analysis resulted in regions under potential selection characterizing heritage pig breeds, known to harbour a long-established breeding history. In particular, the selection focus on fitness-related traits was underlined by various ROHs harbouring disease resistance or tolerance-associated genes. Moreover, we identified potential selection signatures associated with ear morphology, which confirmed known candidate genes as well as uncovered a missense mutation in the ABCA6 gene potentially supporting ear cartilage formation. Conclusions: The results of this study highlight the strengths and unique features of rule- and model-based approaches as well as demonstrate their potential for ROH analysis in animal populations. We provide a workflow for ROH detection, evaluating the major steps from filtering for high-quality SNP sets to intersecting ROH regions. Formula-based estimations defining ROHs for rule-based method show its limits, particularly for efficient detection of smaller ROHs. Moreover, we emphasize the role of ROH detection for the identification of potential footprints of selec- tion in pigs, displaying their breed-specific characteristics or favourable phenotypes

TADA – a Machine Learning Tool for Functional Annotation based Prioritisation of Putative Pathogenic CNVs

Few methods have been developed to investigate copy number variants (CNVs) based on their predicted pathogenicity. We introduce TADA, a method to prioritise pathogenic CNVs through assisted manual filtering and automated classification, based on an extensive catalogue of functional annotation supported by rigourous enrichment analysis. We demonstrate that our classifiers are able to accurately predict pathogenic CNVs, outperforming current alternative methods, and produce a well-calibrated pathogenicity score. Our results suggest that functional annotation-based prioritisation of pathogenic CNVs is a promising approach to support clinical diagnostics and to further the understanding of mechanisms controlling the disease impact of larger genomic alterations

A likelihood ratio based method to predict exact pedigrees for complex families from next-generation sequencing data

MOTIVATION: Next generation sequencing (NGS) technology considerably changed the way we screen for pathogenic mutations in rare Mendelian disorders. However, the identification of the disease- causing mutation amongst thousands of variants of partly unknown relevance is still challenging and efficient techniques that reduce the genomic search space play a decisive role. Often segregation- or linkage analysis are used to prioritize candidates, however, these approaches require correct information about the degree of relationship among the sequenced samples. For quality assurance an automated control of pedigree structures and sample assignment is therefore highly desirable in order to detect label mix-ups that might otherwise corrupt downstream analysis. RESULTS: We developed an algorithm based on likelihood ratios that discriminates between different classes of relationship for an arbitrary number of genotyped samples. By identifying the most likely class we are able to reconstruct entire pedigrees iteratively, even for highly consanguineous families. We tested our approach on exome data of different sequencing studies and achieved high precision for all pedigree predictions. By analyzing the precision for varying degrees of relatedness or inbreeding we could show that a prediction is robust down to magnitudes of a few hundred loci. AVAILABILITY: A java standalone application that computes the relationships between multiple samples as well as a Rscript that visualizes the pedigree information is available for download as well as a web service at www.gene-talk.de CONTACT: [email protected]

Dev Dyn

Ror2 is a receptor tyrosine kinase mutated in the human syndromes Brachydactyly type B (BDB) and recessive Robinow syndrome (RS). In this study, we used the chick as a model to investigate the role of Ror2 in skeletogenesis and to elucidate the functional consequences of Ror2 mutations. For this purpose, we cloned chicken Ror2 and analyzed its expression pattern at various embryonic stages by in situ hybridization and immunolabeling. We document expression of cRor2 in several organs, including mesonephros, heart, nervous system, intestine and cartilage. The high conservation of expression when compared with the mouse underlines the validity of the chick as a model system. Using replication-competent retroviral vector-mediated overexpression, we analyzed the functional consequences of truncating BDB and RS mutations in the developing chick limb. Overexpression of Ror2 mutants led to a disturbance of growth plate architecture and a severe block of chondrocyte differentiation, demonstrating the functional importance of Ror2 in skeletogenesis

Eur J Human Genet

Heterozygous missense mutations in the serine-threonine kinase receptor BMPR1B result typically in brachydactyly type A2 (BDA2), whereas mutations in the corresponding ligand GDF5 cause brachydactyly type C (BDC). Mutations in the GDF inhibitor Noggin (NOG) or activating mutations in GDF5 cause proximal symphalangism (SYM1). Here, we describe a novel mutation in BMPR1B (R486Q) that is associated with either BDA2 or a BDC/SYM1-like phenotype. Functional investigations of the R486Q mutation were performed and compared with the previously reported BDA2-causing mutation R486W and WT BMPR1B. Overexpression of the mutant receptors in chicken micromass cultures resulted in a strong inhibition of chondrogenesis with the R486Q mutant, showing a stronger effect than the R486W mutant. To investigate the consequences of the BMPR1B mutations on the intracellular signal transduction, we used stably transfected C2C12 cells and measured the activity of SMAD-dependent and SMAD-independent pathways. SMAD activation after stimulation with GDF5 was suppressed in both mutants. Alkaline phosphatase induction showed an almost complete loss of activation by both mutants. Our data extend the previously known mutational and phenotypic spectrum associated with mutations in BMPR1B. Disturbances of NOG-GDF5-BMPR1B signaling cascade can result in similar clinical manifestations depending on the quantitative effect and mode of action of the specific mutations within the same functional pathway

A novel COL1A2 C-propeptide cleavage site mutation causing high bone mass osteogenesis imperfecta with a regional distribution pattern

Osteogenesis imperfecta (OI) is typically characterized by low bone mass and increased bone fragility caused by heterozygous mutations in the type I procollagen genes (COL1A1/COL1A2). We report two cases of a 56-year-old woman and her 80-year-old mother who suffered from multiple vertebral and non-vertebral fractures with onset in early childhood. A full osteologic assessment including dual-energy X-ray absorptiometry (DXA), high-resolution peripheral quantitative computed tomography (HR-pQCT), and serum analyses pointed to a high bone mineral density (BMD) in the hip (DXA Z-score + 3.7 and + 3.9) but low to normal bone mass in the spine and preserved bone microstructure in the distal tibia. Serum markers of bone formation and bone resorption were elevated. Using whole exome sequencing, we identified a novel mutation in the COL1A2 gene causing a p. (Asp1120Gly) substitution at the protein level and affecting the type I procollagen C-propeptide cleavage site. In line with previously reported cases, our data independently prove the existence of an unusual phenotype of high bone mass OI caused by a mutation in the procollagen C-propeptide cleavage with a clinically persistent phenotype through adulthood

Identification of a molecular defect in a stillborn fetus with perinatal lethal hypophosphatasia using a disease-associated genome sequencing approach

Lethal skeletal disorders represent a heterogeneous and clinically variable group of genetic conditions, usually difficult to diagnose without post-mortem radiological assessment. Here we report on a stillborn patient delivered at 22 weeks of gestation who presented with severe skeletal symptoms comprising limb shortening and intrauterine fractures detected upon prenatal ultrasound and autopsy examination. Since post-mortem X-ray was refused and no phenotypic diagnosis could be attempted, we performed next-generation sequencing (NGS) of 2741 genes associated with all known Mendelian disorders. With this strategy, we were able to demonstrate the diagnosis at a molecular level, which turned out to be perinatal lethal hypophosphatasia (HPP). This severe form of HPP represents an inborn defect of ossification often resulting in stillbirth or postnatal death. The NGS panel revealed compound heterozygous ALPL missense mutations: c.1283G>C(p.Arg428Pro) and c.1363G>A(p.Gly455Ser). Mutations detected in our case, although previously described in other patients, have not been reported to co-occur in a single individual. The diagnosis established in our index using the NGS-based approach could have been successfully reached by standard radiography. Thus, our report points to the importance of X-ray examination in stillborn cases and highlights the emerging role of NGS strategies in the diagnostic process of prenatally manifesting skeletal disorders

Am J Med Genet A

We report on a girl with delayed mental and motor development, ophthalmological abnormalities, and peripheral neuropathy. Chromosome analysis suggested a deletion within chromosome 8p. Further investigation by array-based comparative genomic hybridization (array-CGH) delineated an 8 Mb interstitial deletion on the short arm of chromosome 8. The breakpoints are located at chromosome bands 8p12 and 8p21.2. Forty-two known genes including gonadotropin-releasing hormone 1 (GNRH1), transcription factor EBF2, exostosin-like 3 (EXTL3), glutathione reductase (GSR), and neuregulin 1 (NRG1), are located within the deleted region on chromosome 8p. A comparison of our patient with the cases described in the literature is presented, and we discuss the genotype-phenotype correlation in our patient. This is the first report of array-CGH analysis of an interstitial deletion at chromosome 8p

Xq27.1 palindrome mediated interchromosomal insertion likely causes familial congenital bilateral laryngeal abductor paralysis (Plott syndrome)

Bilateral laryngeal abductor paralysis is a rare entity and the second most common cause of stridor in newborns. So far, no conclusive genetic or chromosomal aberration has been reported for X-linked isolated bilateral vocal cord paralysis, also referred to as Plott syndrome. Via whole genome sequencing (WGS), we identified a complex interchromosomal insertion in a large family with seven affected males. The 404 kb inserted fragment originates from chromosome 10q21.3, contains no genes and is inserted inversionally into the intergenic chromosomal region Xq27.1, 82 kb centromeric to the nearest gene SOX3. The patterns found at the breakpoint junctions resemble typical characteristics that arise in replication-based mechanisms with long-distance template switching. Non protein-coding insertions into the same genomic region have been described to result in different phenotypes, indicating that the phenotypic outcome likely depends on the introduction of regulatory elements. In conclusion, our data adds Plott syndrome as another entity, likely caused by the insertion of non-coding DNA into the intergenic chromosomal region Xq27.1. In this regard, we demonstrate the importance of WGS as a powerful diagnostic test in unsolved genetic diseases, as this genomic rearrangement has not been detected by current first-line diagnostic tests, i.e., exome sequencing and chromosomal microarray analysis

Structure and deformation of the Kermadec forearc in response to subduction of the Pacific oceanic plate

The Tonga-Kermadec forearc is deforming in response to on-going subduction of the Pacific Plate beneath the Indo-Australian Plate. Previous research has focussed on the structural development of the forearc where large bathymetric features such as the Hikurangi Plateau and Louisville Ridge seamount chain are being subducted. Consequently, knowledge of the ‘background’ forearc in regions of normal plate convergence is limited. We report on an ∼250-km-long multichannel seismic reflection profile that was shot perpendicular to the Tonga-Kermadec trench at ∼28°S to determine the lateral and temporal variations in the structure, stratigraphy and deformation of the Kermadec forearc resulting solely from Pacific Plate subduction. Interpretation of the seismic profile, in conjunction with regional swath bathymetry data, shows that the Pacific Plate exhibits horst and graben structures that accommodate bending-induced extensional stresses, generated as the trenchward dip of the crust increases. Trench infill is also much thicker than expected at 1 km which, we propose, results from increased sediment flux into and along the trench. Pervasive normal faulting of the mid-trench slope most likely accommodates the majority of the observed forearc extension in response to basal subduction erosion, and a structural high is located between the mid- and upper-trench slopes. We interpret this high as representing a dense and most likely structurally robust region of crust lying beneath this region. Sediment of the upper-trench slope documents depositional hiatuses and on-going uplift of the arc. Strong along-arc currents appear to erode the Kermadec volcanic arc and distribute this sediment to the surrounding basins, while currents over the forearc redistribute deposits as sediment waves. Minor uplift of the transitional Kermadec forearc, observed just to the north of the profile, appears to relate to an underlying structural trend as well as subduction of the Louisville Ridge seamount chain 250 km to the north. Relative uplift of the Kermadec arc is observed from changes in the tilt of upper-trench slope deposits and extensional faulting of the basement immediately surrounding the Louisville Ridge