Autism spectrum disorder (ASD) is heterogeneous in symptom and etiology. Rare copy number variations (CNVs) are important genetic factors contributing to ASD. Currently chromosomal microarray (CMA) detecting CNVs is recommended as a first-tier diagnostic assay, largely based on research in North America and Europe. The feature of rare CNVs has not been well characterized in ASD cohorts from non-European ancestry. In this study, high resolution CMA was utilized to investigate rare CNVs in a Chinese cohort of ASD (n = 401, including 177 mildly/moderately and 224 severely affected individuals), together with an ancestry-matched control cohort (n = 197). Diagnostic yield was about 4.2%, with 17 clinically significant CNVs identified in ASD individuals, of which 12 CNVs overlapped with recurrent autism risk loci or genes. Autosomal rare CNV burden analysis showed an overrepresentation of rare loss events in ASD cohort, whereas the rate of rare gain events correlated with the phenotypic severity. Further analysis showed rare losses disrupting genes highly intolerant of loss-of-function variants were enriched in the ASD cohort. Among these highly constrained genes disrupted by rare losses, RIMS2 is a promising candidate contributing to ASD risk. This pilot study evaluated clinical utility of CMA and the feature of rare CNVs in Chinese ASD, with candidate genes identified as potential risk factors.

Cleidocranial dysplasia (CCD) is a rare autosomal dominant disorder mainly characterized by hypoplastic or absent clavicles, delayed closure of the fontanelles, multiple dental abnormalities, and short stature. Runt-related transcription factor 2 (RUNX2) gene variants can cause CCD, but are not identified in all CCD patients.

PYGL mutations can cause liver phosphorylase deficiency, resulting in a glycogenolysis disorder, namely, glycogen storage disease (GSD) VI. The disease is rarely reported in the Chinese population. GSD VI is mainly characterized in untreated children by hepatomegaly, growth retardation and elevated liver transaminases.

We generated a knockout mouse for the neuronal-specific β-tubulin isoform Tubb3 to investigate its role in nervous system formation and maintenance. Tubb3-/- mice have no detectable neurobehavioral or neuropathological deficits, and upregulation of mRNA and protein of the remaining β-tubulin isotypes results in equivalent total β-tubulin levels in Tubb3-/- and wild-type mice. Despite similar levels of total β-tubulin, adult dorsal root ganglia lacking TUBB3 have decreased growth cone microtubule dynamics and a decreased neurite outgrowth rate of 22% in vitro and in vivo. The effect of the 22% slower growth rate is exacerbated for sensory recovery, where fibers must reinnervate the full volume of the skin to recover touch function. Overall, these data reveal that, while TUBB3 is not required for formation of the nervous system, it has a specific role in the rate of peripheral axon regeneration that cannot be replaced by other β-tubulins.

Neurodevelopment is a transcriptionally orchestrated process. Cyclin K, a regulator of transcription encoded by CCNK, is thought to play a critical role in the RNA polymerase II-mediated activities. However, dysfunction of CCNK has not been linked to genetic disorders. In this study, we identified three unrelated individuals harboring de novo heterozygous copy number loss of CCNK in an overlapping 14q32.3 region and one individual harboring a de novo nonsynonymous variant c.331A>G (p.Lys111Glu) in CCNK. These four individuals, though from different ethnic backgrounds, shared a common phenotype of developmental delay and intellectual disability (DD/ID), language defects, and distinctive facial dysmorphism including high hairline, hypertelorism, thin eyebrows, dysmorphic ears, broad nasal bridge and tip, and narrow jaw. Functional assay in zebrafish larvae showed that Ccnk knockdown resulted in defective brain development, small eyes, and curly spinal cord. These defects were partially rescued by wild-type mRNA coding CCNK but not the mRNA with the identified likely pathogenic variant c.331A>G, supporting a causal role of CCNK variants in neurodevelopmental disorders. Taken together, we reported a syndromic neurodevelopmental disorder with DD/ID and facial characteristics caused by CCNK variations, possibly through a mechanism of haploinsufficiency.

Though numerous studies have focused on the cell wall disassembly of bananas during the ripening process, the modification of homogalacturonan (HG) during fruit development remains exclusive. To better understand the role of HGs in controlling banana fruit growth and ripening, RNA-Seq, qPCR, immunofluorescence labeling, and biochemical methods were employed to reveal their dynamic changes in banana peels during these processes. Most HG-modifying genes in banana peels showed a decline in expression during fruit development. Four polygalacturonase and three pectin acetylesterases showing higher expression levels at later developmental stages than earlier ones might be related to fruit expansion. Six out of the 10 top genes in the Core Enrichment Gene Set were HG degradation genes, and all were upregulated after softening, paralleled to the significant increase in HG degradation enzyme activities, decline in peel firmness, and the epitope levels of 2F4, CCRC-M38, JIM7, and LM18 antibodies. Most differentially expressed alpha-1,4-galacturonosyltransferases were upregulated by ethylene treatment, suggesting active HG biosynthesis during the fruit softening process. The epitope level of the CCRC-M38 antibody was positively correlated to the firmness of banana peel during fruit development and ripening. These results have provided new insights into the role of cell wall HGs in fruit development and ripening.

Dendritic spines are small actin-rich protrusions essential for the formation of functional circuits in the mammalian brain. During development, spines begin as dynamic filopodia-like protrusions that are then replaced by relatively stable spines containing an expanded head. Remodeling of the actin cytoskeleton plays a key role in the formation and modification of spine morphology, however many of the underlying regulatory mechanisms remain unclear. Capping protein (CP) is a major actin regulating protein that caps the barbed ends of actin filaments, and promotes the formation of dense branched actin networks. Knockdown of CP impairs the formation of mature spines, leading to an increase in the number of filopodia-like protrusions and defects in synaptic transmission. Here, we show that CP promotes the stabilization of dendritic protrusions, leading to the formation of stable mature spines. However, the localization and function of CP in dendritic spines requires interactions with proteins containing a capping protein interaction (CPI) motif. We found that the CPI motif-containing protein Twinfilin-1 (Twf1) also localizes to spines where it plays a role in CP spine enrichment. The knockdown of Twf1 leads to an increase in the density of filopodia-like protrusions and a decrease in the stability of dendritic protrusions, similar to CP knockdown. Finally, we show that CP directly interacts with Shank and regulates its spine accumulation. These results suggest that spatiotemporal regulation of CP in spines not only controls the actin dynamics underlying the formation of stable postsynaptic spine structures, but also plays an important role in the assembly of the postsynaptic apparatus underlying synaptic function.

Noonan syndrome (NS) and Noonan syndrome with multiple lentigines (NSML) are autosomal dominant developmental disorders. NS and NSML are caused by abnormalities in genes that encode proteins related to the RAS-MAPK pathway, including PTPN11, RAF1, BRAF, and MAP2K. In this study, we diagnosed ten NS or NSML patients via targeted sequencing or whole exome sequencing (TS/WES).

Human height is a complex trait with a strong genetic basis. Recently, a significant association between rare copy number variations (CNVs) and short stature has been identified, and candidate genes in these rare CNVs are being explored. This study aims to evaluate the association between mutations in ARID1B gene and short stature, both the syndromic and non-syndromic form.

Disorders of sexual development (DSD) are estimated to occur in 1 of 4500 births. Since the genetic etiology of DSD is highly heterogeneous, obtaining a definitive molecular diagnosis by single gene test is challenging. Utilizing a high-throughput sequencing upfront is proposed as an efficient approach to aid in the diagnosis. This study aimed to examine the diagnostic yield of next-generation sequencing in DSD. 32 DSD patients that previously received clinical examinations and single gene tests were selected, with or without a diagnosis. Prior single gene tests were masked, and then samples went through targeted next-generation sequencing of 80 genes from which the diagnostic yield was assessed. A likely diagnosis, with pathogenic or likely pathogenic variants identified, was obtained from nine of the 32 patients (i.e., 28.1%, versus 10% by single gene tests). In another five patients (15.6%), variants of uncertain significance were found. Among 18 variants identified (i.e., 17 single nucleotide variants and one small deletion), eight had not been previously reported. This study supports the notion that next-generation sequencing can be an efficient tool in the clinical diagnosis and variant discovery in DSD.

Idiopathic short stature (ISS) refers to short stature with no evident etiologies. The custom genome-wide microarray specifically designed to cover height-related genes may be helpful to detect copy number variations (CNVs) in ISS patients, which may be missed by the general microarray. The aim of the study was to validate the applicability of the custom microarray and to analyze CNVs in Chinese ISS children.

Simple copy number variations (CNVs) detected by chromosomal microarray (CMA) can result from complex structural changes. Therefore, it is necessary to characterize potential structural changes that cause pathogenic CNVs. We applied whole-genome low-coverage sequencing (WGLCS) to concurrently detect pathogenic CNVs and their associated chromosomal rearrangements in 15 patients. All the patients had an average of 2-3 pathogenic CNVs involving 1-2 chromosomes. WGLCS identified all the 34 pathogenic CNVs found by microarray. By identifying chimeric read pairs, WGLCS mapped 70 breakpoints in these patients, of which 47 were finely mapped at the nucleotide level and confirmed by subsequent PCR amplification and Sanger sequencing of the junction fragments. In 15 patients, structural rearrangements were defined at molecular level in 13 patients. In 13 patients, WGLCS reveal no additional results in two patients. In another 11 patients, WGLCS revealed new breakpoints or finely mapped the genes disrupted by breakpoints or 1-6 bp microhomology and/or short insertion (4-70 bp) in the breakpoints junctions. However, structural changes in the other two patients still remained unclear after WGLCS was performed. The structural alteration identified in the 13 patients could be divided into the following categories: (1) interstitial inverted duplication with concomitant terminal deletion (inv dup del) (P1,P4,P9,P11); (2) the product of pericentric inversion (P5); (3) ring chromosome (P8); (4) interstitial duplication and/or triplication (P6, P7); and (5) +der(22)t(11;22) (P2,P15); (6) complex structural rearrangements (P3,P12,P14). WGLCS displayed the ability to discover CNVs and define breakpoints and its disrupted genes and its surrounding sequences in one experiment at base-pair-resolution, which help us to learn more about the mechanisms of formation of observed genomic rearrangements, and in which DNA replicative/repair mechanism might contribute to the formation of complex rearrangements in 11 patients. Clear karyotype at molecular level could help provide an accurate evaluation of recurrent risk and guide prenatal diagnosis or reproductive planning.

Autism spectrum disorder (ASD) is a group of clinically and genetically heterogeneous neurodevelopmental disorders. Recent tremendous advances in the whole exome sequencing (WES) enable rapid identification of variants associated with ASD including single nucleotide variations (SNVs) and indels. To further explore genetic etiology of ASD in Chinese children with negative findings of copy number variants (CNVs), we applied WES in 80 simplex families with a single affected offspring with ASD or suspected ASD, and validated variations predicted to be damaging by Sanger sequencing. The results showed that an overall diagnostic yield of 8.8% (9.2% in the group of ASD and 6.7% in the group of suspected ASD) was observed in our cohort. Among patients with diagnosed ASD, developmental delay or intellectual disability (DD/ID) was the most common comorbidity with a diagnostic yield of 13.3%, followed by seizures (50.0%) and craniofacial anomalies (40.0%). All of identified de novo SNVs and indels among patients with ASD were loss of function (LOF) variations and were slightly more frequent among female (male vs. female: 7.3% vs. 8.5%). A total of seven presumed causative genes (CHD8, AFF2, ADNP, POGZ, SHANK3, IL1RAPL1, and PTEN) were identified in this study. In conclusion, WES is an efficient diagnostic tool for diagnosed ASD especially those with negative findings of CNVs and other neurological disorders in clinical practice, enabling early identification of disease related genes and contributing to precision and personalized medicine.

Developmental delay (DD) and intellectual disability (ID) are frequently associated with a broad spectrum of additional phenotypes. Chromosomal microarray analysis (CMA) has been recommended as a first-tier test for DD/ID in general, whereas the diagnostic yield differs significantly among DD/ID patients with different comorbid conditions.

Dendritic spines undergo actin-based growth and shrinkage during synaptic plasticity, in which the actin depolymerizing factor (ADF)/cofilin family of actin-associated proteins are important. Elevated ADF/cofilin activities often lead to reduced spine size and immature spine morphology but can also enhance synaptic potentiation in some cases. Thus, ADF/cofilin may have distinct effects on postsynaptic structure and function. We found that ADF/cofilin-mediated actin dynamics regulated AMPA receptor (AMPAR) trafficking during synaptic potentiation, which was distinct from actin's structural role in spine morphology. Specifically, elevated ADF/cofilin activity markedly enhanced surface addition of AMPARs after chemically induced long-term potentiation (LTP), whereas inhibition of ADF/cofilin abolished AMPAR addition. We found that chemically induced LTP elicited a temporal sequence of ADF/cofilin dephosphorylation and phosphorylation that underlies AMPAR trafficking and spine enlargement. These findings suggest that temporally regulated ADF/cofilin activities function in postsynaptic modifications of receptor number and spine size during synaptic plasticity.

Phenylketonuria (PKU), caused by phenylalanine hydroxylase (PAH) gene variants, is a common autosomal inherited metabolic disease. So far, 1111 PAH variants have been revealed. The residual activity of the PAH variants is the key determinant of the metabolic phenotype and BH4 responsiveness in PKU patients. In this study, the spectrum of PAH variants in 1083 Chinese PKU patients was analyzed. Then 20 variants (p.L52F, p.R86P, p.L128P, p.L142P, p.D163N, p.C203G, p.E214G, p.F260L, p.M276T, p.L311R, p.P314A, p.L364F, p.Q375H, p.F382I, p.A395S, p.V412D, p.E108*, p.C203*, p.C284* and p.E353*) were expressed in COS-7 cells. The residual activities and protein expression levels were detected by isotope-dilution liquid chromatography-electrospray ionization tandem mass spectrometry (LC-ESI-MS/MS) and Western blotting, respectively. We compared the results of the phenotypic prediction based on APV and PAH activity respectively, and further explored the relationship between residual activity and phenotype in PKU patients. We reported 9 newly discovered PAH variants for the first time, thereby expanding the spectrum of PAH variants. Among the 20 variants in our assay, 8 variants showed mild impaired residual activities (48-92%) and approximately normal protein expression levels compared to the wild-type PAH. In contrast, 9 variants showed severely impaired residual activities (0-34%) and reduced protein expression. However, three variants (p.L52F, p.F260L and p.P314A) showed impaired residual activities (5%, 32% and 29%), although the proteins were well expressed. We assigned APV scores for 14 variants, in which the results of the phenotypic prediction were consistent for 12/14 (86%) variants based on APV and residual activity respectively, and the residual activity correctly predicted 17/22 (77%) of the patients. Our study helped to further understand the genotype-phenotype correlation in PKU patients.

Next generation sequencing-based exome sequencing can be used to identify genetic abnormalities in patients believed to be suffering from primary hyperoxaluria. We outline our efforts to improve the diagnostic capacity of exome sequencing for these patients.