Downregulation of genes outside the deleted region in individuals with 22q11.2 deletion syndrome

Next-generation sequencing profiling of miRNAs in individuals with 22q11.2 deletion syndrome revealed altered expression of miR-185-5p

BACKGROUND

The 22q11.2 deletion syndrome (22q11.2DS) is a microdeletion syndrome with highly variable phenotypic manifestations, even though most patients present the typical 3 Mb microdeletion, usually affecting the same ~ 106 genes. One of the genes affected by this deletion is DGCR8, which plays a crucial role in miRNA biogenesis. Therefore, the haploinsufficiency of DGCR8 due to this microdeletion can alter the modulation of the expression of several miRNAs involved in a range of biological processes.

RESULTS

In this study, we used next-generation sequencing to evaluate the miRNAs profiles in the peripheral blood of 12 individuals with typical 22q11DS compared to 12 healthy matched controls. We used the DESeq2 package for differential gene expression analysis and the DIANA-miTED dataset to verify the expression of differentially expressed miRNAs in other tissues. We used miRWalk to predict the target genes of differentially expressed miRNAs. Here, we described two differentially expressed miRNAs in patients compared to controls: hsa-miR-1304-3p, located outside the 22q11.2 region, upregulated in patients, and hsa-miR-185-5p, located in the 22q11.2 region, which showed downregulation. Expression of miR-185-5p is observed in tissues frequently affected in patients with 22q11DS, and previous studies have reported its downregulation in individuals with 22q11DS. hsa-miR-1304-3p has low expression in blood and, thus, needs more validation, though using a sensitive technology allowed us to identify differences in expression between patients and controls.

CONCLUSIONS

Thus, lower expression of miR-185-5p can be related to the 22q11.2 deletion and DGCR8 haploinsufficiency, leading to phenotypic consequences in 22q11.2DS patients, while higher expression of hsa-miR-1304-3p might be related to individual genomic variances due to the heterogeneous background of the Brazilian population.

Variants in Candidate Genes for Phenotype Heterogeneity in Patients with the 22q11.2 Deletion Syndrome

22q11.2 deletion syndrome (22q11.2DS) is a microdeletion syndrome with a broad and heterogeneous phenotype, even though most of the deletions present similar sizes, involving ∼3 Mb of DNA. In a relatively large population of a Brazilian 22q11.2DS cohort (60 patients), we investigated genetic variants that could act as genetic modifiers and contribute to the phenotypic heterogeneity, using a targeted NGS (Next Generation Sequencing) with a specific Ion AmpliSeq panel to sequence nine candidate genes (CRKL, MAPK1, HIRA, TANGO2, PI4KA, HDAC1, ZDHHC8, ZFPM2, and JAM3), mapped in and outside the 22q11.2 hemizygous deleted region. In silico prediction was performed, and the whole-genome sequencing annotation analysis package (WGSA) was used to predict the possible pathogenic effect of single nucleotide variants (SNVs). For the in silico prediction of the indels, we used the genomic variants filtered by a deep learning model in NGS (GARFIELD-NGS). We identified six variants, 4 SNVs and 2 indels, in MAPK1, JAM3, and ZFPM2 genes with possibly synergistic deleterious effects in the context of the 22q11.2 deletion. Our results provide the opportunity for the discovery of the co-occurrence of genetic variants with 22q11.2 deletions, which may influence the patients´ phenotype.

Coexisting Conditions Modifying Phenotypes of Patients with 22q11.2 Deletion Syndrome

22q11.2 deletion syndrome (22q11.2DS) is the most common genomic disorder with an extremely broad phenotypic spectrum. The aim of our study was to investigate how often the additional variants in the genome can affect clinical variation among patients with the recurrent deletion. To examine the presence of additional variants affecting the phenotype, we performed microarray in 82 prenatal and 77 postnatal cases and performed exome sequencing in 86 postnatal patients with 22q11.2DS. Within those 159 patients where array was performed, 5 pathogenic and 5 likely pathogenic CNVs were identified outside of the 22q11.2 region. This indicates that in 6.3% cases, additional CNVs most likely contribute to the clinical presentation. Additionally, exome sequencing in 86 patients revealed 3 pathogenic (3.49%) and 5 likely pathogenic (5.81%) SNVs and small CNV. These results show that the extension of diagnostics with genome-wide methods can reveal other clinically relevant changes in patients with 22q11 deletion syndrome.

The 22q11.2 region regulates presynaptic gene-products linked to schizophrenia

It is unclear how the 22q11.2 deletion predisposes to psychiatric disease. To study this, we generated induced pluripotent stem cells from deletion carriers and controls and utilized CRISPR/Cas9 to introduce the heterozygous deletion into a control cell line. Here, we show that upon differentiation into neural progenitor cells, the deletion acted in trans to alter the abundance of transcripts associated with risk for neurodevelopmental disorders including autism. In excitatory neurons, altered transcripts encoded presynaptic factors and were associated with genetic risk for schizophrenia, including common and rare variants. To understand how the deletion contributed to these changes, we defined the minimal protein-protein interaction network that best explains gene expression alterations. We found that many genes in 22q11.2 interact in presynaptic, proteasome, and JUN/FOS transcriptional pathways. Our findings suggest that the 22q11.2 deletion impacts genes that may converge with psychiatric risk loci to influence disease manifestation in each deletion carrier.

Positively selected genes in the hoary bat (Lasiurus cinereus) lineage: prominence of thymus expression, immune and metabolic function, and regions of ancient synteny

Background

Bats of the genus Lasiurus occur throughout the Americas and have diversified into at least 20 species among three subgenera. The hoary bat (Lasiurus cinereus) is highly migratory and ranges farther across North America than any other wild mammal. Despite the ecological importance of this species as a major insect predator, and the particular susceptibility of lasiurine bats to wind turbine strikes, our understanding of hoary bat ecology, physiology, and behavior remains poor.

Methods

To better understand adaptive evolution in this lineage, we used whole-genome sequencing to identify protein-coding sequence and explore signatures of positive selection. Gene models were predicted with Maker and compared to seven well-annotated and phylogenetically representative species. Evolutionary rate analysis was performed with PAML.

Results

Of 9,447 single-copy orthologous groups that met evaluation criteria, 150 genes had a significant excess of nonsynonymous substitutions along the L. cinereus branch (P < 0.001 after manual review of alignments). Selected genes as a group had biased expression, most strongly in thymus tissue. We identified 23 selected genes with reported immune functions as well as a divergent paralog of Steep1 within suborder Yangochiroptera. Seventeen genes had roles in lipid and glucose metabolic pathways, partially overlapping with 15 mitochondrion-associated genes; these adaptations may reflect the metabolic challenges of hibernation, long-distance migration, and seasonal variation in prey abundance. The genomic distribution of positively selected genes differed significantly from background expectation by discrete Kolmogorov-Smirnov test (P < 0.001). Remarkably, the top three physical clusters all coincided with islands of conserved synteny predating Mammalia, the largest of which shares synteny with the human cat-eye critical region (CECR) on 22q11. This observation coupled with the expansion of a novel Tbx1-like gene family may indicate evolutionary innovation during pharyngeal arch development: both the CECR and Tbx1 cause dosage-dependent congenital abnormalities in thymus, heart, and head, and craniodysmorphy is associated with human orthologs of other positively selected genes as well.

Integration of genetic, transcriptomic, and clinical data provides insight into 16p11.2 and 22q11.2 CNV genes

BACKGROUND

Deletions and duplications of the multigenic 16p11.2 and 22q11.2 copy number variant (CNV) regions are associated with brain-related disorders including schizophrenia, intellectual disability, obesity, bipolar disorder, and autism spectrum disorder (ASD). The contribution of individual CNV genes to each of these identified phenotypes is unknown, as well as the contribution of these CNV genes to other potentially subtler health implications for carriers. Hypothesizing that DNA copy number exerts most effects via impacts on RNA expression, we attempted a novel in silico fine-mapping approach in non-CNV carriers using both GWAS and biobank data.

METHODS

We first asked whether gene expression level in any individual gene in the CNV region alters risk for a known CNV-associated behavioral phenotype(s). Using transcriptomic imputation, we performed association testing for CNV genes within large genotyped cohorts for schizophrenia, IQ, BMI, bipolar disorder, and ASD. Second, we used a biobank containing electronic health data to compare the medical phenome of CNV carriers to controls within 700,000 individuals in order to investigate the full spectrum of health effects of the CNVs. Third, we used genotypes for over 48,000 individuals within the biobank to perform phenome-wide association studies between imputed expressions of individual 16p11.2 and 22q11.2 genes and over 1500 health traits.

RESULTS

Using large genotyped cohorts, we found individual genes within 16p11.2 associated with schizophrenia (TMEM219, INO80E, YPEL3), BMI (TMEM219, SPN, TAOK2, INO80E), and IQ (SPN), using conditional analysis to identify upregulation of INO80E as the driver of schizophrenia, and downregulation of SPN and INO80E as increasing BMI. We identified both novel and previously observed over-represented traits within the electronic health records of 16p11.2 and 22q11.2 CNV carriers. In the phenome-wide association study, we found seventeen significant gene-trait pairs, including psychosis (NPIPB11, SLX1B) and mood disorders (SCARF2), and overall enrichment of mental traits.

CONCLUSIONS

Our results demonstrate how integration of genetic and clinical data aids in understanding CNV gene function and implicates pleiotropy and multigenicity in CNV biology.

Genome-wide study of copy number variation implicates multiple novel loci for schizophrenia risk in Han Chinese family trios

Schizophrenia (SCZ) is a severe neuropsychiatric disorder that affects 1% of the global population. Copy number variations (CNVs) have been shown to play a critical role in its pathophysiology; however, only case-control studies on SCZ susceptibility CNVs have been conducted in Han Chinese. Here, we performed an array comparative genomic hybridization-based genome-wide CNV analysis in 100 Chinese family trios with SCZ. Burden test suggested that the SCZ probands carried more duplications than their healthy parents and unrelated healthy controls. Besides, five CNV loci were firstly reported to be associated with SCZ here, including both unbalanced transmitted CNVs and enriched de novo CNVs. Moreover, two genes (CTDSPL and MGAM) in these CNVs showed significant SCZ relevance in the expression level. Our findings support the crucial role of CNVs in the etiology of SCZ and provide new insights into the underlying mechanism of SCZ pathogenesis.