Combinatorial patterns of gene expression changes contribute to variable expressivity of the developmental delay-associated 16p12.1 deletion

Genome Med. 2021 Oct 18;13(1):163. doi: 10.1186/s13073-021-00982-z.

Abstract

Background: Recent studies have suggested that individual variants do not sufficiently explain the variable expressivity of phenotypes observed in complex disorders. For example, the 16p12.1 deletion is associated with developmental delay and neuropsychiatric features in affected individuals, but is inherited in > 90% of cases from a mildly-affected parent. While children with the deletion are more likely to carry additional "second-hit" variants than their parents, the mechanisms for how these variants contribute to phenotypic variability are unknown.

Methods: We performed detailed clinical assessments, whole-genome sequencing, and RNA sequencing of lymphoblastoid cell lines for 32 individuals in five large families with multiple members carrying the 16p12.1 deletion. We identified contributions of the 16p12.1 deletion and "second-hit" variants towards a range of expression changes in deletion carriers and their family members, including differential expression, outlier expression, alternative splicing, allele-specific expression, and expression quantitative trait loci analyses.

Results: We found that the deletion dysregulates multiple autism and brain development genes such as FOXP1, ANK3, and MEF2. Carrier children also showed an average of 5323 gene expression changes compared with one or both parents, which matched with 33/39 observed developmental phenotypes. We identified significant enrichments for 13/25 classes of "second-hit" variants in genes with expression changes, where 4/25 variant classes were only enriched when inherited from the noncarrier parent, including loss-of-function SNVs and large duplications. In 11 instances, including for ZEB2 and SYNJ1, gene expression was synergistically altered by both the deletion and inherited "second-hits" in carrier children. Finally, brain-specific interaction network analysis showed strong connectivity between genes carrying "second-hits" and genes with transcriptome alterations in deletion carriers.

Conclusions: Our results suggest a potential mechanism for how "second-hit" variants modulate expressivity of complex disorders such as the 16p12.1 deletion through transcriptomic perturbation of gene networks important for early development. Our work further shows that family-based assessments of transcriptome data are highly relevant towards understanding the genetic mechanisms associated with complex disorders.

Keywords: Complex disorders; Copy-number variant; Developmental disorders; Inherited variants; RNA sequencing; Whole-genome sequencing.

Publication types

  • Research Support, N.I.H., Extramural
  • Research Support, Non-U.S. Gov't

MeSH terms

  • Ankyrins / genetics
  • Autistic Disorder / genetics
  • Biological Variation, Population*
  • Brain
  • Chromosome Deletion*
  • Exome Sequencing
  • Family
  • Forkhead Transcription Factors / genetics
  • Gene Expression*
  • Humans
  • Phenotype
  • Phosphoric Monoester Hydrolases / genetics
  • Repressor Proteins / genetics
  • Transcription Factors / genetics
  • Whole Genome Sequencing
  • Zinc Finger E-box Binding Homeobox 2 / genetics

Substances

  • ANK3 protein, human
  • Ankyrins
  • FOXP1 protein, human
  • Forkhead Transcription Factors
  • Repressor Proteins
  • Transcription Factors
  • ZEB2 protein, human
  • Zinc Finger E-box Binding Homeobox 2
  • Phosphoric Monoester Hydrolases
  • phosphoinositide 5-phosphatase