Assortative mating and parental genetic relatedness contribute to the pathogenicity of variably expressive variants.

We examined more than 97,000 families from four neurodevelopmental disease cohorts and the UK Biobank to identify phenotypic and genetic patterns in parents contributing to neurodevelopmental disease risk in children. We identified within- and cross-disorder correlations between six phenotypes in parents and children, such as obsessive-compulsive disorder (R = 0.32-0.38, p < 10-126). We also found that measures of sub-clinical autism features in parents are associated with several autism severity measures in children, including biparental mean Social Responsiveness Scale scores and proband Repetitive Behaviors Scale scores (regression coefficient = 0.14, p = 3.38 × 10-4). We further describe patterns of phenotypic similarity between spouses, where spouses show correlations for six neurological and psychiatric phenotypes, including a within-disorder correlation for depression (R = 0.24-0.68, p < 0.001) and a cross-disorder correlation between anxiety and bipolar disorder (R = 0.09-0.22, p < 10-92). Using a simulated population, we also found that assortative mating can lead to increases in disease liability over generations and the appearance of "genetic anticipation" in families carrying rare variants. We identified several families in a neurodevelopmental disease cohort where the proband inherited multiple rare variants in disease-associated genes from each of their affected parents. We further identified parental relatedness as a risk factor for neurodevelopmental disorders through its inverse relationship with variant pathogenicity and propose that parental relatedness modulates disease risk by increasing genome-wide homozygosity in children (R = 0.05-0.26, p < 0.05). Our results highlight the utility of assessing parent phenotypes and genotypes toward predicting features in children who carry rare variably expressive variants and implicate assortative mating as a risk factor for increased disease severity in these families.

Integrating RNA-Seq into genome sequencing workflow enhances the analysis of structural variants causing neurodevelopmental disorders.

BACKGROUND: Molecular diagnosis of neurodevelopmental disorders (NDDs) is mainly based on exome sequencing (ES), with a diagnostic yield of 31% for isolated and 53% for syndromic NDD. As sequencing costs decrease, genome sequencing (GS) is gradually replacing ES for genome-wide molecular testing. As many variants detected by GS only are in deep intronic or non-coding regions, the interpretation of their impact may be difficult. Here, we showed that integrating RNA-Seq into the GS workflow can enhance the analysis of the molecular causes of NDD, especially structural variants (SVs), by providing valuable complementary information such as aberrant splicing, aberrant expression and monoallelic expression. METHODS: We performed trio-GS on a cohort of 33 individuals with NDD for whom ES was inconclusive. RNA-Seq on skin fibroblasts was then performed in nine individuals for whom GS was inconclusive and optical genome mapping (OGM) was performed in two individuals with an SV of unknown significance. RESULTS: We identified pathogenic or likely pathogenic variants in 16 individuals (48%) and six variants of uncertain significance. RNA-Seq contributed to the interpretation in three individuals, and OGM helped to characterise two SVs. CONCLUSION: Our study confirmed that GS significantly improves the diagnostic performance of NDDs. However, most variants detectable by GS alone are structural or located in non-coding regions, which can pose challenges for interpretation. Integration of RNA-Seq data overcame this limitation by confirming the impact of variants at the transcriptional or regulatory level. This result paves the way for new routinely applicable diagnostic protocols.

RLIM Is a Candidate Dosage-Sensitive Gene for Individuals with Varying Duplications of Xq13, Intellectual Disability, and Distinct Facial Features.

Interpretation of the significance of maternally inherited X chromosome variants in males with neurocognitive phenotypes continues to present a challenge to clinical geneticists and diagnostic laboratories. Here we report 14 males from 9 families with duplications at the Xq13.2-q13.3 locus with a common facial phenotype, intellectual disability (ID), distinctive behavioral features, and a seizure disorder in two cases. All tested carrier mothers had normal intelligence. The duplication arose de novo in three mothers where grandparental testing was possible. In one family the duplication segregated with ID across three generations. RLIM is the only gene common to our duplications. However, flanking genes duplicated in some but not all the affected individuals included the brain-expressed genes NEXMIF, SLC16A2, and the long non-coding RNA gene FTX. The contribution of the RLIM-flanking genes to the phenotypes of individuals with different size duplications has not been fully resolved. Missense variants in RLIM have recently been identified to cause X-linked ID in males, with heterozygous females typically having normal intelligence and highly skewed X chromosome inactivation. We detected consistent and significant increase of RLIM mRNA and protein levels in cells derived from seven affected males from five families with the duplication. Subsequent analysis of MDM2, one of the targets of the RLIM E3 ligase activity, showed consistent downregulation in cells from the affected males. All the carrier mothers displayed normal RLIM mRNA levels and had highly skewed X chromosome inactivation. We propose that duplications at Xq13.2-13.3 including RLIM cause a recognizable but mild neurocognitive phenotype in hemizygous males.

De Novo SOX6 Variants Cause a Neurodevelopmental Syndrome Associated with ADHD, Craniosynostosis, and Osteochondromas.

SOX6 belongs to a family of 20 SRY-related HMG-box-containing (SOX) genes that encode transcription factors controlling cell fate and differentiation in many developmental and adult processes. For SOX6, these processes include, but are not limited to, neurogenesis and skeletogenesis. Variants in half of the SOX genes have been shown to cause severe developmental and adult syndromes, referred to as SOXopathies. We here provide evidence that SOX6 variants also cause a SOXopathy. Using clinical and genetic data, we identify 19 individuals harboring various types of SOX6 alterations and exhibiting developmental delay and/or intellectual disability; the individuals are from 17 unrelated families. Additional, inconstant features include attention-deficit/hyperactivity disorder (ADHD), autism, mild facial dysmorphism, craniosynostosis, and multiple osteochondromas. All variants are heterozygous. Fourteen are de novo, one is inherited from a mosaic father, and four offspring from two families have a paternally inherited variant. Intragenic microdeletions, balanced structural rearrangements, frameshifts, and nonsense variants are predicted to inactivate the SOX6 variant allele. Four missense variants occur in residues and protein regions highly conserved evolutionarily. These variants are not detected in the gnomAD control cohort, and the amino acid substitutions are predicted to be damaging. Two of these variants are located in the HMG domain and abolish SOX6 transcriptional activity in vitro. No clear genotype-phenotype correlations are found. Taken together, these findings concur that SOX6 haploinsufficiency leads to a neurodevelopmental SOXopathy that often includes ADHD and abnormal skeletal and other features.

Exome sequencing as a first-tier test for copy number variant detection: retrospective evaluation and prospective screening in 2418 cases.

BACKGROUND: Despite the availability of whole exome (WES) and genome sequencing (WGS), chromosomal microarray (CMA) remains the first-line diagnostic test in most rare disorders diagnostic workup, looking for copy number variations (CNVs), with a diagnostic yield of 10%-20%. The question of the equivalence of CMA and WES in CNV calling is an organisational and economic question, especially when ordering a WGS after a negative CMA and/or WES. METHODS: This study measures the equivalence between CMA and GATK4 exome sequencing depth of coverage method in detecting coding CNVs on a retrospective cohort of 615 unrelated individuals. A prospective detection of WES-CNV on a cohort of 2418 unrelated individuals, including the 615 individuals from the validation cohort, was performed. RESULTS: On the retrospective validation cohort, every CNV detectable by the method (ie, a CNV with at least one exon not in a dark zone) was accurately called (64/64 events). In the prospective cohort, 32 diagnoses were performed among the 2418 individuals with CNVs ranging from 704 bp to aneuploidy. An incidental finding was reported. The overall increase in diagnostic yield was of 1.7%, varying from 1.2% in individuals with multiple congenital anomalies to 1.9% in individuals with chronic kidney failure. CONCLUSION: Combining single-nucleotide variant (SNV) and CNV detection increases the suitability of exome sequencing as a first-tier diagnostic test for suspected rare Mendelian disorders. Before considering the prescription of a WGS after a negative WES, a careful reanalysis with updated CNV calling and SNV annotation should be considered.

Terroir Influence on Polyphenol Metabolism from Grape Canes: A Spatial Metabolomic Study at Parcel Scale.

The composition of bioactive polyphenols from grape canes, an important viticultural byproduct, was shown to be varietal-dependent; however, the influence of soil-related terroir factors remains unexplored. Using spatial metabolomics and correlation-based networks, we investigated how continuous changes in soil features and topography may impact the polyphenol composition in grape canes. Soil properties, topography, and grape cane extracts were analyzed at georeferenced points over 3 consecutive years, followed by UPLC-DAD-MS-based metabolomic analysis targeting 42 metabolites. Principal component analyses on intra-vintage metabolomic data presented a good reproducibility in relation to geographic coordinates. A correlation-driven approach was used to explore the combined influence of soil and topographic variables on metabolomic responses. As a result, a metabolic cluster including flavonoids was correlated with elevation and curvature. Spatial metabolomics driven by correlation-based networks represents a powerful approach to spatialize field-omics data and may serve as new field-phenotyping tool in precision agriculture.

RPL13 Variants Cause Spondyloepimetaphyseal Dysplasia with Severe Short Stature.

Variants in genes encoding ribosomal proteins have thus far been associated with Diamond-Blackfan anemia, a rare inherited bone marrow failure, and isolated congenital asplenia. Here, we report one de novo missense variant and three de novo splice variants in RPL13, which encodes ribosomal protein RPL13 (also called eL13), in four unrelated individuals with a rare bone dysplasia causing severe short stature. The three splice variants (c.477+1G>T, c.477+1G>A, and c.477+2 T>C) result in partial intron retention, which leads to an 18-amino acid insertion. In contrast to observations from Diamond-Blackfan anemia, we detected no evidence of significant pre-rRNA processing disturbance in cells derived from two affected individuals. Consistently, we showed that the insertion-containing protein is stably expressed and incorporated into 60S subunits similar to the wild-type protein. Erythroid proliferation in culture and ribosome profile on sucrose gradient are modified, suggesting a change in translation dynamics. We also provide evidence that RPL13 is present at high levels in chondrocytes and osteoblasts in mouse growth plates. Taken together, we show that the identified RPL13 variants cause a human ribosomopathy defined by a rare skeletal dysplasia, and we highlight the role of this ribosomal protein in bone development.

Complex Compound Inheritance of Lethal Lung Developmental Disorders Due to Disruption of the TBX-FGF Pathway.

Primary defects in lung branching morphogenesis, resulting in neonatal lethal pulmonary hypoplasias, are incompletely understood. To elucidate the pathogenetics of human lung development, we studied a unique collection of samples obtained from deceased individuals with clinically and histopathologically diagnosed interstitial neonatal lung disorders: acinar dysplasia (n = 14), congenital alveolar dysplasia (n = 2), and other lethal lung hypoplasias (n = 10). We identified rare heterozygous copy-number variant deletions or single-nucleotide variants (SNVs) involving TBX4 (n = 8 and n = 2, respectively) or FGF10 (n = 2 and n = 2, respectively) in 16/26 (61%) individuals. In addition to TBX4, the overlapping ∼2 Mb recurrent and nonrecurrent deletions at 17q23.1q23.2 identified in seven individuals with lung hypoplasia also remove a lung-specific enhancer region. Individuals with coding variants involving either TBX4 or FGF10 also harbored at least one non-coding SNV in the predicted lung-specific enhancer region, which was absent in 13 control individuals with the overlapping deletions but without any structural lung anomalies. The occurrence of rare coding variants involving TBX4 or FGF10 with the putative hypomorphic non-coding SNVs implies a complex compound inheritance of these pulmonary hypoplasias. Moreover, they support the importance of TBX4-FGF10-FGFR2 epithelial-mesenchymal signaling in human lung organogenesis and help to explain the histopathological continuum observed in these rare lethal developmental disorders of the lung.

Rare and de novo duplications containing TCF20 are associated with a neurodevelopmental disorder.

Transcriptor co-activator factor 20 gene (TCF20) encodes a nuclear chromatin-binding protein involved in regulation of gene expression. In human pathology, pathogenic variants or deletions in TCF20 were identified in patients with developmental delay, variable intellectual disability and behavioral impairment (OMIM: 618430). The shared core phenotype includes developmental delay, hypotonia, motor delay, autism spectrum disorders, neurobehavioral anomalies, neurological features such as ataxia, seizures, movement disorders, structural brain anomalies, craniofacial features and various congenital anomalies. Most pathogenic variants are loss-of-function variants. Duplication including TCF20 was suspected to cause a neurodevelopmental disorder (NDD) with mirror traits compared to patients with TCF20 deletions. In the present study, we report three patients from three unrelated families with NDD with a de novo duplication at 22q13.2 encompassing TCF20. We propose that the TCF20 duplication could be involved in a new 22q13.2 microduplication syndrome with high penetrance, enlarging the genotype-phenotype knowledge of TCF20-associated NDDs.

Novel interstitial 2q12.3q13 microdeletion predisposes to developmental delay and behavioral problems.

Microarray-based comparative genomic hybridization (aCGH) is being increasingly applied to delineate novel genomic disorders and related syndromes in patients with developmental delay. In this study, detailed clinical and cytogenetic data of three unrelated patients with interstitial 2q12.3q13 microdeletion were described and compared with thirteen 2q12.3q13 microdeletion patients, gathered from the medical literature and public databases. 60 K aCGH analysis revealed three overlapping 2q12.3q13 microdeletions measuring 1.88 Mb in patient 1, 1.25 Mb in patient 2, and 0.41 Mb in patient 3, respectively. Confirmation and segregation studies were performed using fluorescence in situ hybridization (FISH) and quantitative real-time PCR. Variable clinical features of 2q12.3q13 microdeletion including microcephaly, prenatal growth retardation, developmental delay, short stature, behavioral problems, learning difficulties, skeletal anomalies, congenital heart defects, and features of ectodermal dysplasia were observed. The boundaries and sizes of the 2q12.3q13 deletions in the sixteen patients were different, but an overlapping region of 249 kb in 2q12.3 was defined. The SRO (smallest region of overlap) encompasses four genes, including LIMS1, RANBP2, CCDC138, and EDAR. Among these genes, RANBP2 is a strong candidate gene for neurological phenotype and genetic susceptibility to viral infections. To our knowledge, this is the first published report of 2q12.3q13 microdeletion syndrome and our observations strongly suggest that these recurrent CNVs may be a novel risk factor for developmental delay with variable expressivity and incomplete penetrance.

Whole genome copy number analysis in search of new prognostic biomarkers in first line treatment of mantle cell lymphoma. A study by the LYSA group.

Mantle cell lymphoma (MCL) is a lymphoproliferative disorder characterized by the t(11;14)(q13;q32) CCND1/IGH translocation. This lymphoma is however extremely heterogeneous in terms of molecular alterations. Moreover, the course of the disease can vary greatly between indolent forms with slow progression and aggressive conditions rapidly pejorative. The identification of early markers allowing to predict individual patients outcome has however been unsuccessful so far. The LyMa trial treated homogeneously a cohort of young MCL patients. This appeared as a good opportunity to search for biomarkers of response to therapy. DNA extracted from diagnostic paraffin-embedded lymph node biopsies from 100 patients with newly diagnosed MCL, homogeneously treated in this prospective clinical trial, were investigated for copy number alterations and copy neutral loss of heterozygosity using the Oncoscan SNP-array scanning the whole genome. An independent confirmatory cohort was used to strengthen the possibly relevant anomalies observed. Here we describe the recurrent anomalies identified with this technique. Deletions of 17p(TP53) and 9p(CDKN2A) were more frequent in refractory or early relapsing patients (10%), but had no significant impact in univariate analysis on progression-free (PFS) or overall survival (OS). Regardless of the presence of TP53 or CDKN2A deletions, gains in 7p22 (8,5%) were associated with better PFS in univariate but not in multivariate analysis including MCL International Prognostic Index and treatment. Gains of 11q(CCDN1), suggesting gains of the CCND1/IGH fusion, were associated with worse OS and PFS in univariate and multivariate analyses. This worse prognosis impact was confirmed by FISH in an independent confirmatory cohort. This work, using a whole genome approach, confirms the broad genomic landscape of MCL and shows that gains of the CCND1/IGH fusion can be considered as a new prognostic structural variant. Genomic abnormalities of prognostic impact could be useful to strengthen or de-escalate treatment schedules or choosing targeted therapies or CART-cells.

Identification of mobile retrocopies during genetic testing: Consequences for routine diagnosis.

Human retrocopies, that is messenger RNA transcripts benefitting from the long interspersed element 1 machinery for retrotransposition, may have specific consequences for genomic testing. Next genetration sequencing (NGS) techniques allow the detection of such mobile elements but they may be misinterpreted as genomic duplications or be totally overlooked. We report eight observations of retrocopies detected during diagnostic NGS analyses of targeted gene panels, exome, or genome sequencing. For seven cases, while an exons-only copy number gain was called, read alignment inspection revealed a depth of coverage shift at every exon-intron junction where indels were also systematically called. Moreover, aberrant chimeric read pairs spanned entire introns or were paired with another locus for terminal exons. The 8th retrocopy was present in the reference genome and thus showed a normal NGS profile. We emphasize the existence of retrocopies and strategies to accurately detect them at a glance during genetic testing and discuss pitfalls for genetic testing.

Pseudodicentric Chromosome Originating from Autosomes 9 and 21 in a Male Patient with Oligozoospermia.

Genetic factors are responsible for 15% of male infertility conditions. Numerical and structural chromosomal anomalies (related to the Y chromosome or to the autosomes) are validated genetic factors leading to spermatogenic quantitative defects with a frequency depending on the severity of the phenotype. The most frequent structural chromosomal rearrangements of autosomes are translocations and inversions, whereas dicentric chromosomes involving autosomes are rare. We report a man bearing a pseudodicentric chromosome (9;21) and presenting with oligozoospermia. Extensive cytogenetic analyses were necessary to determine the precise nature of the derivative chromosome and to discount the presence of interstitial telomeric sequences. Defects in spermatogenesis and abnormal segregation at meiosis for existing spermatozoa are proposed and are the likely cause of the reproductive phenotype of the patient.

Correction: Variants in MED12L, encoding a subunit of the Mediator kinase module, are responsible for intellectual disability associated with transcriptional defect.

In the Acknowledgements section of the paper the authors neglected to mention that the study was supported by a grant from the National Human Genome Research Institute (NHGRI) UM1HG007301 (S.H., M.L.T.). In addition, the award of MD was associated with the authors Michelle L. Thompson and Susan Hiatt instead of PhD. The PDF and HTML versions of the Article have been modified accordingly.

Variants in MED12L, encoding a subunit of the mediator kinase module, are responsible for intellectual disability associated with transcriptional defect.

PURPOSE: Mediator is a multiprotein complex that allows the transfer of genetic information from DNA binding proteins to the RNA polymerase II during transcription initiation. MED12L is a subunit of the kinase module, which is one of the four subcomplexes of the mediator complex. Other subunits of the kinase module have been already implicated in intellectual disability, namely MED12, MED13L, MED13, and CDK19. METHODS: We describe an international cohort of seven affected individuals harboring variants involving MED12L identified by array CGH, exome or genome sequencing. RESULTS: All affected individuals presented with intellectual disability and/or developmental delay, including speech impairment. Other features included autism spectrum disorder, aggressive behavior, corpus callosum abnormality, and mild facial morphological features. Three individuals had a MED12L deletion or duplication. The other four individuals harbored single-nucleotide variants (one nonsense, one frameshift, and two splicing variants). Functional analysis confirmed a moderate and significant alteration of RNA synthesis in two individuals. CONCLUSION: Overall data suggest that MED12L haploinsufficiency is responsible for intellectual disability and transcriptional defect. Our findings confirm that the integrity of this kinase module is a critical factor for neurological development.

Rare variants in the genetic background modulate cognitive and developmental phenotypes in individuals carrying disease-associated variants.

PURPOSE: To assess the contribution of rare variants in the genetic background toward variability of neurodevelopmental phenotypes in individuals with rare copy-number variants (CNVs) and gene-disruptive variants. METHODS: We analyzed quantitative clinical information, exome sequencing, and microarray data from 757 probands and 233 parents and siblings who carry disease-associated variants. RESULTS: The number of rare likely deleterious variants in functionally intolerant genes ("other hits") correlated with expression of neurodevelopmental phenotypes in probands with 16p12.1 deletion (n=23, p=0.004) and in autism probands carrying gene-disruptive variants (n=184, p=0.03) compared with their carrier family members. Probands with 16p12.1 deletion and a strong family history presented more severe clinical features (p=0.04) and higher burden of other hits compared with those with mild/no family history (p=0.001). The number of other hits also correlated with severity of cognitive impairment in probands carrying pathogenic CNVs (n=53) or de novo pathogenic variants in disease genes (n=290), and negatively correlated with head size among 80 probands with 16p11.2 deletion. These co-occurring hits involved known disease-associated genes such as SETD5, AUTS2, and NRXN1, and were enriched for cellular and developmental processes. CONCLUSION: Accurate genetic diagnosis of complex disorders will require complete evaluation of the genetic background even after a candidate disease-associated variant is identified.

11q24.2q24.3 microdeletion in two families presenting features of Jacobsen syndrome, without intellectual disability: Role of FLI1, ETS1, and SENCR long noncoding RNA.

This report presents two families with interstitial 11q24.2q24.3 deletion, associated with malformations, hematologic features, and typical facial dysmorphism, observed in Jacobsen syndrome (JS), except for intellectual disability (ID). The smallest 700 Kb deletion contains only two genes: FLI1 and ETS1, and a long noncoding RNA, SENCR, narrowing the minimal critical region for some features of JS. Consistent with recent literature, it adds supplemental data to confirm the crucial role of FLI1 and ETS1 in JS, namely FLI1 in thrombocytopenia and ETS1 in cardiopathy and immune deficiency. It also supports that combined ETS1 and FLI1 haploinsufficiency explains dysmorphic features, notably ears, and nose anomalies. Moreover, it raises the possibility that SENCR, a long noncoding RNA, could be responsible for limb defects, because of its early role in endothelial cell commitment and function. Considering ID and autism spectrum disorder, which are some of the main features of JS, a participation of ETS1, FLI1, or SENCR cannot be excluded. But, considering the normal neurodevelopment of our patients, their role would be either minor or with an important variability in penetrance. Furthermore, according to literature, ARHGAP32 and KIRREL3 seem to be the strongest candidate genes in the 11q24 region for other Jacobsen patients.

Risk estimation of uniparental disomy of chromosome 14 or 15 in a fetus with a parent carrying a non-homologous Robertsonian translocation. Should we still perform prenatal diagnosis?

OBJECTIVE: Uniparental disomy (UPD) testing is currently recommended during pregnancy in fetuses carrying a balanced Robertsonian translocation (ROB) involving chromosome 14 or 15, both chromosomes containing imprinted genes. The overall risk that such a fetus presents a UPD has been previously estimated to be around ~0.6-0.8%. However, because UPD are rare events and this estimate has been calculated from a number of studies of limited size, we have reevaluated the risk of UPD in fetuses for whom one of the parents was known to carry a nonhomologous ROB (NHROB). METHOD: We focused our multicentric study on NHROB involving chromosome 14 and/or 15. A total of 1747 UPD testing were performed in fetuses during pregnancy for the presence of UPD(14) and/or UPD(15). RESULT: All fetuses were negative except one with a UPD(14) associated with a maternally inherited rob(13;14). CONCLUSION: Considering these data, the risk of UPD following prenatal diagnosis of an inherited ROB involving chromosome 14 and/or 15 could be estimated to be around 0.06%, far less than the previous estimation. Importantly, the risk of miscarriage following an invasive prenatal sampling is higher than the risk of UPD. Therefore, we do not recommend prenatal testing for UPD for these pregnancies and parents should be reassured.

Familial autosomal dominant severe ankyloglossia with tooth abnormalities.

Ankyloglossia is a congenital oral anomaly characterized by the presence of a hypertrophic and short lingual frenulum. Mutations in the gene encoding the transcription factor TBX22 have been involved in isolated ankyloglossia and X-linked cleft palate. The knockout of Lgr5 in mice results in ankyloglossia. Here, we report a five-generation family including patients with severe ankyloglossia and missing lower central incisors. Two members of this family also exhibited congenital anorectal malformations. In this report, male-to-male transmission was in favor of an autosomal dominant inheritance, which allowed us to exclude the X-linked TBX22 gene. Linkage analysis using short tandem repeat markers located in the vicinity of LGR5 excluded this gene as a potential candidate. These results indicate genetic heterogeneity for ankyloglossia. Further investigations with additional families are required in order to identify novel candidate genes.

Comparison of Tumor- and Bone Marrow-Derived Mesenchymal Stromal/Stem Cells from Patients with High-Grade Osteosarcoma.

Osteosarcoma (OS) is suspected to originate from dysfunctional mesenchymal stromal/stem cells (MSC). We sought to identify OS-derived cells (OSDC) with potential cancer stem cell (CSC) properties by comparing OSDC to MSC derived from bone marrow of patients. This study included in vitro characterization with sphere forming assays, differentiation assays, cytogenetic analysis, and in vivo investigations of their tumorigenicity and tumor supportive capacities. Primary cell lines were isolated from nine high-grade OS samples. All primary cell lines demonstrated stromal cell characteristics. Compared to MSC, OSDC presented a higher ability to form sphere clones, indicating a potential CSC phenotype, and were more efficient at differentiation towards osteoblasts. None of the OSDC displayed the complex chromosome rearrangements typical of high grade OS and none of them induced tumors in immunodeficient mice. However, two OSDC demonstrated focused genomic abnormalities. Three out of seven, and six out of seven OSDC showed a supportive role on local tumor development, and on metastatic progression to the lungs, respectively, when co-injected with OS cells in nude mice. The observation of OS-associated stromal cells with rare genetic abnormalities and with the capacity to sustain tumor progression may have implications for future tumor treatments.