Long-read whole genome sequencing reveals HOXD13 alterations in synpolydactyly.

Synpolydactyly 1, also called syndactyly type II (SDTY2), is a genetic limb malformation characterized by polydactyly with syndactyly involving the webbing of the third and fourth fingers, and the fourth and fifth toes. It is caused by heterozygous alterations in HOXD13 with incomplete penetrance and phenotypic variability. In our study, a five-generation family with an SPD phenotype was enrolled in our Rare Disease Genomics Protocol. A comprehensive examination of three generations using Illumina short-read whole-genome sequencing (WGS) did not identify any causative variants. Subsequent WGS using Pacific Biosciences (PacBio) long-read HiFi Circular Consensus Sequencing (CCS) revealed a heterozygous 27-bp duplication in the polyalanine tract of HOXD13. Sanger sequencing of all available family members confirmed that the variant segregates with affected individuals. Reanalysis of an unrelated family with a similar SPD phenotype uncovered a 21-bp (7-alanine) duplication in the same region of HOXD13. Although ExpansionHunter identified these events in most individuals in a retrospective analysis, low sequence coverage due to high GC content in the HOXD13 polyalanine tract makes detection of these events challenging. Our findings highlight the value of long-read WGS in elucidating the molecular etiology of congenital limb malformation disorders.

Defining the transcriptome of PIK3CA-altered cells in a human capillary malformation using single cell long-read sequencing.

PIK3CA-related overgrowth spectrum (PROS) disorders are caused by somatic mosaic variants that result in constitutive activation of the phosphatidylinositol-3-kinase/AKT/mTOR pathway. Promising responses to molecularly targeted therapy have been reported, although identification of an appropriate agent can be hampered by the mosaic nature and corresponding low variant allele frequency of the causal variant. Moreover, our understanding of the molecular consequences of these variants-for example how they affect gene expression profiles-remains limited. Here we describe in vitro expansion of a human capillary malformation followed by molecular characterization using exome sequencing, single cell gene expression, and targeted long-read single cell RNA-sequencing in a patient with clinical features consistent with Megalencephaly-Capillary Malformation Syndrome (MCAP, a PROS condition). These approaches identified a targetable PIK3CA variant with expression restricted to PAX3+ fibroblast and undifferentiated keratinocyte populations. This study highlights the innovative combination of next-generation single cell sequencing methods to better understand unique transcriptomic profiles and cell types associated with MCAP, revealing molecular intricacies of this genetic syndrome.

Expanding the clinical phenotype of FGFR1 internal tandem duplication.

Closed spinal dysraphism (SD) is a type of neural tube defect originating during early embryonic development whereby the neural tissue of the spinal defect remains covered by skin, often coinciding with markers of cutaneous stigmata. It is hypothesized that these events are caused by multifactorial processes, including genetic and environmental causes. We present an infant with a unique congenital midline lesion associated with a closed SD. Through comprehensive molecular profiling of the intraspinal lesion and contiguous skin lesion, an internal tandem duplication (ITD) of the kinase domain of the fibroblast growth factor receptor 1 (FGFR1) gene was found. This ITD variant is somatic mosaic in nature as supported by a diminished variant allele frequency in the lesional tissue and by its absence in peripheral blood. FGFR1 ITD results in constitutive activation of the receptor tyrosine kinase to promote cell growth, differentiation, and survival through RAS/MAPK signaling. Identification of FGFR1 ITD outside of central nervous system tumors is exceedingly rare, and this report broadens the phenotypic spectrum of somatic mosaic FGFR1-related disease.

Analysis of expressed sequence tags generated from full-length enriched cDNA libraries of melon.

BACKGROUND: Melon (Cucumis melo), an economically important vegetable crop, belongs to the Cucurbitaceae family which includes several other important crops such as watermelon, cucumber, and pumpkin. It has served as a model system for sex determination and vascular biology studies. However, genomic resources currently available for melon are limited. RESULT: We constructed eleven full-length enriched and four standard cDNA libraries from fruits, flowers, leaves, roots, cotyledons, and calluses of four different melon genotypes, and generated 71,577 and 22,179 ESTs from full-length enriched and standard cDNA libraries, respectively. These ESTs, together with ~35,000 ESTs available in public domains, were assembled into 24,444 unigenes, which were extensively annotated by comparing their sequences to different protein and functional domain databases, assigning them Gene Ontology (GO) terms, and mapping them onto metabolic pathways. Comparative analysis of melon unigenes and other plant genomes revealed that 75% to 85% of melon unigenes had homologs in other dicot plants, while approximately 70% had homologs in monocot plants. The analysis also identified 6,972 gene families that were conserved across dicot and monocot plants, and 181, 1,192, and 220 gene families specific to fleshy fruit-bearing plants, the Cucurbitaceae family, and melon, respectively. Digital expression analysis identified a total of 175 tissue-specific genes, which provides a valuable gene sequence resource for future genomics and functional studies. Furthermore, we identified 4,068 simple sequence repeats (SSRs) and 3,073 single nucleotide polymorphisms (SNPs) in the melon EST collection. Finally, we obtained a total of 1,382 melon full-length transcripts through the analysis of full-length enriched cDNA clones that were sequenced from both ends. Analysis of these full-length transcripts indicated that sizes of melon 5' and 3' UTRs were similar to those of tomato, but longer than many other dicot plants. Codon usages of melon full-length transcripts were largely similar to those of Arabidopsis coding sequences. CONCLUSION: The collection of melon ESTs generated from full-length enriched and standard cDNA libraries is expected to play significant roles in annotating the melon genome. The ESTs and associated analysis results will be useful resources for gene discovery, functional analysis, marker-assisted breeding of melon and closely related species, comparative genomic studies and for gaining insights into gene expression patterns.

Molecular classification of a complex structural rearrangement of the RB1 locus in an infant with sporadic, isolated, intracranial, sellar region retinoblastoma.

Retinoblastoma is a childhood cancer of the retina involving germline or somatic alterations of the RB Transcriptional Corepressor 1 gene, RB1. Rare cases of sellar-suprasellar region retinoblastoma without evidence of ocular or pineal tumors have been described. A nine-month-old male presented with a sellar-suprasellar region mass. Histopathology showed an embryonal tumor with focal Flexner-Wintersteiner-like rosettes and loss of retinoblastoma protein (RB1) expression by immunohistochemistry. DNA array-based methylation profiling confidently classified the tumor as pineoblastoma group A/intracranial retinoblastoma. The patient was subsequently enrolled on an institutional translational cancer research protocol and underwent comprehensive molecular profiling, including paired tumor/normal exome and genome sequencing and RNA-sequencing of the tumor. Additionally, Pacific Biosciences (PacBio) Single Molecule Real Time (SMRT) sequencing was performed from comparator normal and disease-involved tissue to resolve complex structural variations. RNA-sequencing revealed multiple fusions clustered within 13q14.1-q21.3, including a novel in-frame fusion of RB1-SIAH3 predicted to prematurely truncate the RB1 protein. SMRT sequencing revealed a complex structural rearrangement spanning 13q14.11-q31.3, including two somatic structural variants within intron 17 of RB1. These events corresponded to the RB1-SIAH3 fusion and a novel RB1 rearrangement expected to correlate with the complete absence of RB1 protein expression. Comprehensive molecular analysis, including DNA array-based methylation profiling and sequencing-based methodologies, were critical for classification and understanding the complex mechanism of RB1 inactivation in this diagnostically challenging tumor.