Direct evaluation of neuroaxonal degeneration with the causative genes of neurodegenerative diseases in Drosophila using the automated axon quantification system, MeDUsA.

Drosophila is an excellent model organism for studying human neurodegenerative diseases (NDs). However, there is still almost no experimental system that could directly observe the degeneration of neurons and automatically quantify axonal degeneration. In this study, we created MeDUsA (a 'method for the quantification of degeneration using fly axons'), a standalone executable computer program based on Python that combines a pre-trained deep-learning masking tool with an axon terminal counting tool. This software automatically quantifies the number of retinal R7 axons in Drosophila from a confocal z-stack image series. Using this software, we were able to directly demonstrate that axons were degenerated by the representative causative genes of NDs for the first time in Drosophila. The fly retinal axon is an excellent experimental system that is capable of mimicking the pathology of axonal degeneration in human NDs. MeDUsA rapidly and accurately quantifies axons in Drosophila photoreceptor neurons. It enables large-scale research into axonal degeneration, including screening to identify genes or drugs that mediate axonal toxicity caused by ND proteins and diagnose the pathological significance of novel variants of human genes in axons.

Truncating variants of the sterol recognition region of SHH cause hypertelorism phenotype rather than hypotelorism-holoprosencephaly.

Sonic hedgehog signaling molecule (SHH) is a key molecule in the cilia-mediated signaling pathway and a critical morphogen in embryogenesis. The association between loss-of-function variants of SHH and holoprosencephaly is well established. In mice experiments, reduced or increased signaling of SHH have been shown to be associated with narrowing or excessive expansion of the facial midline, respectively. Herein, we report two unrelated patients with de novo truncating variants of SHH presenting with hypertelorism rather than hypotelorism. The first patient was a 13-year-old girl. Her facial features included hypertelorism, strabismus, telecanthus, malocclusion, frontal bossing, and wide widow's peak. She had borderline developmental delay and agenesis of the corpus callosum. She had a nonsense variant of SHH: Chr7(GRCh38):g.155802987C > T, NM_000193.4:c.1302G > A, p.(Trp434*). The second patient was a 25-year-old girl. Her facial features included hypertelorism and wide widow's peak. She had developmental delay and agenesis of the corpus callosum. She had a frameshift variant of SHH: Chr7(GRCh38):g.155803072_155803074delCGGinsT, NM_000193.4:c.1215_1217delCCGinsA, p.(Asp405Glufs*92). The hypertelorism phenotype contrasts sharply with the prototypical hypotelorism-holoprosencephaly phenotype associated with loss-of-function of SHH. We concluded that a subset of truncating variants of SHH could be associated with hypertelorism rather than hypotelorism.

BCS1L mutations produce Fanconi syndrome with developmental disability.

Fanconi syndrome is a functional disorder of the proximal tubule, characterized by pan-aminoaciduria, glucosuria, hypophosphatemia, and metabolic acidosis. With the advancements in gene analysis technologies, several causative genes are identified for Fanconi syndrome. Several mitochondrial diseases cause Fanconi syndrome and various systemic symptoms; however, it is rare that the main clinical symptoms in such disorders are Fanconi syndrome without systematic active diseases like encephalomyopathy or cardiomyopathy. In this study, we analyzed two families exhibiting Fanconi syndrome, developmental disability and mildly elevated liver enzyme levels. Whole-exome sequencing (WES) detected compound heterozygous known and novel BCS1L mutations, which affect the assembly of mitochondrial respiratory chain complex III, in both cases. The pathogenicity of these mutations has been established in several mitochondria-related functional analyses in this study. Mitochondrial diseases with isolated renal symptoms are uncommon; however, this study indicates that mitochondrial respiratory chain complex III deficiency due to BCS1L mutations cause Fanconi syndrome with developmental disability as the primary indications.

A patient with compound heterozygosity of SMPD4: Another example of utility of exome-based copy number analysis in autosomal recessive disorders.

For the efficient diagnosis of rare and undiagnosed diseases, the parallel detection of copy number variants (CNVs) and single nucleotide variants using exome analysis is required. Recently, our group reported the usefulness of a program called EXCAVATOR2, which screens for CNVs from aligned exome data in bam format. This method is expected to contribute to the identification of structural variants and to improve the diagnosis rate, especially for the diagnosis of autosomal recessive disease, when a conventional exome analysis identifies a pathogenic variant in one allele but not the other. Here we report a 2-year-old Japanese boy with an undiagnosed disease. He had severe neonatal asphyxia, severe intellectual disability, intractable seizures, cerebellar and brainstem hypoplasia and dysmorphic features including a prominent supraorbital ridge, thin upper lip, and prominent antihelix. An exome analysis reinforced with a copy number analysis using the EXCAVATOR2 method revealed that the patient had a hemizygous variant in chr2(GRCh37):g.130925108G>A, NM_017951.4 c.832C>T, p.(Arg278*) in SMPD4 that was derived from his father and a deletion of SMPD4 derived from his mother. The presence of the deletion spanning SMPD4 was confirmed by short-read and long-read whole-genome sequencing. The successful diagnosis of this reported patient demonstrates the diagnostic utility of EXCAVATOR2 and overcomes the weakness of exome analysis for the detection of autosomal recessive diseases in nonconsanguineous families, significantly impacting genetic counseling for family planning.

The p.Thr395Met missense variant of NFIA found in a patient with intellectual disability is a defective variant.

Nuclear factor one A (NFIA) is a transcription factor that regulates the development of the central nervous system. Haploinsufficiency of the NFIA gene causes NFIA-related disorder, which includes brain abnormalities and intellectual disability, with or without urinary tract defects. Intragenic deletions, nonsense variants, frameshift variants, and missense variants in one allele of the NFIA gene have been reported to cause various neurological and urogenital symptoms. Here we report a 10-year-old male patient with developmental delay, coarctation of the aorta, and distinctive facial features. Exome analysis identified a rare de novo heterozygous missense variant p.Thr395Met in NFIA. We employed zebrafish as a model organism in our NFIA analysis and found that nfia-/- zebrafish initially showed a loss of commissural axons in the brain, and eventually underwent growth retardation resulting in premature death. Impairment of the commissural neurons in nfia-/- zebrafish embryos could be restored by the expression of wild-type human NFIA protein, but not of mutant human protein harboring the p.Thr395Met substitution, indicating that this variant affects the function of NFIA protein. Taken together, we suggest that the p.Thr395Met allele in the NFIA gene is relevant to the pathogenesis of NFIA-related disorder.

The novel and recurrent variants in exon 31 of CREBBP in Japanese patients with Menke-Hennekam syndrome.

Menke-Hennekam syndrome-1 (MKHK1) is a congenital disorder caused by the heterozygous variants in exon 30 or 31 of CREBBP (CREB binding protein) gene mapped on 16p13.3. It is characterized by psychomotor delay, variable impairment of intellectual disability (ID), feeding difficulty, autistic behavior, hearing impairment, short stature, microcephaly, and facial dysmorphisms. The CREBBP loss-of-function variants cause Rubinstein-Taybi syndrome-1 (RSTS1). The function of CREBBP leading to MKHK1 has not been clarified so far, and the phenotype of MKHK1 significantly differs from that of RSTS1. We examined six patients with de novo pathogenic variants affecting the last exon of CREBBP, and they shared the clinical features of MKHK1. This study revealed that one frameshift and three nonsense variants of CREBBP cause MKHK1, and inferred that the nonsense variants of the last exon could further help in the elucidation of the etiology of MKHK1.

CDK19-related disorder results from both loss-of-function and gain-of-function de novo missense variants.

PURPOSE: To expand the recent description of a new neurodevelopmental syndrome related to alterations in CDK19. METHODS: Individuals were identified through international collaboration. Functional studies included autophosphorylation assays for CDK19 Gly28Arg and Tyr32His variants and in vivo zebrafish assays of the CDK19G28R and CDK19Y32H. RESULTS: We describe 11 unrelated individuals (age range: 9 months to 14 years) with de novo missense variants mapped to the kinase domain of CDK19, including two recurrent changes at residues Tyr32 and Gly28. In vitro autophosphorylation and substrate phosphorylation assays revealed that kinase activity of protein was lower for p.Gly28Arg and higher for p.Tyr32His substitutions compared with that of the wild-type protein. Injection of CDK19 messenger RNA (mRNA) with either the Tyr32His or the Gly28Arg variants using in vivo zebrafish model significantly increased fraction of embryos with morphological abnormalities. Overall, the phenotype of the now 14 individuals with CDK19-related disorder includes universal developmental delay and facial dysmorphism, hypotonia (79%), seizures (64%), ophthalmologic anomalies (64%), and autism/autistic traits (56%). CONCLUSION: CDK19 de novo missense variants are responsible for a novel neurodevelopmental disorder. Both kinase assay and zebrafish experiments showed that the pathogenetic mechanism may be more diverse than previously thought.

Haploinsufficiency of PRR12 causes a spectrum of neurodevelopmental, eye, and multisystem abnormalities.

PURPOSE: Proline Rich 12 (PRR12) is a gene of unknown function with suspected DNA-binding activity, expressed in developing mice and human brains. Predicted loss-of-function variants in this gene are extremely rare, indicating high intolerance of haploinsufficiency. METHODS: Three individuals with intellectual disability and iris anomalies and truncating de novo PRR12 variants were described previously. We add 21 individuals with similar PRR12 variants identified via matchmaking platforms, bringing the total number to 24. RESULTS: We observed 12 frameshift, 6 nonsense, 1 splice-site, and 2 missense variants and one patient with a gross deletion involving PRR12. Three individuals had additional genetic findings, possibly confounding the phenotype. All patients had developmental impairment. Variable structural eye defects were observed in 12/24 individuals (50%) including anophthalmia, microphthalmia, colobomas, optic nerve and iris abnormalities. Additional common features included hypotonia (61%), heart defects (52%), growth failure (54%), and kidney anomalies (35%). PrediXcan analysis showed that phecodes most strongly associated with reduced predicted PRR12 expression were enriched for eye- (7/30) and kidney- (4/30) phenotypes, such as wet macular degeneration and chronic kidney disease. CONCLUSION: These findings support PRR12 haploinsufficiency as a cause for a novel disorder with a wide clinical spectrum marked chiefly by neurodevelopmental and eye abnormalities.

Coloboma may be a shared feature in a spectrum of disorders caused by mutations in the WDR37-PACS1-PACS2 axis.

We report a male adult with early infantile-onset epilepsy, facial dysmorphism, and iridal and choroidal coloboma who had a de novo heterozygous mutation in PACS2, that is, c.625G > A p.(Glu209Lys). This specific mutation was previously reported in a patient with PACS2-related disorder (early infantile epileptic encephalopathy 66). De novo heterozygous mutations in WDR37 have been shown to cause a novel human disorder, neurooculocardiogenitourinary syndrome (NOCGUS syndrome) (OMIM #618652), characterized by intellectual disability, facial dysmorphism, and coloboma. According to large-scale interactome data, WDR37 interacts most strongly, by far, with PACS1 and PACS2. Clinically, coloboma has been described as a feature in a WDR37-related disorder and a PACS1-related disorder (Schuurs-Hoeijmakers syndrome), but not in a PACS2-related disorder. Our review of the phenotypes of three human disorders caused by WDR37, PACS1, and PACS2 mutations showed a significant overlap of epilepsy, intellectual disability, cerebellar atrophy, and facial features. The present observation of coloboma as a shared feature among these three disorders suggests that this group of genes may be involved in ocular development. We propose that dysregulation of the WDR37-PACS1-PACS2 axis results in a spectrum that is recognizable by intellectual disability, distinctive facial features, and coloboma.

Biallelic loss of OTUD7A causes severe muscular hypotonia, intellectual disability, and seizures.

The heterozygous deletion of 15q13.3 is a recurrently observed microdeletion syndrome associated with a relatively mild phenotype including learning disability and language impairment. In contrast, the homozygous deletion of 15q13.3 is extremely rare and is associated with a much severer phenotype that includes epileptic encephalopathy, profound intellectual disability, and hypotonia. Which of the genes within the deleted interval is responsible for the more severe features when biallelically deleted is currently unknown. Here, we report a patient with profound hypotonia, severe intellectual disability, and seizures who had biallelic loss-of-function variants in OTUD7A: a 15q13.3 deletion including the OTUD7A locus, and a frameshift OTUD7A variant c.1125del, p.(Glu375Aspfs*11). Unexpectedly, both aberrations occurred de novo. Our experiment using Caenorhabditis elegans showed that worms carrying a corresponding homozygous variant in the homolog OTUB-2 exhibited weakened muscle contraction suggestive of aberrant neuromuscular transmission. We concluded that the biallelic complete loss of OTUD7A in humans represents a presumably new autosomal recessive disorder characterized by profound hypotonia, severe intellectual disability, and seizures.

Establishing intellectual disability as the key feature of patients with biallelic RNPC3 variants.

Some mammalian genes contain both major and minor introns, the splicing of which require distinctive major and minor spliceosomes, respectively; these genes are referred to as minor intron containing-genes. RNPC3 (RNA-binding domain-containing protein 3) is one of the proteins that are unique to the minor spliceosome U11/U12 di-snRNP. Only two families with biallelic pathogenic variants in the RNPC3 gene encoding the protein have been reported so far, and the affected members in both families had proportional short stature. While the affected members of the originally identified family did not have intellectual disability, the patients from the other family exhibited intellectual disability. Here, we report on a patient with severe primordial microcephalic dwarfism and intellectual disability who carried compound heterozygous variants in RNPC3 (NM_017619.3): c.261dup, p.Leu88Thrfs*11 and c.1228T>G, p.Phe410Val. The single nucleotide substitution c.1228T>G had a very high predictive score for pathogenicity: the p.Phe410 residue is highly conserved down to fish. Based on ACMG (American College of Medical Genetics and Genomics) guideline, this non-synonymous variant was scored as likely pathogenic. This documentation of yet another patient with biallelic RNPC3 variants exhibiting intellectual disability lends further support to the notion that intellectual disability is a key feature of the spectrum of RNPC3-related disorders.

Clinical spectrum of individuals with de novo EBF3 variants or deletions.

Hypotonia, ataxia and delayed development syndrome (HADDS) (MIM#617330) is a neurodevelopmental disorder caused by heterozygous pathogenic variants in EBF3 (MIM; 607,407), which is located on chromosome 10q26, and was first reported in 2017. To date, missense, nonsense and frameshift variants have been reported as causes of HADDS, and EBF3 pathogenic variants have been predicted to result in nonsense-mediated mRNA decay and haploinsufficiency. It was also reported that total deletion of EBF3 associated with a 10q26.3 microdeletion also causes HADDS symptoms, supporting the concept that HADDS results from haploinsufficiency of EBF3. Here, we report eight unrelated individuals with heterozygous pathogenic variants of EBF3 or haploinsufficiency of EBF3 due to 10q26 deletion, who exhibit clinical findings including craniofacial features of HADDS. In a detailed examination of clinical manifestations in this study, revealed that neurogenic bladder was diagnosed in infancy (the median 6.5 months), was more frequent than previously reported, and required cystostomy in all but one case. For psychomotor delay, it was also found that their motor/skills values were significantly lower than their cognition/adaptation values (p = 0.0016; paired t-test). Therefore, that HADDS is a recognizable syndrome that shares its characteristic facial features, and that neurogenic bladder diagnosed in infancy and psychomotor delay with marked delay in motor/skills are noteworthy findings in the diagnosis and management of individuals with HADDS.

Fork-shaped mandibular incisors as a novel phenotype of LRP5-associated disorder.

The LRP5 gene encodes a Wnt signaling receptor to which Wnt binds directly. In humans, pathogenic monoallelic variants in LRP5 have been associated with increased bone density and exudative vitreoretinopathy. In mice, LRP5 plays a role in tooth development, including periodontal tissue stability and cementum formation. Here, we report a 14-year-old patient with a de novo non-synonymous variant, p.(Val1245Met), in LRP5 who exhibited mildly reduced bone density and mild exudative vitreoretinopathy together with a previously unreported phenotype consisting of dental abnormalities that included fork-like small incisors with short roots and an anterior open bite, molars with a single root, and severe taurodontism. In that exudative vitreoretinopathy has been reported to be associated with heterozygous loss-of-function variants of LRP5 and that our patient reported here with the p.(Val1245Met) variant had mild exudative vitreoretinopathy, the variant can be considered as an incomplete loss-of-function variant. Alternatively, the p.(Val1245Met) variant can be considered as exerting a dominant-negative effect, as no patients with truncating LRP5 variants and exudative vitreoretinopathy have been reported to exhibit dental anomalies. The documentation of dental anomalies in the presently reported patient strongly supports the notion that LRP5 plays a critical role in odontogenesis in humans, similar to its role in mice.

Recurrent NFIA K125E substitution represents a loss-of-function allele: Sensitive in vitro and in vivo assays for nontruncating alleles.

Nuclear factor I A (NFIA) is a transcription factor that belongs to the NFI family. Truncating variants or intragenic deletion of the NFIA gene are known to cause the human neurodevelopmental disorder known as NFIA-related disorder, but no patient heterozygous for a missense mutation has been reported. Here, we document two unrelated patients with typical phenotypic features of the NFIA-related disorder who shared a missense variant p.Lys125Glu (K125E) in the NFIA gene. Patient 1 was a 6-year-old female with global developmental delay, corpus callosum anomaly, macrocephaly, and dysmorphic facial features. Patient 2 was a 14-month-old male with corpus callosum anomaly and macrocephaly. By using Drosophila and zebrafish models, we functionally evaluated the effect of the K125E substitution. Ectopic expression of wild-type human NFIA in Drosophila caused developmental defects such as eye malformation and premature death, while that of human NFIA K125E variant allele did not. nfia-deficient zebrafish embryos showed defects of midline-crossing axons in the midbrain/hindbrain boundary. This impairment of commissural neurons was rescued by expression of wild-type human NFIA, but not by that of mutant variant harboring K125E substitution. In accordance with these in vivo functional analyses, we showed that the K125E mutation impaired the transcriptional regulation of HES1 promoter in cultured cells. Taken together, we concluded that the K125E variant in the NFIA gene is a loss-of-function mutation.

Variants in KIF2A cause broad clinical presentation; the computational structural analysis of a novel variant in a patient with a cortical dysplasia, complex, with other brain malformations 3.

Cortical dysplasia, complex, with other brain malformations 3 (CDCBM3) is a rare autosomal dominant syndrome caused by Kinesin family Member 2A (KIF2A) gene mutation. Patients with CDCBM3 exhibit posterior dominant agyria/pachygyria with severe motor dysfunction. Here, we report an 8-year-old boy with CDCBM3 showing a typical, but relatively mild, clinical presentation of CDCBM3 features. Whole-exome sequencing identified a heterozygous mutation of NM_001098511.2:c.1298C>A [p.(Ser433Tyr)]. To our knowledge, the mutation has never been reported previously. The variant was located distal to the nucleotide binding domain (NBD), in which previously-reported variants in CDCBM3 patients have been located. The computational structural analysis showed the p.433 forms the pocket with NBD. Variants in KIF2A have been reported in the NBD for CDCBM3, in the kinesin motor 3 domain, but not in the NBD in epilepsy, and outside of the kinesin motor domain in autism spectrum syndrome, respectively. Our patient has a variant, that is not in the NBD but at the pocket with the NBD, resulting in a clinical features of CDCBM3 with mild symptoms. The clinical findings of patients with KIF2A variants appear restricted to the central nervous system and facial anomalies. We can call this spectrum "KIF2A syndrome" with variable severity.

Biallelic Mutations in the LSR Gene Cause a Novel Type of Infantile Intrahepatic Cholestasis.

We identified biallelic pathogenic mutations in the Lipolysis-stimulated lipoprotein receptor (LSR) gene in a patient with infantile intrahepatic cholestasis. We established that mutations in the LSR gene, which encodes a protein which is critical for the formation of tricellular tight junctions in the liver, are a novel cause of pediatric cholestasis.

Parallel detection of single nucleotide variants and copy number variants with exome analysis: Validation in a cohort of 700 undiagnosed patients.

Copy number variants (CNVs) are significant causes of rare and undiagnosed diseases. Parallel detection of single nucleotide variants (SNVs) and CNVs with exome analysis, if feasible, would shorten the diagnostic closure in a timely manner. We validated such "parallel" approach through a cohort study of 791 undiagnosed patients. In addition to routine exome analysis, we applied an innovative algorithm EXCAVATOR2 which enhances sensitivity by paradoxically exploiting read depth data that covers nonexonic regions where baits were not originally intended to hybridize. About 48 patients had copy number variations, 42 deletions, and 6 duplications with a resolution of 0.51-14.7 mega base pairs. Importantly from a clinical standpoint, we identified three patients with "dual diagnosis" due to concurrent pathogenic CNV and SNV. We suggest "hitting two birds with one stone" approach to exome data is an efficient strategy in deciphering undiagnosed patients and may well be considered as a first-tier genetic test.

Consecutive medical exome analysis at a tertiary center: Diagnostic and health-economic outcomes.

The utility of whole exome analysis has been extensively demonstrated in research settings, but its clinical utility as a first-tier genetic test has not been well documented from diagnostic and health economic standpoints in real-life clinical settings. We performed medical exome analyses focusing on a clinically interpretable portion of the genome (4,813 genes) as a first-tier genetic test for 360 consecutive patients visiting a genetics clinic at a tertiary children's hospital in Japan, over a 3-year period. Bioinformatics analyses were conducted using standard software. A molecular diagnosis was made in 171 patients involving a total of 107 causative genes. Among these 107 causative genes, 57 genes were classified as genes with potential organ-specific interventions and management strategies. Clinically relevant results were obtained in 26% of the total cohort and 54% of the patients with a definitive molecular diagnosis. Performing the medical exome analysis at the time of the initial visit to the tertiary center, rather than after visits to pertinent specialists, brain MRI examination, and G-banded chromosome testing, would have reduced the financial cost by 197 euros according to retrospective calculation under multiple assumption. The present study demonstrated a high diagnostic yield (47.5%) for singleton medical exome analysis as a first-tier test in a real-life setting. Medical exome analysis yielded clinically relevant information in a quarter of the total patient cohort. The application of genomic testing during the initial visit to a tertiary medical center could be a rational approach to the diagnosis of patients with suspected genetic disorders.

Protein elongation variant of PUF60: Milder phenotypic end of the Verheij syndrome.

The PUF60 gene encodes a ubiquitously expressed essential splicing factor that is recruited to the U2snRNA complex. The complex binds to the 3' splice site of exons in specific target genes and regulates the inclusion or exclusion of such exons. Recently, pathogenic variants of PUF60 have been shown to cause a relatively specific and potentially recognizable pattern of malformation referred to as Verheij syndrome. Here, we report a 12-year-old female patient with a de novo mutation in PUF60 whose phenotype was representative of the milder end of the phenotypic spectrum of Verheij syndrome; the de novo mutation was a frameshift mutation p.(Ser558Cysfs*21) that resulted in the addition of 21 extra amino acids at the carboxy end of the protein. Among the frequent features of Verheij syndrome, the patient exhibited coloboma, cervical spinal segmentation defects, and borderline intellectual functioning, but lacked cardiac abnormalities, deafness, and urogenital abnormalities. The results of RNA analysis using peripheral blood showed the escape of the mutant allele from nonsense-mediated mRNA decay, possibly accounting for the mild phenotype in the presently reported patient. Based on our clinical observations, we inferred that two embryologic processes, closure of the ocular plate and cervical spinal segmentation, are particularly susceptible to deficient PUF60-mediated splicing regulation, compared with other embryogenetic processes leading to the central nervous system, heart, ear, and kidney.

Shortfall of exome analysis for diagnosis of Shwachman-Diamond syndrome: Mismapping due to the pseudogene SBDSP1.

Shwachman-Diamond syndrome characterized by metaphyseal dysplasia, pancreatic insufficiency, and pancytopenia is caused by biallelic mutations in SBDS. Gene conversion between SBDS and its pseudogene SBDSP1 is the major cause. Here, we report two unrelated patients with Shwachman-Diamond syndrome who were shown to be compound heterozygotes for relatively frequent pathogenic alleles (the 258+2T>C allele and another allele composed of 183-184TA>CT and 201A>G) using an established polymerase chain reaction sequencing assay with SBDS-specific primers. Exome analysis of the patients showed discrepant results: 258+2T>C with variant allele frequency around 0.85, and no variants detected for the 183-184TA>CT allele. Parental exome analysis of the two families further supported this notion. Confronted with two patients with an unexpected segregation pattern, we performed a transcriptome analysis of peripheral blood-derived mRNA to demonstrate that the results were compatible with those obtained using SBDS-specific PCR primers. Both alleles could be accounted for by gene conversion events. The diagnostic discrepancy can be accounted for by a decreased efficiency in the computational mapping of the reads with 183-184TA>CT and 201A>G to the reference sequence of the SBDS locus during exome analysis. This report highlights the pitfall of exome analysis for genes with pseudogenes, such as SBDS and the alternative use of RNA-seq is recommended to circumvent this problem.