Resolving genomic disorder-associated breakpoints within segmental DNA duplications using massively parallel sequencing

Parallelized engineering of mutational models using piggyBac transposon delivery of CRISPR libraries

New technologies and large-cohort studies have enabled novel variant discovery and association at unprecedented scale, yet functional characterization of these variants remains paramount to deciphering disease mechanisms. Approaches that facilitate parallelized genome editing of cells of interest or induced pluripotent stem cells (iPSCs) have become critical tools toward this goal. Here, we developed an approach that incorporates libraries of CRISPR-Cas9 guide RNAs (gRNAs) together with inducible Cas9 into a piggyBac (PB) transposon system to engineer dozens to hundreds of genomic variants in parallel against isogenic cellular backgrounds. This method empowers loss-of-function (LoF) studies through the introduction of insertions or deletions (indels) and copy-number variants (CNVs), though generating specific nucleotide changes is possible with prime editing. The ability to rapidly establish high-quality mutational models at scale will facilitate the development of isogenic cellular collections and catalyze comparative functional genomic studies investigating the roles of hundreds of genes and mutations in development and disease.

Possible association of 16p11.2 copy number variation with altered lymphocyte and neutrophil counts

Recurrent copy-number variations (CNVs) at chromosome 16p11.2 are associated with neurodevelopmental diseases, skeletal system abnormalities, anemia, and genitourinary defects. Among the 40 protein-coding genes encompassed within the rearrangement, some have roles in leukocyte biology and immunodeficiency, like SPN and CORO1A. We therefore investigated leukocyte differential counts and disease in 16p11.2 CNV carriers. In our clinically-recruited cohort, we identified three deletion carriers from two families (out of 32 families assessed) with neutropenia and lymphopenia. They had no deleterious single-nucleotide or indel variant in known cytopenia genes, suggesting a possible causative role of the deletion. Noticeably, all three individuals had the lowest copy number of the human-specific BOLA2 duplicon (copy-number range: 3–8). Consistent with the lymphopenia and in contrast with the neutropenia associations, adult deletion carriers from UK biobank (n = 74) showed lower lymphocyte (Padj = 0.04) and increased neutrophil (Padj = 8.31e-05) counts. Mendelian randomization studies pinpointed to reduced CORO1A, KIF22, and BOLA2-SMG1P6 expressions being causative for the lower lymphocyte counts. In conclusion, our data suggest that 16p11.2 deletion, and possibly also the lowest dosage of the BOLA2 duplicon, are associated with low lymphocyte counts. There is a trend between 16p11.2 deletion with lower copy-number of the BOLA2 duplicon and higher susceptibility to moderate neutropenia. Higher numbers of cases are warranted to confirm the association with neutropenia and to resolve the involvement of the deletion coupled with deleterious variants in other genes and/or with the structure and copy number of segments in the CNV breakpoint regions.

16p13.11p11.2 triplication syndrome: a new recognizable genomic disorder characterized by optical genome mapping and whole genome sequencing

Highly identical segmental duplications (SDs) account for over 5% of the human genome and are enriched in the short arm of the chromosome 16. These SDs are susceptibility factors for recurrent chromosomal rearrangements mediated by non-allelic homologous recombination (NAHR). Chromosomal microarray analysis (CMA) has been widely used as the first-tier test for individuals with developmental disabilities and/or congenital anomalies and several genomic disorders involving the 16p-arm have been identified with this technique. However, the resolution of CMA and the limitations of short-reads whole genome sequencing (WGS) technology do not allow the full characterization of the most complex chromosomal rearrangements. Herein, we report on two unrelated patients with a de novo 16p13.11p11.2 triplication associated with a 16p11.2 duplication, detected by CMA. These patients share a similar phenotype including hypotonia, severe neurodevelopmental delay with profound speech impairment, hyperkinetic behavior, conductive hearing loss, and distinctive facial features. Short-reads WGS could not map precisely any of the rearrangement's breakpoints that lie within SDs. We used optical genome mapping (OGM) to determine the relative orientation of the triplicated and duplicated segments as well as the genomic positions of the breakpoints, allowing us to propose a mechanism involving recombination between allelic SDs and a NAHR event. In conclusion, we report a new clinically recognizable genomic disorder. In addition, the mechanism of these complex chromosomal rearrangements involving SDs could be unraveled by OGM.

A Combined Proteomics and Metabolomics Profiling to Investigate the Genetic Heterogeneity of Autistic Children

Autism spectrum disorder (ASD) has become one of the most common neurological developmental disorders in children. However, the study of ASD diagnostic markers faces significant challenges due to the existence of heterogeneity. In this study, genetic testing was performed on children who were clinically diagnosed with ASD. Children with ASD susceptibility genes and healthy controls were studied. The proteomics of plasma and peripheral blood mononuclear cells (PBMCs) as well as plasma metabolomics were carried out. The results showed that although there was genetic heterogeneity in children with ASD, the differentially expressed proteins (DEPs) in plasma, peripheral blood mononuclear cells, and differential metabolites in plasma could still effectively distinguish autistic children from controls. The mechanism associated with them focuses on several common and previously reported mechanisms of ASD. The biomarkers for ASD diagnosis could be found by taking differentially expressed proteins and differential metabolites into consideration. Integrating omics data, glycerophospholipid metabolism and N-glycan biosynthesis might play a critical role in the pathogenesis of ASD.

Genetic analysis of N6-methyladenosine modification genes in Parkinson's disease

Parkinson's disease (PD) is a common neurodegenerative disease with a relatively unclear etiology. Previous studies have shown that N6-methyladenosine (m6A) is a vital RNA modification enriched in brain tissue, and that the genes involved in m6A modification are implicated in various neurologic diseases. Here, we conducted a comprehensive genetic analysis using targeted sequencing with molecular inversion probes (MIPs) to identify m6A-modification genes (including METTL3, METTL14, WTAP, FTO, ALKBH5, YTHDF1, YTHDF2, YTHDF3, HNRNPC, and ELAVL1) in a total of 1647 sporadic PD patients and 1372 controls of Han Chinese origin. PD patients were divided into early-onset PD (EOPD) and late-onset PD (LOPD) based on whether the onset of motor symptoms occurred before or after 50 years of age. Rare variants were subjected to gene-based burden tests and common variants were subjected to single-variant association analyses. As a result, we identified 214 rare variants in all 10 m6A-modification genes and 16 common variants in 7 genes. Gene-wise association analyses of rare variants in each m6A-modification gene did not achieved a p value of less than 0.05 in either total cohorts or 2 age groups. In fact, p values greater than 0.05 were found when conducting single-variant association analyses on common variants of these genes between PD and control patients. Our comprehensive analyses of m6A-modification genes suggest that there is no significant association between these 10 m6A-modification genes and the risk of sporadic PD.

Genetic analysis of DNA methylation and hydroxymethylation genes in Parkinson's disease

DNA methylation is an important regulatory mechanism of Parkinson's disease (PD). To investigate the relationship between DNA methylation and hydroxymethylation genes and PD, we performed gene-targeted sequencing using molecular inversion probes in a Chinese PD population. We sequenced 12 genes related to DNA methylation and hydroxymethylation in 1657 patients and 1394 control subjects. We conducted genewise association analyses of rare variants detected in the present study and identified the TET1 gene as important in PD (p = 0.0037738, 0.013, 0.019521 (b.collapse test, variable threshold test, and skat-o test, respectively; sex + age as covariates). However, no positive results were observed when conducting association analyses on common variants in these genes. We performed a comprehensive analysis of associations between variants of DNA methylation and hydroxymethylation genes and PD, resulting in determination that TET1 might play a role in PD.

Systematically analyzing rare variants of autosomal-dominant genes for sporadic Parkinson's disease in a Chinese cohort

Studies have shown that rare variants of Mendelian genes for Parkinson's disease (PD) contribute to sporadic PD in the Caucasian population, which lacked confirmation in the Chinese population. Because the autosomal-dominant PD (AD-PD) had a phenotype closely resembling sporadic PD, we performed a systematic analysis of 7 AD-PD genes (SNCA, LRRK2, GIGYF2, VPS35, EIF4G1, DNAJC13, and CHCHD2) in 1456 Chinese sporadic PD patients and 1568 controls. Overall, 72 rare variants were identified, 7 of which were classified as likely pathogenic, 63 of which were categorized as of uncertain significance, and 2 of them were predicted to be likely benign. These AD-PD genes represented a clear enrichment of rare variants in PD patients from a burden analysis (p = 0.003), and significant differences could still be observed when likely pathogenic variants were removed (p = 0.027). The gene-based association testing also reached significance for LRRK2 (p = 0.004) and remained statistically significant after the Bonferroni correction. This report suggested that rare variants of AD-PD genes had a role in the Chinese sporadic PD cohort, especially for those rare variants of LRRK2.

An exploratory study of predisposing genetic factors for DiGeorge/velocardiofacial syndrome

DiGeorge/velocardiofacial syndrome (DGS/VCFS) is a disorder caused by a 22q11.2 deletion mediated by non-allelic homologous recombination (NAHR) between low-copy repeats (LCRs). We have evaluated the role of LCR22 genomic architecture and PRDM9 variants as DGS/VCFS predisposing factors. We applied FISH using fosmid probes on chromatin fibers to analyze the number of tandem repeat blocks in LCR22 in two DGS/VCFS fathers-of-origin with proven 22q11.2 NAHR susceptibility. Results revealed copy number variations (CNVs) of L9 and K3 fosmids in these individuals compared to controls. The total number of L9 and K3 copies was also characterized using droplet digital PCR (ddPCR). Although we were unable to confirm variations, we detected an additional L9 amplicon corresponding to a pseudogene. Moreover, none of the eight DGS/VCFS parents-of-origin was heterozygote for the inv(22)(q11.2) haplotype. PRDM9 sequencing showed equivalent allelic distributions between DGS/VCFS parents-of-origin and controls, although a new PRDM9 allele (L50) was identified in one case. Our results support the hypothesis that LCR22s variations influences 22q11.2 NAHR events, however further studies are needed to confirm this association and clarify the contribution of pseudogenes and rare PDRM9 alleles to NAHR susceptibility.

Emergence of a Homo sapiens-specific gene family and chromosome 16p11.2 CNV susceptibility

Genetic differences that specify unique aspects of human evolution have typically been identified by comparative analyses between the genomes of humans and closely related primates, including more recently the genomes of archaic hominins. Not all regions of the genome, however, are equally amenable to such study. Recurrent copy number variation (CNV) at chromosome 16p11.2 accounts for approximately 1% of cases of autism and is mediated by a complex set of segmental duplications, many of which arose recently during human evolution. Here we reconstruct the evolutionary history of the locus and identify bolA family member 2 (BOLA2) as a gene duplicated exclusively in Homo sapiens. We estimate that a 95-kilobase-pair segment containing BOLA2 duplicated across the critical region approximately 282 thousand years ago (ka), one of the latest among a series of genomic changes that dramatically restructured the locus during hominid evolution. All humans examined carried one or more copies of the duplication, which nearly fixed early in the human lineage--a pattern unlikely to have arisen so rapidly in the absence of selection (P < 0.0097). We show that the duplication of BOLA2 led to a novel, human-specific in-frame fusion transcript and that BOLA2 copy number correlates with both RNA expression (r = 0.36) and protein level (r = 0.65), with the greatest expression difference between human and chimpanzee in experimentally derived stem cells. Analyses of 152 patients carrying a chromosome 16p11. rearrangement show that more than 96% of breakpoints occur within the H. sapiens-specific duplication. In summary, the duplicative transposition of BOLA2 at the root of the H. sapiens lineage about 282 ka simultaneously increased copy number of a gene associated with iron homeostasis and predisposed our species to recurrent rearrangements associated with disease.

Variant discovery and breakpoint region prediction for studying the human 22q11.2 deletion using BAC clone and whole genome sequencing analysis

Velo-cardio-facial syndrome/DiGeorge syndrome/22q11.2 deletion syndrome (22q11.2DS) is caused by meiotic non-allelic homologous recombination events between flanking low copy repeats termed LCR22A and LCR22D, resulting in a 3 million base pair (Mb) deletion. Due to their complex structure, large size and high sequence identity, genetic variation within LCR22s among different individuals has not been well characterized. In this study, we sequenced 13 BAC clones derived from LCR22A/D and aligned them with 15 previously available BAC sequences to create a new genetic variation map. The thousands of variants identified by this analysis were not uniformly distributed in the two LCR22s. Moreover, shared single nucleotide variants between LCR22A and LCR22D were enriched in the Breakpoint Cluster Region pseudogene (BCRP) block, suggesting the existence of a possible recombination hotspot there. Interestingly, breakpoints for atypical 22q11.2 rearrangements have previously been located to BCRPs To further explore this finding, we carried out in-depth analyses of whole genome sequence (WGS) data from two unrelated probands harbouring a de novo 3Mb 22q11.2 deletion and their normal parents. By focusing primarily on WGS reads uniquely mapped to LCR22A, using the variation map from our BAC analysis to help resolve allele ambiguity, and by performing PCR analysis, we infer that the deletion breakpoints were most likely located near or within the BCRP module. In summary, we found a high degree of sequence variation in LCR22A and LCR22D and a potential recombination breakpoint near or within the BCRP block, providing a starting point for future breakpoint mapping using additional trios.

Spherical nucleic acids as a divergent platform for synthesizing RNA-nanoparticle conjugates through enzymatic ligation

Herein, we describe a rapid, divergent method for using spherical nucleic acids (SNAs) as a universal platform for attaching RNA to DNA-modified nanoparticles using enzyme-mediated techniques. This approach provides a sequence-specific method for the covalent attachment of one or more in vitro transcribed RNAs to a universal SNA scaffold, regardless of RNA sequence. The RNA-nanoparticle constructs are shown to effectively knock down two different gene targets using a single, dual-ligated nanoparticle construct.