NCBP2 modulates neurodevelopmental defects of the 3q29 deletion in Drosophila and Xenopus laevis models

3q29 duplications: A cohort of 46 patients and a literature review

Duplications of the 3q29 cytoband are rare chromosomal copy number variations (CNVs) (overlapping or recurrent ~1.6 Mb 3q29 duplications). They have been associated with highly variable neurodevelopmental disorders (NDDs) with various associated features or reported as a susceptibility factor to the development of learning disabilities and neuropsychiatric disorders. The smallest region of overlap and the phenotype of 3q29 duplications remain uncertain. We here report a French cohort of 31 families with a 3q29 duplication identified by chromosomal microarray analysis (CMA), including 14 recurrent 1.6 Mb duplications, eight overlapping duplications (>1 Mb), and nine small duplications (<1 Mb). Additional genetic findings that may be involved in the phenotype were identified in 11 patients. Focusing on apparently isolated 3q29 duplications, patients present mainly mild NDD as suggested by a high rate of learning disabilities in contrast to a low proportion of patients with intellectual disabilities. Although some are de novo, most of the 3q29 duplications are inherited from a parent with a similar mild phenotype. Besides, the study of small 3q29 duplications does not provide evidence for any critical region. Our data suggest that the overlapping and recurrent 3q29 duplications seem to lead to mild NDD and that a severe or syndromic clinical presentation should warrant further genetic analyses.

Cross-species analysis identifies mitochondrial dysregulation as a functional consequence of the schizophrenia-associated 3q29 deletion

The 1.6-megabase deletion at chromosome 3q29 (3q29Del) is the strongest identified genetic risk factor for schizophrenia, but the effects of this variant on neurodevelopment are not well understood. We interrogated the developing neural transcriptome in two experimental model systems with complementary advantages: isogenic human cortical organoids and isocortex from the 3q29Del mouse model. We profiled transcriptomes from isogenic cortical organoids that were aged for 2 and 12 months, as well as perinatal mouse isocortex, all at single-cell resolution. Systematic pathway analysis implicated dysregulation of mitochondrial function and energy metabolism. These molecular signatures were supported by analysis of oxidative phosphorylation protein complex expression in mouse brain and assays of mitochondrial function in engineered cell lines, which revealed a lack of metabolic flexibility and a contribution of the 3q29 gene PAK2. Together, these data indicate that metabolic disruption is associated with 3q29Del and is conserved across species.

Genetics and Gene Editing Methods in Xenopus laevis and Xenopus tropicalis

Our understanding of biological systems has for many years been heavily influenced by experimental approaches that exploit genetic methods. These include gain-of-function experiments that overexpress transgenes or ectopically express injected RNA and loss-of-function experiments that knock out genes or knock down RNAs. Here, we review how these methods have been applied in Xenopus frogs and introduce a variety of protocols for genetic manipulation of Xenopus laevis and Xenopus tropicalis.

Cross-species transcriptomic analysis identifies mitochondrial dysregulation as a functional consequence of the schizophrenia-associated 3q29 deletion

Recent advances in the genetics of schizophrenia (SCZ) have identified rare variants that confer high disease risk, including a 1.6 Mb deletion at chromosome 3q29 with a staggeringly large effect size (O.R. > 40). Understanding the impact of the 3q29 deletion (3q29Del) on the developing CNS may therefore lead to insights about the pathobiology of schizophrenia. To gain clues about the molecular and cellular perturbations caused by the 3q29 deletion, we interrogated transcriptomic effects in two experimental model systems with complementary advantages: isogenic human forebrain cortical organoids and isocortex from the 3q29Del mouse model. We first created isogenic lines by engineering the full 3q29Del into an induced pluripotent stem cell line from a neurotypical individual. We profiled transcriptomes from isogenic cortical organoids that were aged for 2 months and 12 months, as well as day p7 perinatal mouse isocortex, all at single cell resolution. Differential expression analysis by genotype in each cell-type cluster revealed that more than half of the differentially expressed genes identified in mouse cortex were also differentially expressed in human cortical organoids, and strong correlations were observed in mouse-human differential gene expression across most major cell-types. We systematically filtered differentially expressed genes to identify changes occurring in both model systems. Pathway analysis on this filtered gene set implicated dysregulation of mitochondrial function and energy metabolism, although the direction of the effect was dependent on developmental timepoint. Transcriptomic changes were validated at the protein level by analysis of oxidative phosphorylation protein complexes in mouse brain tissue. Assays of mitochondrial function in human heterologous cells further confirmed robust mitochondrial dysregulation in 3q29Del cells, and these effects are partially recapitulated by ablation of the 3q29Del gene PAK2 . Taken together these data indicate that metabolic disruption is associated with 3q29Del and is conserved across species. These results converge with data from other rare SCZ-associated variants as well as idiopathic schizophrenia, suggesting that mitochondrial dysfunction may be a significant but overlooked contributing factor to the development of psychotic disorders. This cross-species scRNA-seq analysis of the SCZ-associated 3q29 deletion reveals that this copy number variant may produce early and persistent changes in cellular metabolism that are relevant to human neurodevelopment.

An artificial neural network model to diagnose non-obstructive azoospermia based on RNA-binding protein-related genes

Non-obstructive azoospermia (NOA) is a severe form of male infertility, but its pathological mechanisms and diagnostic biomarkers remain obscure. Since the dysregulation of RNA-binding proteins (RBPs) had nonnegligible effects on spermatogenesis, we aimed to investigate the functions and diagnosis values of RBPs in NOA. 58 testicular samples (control = 11, NOA = 47) from Gene Expression Omnibus (GEO) were set as the training cohort. Three public datasets, containing GSE45885 (control = 4, NOA = 27), GSE45887 (control = 4, NOA = 16), and GSE145467 (control = 10, NOA = 10), and 44 clinical samples from the local hospital (control = 27, NOA = 17) were used for validation. Through a series of bioinformatical analyses and machine learning algorithms, including genomic difference detection, protein-protein interaction network analysis, LASSO, SVM-RFE, and Boruta, DDX20 and NCBP2 were determined as significant predictors of NOA. Single-cell RNA sequencing of 432 testicular cell samples from NOA patients indicated that DDX20 and NCBP2 were associated with spermatogenesis (false discovery rate < 0.05). Based on the transcriptome expressions of DDX20 and NCBP2, we constructed multiple diagnosis models using logistic regression, random forest, and artificial neural network (ANN). The ANN model exhibited the most reliable predictive performance in the training cohort (AUC = 0.840), GSE45885 (AUC = 0.731), GSE45887 (AUC = 0.781), GSE145467 (AUC = 0.850), and local cohort (AUC = 0.623). Totally, an ANN diagnosis model based on RBP DDX20 and NCBP2 was developed and externally validated in NOA, functioning as a promising tool in clinical practice.

Establishment of a schizophrenia classifier based on peripheral blood signatures and investigation of pathogenic miRNA-mRNA regulation

To date, the diagnosis of schizophrenia (SCZ) mainly relies on patients' or guardians' self-reports and clinical observation, and the pathogenesis of SCZ remains elusive. In this study, we sought to develop a reliable classifier for diagnosing SCZ patients and provide clues to the etiology and pathogenesis of SCZ. Based on the high throughput sequencing analysis of peripheral blood miRNA expression profile and weighted gene co-expression network analysis (WGCNA) in our previous study, we selected eleven hub miRNAs for validation by qRT-PCR in 51 SCZ patients and 51 controls. miR-939-5p, miR-4732-3p let-7d-3p, and miR-142-3p were confirmed to be significantly up-regulated, and miR-30e-3p and miR-23a-3p were down-regulated in SCZ patients. miR-30e-3p with the most considerable fold change and statistically significance was selected for targeting validation. We first performed bioinformatics prediction followed by qRT-PCR and verified the up-regulation of potential target mRNAs (ABI1, NMT1, HMGB1) expression. Next, we found that the expression level of ABI1 was significantly up-regulated in SH-SY5Y cells transfected with miR-30e-3p mimics. Lastly, we conducted a luciferase assay in 293T cells confirming that miR-30e-3p could directly bind with the 3'untranslated region (3'-UTR) of ABI1, revealing that miR-30e-3p might play a role in the polymerization of neuronal actin and the reconstruction of the cytoskeleton via the downstream regulation of ABI1. In addition, we constructed a classifier by a series of bioinformatics algorithms and evaluated its diagnostic performance. It appears that the classifier consists of miRNAs and mRNAs possess a better discrimination performance than individual miRNA or mRNA in SCZ.

Identification of potential pathways and biomarkers linked to progression in ALS

OBJECTIVE

To identify potential diagnostic and prognostic biomarkers for clinical management and clinical trials in amyotrophic lateral sclerosis.

METHODS

We analysed proteomics data of ALS patient-induced pluripotent stem cell-derived motor neurons available through the AnswerALS consortium. After stratifying patients using clinical ALSFRS-R and ALS-CBS scales, we identified differentially expressed proteins indicative of ALS disease severity and progression rate as candidate ALS-related and prognostic biomarkers. Pathway analysis for identified proteins was performed using STITCH. Protein sets were correlated with the effects of drugs using the Connectivity Map tool to identify compounds likely to affect similar pathways. RNAi screening was performed in a Drosophila TDP-43 ALS model to validate pathological relevance. A statistical classification machine learning model was constructed using ridge regression that uses proteomics data to differentiate ALS patients from controls.

RESULTS

We identified 76, 21, 71 and 1 candidate ALS-related biomarkers and 22, 41, 27 and 64 candidate prognostic biomarkers from patients stratified by ALSFRS-R baseline, ALSFRS-R progression slope, ALS-CBS baseline and ALS-CBS progression slope, respectively. Nineteen proteins enhanced or suppressed pathogenic eye phenotypes in the ALS fly model. Nutraceuticals, dopamine pathway modulators, statins, anti-inflammatories and antimicrobials were predicted starting points for drug repurposing using the connectivity map tool. Ten diagnostic biomarker proteins were predicted by machine learning to identify ALS patients with high accuracy and sensitivity.

INTERPRETATION

This study showcases the powerful approach of iPSC-motor neuron proteomics combined with machine learning and biological confirmation in the prediction of novel mechanisms and diagnostic and predictive biomarkers in ALS.

Case report: 7p22.3 deletion and 8q24.3 duplication in a patient with epilepsy and psychomotor delay-Does both possibly act to modulate a candidate gene region for the patient's phenotype?

Background: Psychomotor delay, epilepsy and dysmorphic features are clinical signs which are described in multiple syndromes due to chromosomal imbalances or mutations involving key genes implicated in the stages of Early Embryonic Development. In this context, we report a 10 years old Tunisian patient with these three signs. Our objective is to determine the cause of developmental, behavioral and facial abnormalities in this patient. Methods: We used banding cytogenetics (karyotype) and Array Comparative Genomic Hybridization (Array CGH) to this purpose. Results: The karyotype was in favor of a derivative of chromosome 7 in the patient and Array CGH analysis revealed a loss of genetic material in 7p22.3-p22.1 (4,56 Mb) with a gain at 8q24.23-q24 (9.20 Mb) resulting from maternal 7/8 reciprocal translocation. An in silico analysis of the unbalanced region was carried out and showed that the 7p22.3-p22.1 deletion contains eight genes. Among them, BRAT1 gene, previously described in several neurodevelopmental diseases, may be a candidate gene which absence could be correlated to the patient's phenotype. However, the 8q24.23-q24 duplication could be involved in the phenotype of this patient. Conclusion: In this study, we report for the first time a 7p deletion/8q duplication in a patient with psychomoteur delay, epilepsy and facial dysmorphism. Our study showed that Array CGH still useful for delivering a conclusive genetic diagnosis for patients having neurodevelopmental abnormalities in the era of next-generation sequencing.

Experimental verification and comprehensive analysis of m7G methylation regulators in the subcluster classification of ischemic stroke

Background: Ischemic stroke (IS) is a fatal cerebrovascular disease involving several pathological mechanisms. Modification of 7-methylguanosine (m7G) has multiple regulatory functions. However, the expression pattern and mechanism of m7G in IS remain unknown. Herein, we aimed to explore the effect of m7G modification on IS. Methods: We screened significantly different m7G-regulated genes in Gene Expression Omnibus datasets, GSE58294 and GSE22255. The random forest (RF) algorithm was selected to identify key m7G-regulated genes that were subsequently validated using the middle cerebral artery occlusion (MCAO) model and quantitative polymerase chain reaction (qPCR). A risk model was subsequently generated using key m7G-regulated genes. Then, "ConsensusClusterPlus" package was used to distinguish different m7G clusters of patients with IS. Simultaneously, between two m7G clusters, differentially expressed genes (DEGs) and immune infiltration differences were also explored. Finally, we investigated functional enrichment and the mRNA-miRNA-transcription factor network of DEGs. Results: RF and qPCR confirmed that EIF3D, CYFIP2, NCBP2, DCPS, and NUDT1 were key m7G-related genes in IS that could accurately predict clinical risk (area under the curve = 0.967). NCBP2 was the most significantly associated gene with immune infiltration. Based on the expression profiles of these key m7G-related genes, the IS group could be divided into two clusters. According to the single-sample gene set enrichment analysis algorithm, four types of immune cells (immature dendritic cells, macrophages, natural killer T cells, and TH1 cells) were significantly different in the two m7G clusters. The functional enrichment of 282 DEGs between the two clusters was mainly concentrated in the "regulation of apoptotic signaling pathway," "cellular response to DNA damage stimulus," "adaptive immune system," and "pyroptosis." The miR-214-LTF-FOXJ1 axis may be a key regulatory pathway for IS. Conclusion: Our findings suggest that EIF3D, CYFIP2, NCBP2, DCPS, and NUDT1 may serve as potential diagnostic biomarkers for IS and that the m7G clusters developed by these genes provide more evidence for the regulation of m7G in IS.

Nuclear import receptors are recruited by FG-nucleoporins to rescue hallmarks of TDP-43 proteinopathy

BACKGROUND

Cytoplasmic mislocalization and aggregation of TAR DNA-binding protein-43 (TDP-43) is a hallmark of the amyotrophic lateral sclerosis and frontotemporal dementia (ALS/FTD) disease spectrum, causing both nuclear loss-of-function and cytoplasmic toxic gain-of-function phenotypes. While TDP-43 proteinopathy has been associated with defects in nucleocytoplasmic transport, this process is still poorly understood. Here we study the role of karyopherin-β1 (KPNB1) and other nuclear import receptors in regulating TDP-43 pathology.

METHODS

We used immunostaining, immunoprecipitation, biochemical and toxicity assays in cell lines, primary neuron and organotypic mouse brain slice cultures, to determine the impact of KPNB1 on the solubility, localization, and toxicity of pathological TDP-43 constructs. Postmortem patient brain and spinal cord tissue was stained to assess KPNB1 colocalization with TDP-43 inclusions. Turbidity assays were employed to study the dissolution and prevention of aggregation of recombinant TDP-43 fibrils in vitro. Fly models of TDP-43 proteinopathy were used to determine the effect of KPNB1 on their neurodegenerative phenotype in vivo.

RESULTS

We discovered that several members of the nuclear import receptor protein family can reduce the formation of pathological TDP-43 aggregates. Using KPNB1 as a model, we found that its activity depends on the prion-like C-terminal region of TDP-43, which mediates the co-aggregation with phenylalanine and glycine-rich nucleoporins (FG-Nups) such as Nup62. KPNB1 is recruited into these co-aggregates where it acts as a molecular chaperone that reverses aberrant phase transition of Nup62 and TDP-43. These findings are supported by the discovery that Nup62 and KPNB1 are also sequestered into pathological TDP-43 aggregates in ALS/FTD postmortem CNS tissue, and by the identification of the fly ortholog of KPNB1 as a strong protective modifier in Drosophila models of TDP-43 proteinopathy. Our results show that KPNB1 can rescue all hallmarks of TDP-43 pathology, by restoring its solubility and nuclear localization, and reducing neurodegeneration in cellular and animal models of ALS/FTD.

CONCLUSION

Our findings suggest a novel NLS-independent mechanism where, analogous to its canonical role in dissolving the diffusion barrier formed by FG-Nups in the nuclear pore, KPNB1 is recruited into TDP-43/FG-Nup co-aggregates present in TDP-43 proteinopathies and therapeutically reverses their deleterious phase transition and mislocalization, mitigating neurodegeneration.

Integrated Analysis of Tumor Mutation Burden and Immune Infiltrates in Hepatocellular Carcinoma

Tumor mutation burdens (TMBs) act as an indicator of immunotherapeutic responsiveness in various tumors. However, the relationship between TMBs and immune cell infiltrates in hepatocellular carcinoma (HCC) is still obscure. The present study aimed to explore the potential diagnostic markers of TMBs for HCC and analyze the role of immune cell infiltration in this pathology. We used OA datasets from The Cancer Genome Atlas database. First, the “maftools” package was used to screen the highest mutation frequency in all samples. R software was used to identify differentially expressed genes (DEGs) according to mutation frequency and perform functional correlation analysis. Then, the gene ontology (GO) enrichment analysis was performed with “clusterProfiler”, “enrichplot”, and “ggplot2” packages. Finally, the correlations between diagnostic markers and infiltrating immune cells were analyzed, and CIBERSORT was used to evaluate the infiltration of immune cells in HCC tissues. As a result, we identified a total of 359 DEGs in this study. These DEGs may affect HCC prognosis by regulating fatty acid metabolism, hypoxia, and the P53 pathway. The top 15 genes were selected as the hub genes through PPI network analysis. SRSF1, SNRPA1, and SRSF3 showed strong similarities in biological effects, NCBP2 was demonstrated as a diagnostic marker of HCC, and high NCBP2 expression was significantly correlated with poor over survival (OS) in HCC. In addition, NCBP2 expression was correlated with the infiltration of B cells (r = 0.364, p = 3.30 × 10⁻¹²), CD8⁺ T cells (r = 0.295, p = 2.71 × 10⁻⁸), CD4⁺ T cells, (r = 0.484, p = 1.37 × 10⁻²¹), macrophages (r = 0.551, p = 1.97 × 10⁻²⁸), neutrophils (r = 0.457, p = 3.26 × 10⁻¹⁹), and dendritic cells (r = 0.453, p = 1.97 × 10⁻¹⁸). Immune cell infiltration analysis revealed that the degree of central memory T-cell (Tcm) infiltration may be correlated with the HCC process. In conclusion, NCBP2 can be used as diagnostic markers of HCC, and immune cell infiltration plays an important role in the occurrence and progression of HCC.

A cross-disorder dosage sensitivity map of the human genome

Rare copy-number variants (rCNVs) include deletions and duplications that occur infrequently in the global human population and can confer substantial risk for disease. In this study, we aimed to quantify the properties of haploinsufficiency (i.e., deletion intolerance) and triplosensitivity (i.e., duplication intolerance) throughout the human genome. We harmonized and meta-analyzed rCNVs from nearly one million individuals to construct a genome-wide catalog of dosage sensitivity across 54 disorders, which defined 163 dosage sensitive segments associated with at least one disorder. These segments were typically gene dense and often harbored dominant dosage sensitive driver genes, which we were able to prioritize using statistical fine-mapping. Finally, we designed an ensemble machine-learning model to predict probabilities of dosage sensitivity (pHaplo & pTriplo) for all autosomal genes, which identified 2,987 haploinsufficient and 1,559 triplosensitive genes, including 648 that were uniquely triplosensitive. This dosage sensitivity resource will provide broad utility for human disease research and clinical genetics.

16p12.1 Deletion Orthologs are Expressed in Motile Neural Crest Cells and are Important for Regulating Craniofacial Development in Xenopus laevis

Copy number variants (CNVs) associated with neurodevelopmental disorders are characterized by extensive phenotypic heterogeneity. In particular, one CNV was identified in a subset of children clinically diagnosed with intellectual disabilities (ID) that results in a hemizygous deletion of multiple genes at chromosome 16p12.1. In addition to ID, individuals with this deletion display a variety of symptoms including microcephaly, seizures, cardiac defects, and growth retardation. Moreover, patients also manifest severe craniofacial abnormalities, such as micrognathia, cartilage malformation of the ears and nose, and facial asymmetries; however, the function of the genes within the 16p12.1 region have not been studied in the context of vertebrate craniofacial development. The craniofacial tissues affected in patients with this deletion all derive from the same embryonic precursor, the cranial neural crest, leading to the hypothesis that one or more of the 16p12.1 genes may be involved in regulating neural crest cell (NCC)-related processes. To examine this, we characterized the developmental role of the 16p12.1-affected gene orthologs, polr3e, mosmo, uqcrc2, and cdr2, during craniofacial morphogenesis in the vertebrate model system, Xenopus laevis. While the currently-known cellular functions of these genes are diverse, we find that they share similar expression patterns along the neural tube, pharyngeal arches, and later craniofacial structures. As these genes show co-expression in the pharyngeal arches where NCCs reside, we sought to elucidate the effect of individual gene depletion on craniofacial development and NCC migration. We find that reduction of several 16p12.1 genes significantly disrupts craniofacial and cartilage formation, pharyngeal arch migration, as well as NCC specification and motility. Thus, we have determined that some of these genes play an essential role during vertebrate craniofacial patterning by regulating specific processes during NCC development, which may be an underlying mechanism contributing to the craniofacial defects associated with the 16p12.1 deletion.

Ventx Family and Its Functional Similarities with Nanog: Involvement in Embryonic Development and Cancer Progression

The Ventx family is one of the subfamilies of the ANTP (antennapedia) superfamily and belongs to the NK-like (NKL) subclass. Ventx is a homeobox transcription factor and has a DNA-interacting domain that is evolutionarily conserved throughout vertebrates. It has been extensively studied in Xenopus, zebrafish, and humans. The Ventx family contains transcriptional repressors widely involved in embryonic development and tumorigenesis in vertebrates. Several studies have documented that the Ventx family inhibited dorsal mesodermal formation, neural induction, and head formation in Xenopus and zebrafish. Moreover, Ventx2.2 showed functional similarities to Nanog and Barx1, leading to pluripotency and neural-crest migration in vertebrates. Among them, Ventx protein is an orthologue of the Ventx family in humans. Studies have demonstrated that human Ventx was strongly associated with myeloid-cell differentiation and acute myeloid leukemia. The therapeutic potential of Ventx family inhibition in combating cancer progression in humans is discussed. Additionally, we briefly discuss genome evolution, gene duplication, pseudo-allotetraploidy, and the homeobox family in Xenopus.

Mitochondrial Proteostasis Requires Genes Encoded in a Neurodevelopmental Syndrome Locus

Eukaryotic cells maintain proteostasis through mechanisms that require cytoplasmic and mitochondrial translation. Genetic defects affecting cytoplasmic translation perturb synapse development, neurotransmission, and are causative of neurodevelopmental disorders, such as Fragile X syndrome. In contrast, there is little indication that mitochondrial proteostasis, either in the form of mitochondrial protein translation and/or degradation, is required for synapse development and function. Here we focus on two genes deleted in a recurrent copy number variation causing neurodevelopmental disorders, the 22q11.2 microdeletion syndrome. We demonstrate that SLC25A1 and MRPL40, two genes present in the microdeleted segment and whose products localize to mitochondria, interact and are necessary for mitochondrial ribosomal integrity and proteostasis. Our Drosophila studies show that mitochondrial ribosome function is necessary for synapse neurodevelopment, function, and behavior. We propose that mitochondrial proteostasis perturbations, either by genetic or environmental factors, are a pathogenic mechanism for neurodevelopmental disorders.SIGNIFICANCE STATEMENT The balance between cytoplasmic protein synthesis and degradation, or cytoplasmic proteostasis, is required for normal synapse function and neurodevelopment. Cytoplasmic and mitochondrial ribosomes are necessary for two compartmentalized, yet interdependent, forms of proteostasis. Proteostasis dependent on cytoplasmic ribosomes is a well-established target of genetic defects that cause neurodevelopmental disorders, such as autism. Here we show that the mitochondrial ribosome is a neurodevelopmentally regulated organelle whose function is required for synapse development and function. We propose that defective mitochondrial proteostasis is a mechanism with the potential to contribute to neurodevelopmental disease.

Convergent and distributed effects of the 3q29 deletion on the human neural transcriptome

The 3q29 deletion (3q29Del) confers high risk for schizophrenia and other neurodevelopmental and psychiatric disorders. However, no single gene in this interval is definitively associated with disease, prompting the hypothesis that neuropsychiatric sequelae emerge upon loss of multiple functionally-connected genes. 3q29 genes are unevenly annotated and the impact of 3q29Del on the human neural transcriptome is unknown. To systematically formulate unbiased hypotheses about molecular mechanisms linking 3q29Del to neuropsychiatric illness, we conducted a systems-level network analysis of the non-pathological adult human cortical transcriptome and generated evidence-based predictions that relate 3q29 genes to novel functions and disease associations. The 21 protein-coding genes located in the interval segregated into seven clusters of highly co-expressed genes, demonstrating both convergent and distributed effects of 3q29Del across the interrogated transcriptomic landscape. Pathway analysis of these clusters indicated involvement in nervous-system functions, including synaptic signaling and organization, as well as core cellular functions, including transcriptional regulation, posttranslational modifications, chromatin remodeling, and mitochondrial metabolism. Top network-neighbors of 3q29 genes showed significant overlap with known schizophrenia, autism, and intellectual disability-risk genes, suggesting that 3q29Del biology is relevant to idiopathic disease. Leveraging "guilt by association", we propose nine 3q29 genes, including one hub gene, as prioritized drivers of neuropsychiatric risk. These results provide testable hypotheses for experimental analysis on causal drivers and mechanisms of the largest known genetic risk factor for schizophrenia and highlight the study of normal function in non-pathological postmortem tissue to further our understanding of psychiatric genetics, especially for rare syndromes like 3q29Del, where access to neural tissue from carriers is unavailable or limited.

Craniofacial features of 3q29 deletion syndrome: Application of next-generation phenotyping technology

3q29 deletion syndrome (3q29del) is a recurrent deletion syndrome associated with neuropsychiatric disorders and congenital anomalies. Dysmorphic facial features have been described but not systematically characterized. This study aims to detail the 3q29del craniofacial phenotype and use a machine learning approach to categorize individuals with 3q29del through analysis of 2D photos. Detailed dysmorphology exam and 2D facial photos were ascertained from 31 individuals with 3q29del. Photos were used to train the next-generation phenotyping algorithm DeepGestalt (Face2Gene by FDNA, Inc, Boston, MA) to distinguish 3q29del cases from controls and all other recognized syndromes. Area under the curve of receiver operating characteristic curves (AUC-ROC) was used to determine the capacity of Face2Gene to identify 3q29del cases against controls. In this cohort, the most common observed craniofacial features were prominent forehead (48.4%), prominent nose tip (35.5%), and thin upper lip vermillion (25.8%). The FDNA technology showed an ability to distinguish cases from controls with an AUC-ROC value of 0.873 (p = 0.006) and led to the inclusion of 3q29del as one of the supported syndromes. This study found a recognizable facial pattern in 3q29del, as observed by trained clinical geneticists and next-generation phenotyping technology. These results expand the potential application of automated technology such as FDNA in identifying rare genetic syndromes, even when facial dysmorphology is subtle.

Functional assessment of the "two-hit" model for neurodevelopmental defects in Drosophila and X. laevis

We previously identified a deletion on chromosome 16p12.1 that is mostly inherited and associated with multiple neurodevelopmental outcomes, where severely affected probands carried an excess of rare pathogenic variants compared to mildly affected carrier parents. We hypothesized that the 16p12.1 deletion sensitizes the genome for disease, while "second-hits" in the genetic background modulate the phenotypic trajectory. To test this model, we examined how neurodevelopmental defects conferred by knockdown of individual 16p12.1 homologs are modulated by simultaneous knockdown of homologs of "second-hit" genes in Drosophila melanogaster and Xenopus laevis. We observed that knockdown of 16p12.1 homologs affect multiple phenotypic domains, leading to delayed developmental timing, seizure susceptibility, brain alterations, abnormal dendrite and axonal morphology, and cellular proliferation defects. Compared to genes within the 16p11.2 deletion, which has higher de novo occurrence, 16p12.1 homologs were less likely to interact with each other in Drosophila models or a human brain-specific interaction network, suggesting that interactions with "second-hit" genes may confer higher impact towards neurodevelopmental phenotypes. Assessment of 212 pairwise interactions in Drosophila between 16p12.1 homologs and 76 homologs of patient-specific "second-hit" genes (such as ARID1B and CACNA1A), genes within neurodevelopmental pathways (such as PTEN and UBE3A), and transcriptomic targets (such as DSCAM and TRRAP) identified genetic interactions in 63% of the tested pairs. In 11 out of 15 families, patient-specific "second-hits" enhanced or suppressed the phenotypic effects of one or many 16p12.1 homologs in 32/96 pairwise combinations tested. In fact, homologs of SETD5 synergistically interacted with homologs of MOSMO in both Drosophila and X. laevis, leading to modified cellular and brain phenotypes, as well as axon outgrowth defects that were not observed with knockdown of either individual homolog. Our results suggest that several 16p12.1 genes sensitize the genome towards neurodevelopmental defects, and complex interactions with "second-hit" genes determine the ultimate phenotypic manifestation.

Dissecting the complexity of CNV pathogenicity: insights from Drosophila and zebrafish models

Genetic architecture predisposes regions of the human genome to copy-number variants, which confer substantial disease risk, most prominently towards neurodevelopmental disorders. These variants typically contain multiple genes and are often associated with extensive pleiotropy and variable phenotypic expressivity. Despite the expansion of the fidelity of CNV detection, and the study of such lesions at the population level, understanding causal mechanisms for CNV phenotypes will require biological testing of constituent genes and their interactions. In this regard, model systems amenable to high-throughput phenotypic analysis of dosage-sensitive genes (and combinations thereof) are beginning to offer improved granularity of CNV-driven pathology. Here, we review the utility of Drosophila and zebrafish models for pathogenic CNV regions, highlight the advances made in discovery of single gene drivers and genetic interactions that determine specific CNV phenotypes, and argue for their validity in dissecting conserved developmental mechanisms associated with CNVs.

Family-effects in the epigenomic response of red blood cells to a challenge test in the European sea bass (Dicentrarchus labrax, L.)

Abstract

Background

In fish, minimally invasive blood sampling is widely used to monitor physiological stress with blood plasma biomarkers. As fish blood cells are nucleated, they might be a source a potential new markers derived from ‘omics technologies. We modified the epiGBS (epiGenotyping By Sequencing) technique to explore changes in genome-wide cytosine methylation in the red blood cells (RBCs) of challenged European sea bass (Dicentrarchus labrax), a species widely studied in both natural and farmed environments.

Results

We retrieved 501,108,033 sequencing reads after trimming, with a mean mapping efficiency of 73.0% (unique best hits). Minor changes in RBC methylome appeared to manifest after the challenge test and a family-effect was detected. Only fifty-seven differentially methylated cytosines (DMCs) close to 51 distinct genes distributed on 17 of 24 linkage groups (LGs) were detected between RBCs of pre- and post-challenge individuals. Thirty-seven of these genes were previously reported as differentially expressed in the brain of zebrafish, most of them involved in stress coping differences. While further investigation remains necessary, few DMC-related genes associated to the Brain Derived Neurotrophic Factor, a protein that favors stress adaptation and fear memory, appear relevant to integrate a centrally produced stress response in RBCs.

Conclusion

Our modified epiGBS protocol was powerful to analyze patterns of cytosine methylation in RBCs of D. labrax and to evaluate the impact of a challenge using minimally invasive blood samples. This study is the first approximation to identify epigenetic biomarkers of exposure to stress in fish.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07420-9.

Rab11fip5 regulates telencephalon development via ephrinB1 recycling

Rab11 family-interacting protein 5 (Rab11fip5) is an adaptor protein that binds to the small GTPase Rab11, which has an important function in endosome recycling and trafficking of cellular proteins to the plasma membrane. Rab11fip5 is involved in many cellular processes, such as cytoskeleton rearrangement, iron uptake and exocytosis in neuroendocrine cells, and is also known as a candidate gene for autism-spectrum disorder. However, the role of Rab11fip5 during early embryonic development is not clearly understood. In this study, we identified Rab11fip5 as a protein that interacts with ephrinB1, a transmembrane ligand for Eph receptors. The PDZ binding motif in ephrinB1 and the Rab-binding domain in Rab11fip5 are necessary for their interaction in a complex. EphrinB1 and Rab11fip5 display overlapping expression in the telencephalon of developing amphibian embryos. The loss of Rab11fip5 function causes a reduction in telencephalon size and a decrease in the expression level of ephrinB1. Moreover, morpholino oligonucleotide-mediated knockdown of Rab11fip5 decreases cell proliferation in the telencephalon. The overexpression of ephrinB1 rescues these defects, suggesting that ephrinB1 recycling by the Rab11/Rab11fip5 complex is crucial for proper telencephalon development.

Summary: Rab11fip5, genetic mutations of which are associated with autism spectrum disorders in humans, plays a role in regulating telencephalon development via recycling of the crucial cargo ephrinB1.

Drosophila models of pathogenic copy-number variant genes show global and non-neuronal defects during development

While rare pathogenic copy-number variants (CNVs) are associated with both neuronal and non-neuronal phenotypes, functional studies evaluating these regions have focused on the molecular basis of neuronal defects. We report a systematic functional analysis of non-neuronal defects for homologs of 59 genes within ten pathogenic CNVs and 20 neurodevelopmental genes in Drosophila melanogaster. Using wing-specific knockdown of 136 RNA interference lines, we identified qualitative and quantitative phenotypes in 72/79 homologs, including 21 lines with severe wing defects and six lines with lethality. In fact, we found that 10/31 homologs of CNV genes also showed complete or partial lethality at larval or pupal stages with ubiquitous knockdown. Comparisons between eye and wing-specific knockdown of 37/45 homologs showed both neuronal and non-neuronal defects, but with no correlation in the severity of defects. We further observed disruptions in cell proliferation and apoptosis in larval wing discs for 23/27 homologs, and altered Wnt, Hedgehog and Notch signaling for 9/14 homologs, including AATF/Aatf, PPP4C/Pp4-19C, and KIF11/Klp61F. These findings were further supported by tissue-specific differences in expression patterns of human CNV genes, as well as connectivity of CNV genes to signaling pathway genes in brain, heart and kidney-specific networks. Our findings suggest that multiple genes within each CNV differentially affect both global and tissue-specific developmental processes within conserved pathways, and that their roles are not restricted to neuronal functions.

Rare copy-number variants (CNVs), or large deletions and duplications in the genome, are associated with both neuronal and non-neuronal clinical features. Previous functional studies for these disorders have primarily focused on understanding the cellular mechanisms for neurological and behavioral phenotypes. To understand how genes within these CNVs contribute to developmental defects in non-neuronal tissues, we assessed 79 homologs of CNV and known neurodevelopmental genes in Drosophila models. We found that most homologs showed developmental defects when knocked down in the adult fly wing, ranging from mild size changes to severe wrinkled wings or lethality. Although a majority of tested homologs showed defects when knocked down specifically in wings or eyes, we found no correlation in the severity of the observed defects in these two tissues. A subset of the homologs showed disruptions in cellular processes in the developing fly wing, including alterations in cell proliferation, apoptosis, and cellular signaling pathways. Furthermore, human CNV genes also showed differences in gene expression patterns and interactions with signaling pathway genes across multiple human tissues. Our findings suggest that genes within CNV disorders affect global developmental processes in both neuronal and non-neuronal tissues.