Complete androgen insensitivity syndrome caused by a deep intronic pseudoexon-activating mutation in the androgen receptor gene

Advances in the understanding of androgen receptor structure and function and in the development of next-generation AR-targeted therapeutics

Androgen receptor (AR) and its ligand androgens are important for development and physiology of various tissues. AR and its ligands also play critical role in the development of various diseases, making it a valuable therapeutic target. AR ligands, both agonists and antagonists, are being widely used to treat pathological conditions, including prostate cancer and hypogonadism. Despite AR being studied widely over the last five decades, the last decade has seen striking advances in the knowledge on AR and discoveries that have the potential to translate to the clinic. This review provides an overview of the advances in AR biology, AR molecular mechanisms of action, and next generation molecules that are currently in development. Several of the areas described in the review are just unraveling and the next decade will bring more clarity on these developments that will put AR at the forefront of both basic biology and drug development.

LINE1-mediated epigenetic repression of androgen receptor transcription causes androgen insensitivity syndrome

Androgen insensitivity syndrome (AIS) is a difference of sex development (DSD) characterized by different degrees of undervirilization in individuals with a 46,XY karyotype despite normal to high gonadal testosterone production. Classically, AIS is explained by hemizygous mutations in the X-chromosomal androgen receptor (AR) gene. Nevertheless, the majority of individuals with clinically diagnosed AIS do not carry an AR gene mutation. Here, we present a patient with a 46,XY karyotype, born with undervirilized genitalia, age-appropriate testosterone levels and no uterus, characteristic for AIS. Diagnostic whole exome sequencing (WES) showed a maternally inherited LINE1 (L1) retrotransposon insertion in the 5' untranslated region (5'UTR) of the AR gene. Long-read nanopore sequencing confirmed this as an insertion of a truncated L1 element of ≈ 2.7 kb and showed an increased DNA methylation at the L1 insertion site in patient-derived genital skin fibroblasts (GSFs) compared to healthy controls. The insertion coincided with reduced AR transcript and protein levels in patient-derived GSFs confirming the clinical diagnosis AIS. Our results underline the relevance of retrotransposons in human disease, and expand the growing list of human diseases associated with them.

Low androgen signaling rescues genome integrity with innate immune response by reducing fertility in humans

Development of the gonads under complex androgen regulation is critical for germ cells specification. In this work we addressed the relationship between androgens and genomic integrity determining human fertility. We used different study groups: individuals with Differences of Sex Development (DSD), including Complete Androgen Insensitivity Syndrome (CAIS) due to mutated androgen receptor (AR), and men with idiopathic nonobstructive azoospermia. Both showed genome integrity status influenced by androgen signaling via innate immune response activation in blood and gonads. Whole proteome analysis connected low AR to interleukin-specific gene expression, while compromised genome stability and tumorigenesis were also supported by interferons. AR expression was associated with predominant DNA damage phenotype, that eliminated AR-positive Sertoli cells as the degeneration of gonads increased. Low AR contributed to resistance from the inhibition of DNA repair in primary leukocytes. Downregulation of androgen promoted apoptosis and specific innate immune response with higher susceptibility in cells carrying genomic instability.

The genetic spectrum of a Chinese series of patients with 46, XY disorders of the sex development

PURPOSE

The etiology of 46, XY disorders of sex development (46, XY DSD) is complex, and studies have shown that different series of patients with 46, XY DSD has different genetic spectrum. In this study, we aimed to investigate the underlying genetic etiology in a Chinese series of patients with 46, XY DSD by whole exome sequencing (WES).

METHODS

Seventy patients with 46, XY DSD were enrolled from the Peking Union Medical College Hospital (Beijing, China). The detailed clinical characteristics were evaluated, and peripheral blood was collected for WES to find the patients' rare variants (RVs) of genes related to 46, XY DSD. The clinical significance of the RVs was annotated according to American College of Medical Genetics and Genomics (ACMG) guidelines.

RESULTS

A total of 57 RVs from nine genes were identified in 56 patients with 46, XY DSD, which include 21 novel RVs and 36 recurrent RVs. Based on the American ACMG guidelines, 43 variants were classified as pathogenic(P) or likely pathogenic (LP) variants and 14 variants were defined as variants of uncertain significance (VUS). P or LP variants were identified in 64.3% (45/70) patients of the series. Thirty-nine, 14, and 4 RVs were involved in the process of androgen synthesis and action, testicular determination and developmental process, and syndromic 46, XY DSD, respectively. The top three genes most frequently affected to cause 46, XY DSD were AR, SRD5A2, and NR5A1. Seven patients were found harboring RVs of the 46, XY DSD pathogenic genes identified in recent years, namely DHX37 in four patients, MYRF in two patients, and PPP2R3C in one patient.

CONCLUSION

We identified 21 novel RVs of nine genes, which extended the genetic spectrum of 46, XY DSD pathogenic variants. Our study showed that 60% of the patients were caused by AR, SRD5A2 or NR5A1 P/LP variants. Therefore, polymerase chain reaction (PCR) amplification and Sanger sequencing of these three genes could be performed first to identify the pathogeny of the patients. For those patients whose pathogenic variants had not been found, whole-exome sequencing could be helpful in determining the etiology.

Androgen Insensitivity Syndrome DUE to Non-Coding Variation in the Androgen Receptor Gene: Review of the Literature and Case Report of a Patient with Mosaic c.-547C>T Variant

Sexual development (SD) is a complex process with strict spatiotemporal regulation of gene expression. Despite advancements in molecular diagnostics, disorders of sexual development (DSD) have a diagnostic rate of ~50%. Androgen insensitivity syndrome (AIS) represents the most common form of 46,XY DSD, with a spectrum of defects in androgen action. Considering the importance of very strict regulation of the SD, it is reasonable to assume that the genetic cause for proportion of the DSD lies in the non-coding part of the genome that regulates proper gene functioning. Here we present a patient with partial AIS (PAIS) due to a mosaic de novo c.-547C>T pathogenic variant in the 5'UTR of androgen receptor (AR) gene. The same mutation was previously described as inherited, in two unrelated patients with complete AIS (CAIS). Thus, our case further confirms the previous findings that variable gene expressivity could be attributed to mosaicism. Mutations in 5'UTR could create new upstream open reading frames (uORFs) or could disrupt the existing one. A recent systematic genome-wide study identified AR as a member of a subset of genes where modifications of uORFs represents an important disease mechanism. Only a small number of studies are reporting non-coding mutations in the AR gene and our case emphasizes the importance of molecular testing of the entire AR locus in AIS patients. The introduction of new methods for comprehensive molecular testing in routine genetic diagnosis, accompanied with new tools for in sillico analysis could improve the genetic diagnosis of AIS, and DSD in general.

Case report: Identification of a frameshift mutation in GC enrichment and the GCC repeat region of the androgen insensitivity receptor (AR) gene in a patient with complete androgen insensitivity syndrome by whole-exome sequencing (WES) combined with specific PCR and deep sequencing

Background: Androgen insensitivity syndrome (AIS) is an X-linked recessive hereditary disease caused due to a reduced or absent function of the androgen receptor (AR) protein encoded by the AR gene (OMIM-Gene# 313,700). Genetic testing is important in the diagnosis, clinical management, and prevention of AIS (MIM# 300,068). The AR (HGNC: 644) pathogenic variant detection rate ranges from 65% to 95% for patients with complete AIS (CAIS) and 40%-45% for patients with partial androgen insensitivity syndrome (PAIS). Identification of a pathogenic mutation in the AR confirms the diagnosis of AIS, especially in the milder forms that may have a phenotypic overlap with other disorders of sex development. Improvement of the molecular diagnostic rate of AIS is urgently required in clinical practice. We reported the results of the molecular diagnosis of a patient with CAIS who failed previously in either the traditional Sanger sequencing or next-generation sequencing (NGS). Using whole-exome sequencing (WES) combined with a special polymerase chain reaction (PCR) and deep sequencing, we successfully identified a pathogenic variant, a hemizygous mutation (c.1395-1396insGA), in the GC-enriched and unstable GCC repeat regions of the AR gene of the proband. Conclusion: The results may be advantageous for the improvement of the detection rate of AIS, as well as other inherited disorders whose disease-causing genes contain GC-enriched and unstable GCC repeat regions.

Modulation of RNA Splicing by Oligonucleotides: Mechanisms of Action and Therapeutic Implications

Dysregulation of RNA splicing causes many diseases and disorders. Several therapeutic approaches have been developed to correct aberrant alternative splicing events for the treatment of cancers and hereditary diseases, including gene therapy and redirecting splicing, using small molecules or splice switching oligonucleotides (SSO). Significant advances in the chemistry and pharmacology of nucleic acid have led to the development of clinically approved SSO drugs for the treatment of spinal muscular dystrophy and Duchenne muscular dystrophy (DMD). In this review, we discuss the mechanisms of SSO action with emphasis on "less common" approaches to modulate alternative splicing, including bipartite and bifunctional SSO, oligonucleotide decoys for splice factors and SSO-mediated mRNA degradation via AS-NMD and NGD pathways. We briefly discuss the current progress and future perspectives of SSO therapy for rare and ultrarare diseases.

Analysis of Pathogenic Pseudoexons Reveals Novel Mechanisms Driving Cryptic Splicing

Understanding pre-mRNA splicing is crucial to accurately diagnosing and treating genetic diseases. However, mutations that alter splicing can exert highly diverse effects. Of all the known types of splicing mutations, perhaps the rarest and most difficult to predict are those that activate pseudoexons, sometimes also called cryptic exons. Unlike other splicing mutations that either destroy or redirect existing splice events, pseudoexon mutations appear to create entirely new exons within introns. Since exon definition in vertebrates requires coordinated arrangements of numerous RNA motifs, one might expect that pseudoexons would only arise when rearrangements of intronic DNA create novel exons by chance. Surprisingly, although such mutations do occur, a far more common cause of pseudoexons is deep-intronic single nucleotide variants, raising the question of why these latent exon-like tracts near the mutation sites have not already been purged from the genome by the evolutionary advantage of more efficient splicing. Possible answers may lie in deep intronic splicing processes such as recursive splicing or poison exon splicing. Because these processes utilize intronic motifs that benignly engage with the spliceosome, the regions involved may be more susceptible to exonization than other intronic regions would be. We speculated that a comprehensive study of reported pseudoexons might detect alignments with known deep intronic splice sites and could also permit the characterisation of novel pseudoexon categories. In this report, we present and analyse a catalogue of over 400 published pseudoexon splice events. In addition to confirming prior observations of the most common pseudoexon mutation types, the size of this catalogue also enabled us to suggest new categories for some of the rarer types of pseudoexon mutation. By comparing our catalogue against published datasets of non-canonical splice events, we also found that 15.7% of pseudoexons exhibit some splicing activity at one or both of their splice sites in non-mutant cells. Importantly, this included seven examples of experimentally confirmed recursive splice sites, confirming for the first time a long-suspected link between these two splicing phenomena. These findings have the potential to improve the fidelity of genetic diagnostics and reveal new targets for splice-modulating therapies.

Genetic Variation and the Distribution of Variant Types in the Horse

Genetic variation is a key contributor to health and disease. Understanding the link between an individual's genotype and the corresponding phenotype is a major goal of medical genetics. Whole genome sequencing (WGS) within and across populations enables highly efficient variant discovery and elucidation of the molecular nature of virtually all genetic variation. Here, we report the largest catalog of genetic variation for the horse, a species of importance as a model for human athletic and performance related traits, using WGS of 534 horses. We show the extent of agreement between two commonly used variant callers. In data from ten target breeds that represent major breed clusters in the domestic horse, we demonstrate the distribution of variants, their allele frequencies across breeds, and identify variants that are unique to a single breed. We investigate variants with no homozygotes that may be potential embryonic lethal variants, as well as variants present in all individuals that likely represent regions of the genome with errors, poor annotation or where the reference genome carries a variant. Finally, we show regions of the genome that have higher or lower levels of genetic variation compared to the genome average. This catalog can be used for variant prioritization for important equine diseases and traits, and to provide key information about regions of the genome where the assembly and/or annotation need to be improved.

Primary Amenorrhea Due to Anatomical Abnormalities of the Reproductive Tract: Molecular Insight

Congenital anomalies of the female reproductive tract that present with primary amenorrhea involve Müllerian aplasia, also known as Mayer-Rokitansky-Küster-Hauser syndrome (MRKHS), and cervical and vaginal anomalies that completely obstruct the reproductive tract. Karyotype abnormalities do not exclude the diagnosis of MRKHS. Familial cases of Müllerian anomalies and associated malformations of the urinary and skeletal systems strongly suggest a complex genetic etiology, but so far, the molecular mechanism in the vast majority of cases remains unknown. Primary amenorrhea may also be the first presentation of complete androgen insensitivity syndrome, steroid 5α-reductase type 2 deficiency, 17β-hydroxysteroid dehydrogenase type 3 deficiency, and Leydig cells hypoplasia type 1; therefore, these disorders should be considered in the differential diagnosis of the congenital absence of the uterus and vagina. The molecular diagnosis in the majority of these cases can be established.

Clinical, Biochemical, and Molecular Characterization of Indian Children with Clinically Suspected Androgen Insensitivity Syndrome

This study describes the clinical, biochemical, and molecular characteristics of Indian children with 46,XY DSD and suspected androgen insensitivity syndrome (AIS). Fifty children (median age 3.0 years, range 0-16.5 years) with 46,XY DSD and a suspected diagnosis of AIS were enrolled. Sanger sequencing was performed to identify pathogenic variants in the androgen receptor (AR) gene and to study genotype-phenotype correlations. All 5 (100%) patients with CAIS and 14/45 (31%) patients with PAIS had pathogenic/likely pathogenic variants in the AR gene (overall, 14 different variants in 19 patients; 38.8%). There was no significant difference in clinical (cryptorchidism, hypospadias, or external masculinizing score) or biochemical parameters (gonadotropins and testosterone) between patients with or without pathogenic variants. However, patients with AIS were more likely to have a positive family history, be assigned female gender at birth, and present with gynaecomastia at puberty. Three novel pathogenic/likely pathogenic variants, including one splice donor site variant c.2318+1G>A, one frameshift variant p.H790Lfs*40, and one missense variant p.G821E, were identified in 3 patients with CAIS. The missense variant p.G821E was predicted as deleterious, damaging, disease-causing, and likely functionally inactive by in silico analysis and protein modelling study. Two previously not reported pathogenic/likely pathogenic variants, including p.R386H and p.G396R, were identified in patients with PAIS. This study contributes in expanding the spectrum of pathogenic variants in the AR gene in patients with AIS. Only 31% patients with a provisional diagnosis of PAIS had pathogenic variants in the AR gene, suggesting other possible mechanisms or candidate genes may be responsible for such a phenotypic presentation.

Targeting the Non-Coding Genome for the Diagnosis of Disorders of Sex Development

Disorders of sex development (DSD) are a complex group of conditions with highly variable clinical phenotypes, most often caused by failure of gonadal development. DSD are estimated to occur in around 1.7% of all live births. Whilst the understanding of genes involved in gonad development has increased exponentially, approximately 50% of patients with a DSD remain without a genetic diagnosis, possibly implicating non-coding genomic regions instead. Here, we review how variants in the non-coding genome of DSD patients can be identified using techniques such as array comparative genomic hybridization (CGH) to detect copy number variants (CNVs), and more recently, whole genome sequencing (WGS). Once a CNV in a patient's non-coding genome is identified, putative regulatory elements such as enhancers need to be determined within these vast genomic regions. We will review the available online tools and databases that can be used to refine regions with potential enhancer activity based on chromosomal accessibility, histone modifications, transcription factor binding site analysis, chromatin conformation, and disease association. We will also review the current in vitro and in vivo techniques available to demonstrate the functionality of the identified enhancers. The review concludes with a clinical update on the enhancers linked to DSD.

[A deep intronic mutation in AR gene causing androgen insensitivity syndrome: difficulties of diagnostics]

Partial androgen resistance syndrome (PAIS) is the most difficult form of disorders/differences of sex development 46,XY (DSD 46,XY) for choosing of patient management. To date, there are no clear biochemical criteria, especially before puberty, that allow differentiating PAIS from other PAIS-like forms of DSD 46, XY, and genetic verification of the partial form of AIS plays an important role. Meanwhile, according to the literature, mutations in the coding region of AR gene have not been identified in more than 50% of patients with suspected AIS. We performed an extensive analysis of the AR gene in a patient with clinical and laboratory signs of AIS and found a deep intron mutation in the AR gene (p. 2450-42G>A). This variant creates an alternative splice acceptor site resulted a disturbance of the AR function. These findings indicate the need for extensive genetic analysis in a cohort of patients with suspected CPA in the absence of mutations in the AR gene using standard methods of genetic diagnosis.

A novel de novo androgen receptor nonsense mutation in a sex-reversed 46,XY infant

An infant with 46,XY karyotype, and unambiguous female phenotype was found to have testes in the inguinal regions. Capillary sequencing of the androgen receptor (AR) gene identified a hemizygous de novo mutation (NM_000044.6:c.1621G > T) in exon 2 resulting in a termination codon p.(Glu541*) at the DNA binding domain (DBD). This novel nonsense mutation adds to the compendium of AR mutations which result in complete androgen insensitivity syndrome (AIS).

Clinical, Hormonal, Genetic, and Molecular Characteristics in Androgen Insensitivity Syndrome in an Asian Indian Cohort from a Single Centre in Western India

The study aimed to analyze clinical and hormonal phenotype,and genotype in patients with genetically proven androgen insensitivity syndrome (AIS) from Western India. Index patients with pathogenic variants in the androgen receptor (AR) gene were identified from a consecutive 46,XY DSD cohort (n = 150) evaluated with clinical exome sequencing, and their genetically-proven affected relatives were also included. In sum, 15 index cases (9 complete AIS [CAIS] and 6 partial AIS [PAIS]) were identified making AIS the second most common (10%) cause of 46,XY DSD, next to 5α-reductase 2 deficiency (n = 26; 17.3%). Most patients presented late in the postpubertal period with primary amenorrhoea in CAIS (89%) and atypical genitalia with gynecomastia in PAIS (71.4%). All CAIS were reared as females and 83.3% of PAIS as males with no gender dysphoria. Four of 6 patients with available testosterone to dihydrotestosterone ratio had a false elevation (>10). Metastatic dysgerminoma was seen in 1 patient in CAIS, while none in the PAIS group had malignancy. Fifteen different (including 6 novel) pathogenic/likely pathogenic variants in AR were found. Nonsense and frameshift variants exclusively led to CAIS phenotype, whereas missense variants led to variable phenotypes. In this largest, monocentric study from the Asian Indian subcontinent, AIS was the second most common cause of 46,XY DSD with similar phenotype but later presentation when compared to cases in the rest of the world. The study reports 6 novel pathogenic variants in AR.

In silico identification of pseudo-exon activation events in personal genome and transcriptome data

Causative mutations for human genetic disorders have mainly been identified in exonic regions that code for amino acid sequences. Recently, however, it has been reported that mutations in deep intronic regions can also cause certain human genetic disorders by creating novel splice sites, leading to pseudo-exon activation. To investigate how frequently pseudo-exon activation events occur in normal individuals, we conducted in silico identification of such events using personal genome data and corresponding high-quality transcriptome data. With rather stringent conditions, on average, 2.6 pseudo-exon activation events per individual were identified. More pseudo-exon activation events were found in 5' donor splice sites than in 3' acceptor splice sites. Although pseudo-exon activation events have sporadically been reported as causative mutations in genetic disorders, it is revealed in this study that such events can be observed in normal individuals at a certain frequency. We estimate that human genomes typically contain on average at least 10 pseudo-exon activation events. The actual number should be higher than this, because we used stringent criteria to identify pseudo-exon activation events. This suggests that it is worth considering the possibility of pseudo-exon activation when searching for causative mutations of genetic disorders if candidate mutations are not identified in coding regions or RNA splice sites.

Genotype-phenotype correlation and molecular heterogeneity in pyruvate kinase deficiency

Pyruvate kinase (PK) deficiency is a rare recessive congenital hemolytic anemia caused by mutations in the PKLR gene. This study reports the molecular features of 257 patients enrolled in the PKD Natural History Study. Of the 127 different pathogenic variants detected, 84 were missense and 43 non-missense, including 20 stop-gain, 11 affecting splicing, five large deletions, four in-frame indels, and three promoter variants. Within the 177 unrelated patients, 35 were homozygous and 142 compound heterozygous (77 for two missense, 48 for one missense and one non-missense, and 17 for two non-missense variants); the two most frequent mutations were p.R510Q in 23% and p.R486W in 9% of mutated alleles. Fifty-five (21%) patients were found to have at least one previously unreported variant with 45 newly described mutations. Patients with two non-missense mutations had lower hemoglobin levels, higher numbers of lifetime transfusions, and higher rates of complications including iron overload, extramedullary hematopoiesis, and pulmonary hypertension. Rare severe complications, including lower extremity ulcerations and hepatic failure, were seen more frequently in patients with non-missense mutations or with missense mutations characterized by severe protein instability. The PKLR genotype did not correlate with the frequency of complications in utero or in the newborn period. With ICCs ranging from 0.4 to 0.61, about the same degree of clinical similarity exists within siblings as it does between siblings, in terms of hemoglobin, total bilirubin, splenectomy status, and cholecystectomy status. Pregnancy outcomes were similar across genotypes in PK deficient women. This report confirms the wide genetic heterogeneity of PK deficiency.

Deep Intronic Mutation in SERPING1 Caused Hereditary Angioedema Through Pseudoexon Activation

PURPOSE

Hereditary angioedema (HAE) is a rare autosomal dominant life-threatening disease characterized by low levels of C1 inhibitor (type I HAE) or normal levels of ineffective C1 inhibitor (type II HAE), typically occurring as a consequence of a SERPING1 mutation. In some cases, a causal mutation remains undetected after using a standard molecular genetic analysis.

RESULTS

Here we show a long methodological way to the final discovery of c.1029 + 384A > G, a novel deep intronic mutation in intron 6 which is responsible for HAE type I in a large family and has not been identified by a conventional diagnostic approach. This mutation results in de novo donor splice site creation and subsequent pseudoexon inclusion, the mechanism firstly described to occur in SERPING1 in this study. We additionally discovered that the proximal part of intron 6 is a region potentially prone to pseudoexon-activating mutations, since natural alternative exons and additional cryptic sites occur therein. Indeed, we confirmed the existence of at least two different alternative exons in this region not described previously.

CONCLUSIONS

In conclusion, our results suggest that detecting aberrant transcripts, which are often low abundant because of nonsense-mediated decay, requires a modified methodological approach. We suggest SERPING1 intron 6 sequencing and/or tailored mRNA analysis to be routinely used in HAE patients with no mutation identified in the coding sequence.

Features of the fetal gonad in androgen synthesis in the postpubertal testis are preserved in complete androgen insensitivity syndrome due to a novel genetic splice site donor variant in androgen receptor gene intron 1

Mutations in the X-linked androgen receptor (AR) gene cause complete androgen insensitivity syndrome (CAIS). CAIS may cause congenital sexual development disorder, which frequently develops into testicular tumors. Here, we describe a novel splice-site intron 1 mutation in AR leading to improper splicing and AR protein absence in CAIS gonads. We characterized a patient's postpubertal gonadal steroidogenic enzyme expression profile. Localization of both CYP11A1 and CYP17A1 enzymes was restricted to both Leydig tumor cells and adjacent to tumor gonadal tissues. Sertoli cells of the CAIS gonad showed abundant HSD17B3 protein, which is an adult Leydig cell marker that enables the conversion of androstenedione to testosterone. Such HSD17B3 expression is typical for fetal-type Sertoli cells in rodents. The postpubertal CAIS gonad of our patient was completely devoid of androgen signaling pathway activity. Plausibly, the postpubertal Leydig cells consisted of two distinct cell populations: postpubertal fetal-type Leydig cells that persisted as androgen-independent cells and immature adult Leydig cells that failed to differentiate. Taken together, in this CAIS postpubertal testis, both Leydig and fetal-type Sertoli cells participated in testosterone production. Our results indicate the importance of molecular analysis as well as the characterization of steroidogenic enzyme profiling in the CAIS patient's gonad.

Testis Development

Production of sperm and androgens is the main function of the testis. This depends on normal development of both testicular somatic cells and germ cells. A genetic program initiated from the Y chromosome gene sex-determining region Y (SRY) directs somatic cell specification to Sertoli cells that orchestrate further development. They first guide fetal germ cell differentiation toward spermatogenic destiny and then take care of the full service to spermatogenic cells during spermatogenesis. The number of Sertoli cells sets the limits of sperm production. Leydig cells secrete androgens that determine masculine development. Testis development does not depend on germ cells; that is, testicular somatic cells also develop in the absence of germ cells, and the testis can produce testosterone normally to induce full masculinization in these men. In contrast, spermatogenic cell development is totally dependent on somatic cells. We herein review germ cell differentiation from primordial germ cells to spermatogonia and development of the supporting somatic cells. Testicular descent to scrota is necessary for normal spermatogenesis, and cryptorchidism is the most common male birth defect. This is a mild form of a disorder of sex differentiation. Multiple genetic reasons for more severe forms of disorders of sex differentiation have been revealed during the last decades, and these are described along with the description of molecular regulation of testis development.

Systematic analysis of dark and camouflaged genes reveals disease-relevant genes hiding in plain sight

BACKGROUND

The human genome contains "dark" gene regions that cannot be adequately assembled or aligned using standard short-read sequencing technologies, preventing researchers from identifying mutations within these gene regions that may be relevant to human disease. Here, we identify regions with few mappable reads that we call dark by depth, and others that have ambiguous alignment, called camouflaged. We assess how well long-read or linked-read technologies resolve these regions.

RESULTS

Based on standard whole-genome Illumina sequencing data, we identify 36,794 dark regions in 6054 gene bodies from pathways important to human health, development, and reproduction. Of these gene bodies, 8.7% are completely dark and 35.2% are ≥ 5% dark. We identify dark regions that are present in protein-coding exons across 748 genes. Linked-read or long-read sequencing technologies from 10x Genomics, PacBio, and Oxford Nanopore Technologies reduce dark protein-coding regions to approximately 50.5%, 35.6%, and 9.6%, respectively. We present an algorithm to resolve most camouflaged regions and apply it to the Alzheimer's Disease Sequencing Project. We rescue a rare ten-nucleotide frameshift deletion in CR1, a top Alzheimer's disease gene, found in disease cases but not in controls.

CONCLUSIONS

While we could not formally assess the association of the CR1 frameshift mutation with Alzheimer's disease due to insufficient sample-size, we believe it merits investigating in a larger cohort. There remain thousands of potentially important genomic regions overlooked by short-read sequencing that are largely resolved by long-read technologies.

Androgen Insensitivity Syndrome: Clinical Phenotype and Molecular Analysis in a Single Tertiary Center Cohort

OBJECTIVE

The aim of this study was the molecular characterization of the AR gene as the cause of 46,XY disorder in our population.

METHODS

We studied 41, non related, 46,XY disorder of sexual differentiation index cases, having characteristics consistent with androgen insensivity syndrome (AIS). Genomic DNA was isolated from peripheral blood leukocytes of all patients and 25 family members from 17 non-related families.

RESULTS

The AR gene analysis revealed an abnormal sequence in 58.5% of the index patients. All of the complete AIS (CAIS) cases were genetically confirmed, while in the partial form (PAIS) a mutation in AR was detected in only 13 (43.3%). Molecular studies revealed other affected or carrier relatives in 87% of the index cases. The AR mutations were found spread along the whole coding sequence, with a higher prevalence in the ligand binding domain. Nine out of 23 (39%) AR mutations were novel. In 17% of patients with detected AR mutations, somatic mosaicism was detected in leucocyte DNA. In our cohort, long-term follow up gender dysphoria, raised as male or female, was not found. Finally, in suspected PAIS, the identification of AR mutation occurred significantly less than in CAIS patients.

CONCLUSION

Improved knowledge of the components of the AR complex and signaling network might contribute to long term outcome and genetic counseling in AIS patients.

Loss of function, missense, and intronic variants in NOTCH1 confer different risks for left ventricular outflow tract obstructive heart defects in two European cohorts

Loss of function variants in NOTCH1 cause left ventricular outflow tract obstructive defects (LVOTO). However, the risk conferred by rare and noncoding variants in NOTCH1 for LVOTO remains largely uncharacterized. In a cohort of 49 families affected by hypoplastic left heart syndrome, a severe form of LVOTO, we discovered predicted loss of function NOTCH1 variants in 6% of individuals. Rare or low-frequency missense variants were found in 16% of families. To make a quantitative estimate of the genetic risk posed by variants in NOTCH1 for LVOTO, we studied associations of 400 coding and noncoding variants in NOTCH1 in 1,085 cases and 332,788 controls from the UK Biobank. Two rare intronic variants in strong linkage disequilibrium displayed significant association with risk for LVOTO amongst European-ancestry individuals. This result was replicated in an independent analysis of 210 cases and 68,762 controls of non-European and mixed ancestry. In conclusion, carrying rare predicted loss of function variants in NOTCH1 confer significant risk for LVOTO. In addition, the two intronic variants seem to be associated with an increased risk for these defects. Our approach demonstrates the utility of population-based data sets in quantifying the specific risk of individual variants for disease-related phenotypes.

Understanding human DNA variants affecting pre-mRNA splicing in the NGS era

Pre-mRNA splicing, an essential step in eukaryotic gene expression, relies on recognition of short sequences on the primary transcript intron ends and takes place along transcription by RNA polymerase II. Exonic and intronic auxiliary elements may modify the strength of exon definition and intron recognition. Splicing DNA variants (SV) have been associated with human genetic diseases at canonical intron sites, as well as exonic substitutions putatively classified as nonsense, missense or synonymous variants. Their effects on mRNA may be modulated by cryptic splice sites associated to the SV allele, comprehending exon skipping or shortening, and partial or complete intron retention. As splicing mRNA outputs result from combinatorial effects of both intrinsic and extrinsic factors, in vitro functional assays supported by computational analyses are recommended to assist SV pathogenicity assessment for human Mendelian inheritance diseases. The increasing use of next-generating sequencing (NGS) targeting full genomic gene sequence has raised awareness of the relevance of deep intronic SV in genetic diseases and inclusion of pseudo-exons into mRNA. Finally, we take advantage of recent advances in sequencing and computational technologies to analyze alternative splicing in cancer. We explore the Catalog of Somatic Mutations in Cancer (COSMIC) to describe the proportion of splice-site mutations in cis and trans regulatory elements. Genomic data from large cohorts of different cancer types are increasingly available, in addition to repositories of normal and somatic genetic variations. These are likely to bring new insights to understanding the genetic control of alternative splicing by mapping splicing quantitative trait loci in tumors.

RNA splicing analysis in genomic medicine

High-throughput next-generation sequencing technologies have led to a rapid increase in the number of sequence variants identified in clinical practice via diagnostic genetic tests. Current bioinformatic analysis pipelines fail to take adequate account of the possible splicing effects of such variants, particularly where variants fall outwith canonical splice site sequences, and consequently the pathogenicity of such variants may often be missed. The regulation of splicing is highly complex and as a result, in silico prediction tools lack sufficient sensitivity and specificity for reliable use. Variants of all kinds can be linked to aberrant splicing in disease and the need for correct identification and diagnosis grows ever more crucial as novel splice-switching antisense oligonucleotide therapies start to enter clinical usage. RT-PCR provides a useful targeted assay of the splicing effects of identified variants, while minigene assays, massive parallel reporter assays and animal models can also be used for more detailed study of a particular splicing system, given enough time and resources. However, RNA-sequencing (RNA-seq) has the potential to be used as a rapid diagnostic tool in genomic medicine. By utilising data science approaches and machine learning, it may prove possible to finally understand and interpret the 'splicing code' and apply this knowledge in human disease diagnostics.

Meta-analysis of the diagnostic and clinical utility of genome and exome sequencing and chromosomal microarray in children with suspected genetic diseases

Genetic diseases are leading causes of childhood mortality. Whole-genome sequencing (WGS) and whole-exome sequencing (WES) are relatively new methods for diagnosing genetic diseases, whereas chromosomal microarray (CMA) is well established. Here we compared the diagnostic utility (rate of causative, pathogenic, or likely pathogenic genotypes in known disease genes) and clinical utility (proportion in whom medical or surgical management was changed by diagnosis) of WGS, WES, and CMA in children with suspected genetic diseases by systematic review of the literature (January 2011–August 2017) and meta-analysis, following MOOSE/PRISMA guidelines. In 37 studies, comprising 20,068 children, diagnostic utility of WGS (0.41, 95% CI 0.34–0.48, I² = 44%) and WES (0.36, 95% CI 0.33–0.40, I² = 83%) were qualitatively greater than CMA (0.10, 95% CI 0.08–0.12, I² = 81%). Among studies published in 2017, the diagnostic utility of WGS was significantly greater than CMA (P < 0.0001, I² = 13% and I² = 40%, respectively). Among studies featuring within-cohort comparisons, the diagnostic utility of WES was significantly greater than CMA (P < 0.001, I² = 36%). The diagnostic utility of WGS and WES were not significantly different. In studies featuring within-cohort comparisons of WGS/WES, the likelihood of diagnosis was significantly greater for trios than singletons (odds ratio 2.04, 95% CI 1.62–2.56, I² = 12%; P < 0.0001). Diagnostic utility of WGS/WES with hospital-based interpretation (0.42, 95% CI 0.38–0.45, I² = 48%) was qualitatively higher than that of reference laboratories (0.29, 95% CI 0.27–0.31, I² = 49%); this difference was significant among studies published in 2017 (P < .0001, I² = 22% and I² = 26%, respectively). The clinical utility of WGS (0.27, 95% CI 0.17–0.40, I² = 54%) and WES (0.17, 95% CI 0.12–0.24, I² = 76%) were higher than CMA (0.06, 95% CI 0.05–0.07, I² = 42%); this difference was significant for WGS vs CMA (P < 0.0001). In conclusion, in children with suspected genetic diseases, the diagnostic and clinical utility of WGS/WES were greater than CMA. Subgroups with higher WGS/WES diagnostic utility were trios and those receiving hospital-based interpretation. WGS/WES should be considered a first-line genomic test for children with suspected genetic diseases.

Children with suspected genetic disease are more likely to receive a diagnosis when gene sequencing technologies are employed. At present, a microarray test for short pieces of missing or extra chromosomes (chromosomal microarray, CMA) is the recommended test for diagnosing genetic disease in children. The recent use of whole-genome sequencing (WGS) and whole-exome sequencing (WES) for diagnosing genetic disease prompted Stephen Kingsmore at Rady Children’s Institute for Genomic Medicine in San Diego, California, USA, and colleagues to examine the effectiveness of these methods compared with CMA. Analyses of studies comprising over 20,000 children showed that using WGS and WES increased the likelihood of diagnosis and led to a change in inpatient management to prevent or ameliorate adverse health outcomes. The authors advocate the use of WGS or WES as a first-line genomic test for genetic disease in children.

Splicing mutations in human genetic disorders: examples, detection, and confirmation

Precise pre-mRNA splicing, essential for appropriate protein translation, depends on the presence of consensus "cis" sequences that define exon-intron boundaries and regulatory sequences recognized by splicing machinery. Point mutations at these consensus sequences can cause improper exon and intron recognition and may result in the formation of an aberrant transcript of the mutated gene. The splicing mutation may occur in both introns and exons and disrupt existing splice sites or splicing regulatory sequences (intronic and exonic splicing silencers and enhancers), create new ones, or activate the cryptic ones. Usually such mutations result in errors during the splicing process and may lead to improper intron removal and thus cause alterations of the open reading frame. Recent research has underlined the abundance and importance of splicing mutations in the etiology of inherited diseases. The application of modern techniques allowed to identify synonymous and nonsynonymous variants as well as deep intronic mutations that affected pre-mRNA splicing. The bioinformatic algorithms can be applied as a tool to assess the possible effect of the identified changes. However, it should be underlined that the results of such tests are only predictive, and the exact effect of the specific mutation should be verified in functional studies. This article summarizes the current knowledge about the "splicing mutations" and methods that help to identify such changes in clinical diagnosis.

Disorders of sex development

Normal sex development depends on the precise spatio-temporal sequence and coordination of mutually antagonistic activating and repressing factors. These factors regulate the commitment of the unipotential gonad into the binary pathways governing normal sex development. Typically, the presence of the SRY gene on the Y chromosome triggers the cascade of molecular events that lead to male sex development. Disorders of sex development comprise a heterogeneous group of congenital conditions associated with atypical development of internal and external genitalia. These disorders are generally attributed to deviations from the typical progression of sex development. Disorders of sex development can be classified into several categories including chromosomal, gonadal, and anatomic abnormalities. Genetic tools such as microarray analyses and next-generation sequencing techniques have identified novel genetic variants among patients with disorders of sexual development. Most importantly, patient management needs to be individualized, especially for decisions related to sex of rearing, surgical interventions, hormone treatment, and potential for fertility preservation.

The NSIGHT1-randomized controlled trial: rapid whole-genome sequencing for accelerated etiologic diagnosis in critically ill infants

Genetic disorders are a leading cause of morbidity and mortality in infants in neonatal and pediatric intensive care units (NICU/PICU). While genomic sequencing is useful for genetic disease diagnosis, results are usually reported too late to guide inpatient management. We performed an investigator-initiated, partially blinded, pragmatic, randomized, controlled trial to test the hypothesis that rapid whole-genome sequencing (rWGS) increased the proportion of NICU/PICU infants receiving a genetic diagnosis within 28 days. The participants were families with infants aged <4 months in a regional NICU and PICU, with illnesses of unknown etiology. The intervention was trio rWGS. Enrollment from October 2014 to June 2016, and follow-up until November 2016. Of all, 26 female infants, 37 male infants, and 2 infants of undetermined sex were randomized to receive rWGS plus standard genetic tests (n = 32, cases) or standard genetic tests alone (n = 33, controls). The study was terminated early due to loss of equipoise: 73% (24) controls received genomic sequencing as standard tests, and 15% (five) controls underwent compassionate cross-over to receive rWGS. Nevertheless, intention to treat analysis showed the rate of genetic diagnosis within 28 days of enrollment (the primary end-point) to be higher in cases (31%, 10 of 32) than controls (3%, 1 of 33; difference, 28% [95% CI, 10-46%]; p = 0.003). Among infants enrolled in the first 25 days of life, the rate of neonatal diagnosis was higher in cases (32%, 7 of 22) than controls (0%, 0 of 23; difference, 32% [95% CI, 11-53%];p = 0.004). Median age at diagnosis (25 days [range 14-90] in cases vs. 130 days [range 37-451] in controls) and median time to diagnosis (13 days [range 1-84] in cases, vs. 107 days [range 21-429] in controls) were significantly less in cases than controls (p = 0.04). In conclusion, rWGS increased the proportion of NICU/PICU infants who received timely diagnoses of genetic diseases.

Partial androgen insensitivity syndrome caused by a deep intronic mutation creating an alternative splice acceptor site of the AR gene

Although partial androgen insensitivity syndrome (PAIS) is caused by attenuated responsiveness to androgens, androgen receptor gene (AR) mutations on the coding regions and their splice sites have been identified only in <25% of patients with a diagnosis of PAIS. We performed extensive molecular studies including whole exome sequencing in a Japanese family with PAIS, identifying a deep intronic variant beyond the branch site at intron 6 of AR (NM_000044.4:c.2450−42 G > A). This variant created the splice acceptor motif that was accompanied by pyrimidine-rich sequence and two candidate branch sites. Consistent with this, reverse transcriptase (RT)-PCR experiments for cycloheximide-treated lymphoblastoid cell lines revealed a relatively large amount of aberrant mRNA produced by the newly created splice acceptor site and a relatively small amount of wildtype mRNA produced by the normal splice acceptor site. Furthermore, most of the aberrant mRNA was shown to undergo nonsense mediated decay (NMD) and, if a small amount of aberrant mRNA may have escaped NMD, such mRNA was predicted to generate a truncated AR protein missing some functional domains. These findings imply that the deep intronic mutation creating an alternative splice acceptor site resulted in the production of a relatively small amount of wildtype AR mRNA, leading to PAIS.

Novel Variant of the Androgen Receptor Gene in a Patient With Complete Androgen Insensitivity Syndrome and Polyorchidism

Introduction: Complete androgen insensitivity (CAIS) in 65-95% cases is caused by pathogenic allelic variants (mutations) in the gene encoding androgen receptor (AR gene) and is characterized by female phenotype development with a male karyotype (46, XY). Patients are usually diagnosed during puberty and undergo gonadectomy due to increased testicular germ cell tumor risk. Only a few outcomes have been reported in older individuals with postponed gonadectomy. Case presentation: A 48-year-old CAIS patient presented with polyorchidism (four testes) without gonadal malignancies. Genetic testing identified a novel allelic variant in the AR gene [c.2141T>G (p.Phe805Cys)] causing the clinical symptoms. Conclusion: We have described a unique patient with CAIS and polyorchidism without malignancies in her late 40's bearing a novel likely pathogenic variant in the AR gene.

Meta-analysis of the diagnostic and clinical utility of genome and exome sequencing and chromosomal microarray in children with suspected genetic diseases

Genetic diseases are leading causes of childhood mortality. Whole-genome sequencing (WGS) and whole-exome sequencing (WES) are relatively new methods for diagnosing genetic diseases, whereas chromosomal microarray (CMA) is well established. Here we compared the diagnostic utility (rate of causative, pathogenic, or likely pathogenic genotypes in known disease genes) and clinical utility (proportion in whom medical or surgical management was changed by diagnosis) of WGS, WES, and CMA in children with suspected genetic diseases by systematic review of the literature (January 2011-August 2017) and meta-analysis, following MOOSE/PRISMA guidelines. In 37 studies, comprising 20,068 children, diagnostic utility of WGS (0.41, 95% CI 0.34-0.48, I² = 44%) and WES (0.36, 95% CI 0.33-0.40, I² = 83%) were qualitatively greater than CMA (0.10, 95% CI 0.08-0.12, I² = 81%). Among studies published in 2017, the diagnostic utility of WGS was significantly greater than CMA (P < 0.0001, I² = 13% and I² = 40%, respectively). Among studies featuring within-cohort comparisons, the diagnostic utility of WES was significantly greater than CMA (P < 0.001, I² = 36%). The diagnostic utility of WGS and WES were not significantly different. In studies featuring within-cohort comparisons of WGS/WES, the likelihood of diagnosis was significantly greater for trios than singletons (odds ratio 2.04, 95% CI 1.62-2.56, I² = 12%; P < 0.0001). Diagnostic utility of WGS/WES with hospital-based interpretation (0.42, 95% CI 0.38-0.45, I² = 48%) was qualitatively higher than that of reference laboratories (0.29, 95% CI 0.27-0.31, I² = 49%); this difference was significant among studies published in 2017 (P < .0001, I² = 22% and I² = 26%, respectively). The clinical utility of WGS (0.27, 95% CI 0.17-0.40, I² = 54%) and WES (0.17, 95% CI 0.12-0.24, I² = 76%) were higher than CMA (0.06, 95% CI 0.05-0.07, I² = 42%); this difference was significant for WGS vs CMA (P < 0.0001). In conclusion, in children with suspected genetic diseases, the diagnostic and clinical utility of WGS/WES were greater than CMA. Subgroups with higher WGS/WES diagnostic utility were trios and those receiving hospital-based interpretation. WGS/WES should be considered a first-line genomic test for children with suspected genetic diseases.

Identification of the Underlying Androgen Receptor Defect in the Dallas Reifenstein Family

CONTEXT

The Dallas Reifenstein family - first described in 1965 - includes 14 known members with partial androgen insensitivity syndrome (PAIS). However, the underlying molecular defect was never identified.

OBJECTIVE

To identify the underlying genetic defect for PAIS in the Dallas Reifenstein family.

DESIGN

DNA was purified from scrotal skin fibroblasts, and whole exome sequencing was then performed in four affected men in the family. Additional family members - both affected and unaffected - were subsequently recruited to confirm segregation of the candidate mutations with the PAIS phenotype.

PATIENTS

The affected men have PAIS with infertility associated with azoospermia, hypospadias, and gynecomastia.

RESULTS

All four men harbored an intronic variant NC_000023.10:g.66788676A>C between exon 1 and exon 2 of the androgen receptor (AR) canonical transcript NM_000044 (complementary DNA position NM_000044: c.1616+22072A>C) predicted to cause an alternatively spliced AR transcript. Reverse transcription (RT) polymerase chain (PCR) experiments detected the predicted PCR product of the alternatively spliced AR transcript, and the mutation segregated with the PAIS phenotype in this family. The transcript includes the insertion of 185 nucleotides with a premature stop codon at chrX:66863131-66863133, likely resulting in a reduction in AR protein expression due to nonsense-mediated decay.

CONCLUSIONS

An intronic AR mutation was identified in the Dallas Reifenstein family. The findings suggest that in cases of PAIS without identifiable AR mutations in coding regions, intronic AR mutations should be considered.

Non-coding variation in disorders of sex development

Genetic studies in Disorders of Sex Development (DSD), representing a wide spectrum of developmental or functional conditions of the gonad, have mainly been oriented towards the coding genome. Application of genomic technologies, such as whole-exome sequencing, result in a molecular genetic diagnosis in ∼50% of cases with DSD. Many of the genes mutated in DSD encode transcription factors such as SRY, SOX9, NR5A1, and FOXL2, characterized by a strictly regulated spatiotemporal expression. Hence, it can be hypothesized that at least part of the missing genetic variation in DSD can be explained by non-coding mutations in regulatory elements that alter gene expression, either by reduced, mis- or overexpression of their target genes. In addition, structural variations such as translocations, deletions, duplications or inversions can affect the normal chromatin conformation by different mechanisms. Here, we review non-coding defects in human DSD phenotypes and in animal models. The wide variety of non-coding defects found in DSD emphasizes that the regulatory landscape of known and to be discovered DSD genes has to be taken into consideration when investigating the molecular pathogenesis of DSD.

Clinically Focused Molecular Investigation of 1000 Consecutive Families with Inherited Retinal Disease

PURPOSE

To devise a comprehensive multiplatform genetic testing strategy for inherited retinal disease and to describe its performance in 1000 consecutive families seen by a single clinician.

DESIGN

Retrospective series.

PARTICIPANTS

One thousand consecutive families seen by a single clinician.

METHODS

The clinical records of all patients seen by a single retina specialist between January 2010 and June 2016 were reviewed, and all patients who met the clinical criteria for a diagnosis of inherited retinal disease were included in the study. Each patient was assigned to 1 of 62 diagnostic categories, and this clinical diagnosis was used to define the scope and order of the molecular investigations that were performed. The number of nucleotides evaluated in a given subject ranged from 2 to nearly 900 000.

MAIN OUTCOME MEASURES

Sensitivity and false genotype rate.

RESULTS

Disease-causing genotypes were identified in 760 families (76%). These genotypes were distributed across 104 different genes. More than 75% of these 104 genes have coding sequences small enough to be packaged efficiently into an adeno-associated virus. Mutations in ABCA4 were the most common cause of disease in this cohort (173 families), whereas mutations in 80 genes caused disease in 5 or fewer families (i.e., 0.5% or less). Disease-causing genotypes were identified in 576 of the families without next-generation sequencing (NGS). This included 23 families with mutations in the repetitive region of RPGR exon 15 that would have been missed by NGS. Whole-exome sequencing of the remaining 424 families revealed mutations in an additional 182 families, and whole-genome sequencing of 4 of the remaining 242 families revealed 2 additional genotypes that were invisible by the other methods. Performing the testing in a clinically focused tiered fashion would be 6.1% more sensitive and 17.7% less expensive and would have a significantly lower average false genotype rate than using whole-exome sequencing to assess more than 300 genes in all patients (7.1% vs. 128%; P < 0.001).

CONCLUSIONS

Genetic testing for inherited retinal disease is now more than 75% sensitive. A clinically directed tiered testing strategy can increase sensitivity and improve statistical significance without increasing cost.