Project description:Dystrophinopathies caused by variants in the DMD gene are a well-studied muscle disease. The most common type of variant in DMD are large deletions. Very rarely reported forms of variants are chromosomal translocations, inversions and deep intronic variants (DIVs) because they are not detectable by standard diagnostic techniques (sequencing of coding sequence, copy number variant detection). This might be the reason that some clinically and histologically proven dystrophinopathy cases remain unsolved. We used whole genome sequencing (WGS) to screen the entire DMD gene for variants in one of two brothers suffering from typical muscular dystrophy with strongly elevated creatine kinase levels. Although a pathogenic DIV could not be detected, we were able to identify a pericentric inversion with breakpoints in DMD intron 44 and Xq13.3, which could be confirmed by Sanger sequencing in the index as well as in his brother and mother. As this variation affects a major part of DMD it is most likely disease causing. Our findings elucidate that WGS is capable of detecting large structural rearrangements and might be suitable for the genetic diagnostics of dystrophinopathies in the future. In particular, inversions might be a more frequent cause for dystrophinopathies as anticipated and should be considered in genetically unsolved dystrophinopathy cases.
Project description:Next generation sequencing methodologies are facilitating the rapid characterisation of novel structural variants at nucleotide resolution. These approaches are particularly applicable to variants initially identified using alternative molecular methods. We report a child born with bilateral postaxial syndactyly of the feet and bilateral fifth finger clinodactyly. This was presumed to be an autosomal recessive syndrome, due to the family history of consanguinity. Karyotype analysis revealed a homozygous pericentric inversion of chromosome 7 (46,XX,inv(7)(p15q21)x2) which was confirmed to be heterozygous in both unaffected parents. Since the resolution of the karyotype was insufficient to identify any putatively causative gene, we undertook medium-coverage whole genome sequencing using paired-end reads, in order to elucidate the molecular breakpoints. In a two-step analysis, we first narrowed down the region by identifying discordant read-pairs, and then determined the precise molecular breakpoint by analysing the mapping locations of "soft-clipped" breakpoint-spanning reads. PCR and Sanger sequencing confirmed the identified breakpoints, both of which were located in intergenic regions. Significantly, the 7p15 breakpoint was located 523 kb upstream of HOXA13, the locus for hand-foot-genital syndrome. By inference from studies of HOXA locus control in the mouse, we suggest that the inversion has delocalised a HOXA13 enhancer to produce the phenotype observed in our patient. This study demonstrates how modern genetic diagnostic approach can characterise structural variants at nucleotide resolution and provide potential insights into functional regulation.
Project description:ObjectiveAbsence of homozygosity (AOH) is a genetic characteristic known to cause human diseases mainly through autosomal recessive or imprinting mechanisms. The importance and necessity of accurate AOH detection has become more clinically significant in recent years. However, it remains a challenging task for sequencing-based methods thus far.MethodsIn this study, we developed and optimized a new bioinformatic algorithm based on the assessment of minimum sequencing coverage, optimal bin size, the Z-score threshold of four types of allele count and the frequency for accurate genotyping using 28 AOH negative samples, and redefined the AOH detection cutoff value. We showed the performance of chromosome analysis by five-fold coverage whole genome sequencing (CMA-seq) for AOH identification in 27 typical prenatal/postnatal AOH positive samples, which were previously confirmed by chromosomal microarray analysis with single nucleotide polymorphism array (CMA/SNP array).ResultsThe blinded study indicated that for all three forms of AOH, including whole genomic AOH, single chromosomal AOH and segmental AOH, and all kinds of sample types, including chorionic villus sampling, amniotic fluid, cord blood, peripheral blood and abortive tissue, CMA-seq showed equivalent detection power to that of routine CMA/SNP arrays (750K). The subtle difference between the two methods is that CMA-seq is prone to detect small inconsecutive AOHs, while CMA/SNP array reports it as a whole.ConclusionBased on our newly developed bioinformatic algorithm, it is feasible to detect clinically significant AOH using CMA-seq in prenatal diagnosis.
Project description:A comparison of the human genome with that of the chimpanzee is an attractive approach to attempts to understand the specificity of a certain phenotype's development. The two karyotypes differ by one chromosome fusion, nine pericentric inversions, and various additions of heterochromatin to chromosomal telomeres. Only the fusion, which gave rise to human chromosome 2, has been characterized at the sequence level. During the present study, we investigated the pericentric inversion by which chimpanzee chromosome 19 differs from human chromosome 17. Fluorescence in situ hybridization was used to identify breakpoint-spanning bacterial artificial chromosomes (BACs) and plasmid artificial chromosomes (PACs). By sequencing the junction fragments, we localized breakpoints in intergenic regions rich in repetitive elements. Our findings suggest that repeat-mediated nonhomologous recombination has facilitated inversion formation. No addition or deletion of any sequence element was detected at the breakpoints or in the surrounding sequences. Next to the break, at a distance of 10.2-39.1 kb, the following genes were found: NGFR and NXPH3 (on human chromosome 17q21.3) and GUC2D and ALOX15B (on human chromosome 17p13). The inversion affects neither the genomic structure nor the gene-activity state with regard to replication timing of these genes.
Project description:Apparently, balanced chromosomal rearrangements usually have no phenotypic consequences for the carrier. However, in some cases, they may be associated with an abnormal phenotype. We report herein the case of a 4-year-old boy presenting with clinically isolated supravalvular aortic stenosis (SVAS). No chromosomal imbalance was detected by array CGH. The karyotype showed a balanced paracentric chromosome 7 inversion. Breakpoint characterization using paired-end whole-genome sequencing (WGS) revealed an ELN gene disruption in intron 1, accounting for the phenotype. Family study showed that the inversion was inherited, with incomplete penetrance. To our knowledge, this is the first case of a disruption of the ELN gene characterized by WGS. It contributes to refine the genotype-phenotype correlation in ELN disruption. Although this disruption is a rare etiology of SVAS, it cannot be detected by the diagnostic tests usually performed, such as array CGH or sequencing methods (Sanger, panel, or exome sequencing). With the future perspective of WGS as a diagnostic tool, it will be important to include a structural variation analysis in order to detect balanced rearrangements and gene disruption.
Project description:The cost of whole-genome bisulfite sequencing (WGBS) remains a bottleneck for many studies and it is therefore imperative to extract as much information as possible from a given dataset. This is particularly important because even at the recommend 30X coverage for reference methylomes, up to 50% of high-resolution features such as differentially methylated positions (DMPs) cannot be called with current methods as determined by saturation analysis. To address this limitation, we have developed a tool that dynamically segments WGBS methylomes into blocks of comethylation (COMETs) from which lost information can be recovered in the form of differentially methylated COMETs (DMCs). Using this tool, we demonstrate recovery of ?30% of the lost DMP information content as DMCs even at very low (5X) coverage. This constitutes twice the amount that can be recovered using an existing method based on differentially methylated regions (DMRs). In addition, we explored the relationship between COMETs and haplotypes in lymphoblastoid cell lines of African and European origin. Using best fit analysis, we show COMETs to be correlated in a population-specific manner, suggesting that this type of dynamic segmentation may be useful for integrated (epi)genome-wide association studies in the future.
Project description:The BAM and CRAM formats provide a supplementary linear index that facilitates rapid access to sequence alignments in arbitrary genomic regions. Comparing consecutive entries in a BAM or CRAM index allows one to infer the number of alignment records per genomic region for use as an effective proxy of sequence depth in each genomic region. Based on these properties, we have developed indexcov, an efficient estimator of whole-genome sequencing coverage to rapidly identify samples with aberrant coverage profiles, reveal large-scale chromosomal anomalies, recognize potential batch effects, and infer the sex of a sample. Indexcov is available at https://github.com/brentp/goleft under the MIT license.
Project description:Maternal uniparental disomy of chromosome 7 is present in 5-10% of patients with Silver-Russell syndrome (SRS), and duplication of 7p including GRB10 (Growth Factor Receptor-Bound Protein 10), an imprinted gene that affects pre-and postnatal growth retardation, has been associated with the SRS phenotype. Here, we report on a 17 year old girl referred to array-CGH analysis for short stature, psychomotor delay, and relative macrocephaly. Array-CGH analysis showed two copy number variants (CNVs): a ~12.7 Mb gain in 7p13-p11.2, involving GRB10 and an ~9 Mb loss in 7q11.21-q11.23. FISH experiments performed on the proband's mother showed a chromosome 7 pericentric inversion that might have mediated the complex rearrangement harbored by the daughter. Indeed, we found that segmental duplications, of which chromosome 7 is highly enriched, mapped at the breakpoints of both the mother's inversion and the daughter's CNVs. We postulate that pairing of highly homologous sequences might have perturbed the correct meiotic chromosome segregation, leading to unbalanced outcomes and acting as the putative meiotic mechanism that was causative of the proband's rearrangement. Comparison of the girl's phenotype to those of patients with similar CNVs supports the presence of 7p in a locus associated with features of SRS syndrome.
Project description:Spiny lizards in the genus Sceloporus are a model system among squamate reptiles for studies of chromosomal evolution. While most pleurodont iguanians retain an ancestral karyotype formula of 2n = 36 chromosomes, Sceloporus exhibits substantial karyotype variation ranging from 2n = 22 to 46 chromosomes. We present two annotated chromosome-scale genome assemblies for the Plateau Fence Lizard (Sceloporus tristichus) to facilitate research on the role of pericentric inversion polymorphisms on adaptation and speciation. Based on previous karyotype work using conventional staining, the S. tristichus genome is characterized as 2n = 22 with six pairs of macrochromosomes and five pairs of microchromosomes and a pericentric inversion polymorphism on chromosome 7 that is geographically variable. We provide annotated, chromosome-scale genomes for two lizards located at opposite ends of a dynamic hybrid zone that are each fixed for different inversion polymorphisms. The assembled genomes are 1.84-1.87 Gb (1.72 Gb for scaffolds mapping to chromosomes) with a scaffold N50 of 267.5 Mb. Functional annotation of the genomes resulted in ∼15K predicted gene models. Our assemblies confirmed the presence of a 4.62-Mb pericentric inversion on chromosome 7, which contains 62 annotated coding genes with known functions. In addition, we collected population genomics data using double digest RAD-sequencing for 44 S. tristichus to estimate population structure and phylogeny across the Colorado Plateau. These new genomic resources provide opportunities to perform genomic scans and investigate the formation and spread of pericentric inversions in a naturally occurring hybrid zone.
Project description:Whole-genome bisulfite sequencing (WGBS) allows genome-wide DNA methylation profiling, but the associated high sequencing costs continue to limit its widespread application. We used several high-coverage reference data sets to experimentally determine minimal sequencing requirements. We present data-derived recommendations for minimum sequencing depth for WGBS libraries, highlight what is gained with increasing coverage and discuss the trade-off between sequencing depth and number of assayed replicates.