Search

Therapy with mesenchymal stromal cells (MSCs) has shown promise in inflammatory bowel disease—leveraging their immunosuppressive and regenerative properties. However, the potential immunogenic complications of allogenic MSCs sourced from different tissues raise concern. Thus, we assessed the fitness and functionality of autologous intestinal MSCs as a potential platform for cellular therapy. Mucosal biopsy-derived MSCs from Crohn’s disease (n = 11), ulcerative colitis (n = 12), and controls (n = 14) were analyzed by microscopy and flow cytometry for doubling-time, morphology, differentiation potential, and immunophenotype. Gene expression, cell-subtype composition, along with surface marker and secretome changes after IFN-γ priming were measured by bulk and single-cell RNA sequencing coupled with a 30-plex Luminex panel. MSCs expanded ex vivo demonstrate canonical MSC markers, similar growth kinetics, and tripotency regardless of the patient phenotype. Global transcription patterns were similar at baseline though inflammatory bowel disease (IBD) rectal MSCs showed changes in select immunomodulatory genes. IFN-γ priming resulted in upregulation of shared immunoregulatory genes (particularly in PD-1 signaling) and overrode the transcriptional differences observed at baseline. Furthermore, MSCs secrete key immunomodulatory molecules at baseline and in response to IFN-γ including CXCL10, CXCL9, and MCP-1. Overall, MSCs from IBD patients have normal transcriptional and immunomodulatory properties with therapeutic potential and can be sufficiently expanded.

Identifying causative gene(s) within disease-associated large genomic regions of copy-number variants (CNVs) is challenging. Here, by targeted sequencing of genes within schizophrenia (SZ)-associated CNVs in 1,779 SZ cases and 1,418 controls, we identified three rare putative loss-of-function (LoF) mutations in OTU deubiquitinase 7A (OTUD7A) within the 15q13.3 deletion in cases but none in controls. To tie OTUD7A LoF with any SZ-relevant cellular phenotypes, we modeled the OTUD7A LoF mutation, rs757148409, in human induced pluripotent stem cell (hiPSC)-derived induced excitatory neurons (iNs) by CRISPR-Cas9 engineering. The mutant iNs showed a ∼50% decrease in OTUD7A expression without undergoing nonsense-mediated mRNA decay. The mutant iNs also exhibited marked reduction of dendritic complexity, density of synaptic proteins GluA1 and PSD-95, and neuronal network activity. Congruent with the neuronal phenotypes in mutant iNs, our transcriptomic analysis showed that the set of OTUD7A LoF-downregulated genes was enriched for those relating to synapse development and function and was associated with SZ and other neuropsychiatric disorders. These results suggest that OTUD7A LoF impairs synapse development and neuronal function in human neurons, providing mechanistic insight into the possible role of OTUD7A in driving neuropsychiatric phenotypes associated with the 15q13.3 deletion.

Motivation There is widespread interest in identifying genetic variants that exhibit parent-of-origin effects (POEs) wherein the effect of an allele on phenotype expression depends on its parental origin. POEs can arise from different phenomena including genomic imprinting and have been documented for many complex traits. Traditional tests for POEs require family data to determine parental origins of transmitted alleles. As most genome-wide association studies (GWAS) sample unrelated individuals (where allelic parental origin is unknown), the study of POEs in such datasets requires sophisticated statistical methods that exploit genetic patterns we anticipate observing when POEs exist. We propose a method to improve discovery of POE variants in large-scale GWAS samples that leverages potential pleiotropy among multiple correlated traits often collected in such studies. Our method compares the phenotypic covariance matrix of heterozygotes to homozygotes based on a Robust Omnibus Test. We refer to our method as the Parent of Origin Inference using Robust Omnibus Test (POIROT) of multiple quantitative traits. Results Through simulation studies, we compared POIROT to a competing univariate variance-based method which considers separate analysis of each phenotype. We observed POIROT to be well-calibrated with improved power to detect POEs compared to univariate methods. POIROT is robust to non-normality of phenotypes and can adjust for population stratification and other confounders. Finally, we applied POIROT to GWAS data from the UK Biobank using BMI and two cholesterol phenotypes. We identified 338 genome-wide significant loci for follow-up investigation. Availability and implementation The code for this method is available at https://github.com/staylorhead/POIROT-POE.

Background Crohn’s disease (CD) is a heritable chronic inflammatory disorder. Non-coding RNAs (ncRNAs) play an important role in epigenetic regulation by affecting gene expression, but can also directly affect protein function, thus having a substantial impact on biological processes. We investigated whether non-coding RNAs (ncRNA) at diagnosis are dysregulated during CD at different CD locations and future disease behaviors to determine if ncRNA signatures can serve as an index to outcomes. Methods Using subjects belonging to the RISK cohort, we analyzed ncRNA from the ileal biopsies of 345 CD and 71 non-IBD controls, and ncRNA from rectal biopsies of 329 CD and 61 non-IBD controls. Sequence alignment was done (STAR package) using Human Genome version 38 (hg38) as reference panel. The differential expression (DE) analysis was performed with EdgeR package and DE ncRNAs were identified with a threshold of fold change (FC) > 2 and FDR Results In total, we identified 130 CD specific DE ncRNAs (89 in ileum and 41 in rectum) when compared to non-IBD controls. Similarly, 35 DE ncRNAs were identified between B1 and B2 in ileum, whereas no differences among CD disease behaviors were noticed in rectum. We also found inflammation specific ncRNAs between inflamed and non-inflamed groups in ileal biopsies. Overall, we observed that expression of mir1244-2, mir1244-3, mir1244-4, and RN7SL2 were increased during CD, regardless of disease behavior, location, or inflammatory status. Lastly, we tested ncRNA expression at baseline as potential tool to predict the disease status, disease behaviors and disease inflammation at 3-year follow up. Conclusions We have identified ncRNAs that are specific to disease location, disease behavior, and disease inflammation in CD. Both ileal and rectal specific ncRNA are changing over the course of CD, specifically during the disease progression in the intestinal mucosa. Collectively, our findings show changes in ncRNA during CD and may have a clinical utility in early identification and characterization of disease progression.

PurposeOrofacial clefts (OFCs) are common birth defects including cleft lip (CL), cleft lip and palate (CLP), and cleft palate (CP). OFCs have heterogeneous etiologies, complicating clinical diagnostics as it is not always apparent if the cause is Mendelian, environmental, or multifactorial. Sequencing is not currently performed for isolated or sporadic OFCs, so we estimated the diagnostic yield for 418 genes in 841 cases and 294 controls.MethodsWe evaluated 418 genes using genome sequencing and curated variants to assess their pathogenicity using American College of Medical Genetics criteria.Results9.04% of cases and 1.02% of controls had ‘likely pathogenic’ (LP) variants (pCTNND1andIRF6, accounted for more than half of the yield (4.64% of cases). Most variants (61.8%) were ‘variants of uncertain significance’ (VUSs), occurring more frequently in cases (p=0.004), but no individual gene showed a significant excess of VUSs.ConclusionThese results underscore the etiological heterogeneity of OFCs and suggest sequencing could reduce the diagnostic gap in OFCs.

Background The gut and oral microbiome have independently been shown to be associated with inflammatory bowel disease (IBD). However, it is not known to what extent gut and oral microbial disease markers converge in terms of their composition in IBD. Further, the spatial and temporal variation within the oral microenvironments of IBD remain to be elucidated. Patients and Methods We used a prospectively recruited cohort of patients with IBD (n = 47) and unrelated healthy control patients (n = 18) to examine the spatial and temporal distribution of microbiota within the various oral microenvironments, represented by saliva, tongue, buccal mucosa, and plaque, and compared them with stool. Microbiome characterization was performed using 16S rRNA gene sequencing. Results The oral microbiome displayed IBD-associated dysbiosis, in a site- and taxa-specific manner. Plaque samples depicted a relatively severe degree of dysbiosis, and the disease-associated dysbiotic bacterial groups were predominantly the members of the phylum Firmicutes. Our 16S rRNA gene analyses show that oral microbiota can distinguish patients with IBD from healthy control patients, with salivary microbiota performing the best, closely matched by stool and other oral sites. Longitudinal profiles of microbial composition suggest that some taxa are more consistently perturbed than others, preferentially in a site-dependent fashion. Conclusions Collectively, these data indicate the potential of using oral microbial profiles in screening and monitoring patients with IBD. Furthermore, these results support the importance of spatial and longitudinal microbiome sampling to interpret disease-associated dysbiotic states and eventually to gain insights into disease pathogenesis.

Significance Expansion of 55-200 CGG repeats in the 5′ untranslated region of FMR1 predisposes carriers to fragile X–associated tremor/ataxia syndrome (FXTAS), a late-onset neurodegenerative disorder. FXTAS demonstrates incomplete penetrance, which strongly suggests the presence of genetic modifiers. We performed whole-genome sequencing (WGS) on male premutation carriers (CGG 55–200 ) followed by a functional screen in Drosophila and identified PSMB5 as a strong suppressor of CGG-associated neurodegeneration, thereby presenting a therapeutic strategy for FXTAS.

Fragile X–associated tremor/ataxia syndrome (FXTAS) is a debilitating late-onset neurodegenerative disease in premutation carriers of the expanded CGG repeat in FMR1 that presents with a spectrum of neurological manifestations, such as gait ataxia, intention tremor, and parkinsonism [P. J. Hagerman, R. J. Hagerman, Ann. N. Y. Acad. Sci. 1338, 58–70 (2015); S. Jacquemont et al., JAMA 291, 460–469 (2004)]. Here, we performed whole-genome sequencing (WGS) on male premutation carriers (CGG55–200) and prioritized candidate variants to screen for candidate genetic modifiers using a Drosophila model of FXTAS. We found 18 genes that genetically modulate CGG-associated neurotoxicity in Drosophila, such as Prosbeta5 (PSMB5), pAbp (PABPC1L), e(y)1 (TAF9), and CG14231 (OSGEPL1). Among them, knockdown of Prosbeta5 (PSMB5) suppressed CGG-associated neurodegeneration in the fly as well as in N2A cells. Interestingly, an expression quantitative trait locus variant in PSMB5, PSMB5rs11543947-A, was found to be associated with decreased expression of PSMB5 and delayed onset of FXTAS in human FMR1 premutation carriers. Finally, we demonstrate evidence that PSMB5 knockdown results in suppression of CGG neurotoxicity via both the RAN translation and RNA-mediated toxicity mechanisms, thereby presenting a therapeutic strategy for FXTAS.

We report the genetic analysis of two consanguineous pedigrees of Pakistani ancestry in which two siblings in each family exhibited developmental delay, epilepsy, intellectual disability and aggressive behavior. Whole-genome sequencing was performed in Family 1, and we identified ~80,000 variants located in regions of homozygosity. Of these, 615 variants had a minor allele frequency ≤ 0.001, and 21 variants had CADD scores ≥ 15. Four homozygous exonic variants were identified in both affected siblings

In peripheral blood, DNA methylation (DNAm) patterns in inflammatory bowel disease patients reflect inflammatory status rather than disease status. Here, we studied DNAm in diseased rectal mucosa from ulcerative colitis (UC) patients, conducting a cell-type-specific EWAS in epithelial, immune and fibroblast cells to understand DNAm changes across disease states, course, and clinical outcomes. At diagnosis, rectal mucosa in UC exhibited a lower proportion of epithelial cells and fibroblasts, and higher proportion of immune cells, along with 3504, 910, and 2279 altered DNAm sites as detected in our cell-specific EWAS, respectively. While treatment had significant effects on DNAm of immune cells, its effects on fibroblasts and epithelial cells were attenuated. Those requiring colectomy exhibited cell composition and DNAm patterns at follow-up more like disease onset than patients who did not require colectomy. Collectively DNAm and gene expression analysis suggest that targeting epithelial genes involved in barrier function may improve clinical outcomes.

ObjectiveAPOEis a strong risk factor for Alzheimer’s disease (AD) and associated with higher low-density lipoprotein cholesterol (LDL-C) levels. Moreover, LDL-C is associated with the development of clinically ascertained AD; however, whether this association is present with the underlying neuropathological manifestations of AD or whether it is independent of the effect ofAPOEis unknown and is the focus of this paper.MethodsIndividuals in the Religious Orders Study/Memory and Ageing Project cohorts with longitudinal measures of blood lipids and detailed autopsies were studied. We modelled the relationship between blood lipids and 12 age-related brain pathologies using a linear mixed model adjusted for potential confounding factors and stratified byAPOEgenotype with overall significance determined by meta-analysis. Blood lipids considered were LDL-C, high-density lipoprotein cholesterol and triglycerides. Brain pathologies included AD pathology measured by silver staining (Braak stage, a modified Consortium to Establish a Registry for Alzheimer’s Disease [CERAD] score and global AD pathology) and immunohistochemistry (beta-amyloid and neurofibrillary tangles) as well as cerebral microinfarct, cerebral macroinfarct, cerebral amyloid angiopathy, cerebral atherosclerosis, hippocampal sclerosis, TDP-43 cytoplasmic inclusions and Lewy bodies.Results559 participants (69.1% female) had complete data for analysis. They were followed for a median of 7 years and a median of 3 years prior to dementia onset. LDL-C was associated with all measures of AD neuropathology (neurofibrillary tangles, beta-amyloid, Braak stage, modified CERAD score and global AD pathology) and cerebral amyloid angiopathy independent ofAPOEafter adjusting for age, sex, cholesterol-lowering medication use, body mass index, smoking and education at false discovery rate (FDR) p-value ConclusionsThese findings implicate LDL-C in the pathophysiology of AD independent ofAPOEand suggest LDL-C is a modifiable risk factor for AD.

Genome-wide association studies (GWAS) have identified hundreds of loci associated with Crohns disease (CD), however, as with all complex diseases, deriving pathogenic mechanisms from these non-coding GWAS discoveries has been challenging. To complement GWAS and better define actionable biological targets, we analysed sequenced data from more than 30,000 CD patients and 80,000 population controls. We observe rare coding variants in established CD susceptibility genes as well as ten genes where coding variation directly implicates the gene in disease risk for the first time.

Background Structural rearrangements of the genome, which generally occur during meiosis and result in large-scale (> 1 kb) copy number variants (CNV; deletions or duplications ≥ 1 kb), underlie genomic disorders. Recurrent pathogenic CNVs harbor similar breakpoints in multiple unrelated individuals and are primarily formed via non-allelic homologous recombination (NAHR). Several pathogenic NAHR-mediated recurrent CNV loci demonstrate biases for parental origin of de novo CNVs. However, the mechanism underlying these biases is not well understood. Methods We performed a systematic, comprehensive literature search to curate parent of origin data for multiple pathogenic CNV loci. Using a regression framework, we assessed the relationship between parental CNV origin and the male to female recombination rate ratio. Results We demonstrate significant association between sex-specific differences in meiotic recombination and parental origin biases at these loci (p = 1.07 × 10–14). Conclusions Our results suggest that parental origin of CNVs is largely influenced by sex-specific recombination rates and highlight the need to consider these differences when investigating mechanisms that cause structural variation.

Biliary atresia (BA) is the most common cause of end-stage liver disease in children and the primary indication for pediatric liver transplantation, yet underlying etiologies remain unknown. Approximately 10% of infants affected by BA exhibit various laterality defects (heterotaxy) including splenic abnormalities and complex cardiac malformations-a distinctive subgroup commonly referred to as the biliary atresia splenic malformation (BASM) syndrome. We hypothesized that genetic factors linking laterality features with the etiopathogenesis of BA in BASM patients could be identified through whole-exome sequencing (WES) of an affected cohort. DNA specimens from 67 BASM subjects, including 58 patient-parent trios, from the National Institute of Diabetes and Digestive and Kidney Diseases-supported Childhood Liver Disease Research Network (ChiLDReN) underwent WES. Candidate gene variants derived from a prespecified set of 2,016 genes associated with ciliary dysgenesis and/or dysfunction or cholestasis were prioritized according to pathogenicity, population frequency, and mode of inheritance. Five BASM subjects harbored rare and potentially deleterious biallelic variants in polycystic kidney disease 1 like 1 (PKD1L1), a gene associated with ciliary calcium signaling and embryonic laterality determination in fish, mice, and humans. Heterozygous PKD1L1 variants were found in 3 additional subjects. Immunohistochemical analysis of liver from the one BASM subject available revealed decreased PKD1L1 expression in bile duct epithelium when compared to normal livers and livers affected by other noncholestatic diseases. Conclusion: WES identified biallelic and heterozygous PKD1L1 variants of interest in 8 BASM subjects from the ChiLDReN data set; the dual roles for PKD1L1 in laterality determination and ciliary function suggest that PKD1L1 is a biologically plausible, cholangiocyte-expressed candidate gene for the BASM syndrome.

Nonsyndromic orofacial clefts (OFCs) are a common birth defect and are phenotypically heterogenous in the structure affected by the cleft - cleft lip (CL) and cleft lip and palate (CLP) - as well as other features, such as the severity of the cleft. Here, we focus on bilateral and unilateral clefts as one dimension of OFC severity, because the genetic architecture of these subtypes is not well understood. We tested for subtype-specific genetic associations in 44 bilateral CL (BCL) cases, 434 unilateral CL (UCL) cases, 530 bilateral CLP cases (BCLP), 1123 unilateral CLP (UCLP) cases, and unrelated controls (N = 1626), using a mixed-model approach. While no novel loci were found, the genetic architecture of UCL was distinct compared to BCL, with 44.03% of suggestive loci having different effects between the two subtypes. To further understand the subtype-specific genetic risk factors, we performed a genome-wide scan for modifiers and found a significant modifier locus on 20p11 (p=7.53×10-9), 300kb downstream of PAX1, that associated with higher odds of BCL vs. UCL, and replicated in an independent cohort (p=0.0018) with no effect in BCLP (p>0.05). We further found that this locus was associated with normal human nasal shape. Taken together, these results suggest bilateral and unilateral clefts may have different genetic architectures. Moreover, our results suggest BCL, the rarest form of OFC, may be genetically distinct from the other OFC subtypes. This expands our understanding of modifiers for OFC subtypes and further elucidates the genetic mechanisms behind the phenotypic heterogeneity in OFCs. ispartof: Human Genetics and Genomics Advances vol:2 issue:2 ispartof: location:United States status: published

We generated an online brain pQTL resource for 7,376 proteins through the analysis of genetic and proteomic data derived from post-mortem samples of the dorsolateral prefrontal cortex of 330 older adults. The identified pQTLs tend to be non-synonymous variation, are over-represented among variants associated with brain diseases, and replicate well (77%) in an independent brain dataset. Comparison to a large study of brain eQTLs revealed that about 75% of pQTLs are also eQTLs. In contrast, about 40% of eQTLs were identified as pQTLs. These results are consistent with lower pQTL mapping power and greater evolutionary constraint on protein abundance. The latter is additionally supported by observations of pQTLs with large effects' tending to be rare, deleterious, and associated with proteins that have evidence for fewer protein-protein interactions. Mediation analyses using matched transcriptomic and proteomic data provided additional evidence that pQTL effects are often, but not always, mediated by mRNA. Specifically, we identified roughly 1.6 times more mRNA-mediated pQTLs than mRNA-independent pQTLs (550 versus 341). Our pQTL resource provides insight into the functional consequences of genetic variation in the human brain and a basis for novel investigations of genetics and disease.

Nonsyndromic orofacial clefts (OFCs) are the most common craniofacial birth defect in humans and, like many complex traits, OFCs are phenotypically and etiologically heterogenous. The phenotypic heterogeneity of OFCs extends beyond the structures affected by the cleft (e.g., cleft lip (CL) and cleft lip and palate (CLP) to other features, such as the severity of the cleft. Here, we focus on bilateral and unilateral clefts as one dimension of OFC severity. Unilateral clefts are more frequent than bilateral clefts for both CL and CLP, but the genetic architecture of these subtypes is not well understood, and it is not known if genetic variants predispose for the formation of one subtype over another. Therefore, we tested for subtype-specific genetic associations in 44 bilateral CL (BCL) cases, 434 unilateral CL (UCL) cases, 530 bilateral CLP cases (BCLP), 1123 unilateral CLP (UCLP) cases, and unrelated controls (N = 1626), using the mixed-model approach implemented in GENESIS. While no novel loci were found in subtype-specific analyses comparing cases to controls, the genetic architecture of UCL was distinct compared to BCL, with 43.8% of suggestive loci (p < 1.0×10−5) having non-overlapping confidence intervals between the two subtypes. To further understand the genetic risk factors for severity differences, we then performed a genome-wide scan for modifiers using a similar mixed-model approach and found one genome-wide significant modifier locus on 20p11 (p = 7.53×10−9), 300kb downstream of PAX1, associated with higher odds of BCL compared to UCL, which also replicated in an independent cohort (p = 0.0018) and showed no effect in BCLP (p>0.05). We further found that SNPs at this locus were associated with normal human nasal shape. Taken together, these results suggest bilateral and unilateral clefts may have differences in their genetic architecture, especially between CL and CLP. Moreover, our results suggest BCL, the rarest form of OFC, may be genetically distinct from the other OFC subtypes. This expands our understanding of genetic modifiers for subtypes of OFCs and further elucidates the genetic mechanisms behind the phenotypic heterogeneity in OFCs.

BackgroundIn a recent study, we identified 1189 CpG sites whose DNA methylation (DNAm) level in blood distinguished Crohn’s disease (CD) cases from controls. We also demonstrated that the vast majority of these differences were a consequence of disease, rather than a cause of CD. Since methylation can be influenced by both genetic and environmental factors, here we focus on CpGs under demonstrable genetic control (methylation quantitative trait loci, or mQTLs). By comparing mQTL patterns across disease states and tissue (blood vs. ileum), we may distinguish patterns unique to CD. Such DNAm patterns may be relevant for the developmental origins of CD.MethodsWe investigated three datasets: (i) 402 blood samples from 164 newly diagnosed pediatric CD patients taken at two time points, and 74 non-IBD controls (ii) 780 blood samples from a non-CD adult population and (iii) 40 ileal biopsies (17 CD cases and 23 non-IBD controls) from group (i). Genome-wide DNAm profiling and genotyping were performed using the Illumina MethylationEPIC and Illumina Multi-Ethnic arrays. SNP-CpG associations were tested via linear models adjusted for age, gender, disease status, disease subtype, estimated cell type and three genotype-based principal components. We used a Bonferroni-adjusted significance threshold to identify significantly associated SNP-CpG pairs, but also considered larger sets identified by a false discovery rate criterionResultsIn total, we observed 535,448 SNP-CpG associations between 287,881 SNPs and 12,843 CpG sites (P−14). These associations and their effects are highly consistent across different ages, races, disease states, and tissue types, suggesting that the vast majority of these mQTLs participate in common gene regulation. However, genes near CpGs associated with IBD SNPs were enriched for 18 KEGG pathways relevant to IBD-linked immune function and inflammatory responses. We observed suggestive evidence for a small number of tissue-specific associations and disease-specific ileal associations in ileum, though larger studies will be needed to confirm these results.ConclusionThe vast majority of blood derived mQTLs are commonly shared across individuals. However, we have identified a subset of such, which may be involved in processes related to CD. Independent cohort studies will be required to validate these findings.

Atrioventricular septal defects (AVSD) are a severe congenital heart defect present in individuals with Down syndrome (DS) at a > 2000-fold increased prevalence compared to the general population. This study aimed to identify risk-associated genes and pathways and to examine a potential polygenic contribution to AVSD in DS. We analyzed a total cohort of 702 individuals with DS with or without AVSD, with genomic data from whole exome sequencing, whole genome sequencing, and/or array-based imputation. We utilized sequence kernel association testing and polygenic risk score (PRS) methods to examine rare and common variants. Our findings suggest that the Notch pathway, particularly NOTCH4, as well as genes involved in the ciliome including CEP290 may play a role in AVSD in DS. These pathways have also been implicated in DS-associated AVSD in prior studies. A polygenic component for AVSD in DS has not been examined previously. Using weights based on the largest genome-wide association study of congenital heart defects available (2594 cases and 5159 controls; all general population samples), we found PRS to be associated with AVSD with odds ratios ranging from 1.2 to 1.3 per standard deviation increase in PRS and corresponding liability r2 values of approximately 1%, suggesting at least a small polygenic contribution to DS-associated AVSD. Future studies with larger sample sizes will improve identification and quantification of genetic contributions to AVSD in DS.

SUMMARYGenomic disorders are caused by structural rearrangements of the genome that generally occur during meiosis1. Often the rearrangements result in large-scale (> 1 kb) copy number variants (CNV; deletions or duplications ≥ 1 kb)2,3. Recurrent pathogenic CNVs harbor similar breakpoints in multiple unrelated individuals and are primarily formed via non-allelic homologous recombination (NAHR)3,4. Several pathogenic NAHR-mediated recurrent CNV loci demonstrate biases for parental origin of de novo CNVs5–9. However, the mechanism underlying these biases is not well understood. Here we have curated parent of origin data for multiple pathogenic CNV loci and demonstrate a significant association between sex-specific differences in meiotic recombination and parental origin biases at these loci. Our results suggest that parental-origin of CNVs is largely controlled by sex-specific recombination rates and bring into light the need to consider these differences when seeking to determine the factors underlying risk for structural variation.

While de novo mutations (DNMs) are known to increase risk of congenital defects, DNMs have not been fully explored regarding orofacial clefts (OFCs), one of the most common human birth defects. Therefore, whole-genome sequencing of 756 case-parent trios of European, Colombian, and Taiwanese ancestry was performed to determine the contributions of coding DNMs to OFC risk. Overall, we identified a significant excess of loss-of-function DNMs in genes highly expressed in craniofacial tissues, as well as genes associated with known autosomal dominant OFC syndromes. This analysis also revealed roles for zinc-finger homeobox domain and SOX2-interacting genes in OFC etiology.

Although de novo mutations (DNMs) are known to increase an individual's risk of congenital defects, DNMs have not been fully explored regarding orofacial clefts (OFCs), one of the most common human birth defects. Therefore, whole-genome sequencing of 756 child-parent trios of European, Colombian, and Taiwanese ancestry was performed to determine the contributions of coding DNMs to an individual's OFC risk. Overall, we identified a significant excess of loss-of-function DNMs in genes highly expressed in craniofacial tissues, as well as genes associated with known autosomal dominant OFC syndromes. This analysis also revealed roles for zinc-finger homeobox domain and SOX2-interacting genes in OFC etiology.

In our recent Autism Sequencing Consortium (ASC) study (Satterstrom et al 2020), we show that cohorts ascertained for autism spectrum disorder (ASD) and those ascertained for a broader array of neurodevelopmental disorders (NDDs) differ both in the rate and in the relative distribution of mutations across genes. In secondary analyses, we used a binary classifier to separate the 102 ASD genes that we identified into those more often disrupted in individuals ascertained for ASD versus those more often disrupted in individuals ascertained for NDDs; we refer to these as ASD-predominant (ASDP) and ASDNDD genes, respectively. Protein truncating variants in ASDP genes are more likely to be found in parents, and are more likely to be inherited, consistent with reduced selective pressure acting against disruptive variants in ASDP genes. Moreover, in individuals ascertained for ASD, those harboring ASDNDD gene mutations show, on average, greater delays in walking and lower IQ, when compared to those harboring ASDP gene mutations. The results demonstrate that there will be genes that are more clearly associated with the extremes of a clinical phenotypic spectrum. We suggest that the discovery of genes associated with the extremes of the spectrum will provide further insights into pathways disrupted in ASD, and that such discovery is an important direction for future research. Here, in response to a recent commentary by our esteemed colleagues (Myers et al. 2020), we discuss the ASC findings in greater detail and place our results in the context of future research in ASD.

Schizophrenia occurs in about one in four individuals with 22q11.2 deletion syndrome (22q11.2DS). The aim of this International Brain and Behavior 22q11.2DS Consortium (IBBC) study was to identify genetic factors that contribute to schizophrenia, in addition to the ~20-fold increased risk conveyed by the 22q11.2 deletion. Using whole-genome sequencing data from 519 unrelated individuals with 22q11.2DS, we conducted genome-wide comparisons of common and rare variants between those with schizophrenia and those with no psychotic disorder at age ≥25 years. Available microarray data enabled direct comparison of polygenic risk for schizophrenia between 22q11.2DS and independent population samples with no 22q11.2 deletion, with and without schizophrenia (total n = 35,182). Polygenic risk for schizophrenia within 22q11.2DS was significantly greater for those with schizophrenia (padj = 6.73 × 10−6). Novel reciprocal case–control comparisons between the 22q11.2DS and population-based cohorts showed that polygenic risk score was significantly greater in individuals with psychotic illness, regardless of the presence of the 22q11.2 deletion. Within the 22q11.2DS cohort, results of gene-set analyses showed some support for rare variants affecting synaptic genes. No common or rare variants within the 22q11.2 deletion region were significantly associated with schizophrenia. These findings suggest that in addition to the deletion conferring a greatly increased risk to schizophrenia, the risk is higher when the 22q11.2 deletion and common polygenic risk factors that contribute to schizophrenia in the general population are both present. ispartof: Molecular Psychiatry vol:26 issue:8 pages:1-15 ispartof: location:England status: published

International audience; We present the largest exome sequencing study of autism spectrum disorder (ASD) to date (n = 35,584 total samples, 11,986 with ASD). Using an enhanced analytical framework to integrate de novo and case-control rare variation, we identify 102 risk genes at a false discovery rate of 0.1 or less. Of these genes, 49 show higher frequencies of disruptive de novo variants in individuals ascertained to have severe neurodevelopmental delay, whereas 53 show higher frequencies in individuals ascertained to have ASD; comparing ASD cases with mutations in these groups reveals phenotypic differences. Expressed early in brain development, most risk genes have roles in regulation of gene expression or neuronal communication (i.e., mutations effect neurodevelopmental and neurophysiological changes), and 13 fall within loci recurrently hit by copy number variants. In cells from the human cortex, expression of risk genes is enriched in excitatory and inhibitory neuronal lineages, consistent with multiple paths to an excitatory-inhibitory imbalance underlying ASD.

Author(s): Denson, Lee A; Jurickova, Ingrid; Karns, Rebekah; Shaw, Kelly A; Cutler, David J; Okou, David; Valencia, C Alexander; Dodd, Anne; Mondal, Kajari; Aronow, Bruce J; Haberman, Yael; Linn, Aaron; Price, Adam; Bezold, Ramona; Lake, Kathleen; Jackson, Kimberly; Walters, Thomas D; Griffiths, Anne; Baldassano, Robert N; Noe, Joshua D; Hyams, Jeffrey S; Crandall, Wallace V; Kirschner, Barbara S; Heyman, Melvin B; Snapper, Scott; Guthery, Stephen L; Dubinsky, Marla C; Leleiko, Neal S; Otley, Anthony R; Xavier, Ramnik J; Stevens, Christine; Daly, Mark J; Zwick, Michael E; Kugathasan, Subra | Abstract: BACKGROUND:Granulocyte-macrophage colony-stimulating factor auto-antibodies (GMAbs) suppress neutrophil-extrinsic GM-CSF signaling and increase risk for stricturing behavior in Crohn's disease (CD). We aimed to define clinical, genomic, and functional associations with neutrophil-intrinsic GM-CSF signaling. METHODS:Missense mutations in CSF2RA, CSF2RB, JAK2, STAT5A, and STAT5B were identified using whole-exome sequencing in 543 pediatric inflammatory bowel disease (IBD) patients. Neutrophil-intrinsic GM-CSF signaling was defined using the GM-CSF-induced STAT5 stimulation index (GMSI) in 180 pediatric IBD patients and 26 non-IBD controls. Reduced GM-CSF signaling (GMSI-Lo) was defined as the 20th percentile within the control group. Variation in neutrophil phospho-protein abundance, bacterial killing, and the global pattern of gene expression with the GMSI was determined. RESULTS:We validated 18 potentially damaging missense mutations in CSF2RA and CSF2RB. CSF2RA A17G carriage increased from 10% in those with intact neutrophil GMSI to 32% in those with low GMSI (P = 0.02). The frequency of reduced Staphylococcus aureus killing increased from 17% in those with intact neutrophil GMSI to 35% in GMSI-Lo neutrophils (P = 0.043). Crohn's disease neutrophils with low GMSI exhibited specific alterations in phospho-protein networks and genes regulating cytokine production, wound healing, and cell survival and proliferation. Stricturing behavior increased from 7% in patients with both low GMAb and intact GMSI to 64% in patients with both elevated GMAb and low GMSI (P l 0.0001). CONCLUSIONS:Low/normal neutrophil-intrinsic GM-CSF signaling is associated with CSF2RA missense mutations, alterations in gene expression networks, and higher rates of disease complications in pediatric CD.

Background & aimsIndividuals with monogenic disorders of phagocyte function develop chronic colitis that resembles Crohn's disease (CD). We tested for associations between mutations in genes encoding reduced nicotinamide adenine dinucleotide phosphate (NADPH) oxidases, neutrophil function, and phenotypes of CD in pediatric patients.MethodsWe performed whole-exome sequence analysis to identify mutations in genes encoding NADPH oxidases (such as CYBA, CYBB, NCF1, NCF2, NCF4, RAC1, and RAC2) using DNA from 543 pediatric patients with inflammatory bowel diseases. Blood samples were collected from an additional 129 pediatric patients with CD and 26 children without IBD (controls); we performed assays for neutrophil activation, reactive oxygen species (ROS) production, and bacteria uptake and killing. Whole-exome sequence analysis was performed using DNA from 46 of the children with CD to examine associations with NADPH gene mutations; RNA sequence analyses were performed using blood cells from 46 children with CD to test for variations in neutrophil gene expression associated with ROS production.ResultsWe identified 26 missense mutations in CYBA, CYBB, NCF1, NCF2, and NCF4. Patients with CD who carried mutations in these genes were 3-fold more likely to have perianal disease (P= .0008) and stricturing complications (P= .002) than children with CD without these mutations. Among patients withCD with none of these mutations, 9% had undergone abdominal surgery; among patients with mutations in these NADPH oxidase genes, 31% had undergone abdominal surgery (P=.0004). A higher proportion of neutrophils from children with CD had low ROS production (47%) than from controls (15%) among the 129 patients tested for ROS (P= .002). Minor alleles of the NADPH genes were detected in 7% of children with CD whose neutrophils produced normal levels of ROS vs 38% of children whose neutrophils produced low levels of ROS (P= .009). Neutrophils that produced low levels of ROS had specific alterations in genes that regulate glucose metabolism and antimicrobial responses.ConclusionsWe identified missense mutations in genes that encode NADPH oxidases in children with CD; these were associated with a more aggressive disease course and reduced ROS production by neutrophils from the patients.