Microdeletions and microinsertions causing human genetic disease: common mechanisms of mutagenesis and the role of local DNA sequence complexity

Assisting the analysis of insertions and deletions using regional allele frequencies

Accurate estimation of population allele frequency (AF) is crucial for gene discovery and genetic diagnostics. However, determining AF for frameshift-inducing small insertions and deletions (indels) faces challenges due to discrepancies in mapping and variant calling methods. Here, we propose an innovative approach to assess indel AF. We developed CRAFTS-indels (Calculating Regional Allele Frequency Targeting Small indels), an algorithm that combines AF of distinct indels within a given region and provides "regional AF" (rAF). We tested and validated CRAFTS-indels using three independent datasets: gnomAD v2 (n=125,748 samples), an internal dataset (IGM; n=39,367), and the UK BioBank (UKBB; n=469,835). By comparing rAF against standard AF, we identified rare indels with rAF exceeding standard AF (sAF≤10-4 and rAF>10-4) as "rAF-hi" indels. Notably, a high percentage of rare indels were "rAF-hi", with a higher proportion in gnomAD v2 (11-20%) and IGM (11-22%) compared to the UKBB (5-9% depending on the CRAFTS-indels' parameters). Analysis of the overlap of regions based on their rAF with low complexity regions and with ClinVar classification supported the pertinence of rAF. Using the internal dataset, we illustrated the utility of CRAFTS-indel in the analysis of de novo variants and the potential negative impact of rAF-hi indels in gene discovery. In summary, annotation of indels with cohort specific rAF can be used to handle some of the limitations of current annotation pipelines and facilitate detection of novel gene disease associations. CRAFTS-indels offers a user-friendly approach to providing rAF annotation. It can be integrated into public databases such as gnomAD, UKBB and used by ClinVar to revise indel classifications.

Characterization of novel mutations in the TEL-patch domain of the telomeric factor TPP1 associated with telomere biology disorders

Telomeres are nucleoprotein structures that protect the chromosome ends from degradation and fusion. Telomerase is a ribonucleoprotein complex essential to maintain the length of telomeres. Germline defects that lead to short and/or dysfunctional telomeres cause telomere biology disorders (TBDs), a group of rare and heterogeneous Mendelian diseases including pulmonary fibrosis, dyskeratosis congenita, and Høyeraal-Hreidarsson syndrome. TPP1, a telomeric factor encoded by the gene ACD, recruits telomerase at telomere and stimulates its activity via its TEL-patch domain that directly interacts with TERT, the catalytic subunit of telomerase. TBDs due to TPP1 deficiency have been reported only in 11 individuals. We here report four unrelated individuals with a wide spectrum of TBD manifestations carrying either heterozygous or homozygous ACD variants consisting in the recurrent and previously described in-frame deletion of K170 (K170∆) and three novel missense mutations G179D, L184R, and E215V. Structural and functional analyses demonstrated that the four variants affect the TEL-patch domain of TPP1 and impair telomerase activity. In addition, we identified in the ACD gene several motifs associated with small deletion hotspots that could explain the recurrence of the K170∆ mutation. Finally, we detected in a subset of blood cells from one patient, a somatic TERT promoter-activating mutation that likely provides a selective advantage over non-modified cells, a phenomenon known as indirect somatic genetic rescue. Together, our results broaden the genetic and clinical spectrum of TPP1 deficiency and specify new residues in the TEL-patch domain that are crucial for length maintenance and stability of human telomeres in vivo.

Identification of New Genetic Determinants in Pediatric Patients with Familial Hypercholesterolemia Using a Custom NGS Panel

The most common form of inherited lipid disorders is familial hypercholesterolemia (FH). It is characterized primarily by high concentrations of the clinical triad of low-density lipoprotein cholesterol, tendon xanthomas and premature CVD. The well-known genetic background are mutations in LDLR, APOB and PCSK9 gene. Causative mutations can be found in 60−80% of definite FH patients and 20−30% of those with possible FH. Their occurrence could be attributed to the activity of minor candidate genes, whose causal mechanism has not been fully discovered. The aim of the conducted study was to identify disease-causing mutations in FH-related and candidate genes in pediatric patients from Poland using next generation sequencing (NGS). An NGS custom panel was designed to cover 21 causative and candidate genes linked to primary dyslipidemia. Recruitment was performed using Simon Broome diagnostic criteria. Targeted next generation sequencing was performed on a MiniSeq sequencer (Illumina, San Diego, CA, USA) using a 2 × 150 bp paired-end read module. Sequencing data analysis revealed pathogenic and possibly pathogenic variants in 33 out of 57 studied children. The affected genes were LDLR, APOB, ABCG5 and LPL. A novel pathogenic 7bp frameshift deletion c.373_379delCAGTTCG in the exon 4 of the LDLR gene was found. Our findings are the first to identify the c.373_379delCAGTTCG mutation in the LDLR gene. Furthermore, the double heterozygous carrier of frameshift insertion c.2416dupG in the LDLR gene and missense variant c.10708C>T in the APOB gene was identified. The c.2416dupG variant was defined as pathogenic, as confirmed by its cosegregation with hypercholesterolemia in the proband’s family. Although the APOB c.10708C>T variant was previously detected in hypercholesterolemic patients, our data seem to demonstrate no clinical impact. Two missense variants in the LPL gene associated with elevated triglyceride plasma level (c.106G>A and c.953A>G) were also identified. The custom NGS panel proved to be an effective research tool for identifying new causative aberrations in a genetically heterogeneous disease as familial hypercholesterolemia (FH). Our findings expand the spectrum of variants associated with the FH loci and will be of value in genetic counseling among patients with the disease.

Distinct sequence features underlie microdeletions and gross deletions in the human genome

Microdeletions and gross deletions are important causes (~20%) of human inherited disease and their genomic locations are strongly influenced by the local DNA sequence environment. This notwithstanding, no study has systematically examined their underlying generative mechanisms. Here, we obtained 42,098 pathogenic microdeletions and gross deletions from the Human Gene Mutation Database (HGMD) that together form a continuum of germline deletions ranging in size from 1 to 28,394,429 bp. We analyzed the DNA sequence within 1 kb of the breakpoint junctions and found that the frequencies of non‐B DNA‐forming repeats, GC‐content, and the presence of seven of 78 specific sequence motifs in the vicinity of pathogenic deletions correlated with deletion length for deletions of length ≤30 bp. Further, we found that the presence of DR, GQ, and STR repeats is important for the formation of longer deletions (>30 bp) but not for the formation of shorter deletions (≤30 bp) while significantly (χ², p < 2E−16) more microhomologies were identified flanking short deletions than long deletions (length >30 bp). We provide evidence to support a functional distinction between microdeletions and gross deletions. Finally, we propose that a deletion length cut‐off of 25–30 bp may serve as an objective means to functionally distinguish microdeletions from gross deletions.

FLNC-Associated Myofibrillar Myopathy: New Clinical, Functional, and Proteomic Data

Objective

To determine whether a new indel mutation in the dimerization domain of filamin C (FLNc) causes a hereditary myopathy with protein aggregation in muscle fibers, we clinically and molecularly studied a German family with autosomal dominant myofibrillar myopathy (MFM).

Methods

We performed mutational analysis in 3 generations, muscle histopathology, and proteomic studies of IM protein aggregates. Functional consequences of the FLNC mutation were investigated with interaction and transfection studies and biophysics molecular analysis.

Results

Eight patients revealed clinical features of slowly progressive proximal weakness associated with a heterozygous c.8025_8030delCAAGACinsA (p.K2676Pfs*3) mutation in FLNC. Two patients exhibited a mild cardiomyopathy. MRI of skeletal muscle revealed lipomatous changes typical for MFM with FLNC mutations. Muscle biopsies showed characteristic MFM findings with protein aggregation and lesion formation. The proteomic profile of aggregates was specific for MFM-filaminopathy and indicated activation of the ubiquitin-proteasome system (UPS) and autophagic pathways. Functional studies revealed that mutant FLNc is misfolded, unstable, and incapable of forming homodimers and heterodimers with wild-type FLNc.

Conclusions

This new MFM-filaminopathy family confirms that expression of mutant FLNC leads to an adult-onset muscle phenotype with intracellular protein accumulation. Mutant FLNc protein is biochemically compromised and leads to dysregulation of protein quality control mechanisms. Proteomic analysis of MFM protein aggregates is a potent method to identify disease-relevant proteins, differentiate MFM subtypes, evaluate the relevance of gene variants, and identify novel MFM candidate genes.

Spectrum of Disease Severity in Patients With X-Linked Retinitis Pigmentosa Due to RPGR Mutations

Purpose

The purpose of this study was to perform a detailed longitudinal phenotyping of X-linked retinitis pigmentosa (RP) caused by mutations in the RPGR gene during a long follow-up period.

Methods

An Italian cohort of 48 male patients (from 31 unrelated families) with RPGR-associated RP was clinically assessed at a single center (mean follow-up = 6.5 years), including measurements of best-corrected visual acuity (BCVA), Goldmann visual field (GVF), optical coherence tomography (OCT), fundus autofluorescence (FAF), microperimetry, and full-field electroretinography (ERG).

Results

Patients (29.6 ± 15.2 years) showed a mean BCVA of 0.6 ± 0.7 logMAR, mostly with myopic refraction (79.2%). Thirty patients (62.5%) presented a typical RP fundus, while the remaining sine pigmento RP. Over the follow-up, BCVA significantly declined at a mean rate of 0.025 logMAR/year. Typical RP and high myopia were associated with a significantly faster decline of BCVA. Blindness was driven primarily by GVF loss. ERG responses with a rod-cone pattern of dysfunction were detectable in patients (50%) that were significantly younger and more frequently presented sine pigmento RP. Thirteen patients (27.1%) had macular abnormalities without cystoid macular edema. Patients (50%) with a perimacular hyper-FAF ring were significantly younger, had a higher BCVA and a better-preserved ellipsoid zone band than those with markedly decreased FAF. Patients harboring pathogenic variants in exons 1 to 14 showed a milder phenotype compared to those with ORF15 mutations.

Conclusions

Our monocentric, longitudinal retrospective study revealed a spectrum disease progression in male patients with RPGR-associated RP. Slow disease progression correlated with sine pigmento RP, absence of high myopia, and mutations in RPGR exons 1 to 14.

On the complexity measures of mutation hotspots in human TP53 protein

The role of sequence complexity in 23 051 somatic missense mutations including 73 well-known mutation hotspots across 22 major cancers was studied in human TP53 proteins. A role for sequence complexity in TP53 protein mutations is suggested since (i) the mutation rate significantly increases in low amino acid pair bias complexity; (ii) probability distribution complexity increases following single point substitution mutations and strikingly increases after mutation at the mutation hotspots including six detectable hotspot mutations (R175, G245, R248, R249, R273, and R282); and (iii) the degree of increase in distribution complexity is significantly correlated with the frequency of missense mutations (r = -0.5758, P < 0.0001) across 20 major types of solid tumors. These results are consistent with the hypothesis that amino acid pair bias and distribution probability may be used as novel measures for protein sequence complexity, and the degree of complexity is related to its susceptibility to mutation, as such, it may be used as a predictor for modeling protein mutations in human cancers.

First comprehensive TSC1/TSC2 mutational analysis in Mexican patients with Tuberous Sclerosis Complex reveals numerous novel pathogenic variants

The aim of this study was to improve knowledge of the mutational spectrum causing tuberous sclerosis complex (TSC) in a sample of Mexican patients, given the limited information available regarding this disease in Mexico and Latin America. Four different molecular techniques were implemented to identify from single nucleotide variants to large rearrangements in the TSC1 and TSC2 genes of 66 unrelated Mexican-descent patients that clinically fulfilled the criteria for a definitive TSC diagnosis. The mutation detection rate was 94%, TSC2 pathogenic variants (PV) prevailed over TSC1 PV (77% vs. 23%) and a recurrent mutation site (hotspot) was observed in TSC1 exon 15. Interestingly, 40% of the identified mutations had not been previously reported. The wide range of novels PV made it difficult to establish any genotype-phenotype correlation, but most of the PV conditioned neurological involvement (intellectual disability and epilepsy). Our 3D protein modeling of two variants classified as likely pathogenic demonstrated that they could alter the structure and function of the hamartin (TSC1) or tuberin (TSC2) proteins. Molecular analyses of parents and first-degree affected family members of the index cases enabled us to distinguish familial (18%) from sporadic (82%) cases and to identify one case of apparent gonadal mosaicism.

A recurrent 8 bp frameshifting indel in FOXF1 defines a novel mutation hotspot associated with alveolar capillary dysplasia with misalignment of pulmonary veins

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare lethal lung developmental disease. Affected infants manifest with severe respiratory distress and refractory pulmonary hypertension and uniformly die in the first month of life. Heterozygous point mutations or copy-number variant deletions involving FOXF1 and/or its upstream lung-specific enhancer on 16q24.1 have been identified in the vast majority of ACDMPV patients. We have previously described two unrelated families with a de novo pathogenic frameshift variant c.691_698del (p.Ala231Argfs*61) in the exon 1 of FOXF1. Here, we present a third unrelated ACDMPV family with the same de novo variant and propose that a direct tandem repeat of eight consecutive nucleotides GCGGCGGC within the ~4 kb CpG island in FOXF1 exon 1 is a novel mutation hotspot causative for ACDMPV.

dbCID: a manually curated resource for exploring the driver indels in human cancer

While recent advances in next-generation sequencing technologies have enabled the creation of a multitude of databases in cancer genomic research, there is no comprehensive database focusing on the annotation of driver indels (insertions and deletions) yet. Therefore, we have developed the database of Cancer driver InDels (dbCID), which is a collection of known coding indels that likely to be engaged in cancer development, progression or therapy. dbCID contains experimentally supported and putative driver indels derived from manual curation of literature and is freely available online at http://bioinfo.ahu.edu.cn:8080/dbCID. Using the data deposited in dbCID, we summarized features of driver indels in four levels (gene, DNA, transcript and protein) through comparing with putative neutral indels. We found that most of the genes containing driver indels in dbCID are known cancer genes playing a role in tumorigenesis. Contrary to the expectation, the sequences affected by driver frameshift indels are not larger than those by neutral ones. In addition, the frameshift and inframe driver indels prefer to disrupt high-conservative regions both in DNA sequences and protein domains. Finally, we developed a computational method for discriminating cancer driver from neutral frameshift indels based on the deposited data in dbCID. The proposed method outperformed other widely used non-cancer-specific predictors on an external test set, which demonstrated the usefulness of the data deposited in dbCID. We hope dbCID will be a benchmark for improving and evaluating prediction algorithms, and the characteristics summarized here may assist with investigating the mechanism of indel–cancer association.

Genome-wide characterization of genetic variants and putative regions under selection in meat and egg-type chicken lines

Background

Meat and egg-type chickens have been selected for several generations for different traits. Artificial and natural selection for different phenotypes can change frequency of genetic variants, leaving particular genomic footprints throghtout the genome. Thus, the aims of this study were to sequence 28 chickens from two Brazilian lines (meat and white egg-type) and use this information to characterize genome-wide genetic variations, identify putative regions under selection using Fst method, and find putative pathways under selection.

Results

A total of 13.93 million SNPs and 1.36 million INDELs were identified, with more variants detected from the broiler (meat-type) line. Although most were located in non-coding regions, we identified 7255 intolerant non-synonymous SNPs, 512 stopgain/loss SNPs, 1381 frameshift and 1094 non-frameshift INDELs that may alter protein functions. Genes harboring intolerant non-synonymous SNPs affected metabolic pathways related mainly to reproduction and endocrine systems in the white-egg layer line, and lipid metabolism and metabolic diseases in the broiler line. Fst analysis in sliding windows, using SNPs and INDELs separately, identified over 300 putative regions of selection overlapping with more than 250 genes. For the first time in chicken, INDEL variants were considered for selection signature analysis, showing high level of correlation in results between SNP and INDEL data. The putative regions of selection signatures revealed interesting candidate genes and pathways related to important phenotypic traits in chicken, such as lipid metabolism, growth, reproduction, and cardiac development.

Conclusions

In this study, Fst method was applied to identify high confidence putative regions under selection, providing novel insights into selection footprints that can help elucidate the functional mechanisms underlying different phenotypic traits relevant to meat and egg-type chicken lines. In addition, we generated a large catalog of line-specific and common genetic variants from a Brazilian broiler and a white egg layer line that can be used for genomic studies involving association analysis with phenotypes of economic interest to the poultry industry.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4444-0) contains supplementary material, which is available to authorized users.

Mutation analysis of the CTNS gene in Iranian patients with infantile nephropathic cystinosis: identification of two novel mutations

Nephropathic cystinosis is an inherited lysosomal transport disorder caused by mutations in the CTNS gene that encodes for a lysosomal membrane transporter, cystinosin. Dysfunction in this protein leads to cystine accumulation in the cells of different organs. The accumulation of cystine in the kidneys becomes apparent with renal tubular Fanconi syndrome between 6 and 12 months of age and leads to renal failure in the first decade of life. The aim of this study was to analyze the CTNS mutations in 20 Iranian patients, from 20 unrelated families, all of whom were afflicted with infantile nephropathic cystinosis. In these patients, seven different mutant alleles were found, including two new mutations, c.517T>C; p.Y173H and c.492_515del, that have not been previously reported. In addition, we observed that c.681G>A, the common Middle Eastern mutation, was the most common mutation in our patients. Moreover, a new minisatellite or variable number of tandem repeat marker (KX499495) was identified at the CTNS gene. Seven different alleles were found for this marker, and its allele frequency and heterozygosity degree were calculated in cystinosis patients and healthy individuals.

The first Japanese patient with mandibular hypoplasia, deafness, progeroid features and lipodystrophy diagnosed via POLD1 mutation detection

Mandibular hypoplasia, deafness, progeroid features and lipodystrophy (MDPL) syndrome is a rare autosomal dominant disorder caused by heterozygous POLD1 mutations. To date, 13 patients affected by POLD1 mutation-caused MDPL have been described. We report a clinically undiagnosed 11-year-old male who noted joint contractures at 6 years of age. Targeted exome sequencing identified a known POLD1 mutation [NM_002691.3:c.1812_1814del, p.(Ser605del)] that diagnosed him as the first Japanese/East Asian MDPL case.

Improved Diagnosis of Inherited Retinal Dystrophies by High-Fidelity PCR of ORF15 followed by Next-Generation Sequencing

Retinitis pigmentosa (RP) is the most common form of retinal dystrophy. The disease is characterized by the progressive degeneration of photoreceptors, ultimately leading to blindness. The exon ORF15 of RP GTPase regulator (RPGR) is a mutation hot spot for X-linked RP and one form of cone dystrophy. However, accurate molecular testing of ORF15 is challenging because of a large segment of highly repetitive purine-rich sequence in this exon. ORF15 performs poorly in next-generation sequencing-based panels or whole exome sequencing analysis, whereas Sanger sequencing of ORF15 requires special reagents and PCR conditions with multiple pairs of overlapping primers that often do not provide a clean sequence. Because of these technical difficulties, molecular analysis of ORF15 is performed mostly in research laboratories without validation for clinical application. Herein, we report the development of a single step of high-fidelity PCR followed by next-generation sequencing for accurate mutation detection, which is easily integrated into routine clinical practice. Our approach has improved coverage depth of ORF15 with the ability to detect single-nucleotide variants and deletions/duplications. Using this method, we were able to identify ORF15 pathogenic variants in approximately 31% of undiagnosed RP patients. Our results underline the clinical importance of complete and accurate sequence analysis of ORF15 for patients with retinal dystrophies.

SPAI: an interactive platform for indel analysis

BACKGROUND

Insertions and Deletions (Indels) are the most common form of structural variation in human genome. Indels not only contribute to genetic diversity but also cause diseases. Therefore assessing indels in human genome has become an interesting topic to the research community. This increasing interest on indel calling research has resulted into the development of a good number of indel calling tools. However, all of these tools are command line based and require expertise from Computer Science (CS) to execute them which makes it challenging for researchers from non-CS background.

METHODS

In this paper, we describe an interactive platform named SPAI which stands for Single Platform for Analyzing Indels.

RESULTS

Being a Graphical User Interface (GUI) tool, SPAI facilitates users to run several popular indel calling tools and perform several analyses on the indel calling results without knowing any command line programming.

CONCLUSIONS

SPAI is written in Java and tested in Linux operating system.

A Role for Non-B DNA Forming Sequences in Mediating Microlesions Causing Human Inherited Disease

Missense/nonsense mutations and microdeletions/microinsertions (<21 bp) represent ∼ 76% of all mutations causing human inherited disease, and their occurrence has been associated with sequence motifs (direct, inverted, and mirror repeats; G-quartets) capable of adopting non-B DNA structures. We found that a significant proportion (∼ 21%) of both microdeletions and microinsertions occur within direct repeats, and are explicable by slipped misalignment. A novel mutational mechanism, DNA triplex formation followed by DNA repair, may explain ∼ 5% of microdeletions and microinsertions at mirror repeats. Further, G-quartets, direct, and inverted repeats also appear to play a prominent role in mediating missense mutations, whereas only direct and inverted repeats mediate nonsense mutations. We suggest a mutational mechanism involving slipped strand mispairing, slipped structure formation, and DNA repair, to explain ∼ 15% of missense and ∼ 12% of nonsense mutations yielding perfect direct repeats from imperfect repeats, or the extension of existing direct repeats. Similar proportions of missense and nonsense mutations were explicable by hairpin/loop formation and DNA repair, yielding perfect inverted repeats from imperfect repeats. We also propose a model for single base-pair substitution based on one-electron oxidation reactions at G-quadruplex DNA. Overall, the proposed mechanisms provide support for a role for non-B DNA structures in human gene mutagenesis.

Detection and characterization of small insertion and deletion genetic variants in modern layer chicken genomes

BACKGROUND

Small insertions and deletions (InDels) constitute the second most abundant class of genetic variants and have been found to be associated with many traits and diseases. The present study reports on the detection and characterisation of about 883 K high quality InDels from the whole-genome analysis of several modern layer chicken lines from diverse breeds.

RESULTS

To reduce the error rates seen in InDel detection, this study used the consensus set from two InDel-calling packages: SAMtools and Dindel, as well as stringent post-filtering criteria. By analysing sequence data from 163 chickens from 11 commercial and 5 experimental layer lines, this study detected about 883 K high quality consensus InDels with 93% validation rate and an average density of 0.78 InDels/kb over the genome. Certain chromosomes, viz, GGAZ, 16, 22 and 25 showed very low densities of InDels whereas the highest rate was observed on GGA6. In spite of the higher recombination rates on microchromosomes, the InDel density on these chromosomes was generally lower relative to macrochromosomes possibly due to their higher gene density. About 43-87% of the InDels were found to be fixed within each line. The majority of detected InDels (86%) were 1-5 bases and about 63% were non-repetitive in nature while the rest were tandem repeats of various motif types. Functional annotation identified 613 frameshift, 465 non-frameshift and 10 stop-gain/loss InDels. Apart from the frameshift and stopgain/loss InDels that are expected to affect the translation of protein sequences and their biological activity, 33% of the non-frameshift were predicted as evolutionary intolerant with potential impact on protein functions. Moreover, about 2.5% of the InDels coincided with the most-conserved elements previously mapped on the chicken genome and are likely to define functional elements. InDels potentially affecting protein function were found to be enriched for certain gene-classes e.g. those associated with cell proliferation, chromosome and Golgi organization, spermatogenesis, and muscle contraction.

CONCLUSIONS

The large catalogue of InDels presented in this study along with their associated information such as functional annotation, estimated allele frequency, etc. are expected to serve as a rich resource for application in future research and breeding in the chicken.

Characterization and identification of hidden rare variants in the human genome

Background

By examining the genotype calls generated by the 1000 Genomes Project we discovered that the human reference genome GRCh37 contains almost 20,000 loci in which the reference allele has never been observed in healthy individuals and around 70,000 loci in which it has been observed only in the heterozygous state.

Results

We show that a large fraction of this rare reference allele (RRA) loci belongs to coding, functional and regulatory elements of the genome and could be linked to rare Mendelian disorders as well as cancer. We also demonstrate that classical germline and somatic variant calling tools are not capable to recognize the rare allele when present in these loci. To overcome such limitations, we developed a novel tool, named RAREVATOR, that is able to identify and call the rare allele in these genomic positions. By using a small cancer dataset we compared our tool with two state-of-the-art callers and we found that RAREVATOR identified more than 1,500 germline and 22 somatic RRA variants missed by the two methods and which belong to significantly mutated pathways.

Conclusions

These results show that, to date, the investigation of around 100,000 loci of the human genome has been missed by re-sequencing experiments based on the GRCh37 assembly and that our tool can fill the gap left by other methods. Moreover, the investigation of the latest version of the human reference genome, GRCh38, showed that although the GRC corrected almost all insertions and a small part of SNVs and deletions, a large number of functionally relevant RRAs still remain unchanged. For this reason, also future resequencing experiments, based on GRCh38, will benefit from RAREVATOR analysis results. RAREVATOR is freely available at http://sourceforge.net/projects/rarevator.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1481-9) contains supplementary material, which is available to authorized users.

Hellbender genome sequences shed light on genomic expansion at the base of crown salamanders

Among animals, genome sizes range from 20 Mb to 130 Gb, with 380-fold variation across vertebrates. Most of the largest vertebrate genomes are found in salamanders, an amphibian clade of 660 species. Thus, salamanders are an important system for studying causes and consequences of genomic gigantism. Previously, we showed that plethodontid salamander genomes accumulate higher levels of long terminal repeat (LTR) retrotransposons than do other vertebrates, although the evolutionary origins of such sequences remained unexplored. We also showed that some salamanders in the family Plethodontidae have relatively slow rates of DNA loss through small insertions and deletions. Here, we present new data from Cryptobranchus alleganiensis, the hellbender. Cryptobranchus and Plethodontidae span the basal phylogenetic split within salamanders; thus, analyses incorporating these taxa can shed light on the genome of the ancestral crown salamander lineage, which underwent expansion. We show that high levels of LTR retrotransposons likely characterize all crown salamanders, suggesting that disproportionate expansion of this transposable element (TE) class contributed to genomic expansion. Phylogenetic and age distribution analyses of salamander LTR retrotransposons indicate that salamanders' high TE levels reflect persistence and diversification of ancestral TEs rather than horizontal transfer events. Finally, we show that relatively slow DNA loss rates through small indels likely characterize all crown salamanders, suggesting that a decreased DNA loss rate contributed to genomic expansion at the clade's base. Our identification of shared genomic features across phylogenetically distant salamanders is a first step toward identifying the evolutionary processes underlying accumulation and persistence of high levels of repetitive sequence in salamander genomes.

Are differences in genomic data sets due to true biological variants or errors in genome assembly: an example from two chloroplast genomes

DNA sequencing has been revolutionized by the development of high-throughput sequencing technologies. Plummeting costs and the massive throughput capacities of second and third generation sequencing platforms have transformed many fields of biological research. Concurrently, new data processing pipelines made rapid de novo genome assemblies possible. However, high quality data are critically important for all investigations in the genomic era. We used chloroplast genomes of one Oryza species (O. australiensis) to compare differences in sequence quality: one genome (GU592209) was obtained through Illumina sequencing and reference-guided assembly and the other genome (KJ830774) was obtained via target enrichment libraries and shotgun sequencing. Based on the whole genome alignment, GU592209 was more similar to the reference genome (O. sativa: AY522330) with 99.2% sequence identity (SI value) compared with the 98.8% SI values in the KJ830774 genome; whereas the opposite result was obtained when the SI values in coding and noncoding regions of GU592209 and KJ830774 were compared. Additionally, the junctions of two single copies and repeat copies in the chloroplast genome exhibited differences. Phylogenetic analyses were conducted using these sequences, and the different data sets yielded dissimilar topologies: phylogenetic replacements of the two individuals were remarkably different based on whole genome sequencing or SNP data and insertions and deletions (indels) data. Thus, we concluded that the genomic composition of GU592209 was heterogeneous in coding and non-coding regions. These findings should impel biologists to carefully consider the quality of sequencing and assembly when working with next-generation data.

DDIG-in: detecting disease-causing genetic variations due to frameshifting indels and nonsense mutations employing sequence and structural properties at nucleotide and protein levels

MOTIVATION

Frameshifting (FS) indels and nonsense (NS) variants disrupt the protein-coding sequence downstream of the mutation site by changing the reading frame or introducing a premature termination codon, respectively. Despite such drastic changes to the protein sequence, FS indels and NS variants have been discovered in healthy individuals. How to discriminate disease-causing from neutral FS indels and NS variants is an understudied problem.

RESULTS

We have built a machine learning method called DDIG-in (FS) based on real human genetic variations from the Human Gene Mutation Database (inherited disease-causing) and the 1000 Genomes Project (GP) (putatively neutral). The method incorporates both sequence and predicted structural features and yields a robust performance by 10-fold cross-validation and independent tests on both FS indels and NS variants. We showed that human-derived NS variants and FS indels derived from animal orthologs can be effectively employed for independent testing of our method trained on human-derived FS indels. DDIG-in (FS) achieves a Matthews correlation coefficient (MCC) of 0.59, a sensitivity of 86%, and a specificity of 72% for FS indels. Application of DDIG-in (FS) to NS variants yields essentially the same performance (MCC of 0.43) as a method that was specifically trained for NS variants. DDIG-in (FS) was shown to make a significant improvement over existing techniques.

The tandem repeats enabling reversible switching between the two phases of β-lactamase substrate spectrum

Expansion or shrinkage of existing tandem repeats (TRs) associated with various biological processes has been actively studied in both prokaryotic and eukaryotic genomes, while their origin and biological implications remain mostly unknown. Here we describe various duplications (de novo TRs) that occurred in the coding region of a β-lactamase gene, where a conserved structure called the omega loop is encoded. These duplications that occurred under selection using ceftazidime conferred substrate spectrum extension to include the antibiotic. Under selective pressure with one of the original substrates (amoxicillin), a high level of reversion occurred in the mutant β-lactamase genes completing a cycle back to the original substrate spectrum. The de novo TRs coupled with reversion makes a genetic toggling mechanism enabling reversible switching between the two phases of the substrate spectrum of β-lactamases. This toggle exemplifies the effective adaptation of de novo TRs for enhanced bacterial survival. We found pairs of direct repeats that mediated the DNA duplication (TR formation). In addition, we found different duos of sequences that mediated the DNA duplication. These novel elements—that we named SCSs (same-strand complementary sequences)—were also found associated with β-lactamase TR mutations from clinical isolates. Both direct repeats and SCSs had a high correlation with TRs in diverse bacterial genomes throughout the major phylogenetic lineages, suggesting that they comprise a fundamental mechanism shaping the bacterial evolution.

β-lactamases can adapt to new antibiotics by mutations in their genes. The original and the extended substrate spectrums of β-lactamases define two phases of catalytic activity, and the conversion by point mutations is unidirectional from the initial to the new spectrum. We describe duplication mutations that enable reversible switching between the substrate spectrums, increasing the adaptability of the bacterium. We provide evidence supporting that two distinct groups of short sequences mediated the formation of DNA duplications in β-lactamases: direct repeats and novel elements that we named, SCSs (same-strand complementary sequences). Our study suggests that DNA duplication processes mediated by both direct repeats and SCSs are not just limited to the β-lactamase genes but comprise a fundamental mechanism in bacterial genome evolution.

A novel de novo duplication mutation of PAX6 in a Chinese family with aniridia and other ocular abnormalities

Aniridia is a congenital panocular disorder caused by the mutations of the paired box gene-6 (PAX6). To investigate the clinical characterization and the underlying genetic defect in a Chinese family with aniridia and other ocular abnormalities, we recruited the family members who underwent ophthalmic examination. Two patients in this family, the proband and his affected son, both have bilateral aniridia, foveal hypoplasia and nystagmus. Moreover, the proband also had presenile cataracts, but his affected son did not show cataracts at the time of examination. Sequencing PAX6 revealed that a heterozygous duplication mutation c.95_105dup11, predicted to generate non-functional truncated protein at position Gly36 (p.G36X), was found in the affected individuals but not in any of the unaffected family members including the parents of the proband. Haplotype analysis showed that the proband and his affected son shared a common disease-related haplotype, which was arisen from the proband's unaffected father through crossing-over. In conclusion, we identified a novel de novo duplication mutation of PAX6 in the aniridia and other ocular abnormalities family. This mutation has occurred de novo on a paternal chromosome by direct duplication, which presumably results from replication slippage or unequal non-sister chromatids exchange during spermatogenesis.

Molecular genetic epidemiology of human diseases: from patterns to predictions

Databases of disease-associated or disease-causing mutations allow the study, not only of the molecular mechanisms underlying the primary lesions at the DNA level, but also of the functional consequences of mutation at the phenotypic level. The Human Gene Mutation Database (HGMD) and the bioinformatics analyses of its content provide an illustrative example of this indirect approach to molecular genetic epidemiology. In fact, the Bayesian type of reasoning underlying previous scientific analyses of HGMD data is also reflected in current software tools used to predict the likely disease relevance of a newly detected genetic variant. After a brief resume of the past scientific utility of HGMD, we, therefore, shortly review three representative and commonly used examples of these tools, namely SIFT, PolyPhen-2 and NNSplice.

Wavelet analysis in current cancer genome research: a survey

With the rapid development of next generation sequencing technology, the amount of biological sequence data of the cancer genome increases exponentially, which calls for efficient and effective algorithms that may identify patterns hidden underneath the raw data that may distinguish cancer Achilles' heels. From a signal processing point of view, biological units of information, including DNA and protein sequences, have been viewed as one-dimensional signals. Therefore, researchers have been applying signal processing techniques to mine the potentially significant patterns within these sequences. More specifically, in recent years, wavelet transforms have become an important mathematical analysis tool, with a wide and ever increasing range of applications. The versatility of wavelet analytic techniques has forged new interdisciplinary bounds by offering common solutions to apparently diverse problems and providing a new unifying perspective on problems of cancer genome research. In this paper, we provide a survey of how wavelet analysis has been applied to cancer bioinformatics questions. Specifically, we discuss several approaches of representing the biological sequence data numerically and methods of using wavelet analysis on the numerical sequences.

Genetic analysis of West Nile virus isolates from an outbreak in Idaho, United States, 2006-2007

West Nile virus (WNV) appeared in the U.S. in 1999 and has since become endemic, with yearly summer epidemics causing tens of thousands of cases of serious disease over the past 14 years. Analysis of WNV strains isolated during the 2006–2007 epidemic seasons demonstrates that a new genetic variant had emerged coincidentally with an intense outbreak in Idaho during 2006. The isolates belonging to the new variant carry a 13 nt deletion, termed ID-Δ13, located at the variable region of the 3′UTR, and are genetically related. The analysis of deletions and insertions in the 3′UTR of two major lineages of WNV revealed the presence of conserved repeats and two indel motifs in the variable region of the 3′UTR. One human and two bird isolates from the Idaho 2006–2007 outbreaks were sequenced using Illumina technology and within-host variability was analyzed. Continued monitoring of new genetic variants is important for public health as WNV continues to evolve.

Ionizing radiation and genetic risks. XVII. Formation mechanisms underlying naturally occurring DNA deletions in the human genome and their potential relevance for bridging the gap between induced DNA double-strand breaks and deletions in irradiated germ cells

While much is known about radiation-induced DNA double-strand breaks (DSBs) and their repair, the question of how deletions of different sizes arise as a result of the processing of DSBs by the cell's repair systems has not been fully answered. In order to bridge this gap between DSBs and deletions, we critically reviewed published data on mechanisms pertaining to: (a) repair of DNA DSBs (from basic studies in this area); (b) formation of naturally occurring structural variation (SV) - especially of deletions - in the human genome (from genomic studies) and (c) radiation-induced mutations and structural chromosomal aberrations in mammalian somatic cells (from radiation mutagenesis and radiation cytogenetic studies). The specific aim was to assess the relative importance of the postulated mechanisms in generating deletions in the human genome and examine whether empirical data on radiation-induced deletions in mouse germ cells are consistent with predictions of these mechanisms. The mechanisms include (a) NHEJ, a DSB repair process that does not require any homology and which functions in all stages of the cell cycle (and is of particular relevance in G0/G1); (b) MMEJ, also a DSB repair process but which requires microhomology and which presumably functions in all cell cycle stages; (c) NAHR, a recombination-based DSB repair mechanism which operates in prophase I of meiosis in germ cells; (d) MMBIR, a microhomology-mediated, replication-based mechanism which operates in the S phase of the cell cycle, and (e) strand slippage during replication (involved in the origin of small insertions and deletions (INDELs). Our analysis permits the inference that, between them, these five mechanisms can explain nearly all naturally occurring deletions of different sizes identified in the human genome, NAHR and MMBIR being potentially more versatile in this regard. With respect to radiation-induced deletions, the basic studies suggest that those arising as a result of the operation of NHEJ/MMEJ processes, as currently formulated, are expected to be relatively small. However, data on induced mutations in mouse spermatogonial stem cells (irradiation in G0/G1 phase of the cell cycle and DSB repair presumed to be via NHEJ predominantly) show that most are associated with deletions of different sizes, some in the megabase range. There is thus a 'discrepancy' between what the basic studies suggest and the empirical observations in mutagenesis studies. This discrepancy, however, is only an apparent but not a real one. It can be resolved by considering the issue of deletions in the broader context of and in conjunction with the organization of chromatin in chromosomes and nuclear architecture, the conceptual framework for which already exists in studies carried out during the past fifteen years or so. In this paper, we specifically hypothesize that repair of DSBs induced in chromatin loops may offer a basis to explain the induction of deletions of different sizes and suggest an approach to test the hypothesis. We emphasize that the bridging of the gap between induced DSB and resulting deletions of different sizes is critical for current efforts in computational modeling of genetic risks.

Characterization of bud emergence 46 (BEM46) protein: sequence, structural, phylogenetic and subcellular localization analyses

The bud emergence 46 (BEM46) protein from Neurospora crassa belongs to the α/β-hydrolase superfamily. Recently, we have reported that the BEM46 protein is localized in the perinuclear ER and also forms spots close by the plasma membrane. The protein appears to be required for cell type-specific polarity formation in N. crassa. Furthermore, initial studies suggested that the BEM46 amino acid sequence is conserved in eukaryotes and is considered to be one of the widespread conserved "known unknown" eukaryotic genes. This warrants for a comprehensive phylogenetic analysis of this superfamily to unravel origin and molecular evolution of these genes in different eukaryotes. Herein, we observe that all eukaryotes have at least a single copy of a bem46 ortholog. Upon scanning of these proteins in various genomes, we find that there are expansions leading into several paralogs in vertebrates. Usingcomparative genomic analyses, we identified insertion/deletions (indels) in the conserved domain of BEM46 protein, which allow to differentiate fungal classes such as ascomycetes from basidiomycetes. We also find that exonic indels are able to differentiate BEM46 homologs of different eukaryotic lineage. Furthermore, we unravel that BEM46 protein from N. crassa possess a novel endoplasmic-retention signal (PEKK) using GFP-fusion tagging experiments. We propose that three residues namely a serine 188S, a histidine 292H and an aspartic acid 262D are most critical residues, forming a catalytic triad in BEM46 protein from N. crassa. We carried out a comprehensive study on bem46 genes from a molecular evolution perspective with combination of functional analyses. The evolutionary history of BEM46 proteins is characterized by exonic indels in lineage specific manner.

Molecular topography of the MED12-deleted region in smooth muscle tumors: a possible link between non-B DNA structures and hypermutability

Background

Deletions of the gene encoding mediator subcomplex 12 (MED12) in human smooth muscle tumors rank among the most frequent genomic alterations in human tumors at all. In a minority of these cases, small deletions are found. In an attempt to delineate key features of the deletions aimed at a better understanding of the molecular pathogenesis of uterine smooth muscle tumors we have analyzed 70 MED12 deletions including 46 cases from the literature and 24 own unpublished cases.

Results

The average length of the deletions was 18.7 bp ranging between 2 bp and 43 bp. While in general multitudes of 3 clearly dominated leaving the transcript in frame, deletions of 21, 24, 30, and 33 nucleotides were clearly underrepresented. Within the DNA segment affected deletion breakpoints were not randomly distributed. Most breakpoints clustered within the center of the segment where two peaks of breakpoint clusters could be distinguished. Interestingly, one of these clusters coincides with the loop of a putative folded non-B DNA structure whereas a much lower number of breaks noted in the 5′ and 3′ stem of the structure forming an intramolecular B-helix. The second cluster mainly consisting of 3′ breaks was located in a region downstream adjacent to the stem.

Conclusion

The present study describes for the first time main characteristics of MED12 deletions occurring in smooth muscle tumors. Interestingly, the non-random distribution of breakpoints within the deletion hotspot region may point to a role of non-canonical DNA structures for the occurrence of these mutations and the molecular pathogenesis of uterine smooth muscle tumors, respectively.

Slow DNA loss in the gigantic genomes of salamanders

Evolutionary changes in genome size result from the combined effects of mutation, natural selection, and genetic drift. Insertion and deletion mutations (indels) directly impact genome size by adding or removing sequences. Most species lose more DNA through small indels (i.e., ∼1–30 bp) than they gain, which can result in genome reduction over time. Because this rate of DNA loss varies across species, small indel dynamics have been suggested to contribute to genome size evolution. Species with extremely large genomes provide interesting test cases for exploring the link between small indels and genome size; however, most large genomes remain relatively unexplored. Here, we examine rates of DNA loss in the tetrapods with the largest genomes—the salamanders. We used low-coverage genomic shotgun sequence data from four salamander species to examine patterns of insertion, deletion, and substitution in neutrally evolving non-long terminal repeat (LTR) retrotransposon sequences. For comparison, we estimated genome-wide DNA loss rates in non-LTR retrotransposon sequences from five other vertebrate genomes: Anolis carolinensis, Danio rerio, Gallus gallus, Homo sapiens, and Xenopus tropicalis. Our results show that salamanders have significantly lower rates of DNA loss than do other vertebrates. More specifically, salamanders experience lower numbers of deletions relative to insertions, and both deletions and insertions are skewed toward smaller sizes. On the basis of these patterns, we conclude that slow DNA loss contributes to genomic gigantism in salamanders. We also identify candidate molecular mechanisms underlying these differences and suggest that natural variation in indel dynamics provides a unique opportunity to study the basis of genome stability.

DDIG-in: discriminating between disease-associated and neutral non-frameshifting micro-indels

Micro-indels (insertions or deletions shorter than 21 bps) constitute the second most frequent class of human gene mutation after single nucleotide variants. Despite the relative abundance of non-frameshifting indels, their damaging effect on protein structure and function has gone largely unstudied. We have developed a support vector machine-based method named DDIG-in (Detecting disease-causing genetic variations due to indels) to prioritize non-frameshifting indels by comparing disease-associated mutations with putatively neutral mutations from the 1,000 Genomes Project. The final model gives good discrimination for indels and is robust against annotation errors. A webserver implementing DDIG-in is available at http://sparks-lab.org/ddig.

New clinical and molecular insights on Barth syndrome

Background

Barth syndrome (BS) is an X-linked infantile-onset cardioskeletal disease characterized by cardiomyopathy, hypotonia, growth delay, neutropenia and 3-methylglutaconic aciduria. It is caused by mutations in the TAZ gene encoding tafazzin, a protein involved in the metabolism of cardiolipin, a mitochondrial-specific phospholipid involved in mitochondrial energy production.

Methods

Clinical, biochemical and molecular characterization of a group of six male patients suspected of having BS. Three patients presented early with severe metabolic decompensation including respiratory distress, oxygen desaturation and cardiomyopathy and died within the first year of life. The remaining three patients had cardiomyopathy, hypotonia and growth delay and are still alive. Cardiomyopathy was detected during pregnancy through a routine check-up in one patient. All patients exhibited 3-methylglutaconic aciduria and neutropenia, when tested and five of them also had lactic acidosis.

Results

We confirmed the diagnosis of BS with sequence analysis of the TAZ gene, and found five new mutations, c.641A>G p.His214Arg, c.284dupG (p.Thr96Aspfs*37), c.678_691del14 (p.Tyr227Trpfs*79), g.8009_16445del8437 and g.[9777_9814del38; 9911-?_14402del] and the known nonsense mutation c.367C>T (p.Arg123Term). The two gross rearrangements ablated TAZ exons 6 to 11 and probably originated by non-allelic homologous recombination and by Serial Replication Slippage (SRS), respectively. The identification of the breakpoints boundaries of the gross deletions allowed the direct detection of heterozygosity in carrier females.

Conclusions

Lactic acidosis associated with 3-methylglutaconic aciduria is highly suggestive of BS, whilst the severity of the metabolic decompensation at disease onset should be considered for prognostic purposes. Mutation analysis of the TAZ gene is necessary for confirming the clinical and biochemical diagnosis in probands in order to identify heterozygous carriers and supporting prenatal diagnosis and genetic counseling.

Dehydrated hereditary stomatocytosis linked to gain-of-function mutations in mechanically activated PIEZO1 ion channels

Dehydrated hereditary stomatocytosis is a genetic condition with defective red blood cell membrane properties that causes an imbalance in intracellular cation concentrations. Recently, two missense mutations in the mechanically activated PIEZO1 (FAM38A) ion channel were associated with dehydrated hereditary stomatocytosis. However, it is not known how these mutations affect PIEZO1 function. Here, by combining linkage analysis and whole-exome sequencing in a large pedigree and Sanger sequencing in two additional kindreds and 11 unrelated dehydrated hereditary stomatocytosis cases, we identify three novel missense mutations and one recurrent duplication in PIEZO1, demonstrating that it is the major gene for dehydrated hereditary stomatocytosis. All the dehydrated hereditary stomatocytosis-associated mutations locate at C-terminal half of PIEZO1. Remarkably, we find that all PIEZO1 mutations give rise to mechanically activated currents that inactivate more slowly than wild-type currents. This gain-of-function PIEZO1 phenotype provides insight that helps to explain the increased permeability of cations in red blood cells of dehydrated hereditary stomatocytosis patients. Our findings also suggest a new role for mechanotransduction in red blood cell biology and pathophysiology.

Molecular characterization of cystinuria in south-eastern European countries

Cystinuria is an autosomal recessive disorder caused by defective transport of cystine and dibasic amino acids in the proximal renal tubules and small intestine. So far, more than 128 mutations in SLC3A1 gene, and 93 in SLC7A9 gene have been described as a cause of cystinuria. We present a molecular characterization of the cystinuria in 47 unrelated south-east European families. The molecular methodology included direct sequencing, single strand conformational polymorphism, and restriction fragment length polymorphism. A total of 93 (94.9 %) out of 98 unrelated cystinuric chromosomes have been characterized. Mutations in SLC3A1 gene account for 64.3 % and in SLC7A9 gene for 30.6 % of the cystinuric chromosomes. Ten different mutations in SLC3A1 gene were found, and two of them were novel (C242R and L573X), while in SLC7A9 gene seven mutations were found, of which three were novel (G73R, V375I and c.1048_1051delACTC). The most common mutations in this study were T216M (24.5 %), M467T (16.3 %) and R365L (11.2 %) in SLC3A1 and G105R (21.4 %) in SLC7A9 gene. A population specificity of cystinuria mutations was observed; T216M mutation was the only mutation present among Gypsies, G105R was the most common mutation among Albanians and Macedonians, and R365L among Serbs. The results of this study allowed introduction of rapid, simple and cost-effective genetic diagnosis of cystinuria that enables an early preventive care of affected patients and a prenatal diagnosis in affected families.

The Human Gene Mutation Database (HGMD) and its exploitation in the fields of personalized genomics and molecular evolution

The Human Gene Mutation Database (HGMD) constitutes a comprehensive core collection of data on germ-line mutations in nuclear genes underlying or associated with human inherited disease (http://www.hgmd.org). Data cataloged include single-base-pair substitutions in coding, regulatory, and splicing-relevant regions, micro-deletions and micro-insertions, indels, and triplet repeat expansions, as well as gross gene deletions, insertions, duplications, and complex rearrangements. Each mutation is entered into HGMD only once, in order to avoid confusion between recurrent and identical-by-descent lesions. By March 2012, the database contained in excess of 123,600 different lesions (HGMD Professional release 2012.1) detected in 4,514 different nuclear genes, with new entries currently accumulating at a rate in excess of 10,000 per annum. ∼6,000 of these entries constitute disease-associated and functional polymorphisms. HGMD also includes cDNA reference sequences for more than 98% of the listed genes.

Identification of recurrent type-2 NF1 microdeletions reveals a mitotic nonallelic homologous recombination hotspot underlying a human genomic disorder

Nonallelic homologous recombination (NAHR) is one of the major mechanisms underlying copy number variation in the human genome. Although several disease-associated meiotic NAHR breakpoints have been analyzed in great detail, hotspots for mitotic NAHR are not well characterized. Type-2 NF1 microdeletions, which are predominantly of postzygotic origin, constitute a highly informative model with which to investigate the features of mitotic NAHR. Here, a custom-designed MLPA- and PCR-based approach was used to identify 23 novel NAHR-mediated type-2 NF1 deletions. Breakpoint analysis of these 23 type-2 deletions, together with 17 NAHR-mediated type-2 deletions identified previously, revealed that the breakpoints are nonuniformly distributed within the paralogous SUZ12 and SUZ12P sequences. Further, the analysis of this large group of type-2 deletions revealed breakpoint recurrence within short segments (ranging in size from 57 to 253-bp) as well as the existence of a novel NAHR hotspot of 1.9-kb (termed PRS4). This hotspot harbored 20% (8/40) of the type-2 deletion breakpoints and contains the 253-bp recurrent breakpoint region BR6 in which four independent type-2 deletion breakpoints were identified. Our findings indicate that a combination of an open chromatin conformation and short non-B DNA-forming repeats may predispose to recurrent mitotic NAHR events between SUZ12 and its pseudogene.

More single-nucleotide mutations surround small insertions than small deletions in primates

Early studies have shown that single-nucleotide mutation rates increase close to insertions and deletions, but it is not fully understood how natural selection shapes genome-wide patterns of indels and their nearby single-nucleotide mutations. In this study, we find that, in primates, more single-nucleotide mutations surround small insertions than small deletions. This pattern affects <150 base pair (bp) sequences close to indels and persists under different genomic properties, such as exon/intron/intergenic contexts, repeated/nonrepeated sequences, replication timing, recombination rates, indel density, and guanine-cytosine (GC) content. We propose two different, but not mutually exclusive, hypothetical mechanisms to explain the pattern. One mechanism is that the sequence context preferring insertion formation may also favor nucleotide substitutions. Another mechanism is related to a hypothesis in which indel heterozygosity tends to increase nearby nucleotide substitution rates. It means that if insertions spend more time in heterozygotes, insertions may accumulate more surrounding single-nucleotide changes. In conclusion, we characterize a special genome-wide evolutionary pattern for indels and nearby single-nucleotide changes. This pattern may be driven by natural selection and bias primates' genome evolution and phenotypic variations.

Revertant mosaicism in a human skin fragility disorder results from slipped mispairing and mitotic recombination

Spontaneous gene repair, also called revertant mosaicism, has been documented in several genetic disorders involving organs that undergo self-regeneration, including the skin. Genetic reversion may occur through different mechanisms, and in a single individual, the mutation can be repaired in various ways. Here we describe a disseminated pattern of revertant mosaicism observed in 6 patients with Kindler syndrome (KS), a genodermatosis caused by loss of kindlin-1 (encoded by FERMT1) and clinically characterized by patchy skin pigmentation and atrophy. All patients presented duplication mutations (c.456dupA and c.676dupC) in FERMT1, and slipped mispairing in direct nucleotide repeats was identified as the reversion mechanism in all investigated revertant skin spots. The sequence around the mutations demonstrated high propensity to mutations, favoring both microinsertions and microdeletions. Additionally, in some revertant patches, mitotic recombination generated areas with homozygous normal keratinocytes. Restoration of kindlin-1 expression led to clinically and structurally normal skin. Since loss of kindlin-1 severely impairs keratinocyte proliferation, we predict that revertant cells have a selective advantage that allows their clonal expansion and, consequently, the improvement of the skin condition.

Interlocus gene conversion events introduce deleterious mutations into at least 1% of human genes associated with inherited disease

Establishing the molecular basis of DNA mutations that cause inherited disease is of fundamental importance to understanding the origin, nature, and clinical sequelae of genetic disorders in humans. The majority of disease-associated mutations constitute single-base substitutions and short deletions and/or insertions resulting from DNA replication errors and the repair of damaged bases. However, pathological mutations can also be introduced by nonreciprocal recombination events between paralogous sequences, a phenomenon known as interlocus gene conversion (IGC). IGC events have thus far been linked to pathology in more than 20 human genes. However, the large number of duplicated gene sequences in the human genome implies that many more disease-associated mutations could originate via IGC. Here, we have used a genome-wide computational approach to identify disease-associated mutations derived from IGC events. Our approach revealed hundreds of known pathological mutations that could have been caused by IGC. Further, we identified several dozen high-confidence cases of inherited disease mutations resulting from IGC in ∼1% of all genes analyzed. About half of the donor sequences associated with such mutations are functional paralogous genes, suggesting that epistatic interactions or differential expression patterns will determine the impact upon fitness of specific substitutions between duplicated genes. In addition, we identified thousands of hitherto undescribed and potentially deleterious mutations that could arise via IGC. Our findings reveal the extent of the impact of interlocus gene conversion upon the spectrum of human inherited disease.

Characterizing linkage disequilibrium and evaluating imputation power of human genomic insertion-deletion polymorphisms

Background

Indels are an important cause of human variation and central to the study of human disease. The 1000 Genomes Project Low-Coverage Pilot identified over 1.3 million indels shorter than 50 bp, of which over 890 were identified as potentially disruptive variants. Yet, despite their ubiquity, the local genomic characteristics of indels remain unexplored.

Results

Herein we describe population- and minor allele frequency-based differences in linkage disequilibrium and imputation characteristics for indels included in the 1000 Genomes Project Low-Coverage Pilot for the CEU, YRI and CHB+JPT populations. Common indels were well tagged by nearby SNPs in all studied populations, and were also tagged at a similar rate to common SNPs. Both neutral and functionally deleterious common indels were imputed with greater than 95% concordance from HapMap Phase 3 and OMNI SNP sites. Further, 38 to 56% of low frequency indels were tagged by low frequency SNPs. We were able to impute heterozygous low frequency indels with over 50% concordance. Lastly, our analysis also revealed evidence of ascertainment bias. This bias prevents us from extending the applicability of our results to highly polymorphic indels that could not be identified in the Low-Coverage Pilot.

Conclusions

Although further scope exists to improve the imputation of low frequency indels, our study demonstrates that there are already ample opportunities to retrospectively impute indels for prior genome-wide association studies and to incorporate indel imputation into future case/control studies.

Characterization of a large genomic deletion in four Irish families with C7 deficiency

Inherited deficiency of the seventh complement component (C7) is associated with increased susceptibility to Neisseria meningitidis infections. The disease is rare in most Western countries. Here we report new investigations of a large, but incompletely characterized genomic deletion of exons 8 and 9 [c.739-?_1093+?del], previously identified in three unrelated Irish families with C7 deficiency. We have analysed DNA from one individual, who is homozygous for the deletion, by PCR using primers progressively proximal to the deleted exons. Thus we were able to map the deletion boundaries. Amplification across the breakpoint and sequencing revealed an indel mutation that included a 6.4kb deletion together with an insertion of a novel 8bp sequence [c.739+1262_1270-2387delinsGCAGGCCA]. We demonstrated the same defect in the C7 deficient patients from each family and developed a duplex PCR method to enable the detection of alleles containing the deletion in heterozygotes. A member of a fourth family was found to be homozygous for the deletion defect. Thus, the deletion defect may be a more commonly distributed cause of C7 deficiency in Ireland.

Characterization of the nonallelic homologous recombination hotspot PRS3 associated with type-3 NF1 deletions

Nonallelic homologous recombination (NAHR) is the major mechanism underlying recurrent genomic rearrangements, including the large deletions at 17q11.2 that cause neurofibromatosis type 1 (NF1). Here, we identify a novel NAHR hotspot, responsible for type-3 NF1 deletions that span 1.0 Mb. Breakpoint clustering within this 1-kb hotspot, termed PRS3, was noted in 10 of 11 known type-3 NF1 deletions. PRS3 is located within the LRRC37B pseudogene of the NF1-REPb and NF1-REPc low-copy repeats. In contrast to other previously characterized NAHR hotspots, PRS3 has not developed on a preexisting allelic homologous recombination hotspot. Furthermore, the variation pattern of PRS3 and its flanking regions is unusual since only NF1-REPc (and not NF1-REPb) is characterized by a high single nucleotide polymorphism (SNP) frequency, suggestive of unidirectional sequence transfer via nonallelic homologous gene conversion (NAHGC). By contrast, the previously described intense NAHR hotspots within the CMT1A-REPs, and the PRS1 and PRS2 hotspots underlying type-1 NF1 deletions, experience frequent bidirectional sequence transfer. PRS3 within NF1-REPc was also found to be involved in NAHGC with the LRRC37B gene, the progenitor locus of the LRRC37B-P duplicons, as indicated by the presence of shared SNPs between these loci. PRS3 therefore represents a weak (and probably evolutionarily rather young) NAHR hotspot with unique properties.

An exonic insertion within Tex14 gene causes spermatogenic arrest in pigs

BACKGROUND

Male infertility is an increasing problem in all domestic species including man. Localization and identification of genes involved in defects causing male infertility provide valuable information of specific events in sperm development. Sperm development is a complex process, where diploid spermatogonia develop into haploid, highly specialized spermatozoa. Correct expression and function of various genes and their protein products are required for production of fertile sperm. We have identified an infertility defect in Finnish Yorkshire boars caused by spermatogenic arrest. The aim of this study was to locate the disease associated region using genome wide screen with the PorcineSNP60 Beadchip and identify the causal mutation by candidate gene approach.

RESULTS

In the Finnish Yorkshire pig population the spermatogenic arrest (SA) defect appears to be of genetic origin and causes severe degeneration of germ cells and total absence of spermatozoa. Genome wide scan with the PorcineSNP60 Beadchip localized the SA defect to porcine chromosome 12 in a 2 Mbp region. Sequencing of a candidate gene Tex14 revealed a 51 bp insertion within exon 27, which caused differential splicing of the exon and created a premature translation stop codon. The expression of Tex14 was markedly down regulated in the testis of a SA affected boar compared to control boars and no protein product was identified by Western blotting. The SA insertion sequence was also found within intron 27 in all analyzed animals, thus the insertion appears to be a possible duplication event.

CONCLUSION

In this study we report the identification of a causal mutation for infertility caused by spermatogenic arrest at an early meiotic phase. Our results highlight the role of TEX14 specifically in spermatogenesis and the importance of specific genomic remodeling events as causes for inherited defects.

A novel application of pattern recognition for accurate SNP and indel discovery from high-throughput data: targeted resequencing of the glucocorticoid receptor co-chaperone FKBP5 in a Caucasian population

The detection of single nucleotide polymorphisms (SNPs) and insertion/deletions (indels) with precision from high-throughput data remains a significant bioinformatics challenge. Accurate detection is necessary before next-generation sequencing can routinely be used in the clinic. In research, scientific advances are inhibited by gaps in data, exemplified by the underrepresented discovery of rare variants, variants in non-coding regions and indels. The continued presence of false positives and false negatives prevents full automation and requires additional manual verification steps. Our methodology presents applications of both pattern recognition and sensitivity analysis to eliminate false positives and aid in the detection of SNP/indel loci and genotypes from high-throughput data. We chose FK506-binding protein 51(FKBP5) (6p21.31) for our clinical target because of its role in modulating pharmacological responses to physiological and synthetic glucocorticoids and because of the complexity of the genomic region. We detected genetic variation across a 160 kb region encompassing FKBP5. 613 SNPs and 57 indels, including a 3.3 kb deletion were discovered. We validated our method using three independent data sets and, with Sanger sequencing and Affymetrix and Illumina microarrays, achieved 99% concordance. Furthermore we were able to detect 267 novel rare variants and assess linkage disequilibrium. Our results showed both a sensitivity and specificity of 98%, indicating near perfect classification between true and false variants. The process is scalable and amenable to automation, with the downstream filters taking only 1.5h to analyze 96 individuals simultaneously. We provide examples of how our level of precision uncovered the interactions of multiple loci, their predicted influences on mRNA stability, perturbations of the hsp90 binding site, and individual variation in FKBP5 expression. Finally we show how our discovery of rare variants may change current conceptions of evolution at this locus.

Sequence shortening in the rodent ancestor

Insertions and deletions (indels), together with nucleotide substitutions, are major drivers of sequence evolution. An excess of deletions over insertions in genomic sequences-the so-called deletional bias-has been reported in a wide range of species, including mammals. However, this bias has not been found in the coding sequences of some mammalian species, such as human and mouse. To determine the strength of the deletional bias in mammals, and the influence of mutation and selection, we have quantified indels in both neutrally evolving noncoding sequences and protein-coding sequences, in six mammalian branches: human, macaque, ancestral primate, mouse, rat, and ancestral rodent. The results obtained with an improved algorithm for the placement of insertions in multiple alignments, Prank(+F), indicate that contrary to previous results, the only mammalian branch with a strong deletional bias is the rodent ancestral branch. We estimate that such a bias has resulted in an ~2.5% sequence loss of mammalian syntenic region in the ancestor of the mouse and rat. Further, a comparison of coding and noncoding sequences shows that negative selection is acting more strongly against mutations generating amino acid insertions than against mutations resulting in amino acid deletions. The strength of selection against indels is found to be higher in the rodent branches than in the primate branches, consistent with the larger effective population sizes of the rodents.