Genome-wide identification of pseudogenes capable of disease-causing gene conversion

Datasets-Based IMPDH1 Revisited: Heterozygous Missense Variants for Dominant Retinitis Pigmentosa While Truncation Variants Are Likely Non-Pathogenic

PURPOSE

Heterozygous variants of IMPDH1 are associated with autosomal dominant retinitis pigmentosa (adRP). The current study aims to investigate the characteristics of the adRP-associated variants.

METHODS

IMPDH1 variants from our exome sequencing dataset were retrieved and systemically evaluated through multiple online prediction tools, comparative genomics (in-house dataset, HGMD, and gnomAD), and phenotypic association. Potential pathogenic variants (PPVs) were further confirmed by Sanger sequencing and segregation analysis.

RESULTS

In total, seven heterozygous PPVs (six missenses and one inframe) were identified in 10 families with RP, in which six of the seven might be classified as pathogenic or likely pathogenic while one others as variants of uncertain significance. IMPDH1 variants contributed to 0.7% (10/1519) of RP families in our cohort, ranking the top four genes implicated in adRP. These adRP-associated variants were located in exons 8-10, a region within or downstream of the CBS domain. All these variants were predicted to be damaged by at least three of the six online prediction tools. Two truncation variants were considered non-pathogenic. Hitherto, 41 heterozygous variants of IMPDH1 were detected in 110 families in published literature, including 33 missenses, two inframes, and six truncations (including a synonymous variant affecting splicing). Of the 35 missense and inframe variants, most were clustered in exons 8-10 (77.1%, 27/35), including 18 (51.4%, 18/35) in exon 10 accounting for 70.9% (78/110) of the families. However, truncation variants were enriched in the general population with a pLI value of 0 (tolerated), and the reported variants in patients with RP did not cluster in specific region.

CONCLUSIONS

Our data together with comprehensive analysis of existing datasets suggest that causative variants of IMPDH1 are usually missense and mostly clustered in exons 8-10. Conversely, most missense variants outside this region and truncation variants should be interpreted with great care in clinical gene test.

4-(Phenylselanyl)-2H-chromen-2-one-Loaded Nanocapsule Suspension-A Promising Breakthrough in Pain Management: Comprehensive Molecular Docking, Formulation Design, and Toxicological and Pharmacological Assessments in Mice

Therapies for the treatment of pain and inflammation continue to pose a global challenge, emphasizing the significant impact of pain on patients' quality of life. Therefore, this study aimed to investigate the effects of 4-(Phenylselanyl)-2H-chromen-2-one (4-PSCO) on pain-associated proteins through computational molecular docking tests. A new pharmaceutical formulation based on polymeric nanocapsules was developed and characterized. The potential toxicity of 4-PSCO was assessed using Caenorhabditis elegans and Swiss mice, and its pharmacological actions through acute nociception and inflammation tests were also assessed. Our results demonstrated that 4-PSCO, in its free form, exhibited high affinity for the selected receptors, including p38 MAP kinase, peptidyl arginine deiminase type 4, phosphoinositide 3-kinase, Janus kinase 2, toll-like receptor 4, and nuclear factor-kappa β. Both free and nanoencapsulated 4-PSCO showed no toxicity in nematodes and mice. Parameters related to oxidative stress and plasma markers showed no significant change. Both treatments demonstrated antinociceptive and anti-edematogenic effects in the glutamate and hot plate tests. The nanoencapsulated form exhibited a more prolonged effect, reducing mechanical hypersensitivity in an inflammatory pain model. These findings underscore the promising potential of 4-PSCO as an alternative for the development of more effective and safer drugs for the treatment of pain and inflammation.

IMPDH1-associated autosomal dominant retinitis pigmentosa: natural history of novel variant Lys314Gln and a comprehensive literature search

BACKGROUND

Inosine monophosphate dehydrogenase (IMPDH) is a key regulatory enzyme in the de novo synthesis of the purine base guanine. Mutations in the inosine monophosphate dehydrogenase 1 gene (IMPDH1) are causative for RP10 autosomal dominant retinitis pigmentosa (adRP). This study reports a novel variant in a family with IMPDH1-associated retinopathy. We also performed a comprehensive review of all reported IMPDH1 disease causing variants with their associated phenotype.

MATERIALS AND METHODS

Multimodal imaging and functional studies documented the phenotype including best-corrected visual acuity (BCVA), fundus photograph, fundus autofluorescence (FAF), full field electroretinogram (ffERG), optical coherence tomography (OCT) and visual field (VF) data were collected. A literature search was performed in the PubMed and LOVD repositories.

RESULTS

We report 3 cases from a 2-generation family with a novel heterozygous likely pathogenic variant p. (Lys314Gln) (exon 10). The ophthalmic phenotype showed diffuse outer retinal atrophy with mild pigmentary changes with sparse pigmentary changes. FAF showed early macular involvement with macular hyperautofluorescence (hyperAF) surrounded by hypoAF. Foveal ellipsoid zone island can be found in the youngest patient but not in the older ones. The literature review identified a further 56 heterozygous, 1 compound heterozygous, and 2 homozygous variant. The heterozygous group included 43 missense, 3 in-frame, 1 nonsense, 2 frameshift, 1 synonymous, and 6 intronic variants. Exon 10 was noted as a hotspot harboring 18 variants.

CONCLUSIONS

We report a novel IMPDH1 variant. IMPDH1-associated retinopathy presents most frequently in the first decade of life with early macular involvement.

The Next Generation of Population-Based DFNB16 Carrier Screening and Diagnosis: STRC Copy-Number Variant Analysis from Genome Sequencing Data

BACKGROUND

Deafness, autosomal recessive 16 (DFNB16) is caused by compound heterozygous or homozygous variants in STRC and is the second most common form of genetic hearing loss. Due to the nearly identical sequences of STRC and the pseudogene STRCP1, analysis of this region is challenging in clinical testing.

METHODS

We developed a method that accurately identifies the copy number of STRC and STRCP1 using standard short-read genome sequencing. Then, we used whole genome sequencing (WGS) data to investigate the population distribution of STRC copy number in 6813 neonates and the correlation between STRC and STRCP1 copy number.

RESULTS

The comparison of WGS results with multiplex ligation-dependent probe amplification demonstrated high sensitivity (100%; 95% CI, 97.5%-100%) and specificity (98.8%; 95% CI, 97.7%-99.5%) in detecting heterozygous deletion of STRC from short-read genome sequencing data. The population analysis revealed that 5.22% of the general population has STRC copy number changes, almost half of which (2.33%; 95% CI, 1.99%-2.72%) were clinically significant, including heterozygous and homozygous STRC deletions. There was a strong inverse correlation between STRC and STRCP1 copy number.

CONCLUSIONS

We developed a novel and reliable method to determine STRC copy number based on standard short-read based WGS data. Incorporating this method into analytic pipelines would improve the clinical utility of WGS in the screening and diagnosis of hearing loss. Finally, we provide population-based evidence of pseudogene-mediated gene conversions between STRC and STRCP1.

PΨFinder: a practical tool for the identification and visualization of novel pseudogenes in DNA sequencing data

BACKGROUND

Processed pseudogenes (PΨgs) are disabled gene copies that are transcribed and may affect expression of paralogous genes. Moreover, their insertion in the genome can disrupt the structure or the regulatory region of a gene, affecting its expression level. These events have been identified as occurring mutations during cancer development, thus being able to identify PΨgs and their location will improve their impact on diagnostic testing, not only in cancer but also in inherited disorders.

RESULTS

We have implemented PΨFinder (P-psy-finder), a tool that identifies PΨgs, annotates known ones and predicts their insertion site(s) in the genome. The tool screens alignment files and provides user-friendly summary reports and visualizations. To demonstrate its applicability, we scanned 218 DNA samples from patients screened for hereditary colorectal cancer. We detected 423 PΨgs distributed in 96% of the samples, comprising 7 different parent genes. Among these, we confirmed the well-known insertion site of the SMAD4-PΨg within the last intron of the SCAI gene in one sample. While for the ubiquitous CBX3-PΨg, present in 82.6% of the samples, we found it reversed inserted in the second intron of the C15ORF57 gene.

CONCLUSIONS

PΨFinder is a tool that can automatically identify novel PΨgs from DNA sequencing data and determine their location in the genome with high sensitivity (95.92%). It generates high quality figures and tables that facilitate the interpretation of the results and can guide the experimental validation. PΨFinder is a complementary analysis to any mutational screening in the identification of disease-causing mutations within cancer and other diseases.

Pseudogenes: Four Decades of Discovery

A pseudogene is defined as a genomic DNA sequence that looks like a mutated or truncated version of a known functional gene. Nearly four decades since their first discovery it has been estimated that between ~12,000 and ~20,000 pseudogenes exist in the human genome. Early efforts to characterize functions for pseudogenes were unsuccessful, thus they were considered functionless relics of evolutionary selection, junk DNA or genetic fossils. Remarkably, an increasing number of pseudogenes have been reported to be expressed as RNA transcripts above and beyond levels considered accidental or spurious transcription. There is emerging evidence that some expressed pseudogene transcripts have biological functions and should be defined as a subclass of functional long noncoding RNAs (lncRNA). In this introductory chapter, I briefly summarize the history and the current knowledge of pseudogenes, and highlight the emerging functions of some pseudogenes in human biology and disease. This second iteration of Pseudogenes in Methods in Molecular Biology highlights new methodological approaches to investigate this intriguing family of lncRNAs and the extent of their biological function.

PB-Motif-A Method for Identifying Gene/Pseudogene Rearrangements With Long Reads: An Application to CYP21A2 Genotyping

Long read sequencing technologies have the potential to accurately detect and phase variation in genomic regions that are difficult to fully characterize with conventional short read methods. These difficult to sequence regions include several clinically relevant genes with highly homologous pseudogenes, many of which are prone to gene conversions or other types of complex structural rearrangements. We present PB-Motif, a new method for identifying rearrangements between two highly homologous genomic regions using PacBio long reads. PB-Motif leverages clustering and filtering techniques to efficiently report rearrangements in the presence of sequencing errors and other systematic artifacts. Supporting reads for each high-confidence rearrangement can then be used for copy number estimation and phased variant calling. First, we demonstrate PB-Motif's accuracy with simulated sequence rearrangements of PMS2 and its pseudogene PMS2CL using simulated reads sweeping over a range of sequencing error rates. We then apply PB-Motif to 26 clinical samples, characterizing CYP21A2 and its pseudogene CYP21A1P as part of a diagnostic assay for congenital adrenal hyperplasia. We successfully identify damaging variation and patient carrier status concordant with clinical diagnosis obtained from multiplex ligation-dependent amplification (MLPA) and Sanger sequencing. The source code is available at: github.com/zstephens/pb-motif.

Identification of key pseudogenes in nasopharyngeal carcinoma based on RNA-Seq analysis

BACKGROUND

Nasopharyngeal carcinoma (NPC) is a malignant head and neck tumor, and more than 70% of new cases are in East and Southeast Asia. However, association between NPC and pseudogenes playing important roles in genesis of multiple tumor types is still not clear and needs to be investigated.

METHODS

Using RNA-Sequencing (RNA-seq) technology, we analyzed pseudogene expression in 13 primary NPC and 6 recurrent NPC samples as well as their paracancerous counterparts. Quantitative PCR was used to validate the differentially expressed pseudogenes.

RESULTS

We found 251 differentially expressed pseudogenes including 73 up-regulated and 178 down-regulated ones between primary NPC and paracancerous tissues. Enrichment analysis of gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were conducted to filter out the key pseudogenes. We reported that pseudogenes from cytochrome P450 (CYP) family, such as CYP2F2P, CYP2G1P, CYP4F24P, CYP2B7P and CYP2G2P were significantly down-regulated in NPC compared to paracancerous tissues, while IGHV1OR15-2, IGHV3-11, FCGR1CP and IGHV3-69-1 belonging to Fc gamma receptors were significantly up-regulated. CYP2B7P, CYP2F2P and CYP4F26P were enriched in arachidonic acid metabolism pathway. The qRT-PCR analysis validated the lower expression of pseudogenes CYP2F2P and CYP2B7P in NPC tissues and cell lines compared to paracancerous tissues and normal human nasopharyngeal epithelial cell line. CYP2B7P overexpression weakened migratory and invasive capacity of NPC cell line. Moreover, the expression pattern of those pseudogenes in recurrent NPC tissues was different from the primary NPC.

CONCLUSION

This study suggested the role of pseudogenes in tumorigenesis and progression, potentially functioning as therapeutic targets to NPC.

Not a piece of junk anymore: Pseudogene T04B2.1 performs non-conventional regulatory role and modulates aggregation of α- synuclein and β-amyloid proteins in C. elegans

The conventional notions of pseudogenes being 'junk DNA' have largely been offset as research studies have established their role in multiple biological processes. Our studies towards identification of genetic modulators employing C. elegans model, that associate reproductive health and age-related neurodegenerative diseases, led us to identification and functional characterization of a pseudogene T04B2.1, which when knocked down, exacerbates the aggregation of α-Synuclein and β-Amyloid proteins, induces lipid deposition and alters morphometric endpoints in worms. Whole transcriptome analysis of worms under knockdown condition of T04B2.1 revealed an altered expression of 187 sequences, most of these being non-coding RNAs, miRNAs and piRNAs modulating the RNAi regulatory processes. Our gene ontology and pathway enrichment analysis demonstrated the role of T04B2.1 in protein quality control, metabolic pathways and development. We further performed a signature motif search and successfully identified a common motif that is present between all piRNA and miRNA molecules, which are significantly altered upon T04B2.1 silencing. This study unveils the non-conventional regulatory role of pseudogene T04B2.1 with respect to effects associated with neurodegenerative diseases and encourages further studies to decipher the regulatory mechanism governed by pseudogenes.

Strategies to Study the Functions of Pseudogenes in Mouse Models of Cancer

Aberrant expression of pseudogenes has been observed in many cancer types. Deregulated pseudogenes engage in a multitude of biological processes at the DNA, RNA, and protein levels and eventually facilitate disease progression. To investigate pseudogene functions in cancer, cell lines and cell line transplantation models have been widely used. However, cancer biology is best studied in the context of an intact organism. Here, we present various strategies to investigate pseudogenes in genetically engineered mouse models and discuss advantages and disadvantages of the different approaches.

Methods to Identify and Study the Evolution of Pseudogenes Using a Phylogenetic Approach

The discovery that pseudogenes are involved in important biological processes has excited enthusiasm and increased the research interest on them. An accurate detection and analysis of pseudogenes can be achieved using comparative methods, but only the use of phylogenetic tools can provide accurate information about their birth, their evolution and their death, hence about the impact that they have on genes and genomes. Here, phylogenetic methods that allow for studying pseudogene history are described.

Pseudogenes as Biomarkers and Therapeutic Targets in Human Cancers

Pseudogenes are commonly labeled as "junk DNA" given their perceived nonfunctional status. However, the advent of large-scale genomics projects prompted a revisit of pseudogene biology, highlighting their key functional and regulatory roles in numerous diseases, including cancers. Integrative analyses of cancer data have shown that pseudogenes can be transcribed and even translated, and that pseudogenic DNA, RNA, and proteins can interfere with the activity and function of key protein coding genes, acting as regulators of oncogenes and tumor suppressors. Capitalizing on the available clinical research, we are able to get an insight into the spread and variety of pseudogene biomarker and therapeutic potential. In this chapter, we describe pseudogenes that fulfill their role as diagnostic or prognostic biomarkers, both as unique elements and in collaboration with other genes or pseudogenes. We also report that the majority of prognostic pseudogenes are overexpressed and exert an oncogenic role in colorectal, liver, lung, and gastric cancers. Finally, we highlight a number of pseudogenes that can establish future therapeutic avenues.

CRISPR-Directed Therapeutic Correction at the NCF1 Locus Is Challenged by Frequent Incidence of Chromosomal Deletions

Resurrection of non-processed pseudogenes may increase the efficacy of therapeutic gene editing, upon simultaneous targeting of a mutated gene and its highly homologous pseudogenes. To investigate the potency of this approach for clinical gene therapy of human diseases, we corrected a pseudogene-associated disorder, the immunodeficiency p47^phox-deficient chronic granulomatous disease (p47^phox CGD), using clustered regularly interspaced short palindromic repeats-associated nuclease Cas9 (CRISPR-Cas9) to target mutated neutrophil cytosolic factor 1 (NCF1). Being separated by less than two million base pairs, NCF1 and two pseudogenes are closely co-localized on chromosome 7. In healthy people, a two-nucleotide GT deletion (ΔGT) is present in the NCF1B and NCF1C pseudogenes only. In the majority of patients with p47^phox CGD, the NCF1 gene is inactivated due to a ΔGT transfer from one of the two non-processed pseudogenes. Here we demonstrate that concurrent targeting and correction of mutated NCF1 and its pseudogenes results in therapeutic CGD phenotype correction, but also causes potentially harmful chromosomal deletions between the targeted loci in a p47^phox-deficient CGD cell line model. Therefore, development of genome-editing-based treatment of pseudogene-related disorders mandates thorough safety examination, as well as technological advances, limiting concurrent induction of multiple double-strand breaks on a single chromosome.

Both male and female gamete generating cells produce processed pseudogenes in the human genome

The human genome contains an unusually large number of processed pseudogenes. The fact that processed pseudogenes are roughly 33% more abundant in our X chromosome than in our autosomes suggests that this overabundance is the result of the fact that human oogenesis is much longer than that of non-mammalian species. Here, we analyze the origins of the processed pseudogenes found on the human Y chromosome to determine whether human spermatogenesis also contribute to this overabundance. Our results show that human processed pseudogenes not only retrotranspose to the Y chromosome, but are also produced by genes on the Y chromosome. Furthermore, the fact that X chromosomes are three times more abundant than Y chromosomes likely explains why the euchromatic density of processed pseudogenes is three times higher in the X chromosome than in the Y chromosome. The large number of processed pseudogenes found in our genome is therefore due to the low substrate specificity of the L1 reverse transcriptase responsible for the reverse transcription of germline mRNA molecules into processed pseudogenes, as well as the life-long production of both male and female gametes.

An analysis of effects of heterozygosity in dairy cattle for bovine tuberculosis resistance

Genetic selection of cattle more resistant to bovine tuberculosis (bTB) may offer a complementary control strategy. Hypothesising underlying non-additive genetic variation, we present an approach using genome-wide high density markers to identify genomic loci with dominance effects on bTB resistance and to test previously published regions with heterozygote advantage in bTB. Our data comprised 1151 Holstein-Friesian cows from Northern Ireland, confirmed bTB cases and controls, genotyped with the 700K Illumina BeadChip. Genome-wide markers were tested for associations between heterozygosity and bTB status using marker-based relationships. Results were tested for robustness against genetic structure, and the genotypic frequencies of a significant locus were tested for departures from Hardy-Weinberg equilibrium. Genomic regions identified in our study and in previous publications were tested for dominance effects. Genotypic effects were estimated through ASReml mixed models. A SNP (rs43032684) on chromosome 6 was significant at the chromosome-wide level, explaining 1.7% of the phenotypic variance. In the controls, there were fewer heterozygotes for rs43032684 (P < 0.01) with the genotypic values suggesting that heterozygosity confers a heterozygote disadvantage. The region surrounding rs43032684 had a significant dominance effect (P < 0.01). SNP rs43032684 resides within a pseudogene with a parental gene involved in macrophage response to infection and within a copy-number-variation region previously associated with nematode resistance. No dominance effect was found for the region on chromosome 11, as indicated by a previous candidate region bTB study. These findings require further validation with large-scale data.

Genomic disorders 20 years on-mechanisms for clinical manifestations

Genomic disorders result from copy-number variants (CNVs) or submicroscopic rearrangements of the genome rather than from single nucleotide variants (SNVs). Diverse technologies, including array comparative genomic hybridization (aCGH) and single nucleotide polymorphism (SNP) microarrays, and more recently, whole genome sequencing and whole-exome sequencing, have enabled robust genome-wide unbiased detection of CNVs in affected individuals and in reportedly healthy controls. Sequencing of breakpoint junctions has allowed for elucidation of upstream mechanisms leading to genomic instability and resultant structural variation, whereas studies of the association between CNVs and specific diseases or susceptibility to morbid traits have enhanced our understanding of the downstream effects. In this review, we discuss the hallmarks of genomic disorders as they were defined during the first decade of the field, including genomic instability and the mechanism for rearrangement defined as nonallelic homologous recombination (NAHR); recurrent vs nonrecurrent rearrangements; and gene dosage sensitivity. Moreover, we highlight the exciting advances of the second decade of this field, including a deeper understanding of genomic instability and the mechanisms underlying complex rearrangements, mechanisms for constitutional and somatic chromosomal rearrangements, structural intra-species polymorphisms and susceptibility to NAHR, the role of CNVs in the context of genome-wide copy number and single nucleotide variation, and the contribution of noncoding CNVs to human disease.

Frequent nonallelic gene conversion on the human lineage and its effect on the divergence of gene duplicates

Gene conversion is the copying of a genetic sequence from a "donor" region to an "acceptor." In nonallelic gene conversion (NAGC), the donor and the acceptor are at distinct genetic loci. Despite the role NAGC plays in various genetic diseases and the concerted evolution of gene families, the parameters that govern NAGC are not well characterized. Here, we survey duplicate gene families and identify converted tracts in 46% of them. These conversions reflect a large GC bias of NAGC. We develop a sequence evolution model that leverages substantially more information in duplicate sequences than used by previous methods and use it to estimate the parameters that govern NAGC in humans: a mean converted tract length of 250 bp and a probability of [Formula: see text] per generation for a nucleotide to be converted (an order of magnitude higher than the point mutation rate). Despite this high baseline rate, we show that NAGC slows down as duplicate sequences diverge-until an eventual "escape" of the sequences from its influence. As a result, NAGC has a small average effect on the sequence divergence of duplicates. This work improves our understanding of the NAGC mechanism and the role that it plays in the evolution of gene duplicates.

snRNP proteins in health and disease

Split gene architecture of most human genes requires removal of intervening sequences by mRNA splicing that occurs on large multiprotein complexes called spliceosomes. Mutations compromising several spliceosomal components have been recorded in degenerative syndromes and haematological neoplasia, thereby highlighting the importance of accurate splicing execution in homeostasis of assorted adult tissues. Moreover, insufficient splicing underlies defective development of craniofacial skeleton and upper extremities. This review summarizes recent advances in the understanding of splicing factor function deduced from cryo-EM structures. We combine these data with the characterization of splicing factors implicated in hereditary or somatic disorders, with a focus on potential functional consequences the mutations may elicit in spliceosome assembly and/or performance. Given aberrant splicing or perturbations in splicing efficiency substantially underpin disease pathogenesis, profound understanding of the mis-splicing principles may open new therapeutic vistas. In three major sections dedicated to retinal dystrophies, hereditary acrofacial syndromes, and haematological malignancies, we delineate the noticeable variety of conditions associated with dysfunctional splicing and accentuate recurrent patterns in splicing defects.

Mutations in spliceosomal proteins and retina degeneration

A majority of human genes contain non-coding intervening sequences – introns that must be precisely excised from the pre-mRNA molecule. This event requires the coordinated action of five major small nuclear ribonucleoprotein particles (snRNPs) along with additional non-snRNP splicing proteins. Introns must be removed with nucleotidal precision, since even a single nucleotide mistake would result in a reading frame shift and production of a non-functional protein. Numerous human inherited diseases are caused by mutations that affect splicing, including mutations in proteins which are directly involved in splicing catalysis. One of the most common hereditary diseases associated with mutations in core splicing proteins is retinitis pigmentosa (RP). So far, mutations in more than 70 genes have been connected to RP. While the majority of mutated genes are expressed specifically in the retina, eight target genes encode for ubiquitous core snRNP proteins (Prpf3, Prpf4, Prpf6, Prpf8, Prpf31, and SNRNP200/Brr2) and splicing factors (RP9 and DHX38). Why mutations in spliceosomal proteins, which are essential in nearly every cell in the body, causes a disease that displays such a tissue-specific phenotype is currently a mystery. In this review, we recapitulate snRNP functions, summarize the missense mutations which are found in spliceosomal proteins as well as their impact on protein functions and discuss specific models which may explain why the retina is sensitive to these mutations.

Integrating multiple immunogenetic data sources for feature extraction and mining somatic hypermutation patterns: the case of "towards analysis" in chronic lymphocytic leukaemia

BACKGROUND

Somatic Hypermutation (SHM) refers to the introduction of mutations within rearranged V(D)J genes, a process that increases the diversity of Immunoglobulins (IGs). The analysis of SHM has offered critical insight into the physiology and pathology of B cells, leading to strong prognostication markers for clinical outcome in chronic lymphocytic leukaemia (CLL), the most frequent adult B-cell malignancy. In this paper we present a methodology for integrating multiple immunogenetic and clinocobiological data sources in order to extract features and create high quality datasets for SHM analysis in IG receptors of CLL patients. This dataset is used as the basis for a higher level integration procedure, inspired form social choice theory. This is applied in the Towards Analysis, our attempt to investigate the potential ontogenetic transformation of genes belonging to specific stereotyped CLL subsets towards other genes or gene families, through SHM.

RESULTS

The data integration process, followed by feature extraction, resulted in the generation of a dataset containing information about mutations occurring through SHM. The Towards analysis performed on the integrated dataset applying voting techniques, revealed the distinct behaviour of subset #201 compared to other subsets, as regards SHM related movements among gene clans, both in allele-conserved and non-conserved gene areas. With respect to movement between genes, a high percentage movement towards pseudo genes was found in all CLL subsets.

CONCLUSIONS

This data integration and feature extraction process can set the basis for exploratory analysis or a fully automated computational data mining approach on many as yet unanswered, clinically relevant biological questions.

Navigating highly homologous genes in a molecular diagnostic setting: a resource for clinical next-generation sequencing

PURPOSE

Next-generation sequencing (NGS) is now routinely used to interrogate large sets of genes in a diagnostic setting. Regions of high sequence homology continue to be a major challenge for short-read technologies and can lead to false-positive and false-negative diagnostic errors. At the scale of whole-exome sequencing (WES), laboratories may be limited in their knowledge of genes and regions that pose technical hurdles due to high homology. We have created an exome-wide resource that catalogs highly homologous regions that is tailored toward diagnostic applications.

METHODS

This resource was developed using a mappability-based approach tailored to current Sanger and NGS protocols.

RESULTS

Gene-level and exon-level lists delineate regions that are difficult or impossible to analyze via standard NGS. These regions are ranked by degree of affectedness, annotated for medical relevance, and classified by the type of homology (within-gene, different functional gene, known pseudogene, uncharacterized noncoding region). Additionally, we provide a list of exons that cannot be analyzed by short-amplicon Sanger sequencing.

CONCLUSION

This resource can help guide clinical test design, supplemental assay implementation, and results interpretation in the context of high homology.Genet Med 18 12, 1282-1289.

Interlocus gene conversion explains at least 2.7% of single nucleotide variants in human segmental duplications

Background

Interlocus gene conversion (IGC) is a recombination-based mechanism that results in the unidirectional transfer of short stretches of sequence between paralogous loci. Although IGC is a well-established mechanism of human disease, the extent to which this mutagenic process has shaped overall patterns of segregating variation in multi-copy regions of the human genome remains unknown. One expected manifestation of IGC in population genomic data is the presence of one-to-one paralogous SNPs that segregate identical alleles.

Results

Here, I use SNP genotype calls from the low-coverage phase 3 release of the 1000 Genomes Project to identify 15,790 parallel, shared SNPs in duplicated regions of the human genome. My approach for identifying these sites accounts for the potential redundancy of short read mapping in multi-copy genomic regions, thereby effectively eliminating false positive SNP calls arising from paralogous sequence variation. I demonstrate that independent mutation events to identical nucleotides at paralogous sites are not a significant source of shared polymorphisms in the human genome, consistent with the interpretation that these sites are the outcome of historical IGC events. These putative signals of IGC are enriched in genomic contexts previously associated with non-allelic homologous recombination, including clear signals in gene families that form tandem intra-chromosomal clusters.

Conclusions

Taken together, my analyses implicate IGC, not point mutation, as the mechanism generating at least 2.7 % of single nucleotide variants in duplicated regions of the human genome.

Electronic supplementary material

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

Pseudogene-derived endogenous siRNAs and their function

Pseudogenes were once considered genomic fossils, but recent studies indicate that they may function as gene regulators through the generation of endogenous small interfering RNAs (esiRNAs), antisense RNAs, and decoys for microRNAs. In this review, we summarize pseudogene study methods, emphasizing relevant publicly available resources, and we describe a systematic pipeline to identify pseudogene-derived esiRNAs and their targets, which can lead to a deeper understanding of pseudogene function.

Long non-coding RNAs in cancer: implications for personalized therapy

Long non-coding RNAs (lncRNAs, pseudogenes and circRNAs) have recently come into light as powerful players in cancer pathogenesis and it is becoming increasingly clear that they have the potential of greatly contributing to the spread and success of personalized cancer medicine. In this concise review, we briefly introduce these three classes of long non-coding RNAs. We then discuss their applications as diagnostic and prognostic biomarkers. Finally, we describe their appeal as targets and as drugs, while pointing out the limitations that still lie ahead of their definitive entry into clinical practice.

[High gene conversion frequency between genes encoding 2-deoxyglucose-6-phosphate phosphatase in 3 Saccharomyces species]

Gene conversions are nonreciprocal sequence exchanges between genes. They are relatively common in Saccharomyces cerevisiae, but few studies have investigated the evolutionary fate of gene conversions or their functional impacts. Here, we analyze the evolution and impact of gene conversions between the two genes encoding 2-deoxyglucose-6-phosphate phosphatase in S. cerevisiae, Saccharomyces paradoxus and Saccharomyces mikatae. Our results demonstrate that the last half of these genes are subject to gene conversions among these three species. The greater similarity and the greater percentage of GC nucleotides in the converted regions, as well as the absence of long regions of adjacent common converted sites, suggest that these gene conversions are frequent and occur independently in all three species. The high frequency of these conversions probably result from the fact that they have little impact on the protein sequences encoded by these genes.

Methods to study the occurrence and the evolution of pseudogenes through a phylogenetic approach

During the last few years, the study of pseudogenes has excited enthusiasm, because it has been proven that at least some of them are involved in important biological processes. An accurate detection and analysis of pseudogenes can be achieved using comparative methods, but only the use of phylogenetic tools can provide accurate information about their birth, their evolution and their death, hence about the impact that they have on genes and genomes. Here, phylogenetic methods that allow studying pseudogene history are described.

Computational methods for pseudogene annotation based on sequence homology

The number of complete genome sequences explodes more and more with each passing year. Thus, methods for genome annotation need to be honed constantly to handle the deluge of information. Annotation of pseudogenes (i.e., gene copies that appear not to make a functional protein) in genomes is a persistent problem; here, we overview pseudogene annotation methods that are based on the detection of sequence homology in genomic DNA.

Pseudogene redux with new biological significance

The study of pseudogenes, originally dismissed as genomic relics of evolutionary selection, has seen a resurgence in scientific literature, in addition to being a peculiar topic of discussion in theological debates. For a long time, pseudogenes have been touted as a beacon of natural selection and a definitive proof of evolution due to the slow mutation rate that differentiated them from their parental genes and ultimately caused their genetic demise as functional genes. It now seems that "creationists" have co-opted some recent reports identifying unheralded biological functions to pseudogens and other noncoding RNAs as evidence to undermine the existence of evolution and supporting intelligent design. This issue of Methods in Molecular Biology focused on pseudogenes will certainly not end, nor enter this debate; however, scientists who are also genomics and pseudogene enthusiasts will certainly appreciate that many scientists are thinking about these particular genetic elements in new and interesting ways. With this new interest in a biological significance and "non-junk" role for pseudogenes and other noncoding RNAs, new methods and approaches are being developed to unlock the mystery of these ancient artifacts we know as pseudogenes. In this brief introductory chapter we highlight the renewed interest in pseudogenes and review a rationale for intensification of pseudogene-related research.

Computational methods of identification of pseudogenes based on functionality: entropy and GC content

Spectral entropy and GC content analyses reveal comprehensive structural features of DNA sequences. To illustrate the significance of these features, we analyze the β-esterase gene cluster, including the Est-6 gene and the ψEst-6 putative pseudogene, in seven species of the Drosophila melanogaster subgroup. The spectral entropies show distinctly lower structural ordering for ψEst-6 than for Est-6 in all species studied. However, entropy accumulation is not a completely random process for either gene and it shows to be nucleotide dependent. Furthermore, GC content in synonymous positions is uniformly higher in Est-6 than in ψEst-6, in agreement with the reduced GC content generally observed in pseudogenes and nonfunctional sequences. The observed differences in entropy and GC content reflect an evolutionary shift associated with the process of pseudogenization and subsequent functional divergence of ψEst-6 and Est-6 after the duplication event. The data obtained show the relevance and significance of entropy and GC content analyses for pseudogene identification and for the comparative study of gene-pseudogene evolution.

Mutational analysis of CYP21A2 gene and CYP21A1P pseudogene: long-range PCR on genomic DNA

CYP21A2, the gene that codes for P450c21 (Steroid 21-hydroxylase), has a duplicated pseudogene called CYP21A1P. The gene and the pseudogene share 98 % and 96 % sequence homology in exons and in noncoding sequences, respectively, and are located 30 kb apart within the HLA class III human histocompatibility complex locus on chromosome 6p21.3. CYP21A1P is inactive due to the presence of 11 deteriorated mutations in its coding region. These mutations can be transferred to the functional CYP21A2 through intergenic recombination during meiosis or mitosis and lead to the congenital adrenal hyperplasia (CAH) resulting from 21-hydroxylase deficiency. Conversely, portions of CYP21A2 sequence can be transferred to CYP21A1P, modifying the haplotype. Here, we describe a well-established protocol that can be used to unambiguously study the mutational profile of CYP21A2 gene and CYP21A1P pseudogene. The protocol is based on long-range PCR amplification with allele-specific primers, followed by DNA sequencing of smaller fragments.

Purifying selection against gene conversions in the folate receptor genes of primates

We characterized the gene conversions between the human folate receptor (FOLR) genes and those of five other primate species. We found 26 gene conversions having an average length of 534 nucleotides. The length of these conversions is correlated with sequence similarity, converted regions have a higher GC-content and the average size of converted regions from a functional donor to another functional donor is significantly smaller than the average size from a functional donor to a pseudogene. Furthermore, the few conversions observed in the FOLR1 and FOLR2 genes did not change any amino acids in their coding regions and did not affect their promoter regions. In contrast, the promoter and coding regions of the FOLR3 gene are frequently converted and these conversions changed many amino acids in marmoset. These results suggest that purifying selection is limiting the functional impact that frequent gene conversions have on functional folate receptor genes.

Pseudogenes and their composers: delving in the 'debris' of human genome

Pseudogenes, the nonfunctional homologs of functional genes and thus exemplified as 'genomic fossils' provide intriguing snapshots of the evolutionary history of human genome. These defunct copies generally arise by retrotransposition or duplication followed by various genetic disablements. In this study, focusing on human pseudogenes and their functional homologues we describe their characteristic features and relevance to protein sequence evolution. We recapitulate that pseudogenes harbor disease-causing degenerative sequence variations in conjunction with the immense disease gene association of their progenitors. Furthermore, we also discuss the issue of functional resurrection and the potentiality observed in some pseudogenes to regulate their functional counterparts.

pseudoMap: an innovative and comprehensive resource for identification of siRNA-mediated mechanisms in human transcribed pseudogenes

RNA interference (RNAi) is a gene silencing process within living cells, which is controlled by the RNA-induced silencing complex with a sequence-specific manner. In flies and mice, the pseudogene transcripts can be processed into short interfering RNAs (siRNAs) that regulate protein-coding genes through the RNAi pathway. Following these findings, we construct an innovative and comprehensive database to elucidate siRNA-mediated mechanism in human transcribed pseudogenes (TPGs). To investigate TPG producing siRNAs that regulate protein-coding genes, we mapped the TPGs to small RNAs (sRNAs) that were supported by publicly deep sequencing data from various sRNA libraries and constructed the TPG-derived siRNA-target interactions. In addition, we also presented that TPGs can act as a target for miRNAs that actually regulate the parental gene. To enable the systematic compilation and updating of these results and additional information, we have developed a database, pseudoMap, capturing various types of information, including sequence data, TPG and cognate annotation, deep sequencing data, RNA-folding structure, gene expression profiles, miRNA annotation and target prediction. As our knowledge, pseudoMap is the first database to demonstrate two mechanisms of human TPGs: encoding siRNAs and decoying miRNAs that target the parental gene. pseudoMap is freely accessible at http://pseudomap.mbc.nctu.edu.tw/.

Database URL:http://pseudomap.mbc.nctu.edu.tw/

Pseudogenes: newly discovered players in human cancer

Because they are generally noncoding and thus considered nonfunctional and unimportant, pseudogenes have long been neglected. Recent advances have established that the DNA of a pseudogene, the RNA transcribed from a pseudogene, or the protein translated from a pseudogene can have multiple, diverse functions and that these functions can affect not only their parental genes but also unrelated genes. Therefore, pseudogenes have emerged as a previously unappreciated class of sophisticated modulators of gene expression, with a multifaceted involvement in the pathogenesis of human cancer.

Molecular trajectories leading to the alternative fates of duplicate genes

Gene duplication generates extra gene copies in which mutations can accumulate without risking the function of pre-existing genes. Such mutations modify duplicates and contribute to evolutionary novelties. However, the vast majority of duplicates appear to be short-lived and experience duplicate silencing within a few million years. Little is known about the molecular mechanisms leading to these alternative fates. Here we delineate differing molecular trajectories of a relatively recent duplication event between humans and chimpanzees by investigating molecular properties of a single duplicate: DNA sequences, gene expression and promoter activities. The inverted duplication of the Glutathione S-transferase Theta 2 (GSTT2) gene had occurred at least 7 million years ago in the common ancestor of African great apes and is preserved in chimpanzees (Pan troglodytes), whereas a deletion polymorphism is prevalent in humans. The alternative fates are associated with expression divergence between these species, and reduced expression in humans is regulated by silencing mutations that have been propagated between duplicates by gene conversion. In contrast, selective constraint preserved duplicate divergence in chimpanzees. The difference in evolutionary processes left a unique DNA footprint in which dying duplicates are significantly more similar to each other (99.4%) than preserved ones. Such molecular trajectories could provide insights for the mechanisms underlying duplicate life and death in extant genomes.

Pseudogenes

Pseudogenes are ubiquitous and abundant in genomes. Pseudogenes were once called “genomic fossils” and treated as “junk DNA” several years. Nevertheless, it has been recognized that some pseudogenes play essential roles in gene regulation of their parent genes, and many pseudogenes are transcribed into RNA. Pseudogene transcripts may also form small interfering RNA or decrease cellular miRNA concentration. Thus, pseudogenes regulate tumor suppressors and oncogenes. Their essential functions draw the attention of our research group in my current work on heat shock protein 90: a chaperone of oncogenes. The paper reviews our current knowledge on pseudogenes and evaluates preliminary results of the chaperone data. Current efforts to understand pseudogenes interactions help to understand the functions of a genome.

Evolutionary conservation and disease gene association of the human genes composing pseudogenes

Pseudogenes, the 'genomic fossils' present portrayal of evolutionary history of human genome. The human genes configuring pseudogenes are also now coming forth as important resources in the study of human protein evolution. In this communication, we explored evolutionary conservation of the genes forming pseudogenes over the genes lacking any pseudogene and delving deeper, we probed an evolutionary rate difference between the disease genes in the two groups. We illustrated this differential evolutionary pattern by gene expressivity, number of regulatory miRNA targeting per gene, abundance of protein complex forming genes and lesser percentage of protein intrinsic disorderness. Furthermore, pseudogenes are observed to harbor sequence variations, over their entirety, those become degenerative disease-causing mutations though the disease involvement of their progenitors is still unexplored. Here, we unveiled an immense association of disease genes in the genes casting pseudogenes in human. We interpreted the issue by disease associated miRNA targeting, genes containing polymorphisms in miRNA target sites, abundance of genes having disease causing non-synonymous mutations, disease gene specific network properties, presence of genes having repeat regions, affluence of dosage sensitive genes and the presence of intrinsically unstructured protein regions.

Identification of gene mutations in autosomal dominant polycystic kidney disease through targeted resequencing

Mutations in two large multi-exon genes, PKD1 and PKD2, cause autosomal dominant polycystic kidney disease (ADPKD). The duplication of PKD1 exons 1-32 as six pseudogenes on chromosome 16, the high level of allelic heterogeneity, and the cost of Sanger sequencing complicate mutation analysis, which can aid diagnostics of ADPKD. We developed and validated a strategy to analyze both the PKD1 and PKD2 genes using next-generation sequencing by pooling long-range PCR amplicons and multiplexing bar-coded libraries. We used this approach to characterize a cohort of 230 patients with ADPKD. This process detected definitely and likely pathogenic variants in 115 (63%) of 183 patients with typical ADPKD. In addition, we identified atypical mutations, a gene conversion, and one missed mutation resulting from allele dropout, and we characterized the pattern of deep intronic variation for both genes. In summary, this strategy involving next-generation sequencing is a model for future genetic characterization of large ADPKD populations.

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.

Genome-wide expression of non-coding RNA and global chromatin modification

Traditionally, we know that genomic DNA will produce transcripts named messenger RNA and then translate into protein following the instruction of genetic central dogma, and RNA works here as a pass-by messenger. Now increasing evidence shows that RNA is a key regulator as well as a message transmitter. It is discovered by next-generation sequencing techniques that most genomic DNA are generally transcribed to non-coding RNA, highly beyond the percentage of coding mRNA. These non-coding RNAs (ncRNAs), belonging to several groups, have critical roles in many cellular processes, expanding our understanding of the RNA world. We review here the different categories of ncRNA according to genome location and how ncRNAs guide and recruit chromatin modification complex to specific loci of genome to modulate gene expression by affecting chromatin state.

Gene conversions in the growth hormone gene family of primates: stronger homogenizing effects in the Hominidae lineage

In humans, the growth hormone/chorionic somatomammotropin gene family is composed of five highly similar genes. We characterized the gene conversions that occurred between the growth hormone genes of 11 primate species. We detected 48 conversions using GENECONV and others were only detected using phylogenetic analyses. Gene conversions were detected in all species analyzed, their average size (±standard deviation) is 197.8±230.4 nucleotides, the size of the conversions is correlated with sequence similarity and converted regions are significantly more GC-rich than non-converted regions. Gene conversions have a stronger homogenizing effect in Hominidae genes than in other primate species. They are also less frequent in conserved gene regions and towards functionally important genes. This suggests that the high degree of sequence similarity observed between the growth hormone genes of primate species is a consequence of frequent gene conversions in gene regions which are under little selective constraints.

Analysis of the CYP21A1P pseudogene: indication of mutational diversity and CYP21A2-like and duplicated CYP21A2 genes

The CYP21A1P gene downstream of the XA gene, carrying 15 deteriorated mutations, is a nonfunctional pseudogene that shares 98% nucleotide sequence homology with CYP21A2 located on chromosome 6p21.3. However, these mutations in the CYP21A1P gene are not totally involved in each individual. From our analysis of 100 healthy ethnic Chinese (i.e., Taiwanese) (n=200 chromosomes) using the polymerase chain reaction (PCR) products combined with an amplification-created restriction site (ACRS) method and DNA sequencing, we found that approximately 10% of CYP21A1P alleles (n=195 chromosomes) presented the CYP21A2 sequence; frequencies of P30, V281, Q318, and R356 in that locus were approximately 24%, 21%, 11%, and 34%, respectively, and approximately 90% of the CYP21A1P alleles had 15 mutated loci. In addition, approximately 2.5% (n=5 chromosomes) showed four haplotypes of the 3.7-kb TaqI-produced fragment of the CYP21A2-like gene and one duplicated CYP21A2 gene. We conclude that the pseudogene of the CYP21A1P mutation presents diverse variants. Moreover, the existence of the CYP21A2-like gene is more abundant than that of the duplicated CYP21A2 gene downstream of the XA gene and could not be distinguished from the CYP21A2-TNXB gene; thus, it may be misdiagnosed by previously established methods for congenital adrenal hyperplasia caused by a 21-hydroxylase deficiency.

Genetic diversity of the allodeterminant alr2 in Hydractinia symbiolongicarpus

Hydractinia symbiolongicarpus, a colonial cnidarian (class Hydrozoa) epibiont on hermit crab shells, is well established as a model for genetic studies of allorecognition. Recently, two linked loci, allorecognition (alr) 1 and alr2, were identified by positional cloning and shown to be major determinants of histocompatibility. Both genes encode putative transmembrane proteins with hypervariable extracellular domains similar to immunoglobulin (Ig)-like domains. We sought to characterize the naturally occurring variation at the alr2 locus and to understand the origins of this molecular diversity. We examined full-length cDNA coding sequences derived from a sample of 21 field-collected colonies, including 18 chosen haphazardly and two laboratory reference strains. Of the 35 alleles recovered from the 18 unbiased samples, 34 encoded unique gene products. We identified two distinct structural classes of alleles that varied over a large central region of the gene but both possessed highly polymorphic extracellular domains I, similar to an Ig-like V-set domain. The discovery of structurally chimeric alleles provided evidence that interallelic recombination may contribute to alr2 variation. Comparisons of the genomic region encompassing alr2 from two field-derived haplotypes and one laboratory reference sequence revealed a history of structural variation at the haplotype level as well. Maintenance of large numbers of equally rare alleles in a natural population is a hallmark of negative frequency-dependent selection and is expected to produce high levels of heterozygosity. The observed alr2 allelic diversity is comparable with that found in immune recognition molecules such as human leukocyte antigens, B cell Igs, or natural killer cell Ig-like receptors.

Segmental duplications in the human genome reveal details of pseudogene formation

Duplicated pseudogenes in the human genome are disabled copies of functioning parent genes. They result from block duplication events occurring throughout evolutionary history. Relatively recent duplications (with sequence similarity ≥90% and length ≥1 kb) are termed segmental duplications (SDs); here, we analyze the interrelationship of SDs and pseudogenes. We present a decision-tree approach to classify pseudogenes based on their (and their parents’) characteristics in relation to SDs. The classification identifies 140 novel pseudogenes and makes possible improved annotation for the 3172 pseudogenes located in SDs. In particular, it reveals that many pseudogenes in SDs likely did not arise directly from parent genes, but are the result of a multi-step process. In these cases, the initial duplication or retrotransposition of a parent gene gives rise to a ‘parent pseudogene’, followed by further duplication creating duplicated–duplicated or duplicated–processed pseudogenes, respectively. Moreover, we can precisely identify these parent pseudogenes by overlap with ancestral SD loci. Finally, a comparison of nucleotide substitutions per site in a pseudogene with its surrounding SD region allows us to estimate the time difference between duplication and disablement events, and this suggests that most duplicated pseudogenes in SDs were likely disabled around the time of the original duplication.