Whole-Transcriptome Sequencing-Based Profiling of the Cutaneous Virome in Patients with Secondary Immunodeficiency.

Most viral infections can be self-limited, with no requirement for medical intervention. However, the same viruses can cause severe diseases in patients with compromised immunity due to single-gene diseases, acquired immune deficiency syndrome, or hematologic malignancies or those receiving immunosuppressive drugs. Occasionally, these immunocompromised patients harbor >1 infectious agent, requiring several concomitant diagnostic tests. We have developed, to our knowledge, a previously unreported whole-transcriptome sequencing-based pipeline that allows virome profiling, quantitation, and expression pattern analysis of 926 distinct viruses by sequencing of RNA isolated from a single lesional skin biopsy. This pipeline can also explore host genetics if there is a Mendelian predisposition to infection. We applied this pipeline to 6 Iranian patients with viral-induced skin lesions associated with immune deficiency secondary to HIV, human T-lymphotropic virus 1, chronic lymphocytic leukemia, and post transplant immunosuppression. In 5 cases, definitive human papillomavirus infections were identified, some caused by multiple viral types. In addition to human papillomavirus, coinfection with other viruses (Merkle cell polyomavirus, cytomegalovirus, and human herpesvirus 4) was detected in some lesions. In 1 case, whole-transcriptome sequencing validated the clinical diagnosis of adult T-cell leukemia/lymphoma in a patient with an initial diagnosis of mycosis fungoides/Sézary syndrome. These findings attest to the power of whole-transcriptome sequencing in profiling the cutaneous virome in the context of compromised immunity.

Monogenic etiologies of persistent human papillomavirus infections: A comprehensive systematic review.

PURPOSE: Persistent human papillomavirus infection (PHPVI) causes cutaneous, anogenital, and mucosal warts. Cutaneous warts include common warts, Treeman syndrome, and epidermodysplasia verruciformis, among others. Although more reports of monogenic predisposition to PHPVI have been published with the development of genomic technologies, genetic testing is rarely incorporated into clinical assessments. To encourage broader molecular testing, we compiled a list of the various monogenic etiologies of PHPVI. METHODS: We conducted a systematic literature review to determine the genetic, immunological, and clinical characteristics of patients with PHPVI. RESULTS: The inclusion criteria were met by 261 of 40,687 articles. In 842 patients, 83 PHPVI-associated genes were identified, including 42, 6, and 35 genes with strong, moderate, and weak evidence for causality, respectively. Autosomal recessive inheritance predominated (69%). PHPVI onset age was 10.8 ± 8.6 years, with an interquartile range of 5 to 14 years. GATA2,IL2RG,DOCK8, CXCR4, TMC6, TMC8, and CIB1 are the most frequently reported PHPVI-associated genes with strong causality. Most genes (74 out of 83) belong to a catalog of 485 inborn errors of immunity-related genes, and 40 genes (54%) are represented in the nonsyndromic and syndromic combined immunodeficiency categories. CONCLUSION: PHPVI has at least 83 monogenic etiologies and a genetic diagnosis is essential for effective management.

High Comorbidity of Pediatric Cancers in Patients with Birth Defects: Insights from Whole Genome Sequencing Analysis of Copy Number Variations.

BACKGROUND: Patients with birth defects (BD) exhibit an elevated risk of cancer. We aimed to investigate the potential link between pediatric cancers and BDs, exploring the hypothesis of shared genetic defects contributing to the coexistence of these conditions. METHODS: This study included 1454 probands with BDs (704 females and 750 males), including 619 (42.3%) with and 845 (57.7%) without co-occurrence of pediatric onset cancers. Whole genome sequencing (WGS) was done at 30X coverage through the Kids First/Gabriella Miller X01 Program. RESULTS: 8211 CNV loci were called from the 1454 unrelated individuals. 191 CNV loci classified as pathogenic/likely pathogenic (P/LP) were identified in 309 (21.3%) patients, with 124 (40.1%) of these patients having pediatric onset cancers. The most common group of CNVs are pathogenic deletions covering the region ChrX:52,863,011-55,652,521, seen in 162 patients including 17 males. Large recurrent P/LP duplications >5MB were detected in 33 patients. CONCLUSIONS: This study revealed that P/LP CNVs were common in a large cohort of BD patients with high rate of pediatric cancers. We present a comprehensive spectrum of P/LP CNVs in patients with BDs and various cancers. Notably, deletions involving E2F target genes and genes implicated in mitotic spindle assembly and G2/M checkpoint were identified, potentially disrupting cell-cycle progression and providing mechanistic insights into the concurrent occurrence of BDs and cancers.

Genomic variants exclusively identified in children with birth defects and concurrent malignant tumors predispose to cancer development.

Children with birth defects (BD) express distinct clinical features that often have various medical consequences, one of which is predisposition to the development of cancers. Identification of the underlying genetic mechanisms related to the development of cancer in BD patients would allow for preventive measures. We performed a whole genome sequencing (WGS) study on blood-derived DNA samples from 1566 individuals without chromosomal anomalies, including 454 BD probands with at least one type of malignant tumors, 767 cancer-free BD probands, and 345 healthy individuals. Exclusive recurrent variants were identified in BD-cancer and BD-only patients and mapped to their corresponding genomic regions. We observed statistically significant overlaps for protein-coding/ncRNA with exclusive variants in exons, introns, ncRNAs, and 3'UTR regions. Exclusive exonic variants, especially synonymous variants, tend to occur in prior exons locus in BD-cancer children. Intronic variants close to splicing site (< 500 bp from exon) have little overlaps in BD-cancer and BD-only patients. Exonic variants in non-coding RNA (ncRNA) tend to occur in different ncRNAs exons regardless of the overlaps. Notably, genes with 5' UTR variants are almost mutually exclusive between the two phenotypes. In conclusion, we conducted the first genomic study to explore the impact of recurrent variants exclusive to the two distinguished clinical phenotypes under study, BD with or without cancer, demonstrating enrichment of selective protein-coding/ncRNAs differentially expressed between these two phenotypes, suggesting that selective genetic factors may underlie the molecular processes of pediatric cancer development in BD children.

Whole transcriptome-based skin virome profiling in typical epidermodysplasia verruciformis reveals α-, ß-, and γ-HPV infections.

HPVs are DNA viruses include approximately 450 types that are classified into 5 genera (α-, ß-, γ-, µ-, and ν-HPV). The γ- and ß-HPVs are present in low copy numbers in healthy individuals; however, in patients with an inborn error of immunity, certain species of ß-HPVs can cause epidermodysplasia verruciformis (EV), manifesting as recalcitrant cutaneous warts and skin cancer. EV presents as either typical or atypical. Manifestations of typical EV are limited to the skin and are caused by abnormal keratinocyte-intrinsic immunity to ß-HPVs due to pathogenic sequence variants in TMC6, TMC8, or CIB1. We applied a transcriptome-based computational pipeline, VirPy, to RNA extracted from normal-appearing skin and wart samples of patients with typical EV to explore the viral and human genetic determinants. In 26 patients, 9 distinct biallelic mutations were detected in TMC6, TMC8, and CIB1, 7 of which are previously unreported to our knowledge. Additionally, 20 different HPV species, including 3 α-HPVs, 16 ß-HPVs, and 1 γ-HPV, were detected, 8 of which are reported here for the first time to our knowledge in patients with EV (ß-HPV-37, -47, -80, -151, and -159; α-HPV-2 and -57; and γ-HPV-128). This study expands the TMC6, TMC8, and CIB1 sequence variant spectrum and implicates new HPV subtypes in the pathogenesis of typical EV.

Identification of risk variants related to malignant tumors in children with birth defects by whole genome sequencing.

BACKGROUND: Children with birth defects (BD) are more likely to develop cancer and the increased risk of cancer persists into adulthood. Prior population-based assessments have demonstrated that even non-chromosomal BDs are associated with at least two-fold increase of cancer risk. Identification of variants that are associated with malignant tumor in BD patients without chromosomal anomalies may improve our understanding of the underlying molecular mechanisms and provide clues for early cancer detection in children with BD. METHODS: In this study, whole genome sequencing (WGS) data of blood-derived DNA for 1653 individuals without chromosomal anomalies were acquired from the Kids First Data Resource Center (DRC), including 541 BD probands with at least one type of malignant tumors, 767 BD probands without malignant tumor, and 345 healthy family members who are the parents or siblings of the probands. Recurrent variants exclusively seen in cancer patients were selected and mapped to their corresponding genomic regions. The targeted genes/non-coding RNAs were further reduced using random forest and forward feature selection (ffs) models. RESULTS: The filtered genes/non-coding RNAs, including variants in non-coding areas, showed enrichment in cancer-related pathways. To further support the validity of these variants, blood WGS data of additional 40 independent BD probands, including 25 patients with at least one type of cancers from unrelated projects, were acquired. The counts of variants of interest identified in the Kid First data showed clear deviation in the validation dataset between BD patients with cancer and without cancer. Furthermore, a deep learning model was built to assess the predictive abilities in the 40 patients using variants of interest identified in the Kids First cohort as feature vectors. The accuracies are ~ 75%, with the noteworthy observation that variants mapped to non-coding regions provided the highest accuracy (31 out of 40 patients were labeled correctly). CONCLUSION: We present for the first time a panorama of genetic variants that are associated with cancers in non-chromosomal BD patients, implying that our approach may potentially serve for the early detection of malignant tumors in patients with BD.

Recalcitrant Cutaneous Warts in a Family with Inherited ICOS Deficiency.

Recalcitrant warts, caused by human papillomaviruses (HPVs), can be a cutaneous manifestation of inborn error of immunity. This study investigated the clinical manifestations, immunodeficiency, single-gene susceptibility, and HPV repertoire in a consanguineous family with severe sinopulmonary infections and recalcitrant warts. Clinical and immunologic evaluations, including FACS and lymphocyte transformation test, provided evidence for immunodeficiency. Combined whole-exome sequencing and genome-wide homozygosity mapping were utilized to disclose candidate sequence variants. Whole-transcriptome sequencing was used to concomitantly investigate the HPV genotypes and the consequences of detected sequence variants in the host. The proband, a male aged 41 years, was found to be homozygous for the c.6delG, p.Lys2Asnfs∗17 variant in ICOS, encoding the inducible T-cell costimulator. This variant was located inside the 5 megabase of runs of homozygosity on 2q33.2. RNA sequencing confirmed the deleteriousness of the ICOS variant in three skin biopsies revealing significant downregulation of ICOS and its ligand, ICOSLG. Reads unaligned to the human genome were applied to 926 different viruses, and α-HPV57, ß-HPV107, ß-HPV14, and ß-HPV17 were detected. Collectively, we describe a previously unrecognized inborn error of T-cell immunity to HPVs, indicating that autosomal recessive ICOS deficiency can underlie recalcitrant warts, emphasizing the immunologic underpinnings of recalcitrant warts at the nexus of human and viral genomic variation.

Whole-transcriptome sequencing-based concomitant detection of viral and human genetic determinants of cutaneous lesions.

Severe viral infections of the skin can occur in patients with inborn errors of immunity (IEI). We report an all-in-one whole-transcriptome sequencing-based method by RNA-Seq on a single skin biopsy for concomitantly identifying the cutaneous virome and the underlying IEI. Skin biopsies were obtained from healthy and lesional skin from patients with cutaneous infections suspected to be of viral origin. RNA-Seq was utilized as the first-tier strategy for unbiased human genome-wide rare variant detection. Reads unaligned to the human genome were utilized for the exploration of 926 viruses in a viral genome catalog. In 9 families studied, the patients carried pathogenic variants in 6 human IEI genes, including IL2RG, WAS, CIB1, STK4, GATA2, and DOCK8. Gene expression profiling also confirmed pathogenicity of the human variants and permitted genome-wide homozygosity mapping, which assisted in identification of candidate genes in consanguineous families. This automated, online, all-in-one computational pipeline, called VirPy, enables simultaneous detection of the viral triggers and the human genetic variants underlying skin lesions in patients with suspected IEI and viral dermatosis.

Recalcitrant Warts, Epidermodysplasia Verruciformis, and the Tree-Man Syndrome: Phenotypic Spectrum of Cutaneous Human Papillomavirus Infections at the Intersection of Genetic Variability of Viral and Human Genomes.

Human papillomavirus (HPV) infections can cause common warts, which usually resolve spontaneously or become recalcitrant, resistant to multiple treatments. In rare cases, they transform into cutaneous giant horns resulting in the tree-man syndrome (TMS). Defective ß-HPVs can cause flat warts in epidermodysplasia verruciformis (EV), a genetic disorder. In typical EV, limited to the skin, the mutated genes are critical for keratinocyte-intrinsic immunity, whereas atypical, syndromic EV involves genes controlling T cells. Inborn errors of immunity due to mutations in distinct genes underlying recalcitrant warts and the α-HPV2‒driven TMS have been identified, all disrupting T-cell immunity. Collectively, these observations attest to the wide phenotypic spectrum of cutaneous infections caused by different HPV types at the intersection of the genetic diversity of the viral and human genomes.

Genetic Predisposition to Numerous Large Ulcerating Basal Cell Carcinomas and Response to Immune Therapy.

OBJECTIVE: Well-defined germ-line mutations in the PTCH1 gene are associated with syndromic multiple basal cell carcinomas (BCCs). Here, we used whole exome sequencing (WES) to identify the role of patched-1 in patients with multiple, unusually large BCCs. METHODS: A 72-year old patient presenting with numerous BCCs progressing to large ulcerating lesions was enrolled. WES was used to identify the pathogenic gene locus. RESULTS: Genetic work-up by WES identified a homozygous PTCH1 nonsense mutation in the tumor tissue but not present in her blood cells or in non-lesional skin. In addition, heterozygous missense mutations were identified in three cancer-associated genes (EPHB2, RET, and GALNT12) in blood cells as well as in lesional and non-lesional skin. We also tested systemic immune therapy as a potentially beneficial approach to treat patients with numerous large BCCs on scatted areas of involvement. A rapid and sustained response to nivolumab was noted, suggesting that it is an efficacious drug for long-term therapeutic outcome. CONCLUSION: PTCH1, EPHB2, RET, and GALNT12 may potentially contribute to the synergistic oncogene driven malignant transformation manifesting as multiple, unusually large BCCs.

Humans with inherited T cell CD28 deficiency are susceptible to skin papillomaviruses but are otherwise healthy.

We study a patient with the human papilloma virus (HPV)-2-driven "tree-man" phenotype and two relatives with unusually severe HPV4-driven warts. The giant horns form an HPV-2-driven multifocal benign epithelial tumor overexpressing viral oncogenes in the epidermis basal layer. The patients are unexpectedly homozygous for a private CD28 variant. They have no detectable CD28 on their T cells, with the exception of a small contingent of revertant memory CD4+ T cells. T cell development is barely affected, and T cells respond to CD3 and CD2, but not CD28, costimulation. Although the patients do not display HPV-2- and HPV-4-reactive CD4+ T cells in vitro, they make antibodies specific for both viruses in vivo. CD28-deficient mice are susceptible to cutaneous infections with the mouse papillomavirus MmuPV1. The control of HPV-2 and HPV-4 in keratinocytes is dependent on the T cell CD28 co-activation pathway. Surprisingly, human CD28-dependent T cell responses are largely redundant for protective immunity.

Kindler epidermolysis bullosa-like skin phenotype and downregulated basement membrane zone gene expression in poikiloderma with neutropenia and a homozygous USB1 mutation.

Epidermolysis bullosa (EB) is a genotypically heterogeneous group of disorders characterized by cutaneous blistering and erosions with a tremendous spectrum of severity. One of the distinct forms of EB, Kindler EB (KEB), manifests with blistering and poikiloderma; this subtype of EB is caused by mutations in the FERMT1 gene encoding kindlin-1. In this study, we investigated a patient clinically diagnosed as KEB with reduced FERMT1 gene expression and intensity of immunostaining for kindlin-1. Transmission electron microscopy showed lamina densa reduplication, frequently observed in KEB. However, no mutations were identified in FERMT1 in this patient with consanguineous parents, and this gene resided outside of genomic regions of homozygosity (ROH). Instead, whole-exome sequencing and homozygosity mapping identified a homozygous sequence variant at the +4 position of intron 2 in the USB1 gene, encoding an exoribonuclease required for processing of U6 snRNA, a critical component of spliceosomes. Examination of the patient's RNA by RNA-Seq confirmed the pathogenicity of this variant, causing aberrant splicing predicted to result in loss of function of USB1. Mutations in this gene have been reported in patients with poikiloderma and neutropenia, with a few reported cases in association with skin fragility, a condition distinct from the KEB phenotype. Transcriptome analysis revealed that several genes, expressed in the cutaneous basement membrane zone and previously associated with different subtypes of EB, were differentially downregulated at the mRNA level. EB-associated mRNA downregulation was confirmed at protein levels by skin immunofluorescence. These observations provide a novel mechanism for blistering and erosions in the skin as a result reduced presence of adhesion complexes critical for stable association of epidermis and dermis at the level of cutaneous basement membrane zone.

Arrhythmogenic right ventricular cardiomyopathy in patients with biallelic JUP-associated skin fragility.

Arrhythmogenic right ventricular cardiomyopathy (ARVC), with skin manifestations, has been associated with mutations in JUP encoding plakoglobin. Genotype-phenotype correlations regarding the penetrance of cardiac involvement, and age of onset have not been well established. We examined a cohort of 362 families with skin fragility to screen for genetic mutations with next-generation sequencing-based methods. In two unrelated families, a previously unreported biallelic mutation, JUP: c.201delC; p.Ser68Alafs*92, was disclosed. The consequences of this mutation were determined by expression profiling both at tissue and ultrastructural levels, and the patients were evaluated by cardiac and cutaneous work-up. Whole-transcriptome sequencing by RNA-Seq revealed JUP as the most down-regulated gene among 21 skin fragility-associated genes. Immunofluorescence showed the lack of plakoglobin in the epidermis. Two probands, 2.5 and 22-year-old, with the same homozygous mutation, allowed us to study the cross-sectional progression of cardiac involvements in relation to age. The older patient had anterior T wave inversions, prolonged terminal activation duration (TAD), and RV enlargement by echocardiogram, and together with JUP mutation met definite ARVC diagnosis. The younger patient had no evidence of cardiac disease, but met possible ARVC diagnosis with one major criterion (the JUP mutation). In conclusion, we identified the same biallelic homozygous JUP mutation in two unrelated families with skin fragility, but cardiac findings highlighted age-dependent penetrance of ARVC. Thus, young, phenotypically normal patients with biallelic JUP mutations should be monitored for development of ARVC.

Linear basal cell nevus with a novel mosaic PTCH1 mutation.

The patched tumor suppressor gene (PTCH1) encodes a receptor, which is a key component of the hedgehog signalling pathway. Mutations in PTCH1 are implicated in the development of sporadic basal cell carcinomas (BCC), as well as those in Gorlin Syndrome. Rarely, BCCs may develop in a linear pattern along lines of Blaschko due to cutaneous mosaicism. In cases in which there are other features of Gorlin syndrome, genomic analysis has demonstrated lesional mutations in the Hedgehog signalling pathway. Causative mutations, however, have not been firmly demonstrated in the cases of linear and segmental BCCs in otherwise healthy individuals. Herein, we report a case of a 31 year-old Caucasian woman with linear development of multiple superficial BCCs in a Blaschkoid distribution without other characteristic findings of Gorlin syndrome. Genomic analysis of lesional skin by whole-exome sequencing identified a novel heterozygous mutation PTCH1: NM_000264.3, Exon 15, c.2336-2337insGGTAGGA, p.Asp779Glufs*13 in PTCH1, shared by two discrete samples within the lesion, while no mutations were found in the non-lesional skin or peripheral blood. Given the young age of our patient and linear distribution of BCCs on non-sun exposed skin, our findings suggest segmental mosaicism. The patient was treated with topical 5% imiquimod with histologically confirmed clearance of BCCs in 2 months.

Inherited Interleukin 2-Inducible T-Cell (ITK) Kinase Deficiency in Siblings With Epidermodysplasia Verruciformis and Hodgkin Lymphoma.

Biallelic mutations in the ITK gene cause a T-cell primary immunodeficiency with Epstein-Barr virus (EBV)-lymphoproliferative disorders. We describe a novel association of a homozygous ITK mutation with ß-human papillomavirus (HPV)-positive epidermodysplasia verruciformis. Thus, loss of function in ITK can result in broad dysregulation of T-cell responses to oncogenic viruses, including ß-HPV and EBV.

A novel autosomal recessive GJB2-associated disorder: Ichthyosis follicularis, bilateral severe sensorineural hearing loss, and punctate palmoplantar keratoderma.

Ichthyosis follicularis, a distinct cutaneous entity reported in combination with atrichia, and photophobia has been associated with mutations in MBTPS2. We sought the genetic cause of a novel syndrome of ichthyosis follicularis, bilateral severe sensorineural hearing loss and punctate palmoplantar keratoderma in two families. We performed whole exome sequencing on three patients from two families. The pathogenicity and consequences of mutations were studied in the Xenopus oocyte expression system and by molecular modeling analysis. Compound heterozygous mutations in the GJB2 gene were discovered: a pathogenic c.526A>G; p.Asn176Asp, and a common frameshift mutation, c.35delG; p.Gly12Valfs*2. The p.Asn176Asp missense mutation was demonstrated to significantly reduce the cell-cell gap junction channel activity and increase the nonjunctional hemichannel activity in the Xenopus oocyte expression system. Molecular modeling analyses of the mutant Cx26 protein revealed significant changes in the structural characteristics and electrostatic potential of the Cx26, either in hemichannel or gap junction conformation. Thus, association of a new syndrome of an autosomal recessive disorder of ichthyosis follicularis, bilateral severe sensorineural hearing loss and punctate palmoplantar keratoderma with mutations in GJB2, expands the phenotypic spectrum of the GJB2-associated disorders. The findings attest to the complexity of the clinical consequences of different mutations in GJB2.

The human CIB1-EVER1-EVER2 complex governs keratinocyte-intrinsic immunity to ß-papillomaviruses.

Patients with epidermodysplasia verruciformis (EV) and biallelic null mutations of TMC6 (encoding EVER1) or TMC8 (EVER2) are selectively prone to disseminated skin lesions due to keratinocyte-tropic human ß-papillomaviruses (ß-HPVs), which lack E5 and E8. We describe EV patients homozygous for null mutations of the CIB1 gene encoding calcium- and integrin-binding protein-1 (CIB1). CIB1 is strongly expressed in the skin and cultured keratinocytes of controls but not in those of patients. CIB1 forms a complex with EVER1 and EVER2, and CIB1 proteins are not expressed in EVER1- or EVER2-deficient cells. The known functions of EVER1 and EVER2 in human keratinocytes are not dependent on CIB1, and CIB1 deficiency does not impair keratinocyte adhesion or migration. In keratinocytes, the CIB1 protein interacts with the HPV E5 and E8 proteins encoded by α-HPV16 and γ-HPV4, respectively, suggesting that this protein acts as a restriction factor against HPVs. Collectively, these findings suggest that the disruption of CIB1-EVER1-EVER2-dependent keratinocyte-intrinsic immunity underlies the selective susceptibility to ß-HPVs of EV patients.