Genomic architecture of autism from comprehensive whole-genome sequence annotation.

Fully understanding autism spectrum disorder (ASD) genetics requires whole-genome sequencing (WGS). We present the latest release of the Autism Speaks MSSNG resource, which includes WGS data from 5,100 individuals with ASD and 6,212 non-ASD parents and siblings (total n = 11,312). Examining a wide variety of genetic variants in MSSNG and the Simons Simplex Collection (SSC; n = 9,205), we identified ASD-associated rare variants in 718/5,100 individuals with ASD from MSSNG (14.1%) and 350/2,419 from SSC (14.5%). Considering genomic architecture, 52% were nuclear sequence-level variants, 46% were nuclear structural variants (including copy-number variants, inversions, large insertions, uniparental isodisomies, and tandem repeat expansions), and 2% were mitochondrial variants. Our study provides a guidebook for exploring genotype-phenotype correlations in families who carry ASD-associated rare variants and serves as an entry point to the expanded studies required to dissect the etiology in the ∼85% of the ASD population that remain idiopathic.

Increasing the diagnostic yield of childhood glaucoma cases recruited into the 100,000 Genomes Project.

Childhood glaucoma (CG) encompasses a heterogeneous group of genetic eye disorders that is responsible for approximately 5% of childhood blindness worldwide. Understanding the molecular aetiology is key to improving diagnosis, prognosis and unlocking the potential for optimising clinical management. In this study, we investigated 86 CG cases from 78 unrelated families of diverse ethnic backgrounds, recruited into the Genomics England 100,000 Genomes Project (GE100KGP) rare disease cohort, to improve the genetic diagnostic yield. Using the Genomics England/Genomic Medicine Centres (GE/GMC) diagnostic pipeline, 13 unrelated families were solved (13/78, 17%). Further interrogation using an expanded gene panel yielded a molecular diagnosis in 7 more unrelated families (7/78, 9%). This analysis effectively raises the total number of solved CG families in the GE100KGP to 26% (20/78 families). Twenty-five percent (5/20) of the solved families had primary congenital glaucoma (PCG), while 75% (15/20) had secondary CG; 53% of this group had non-acquired ocular anomalies (including iris hypoplasia, megalocornea, ectopia pupillae, retinal dystrophy, and refractive errors) and 47% had non-acquired systemic diseases such as cardiac abnormalities, hearing impairment, and developmental delay. CYP1B1 was the most frequently implicated gene, accounting for 55% (11/20) of the solved families. We identified two novel likely pathogenic variants in the TEK gene, in addition to one novel pathogenic copy number variant (CNV) in FOXC1. Variants that passed undetected in the GE100KGP diagnostic pipeline were likely due to limitations of the tiering process, the use of smaller gene panels during analysis, and the prioritisation of coding SNVs and indels over larger structural variants, CNVs, and non-coding variants.

Mapping the genetic landscape of treatable inherited metabolic disorders in a large Middle Eastern biobank.

PURPOSE: To date, approximately 1400 inherited metabolic disorders (IMDs) have been described, some of which are treatable. It is estimated that 2% to 3% of live births worldwide are affected by treatable IMDs. Roughly 80% of IMDs are autosomal recessive, leading to a potentially higher incidence in regions with high consanguinity. METHODOLOGY: The study utilized genome sequencing data from 14,060 adult Qatari participants who were recruited by the Qatar Biobank and sequenced by the Qatar Genome Program. The genome sequencing data were analyzed for 125 nuclear genes known to be associated with 115 treatable IMDs. RESULTS: Our study identified 253 pathogenic/likely pathogenic single-nucleotide variations associated with 69 treatable IMDs, including 211 known and 42 novel predicted loss-of-function variants. We estimated that approximately 1 in 13 unrelated individuals (8%) carry a heterozygous pathogenic variant for at least 1 of 46 treatable IMDs. Notably, phenylketonuria/hyperphenylalaninemia and homocystinuria had among the highest carrier frequencies (1 in 68 and 1 in 85, respectively). CONCLUSION: Population-based studies of treatable IMDs, particularly in globally under-studied populations, can identify high-frequency alleles segregating in the community and inform public health policies, including carrier and newborn screening.

A loss-of-function AGTR1 variant in a critically-ill infant with renal tubular dysgenesis: case presentation and literature review.

BACKGROUND: Renal tubular dysgenesis (RTD) is a severe disorder with poor prognosis significantly impacting the proximal tubules of the kidney while maintaining an anatomically normal gross structure. The genetic origin of RTD, involving variants in the ACE, REN, AGT, and AGTR1 genes, affects various enzymes or receptors within the Renin angiotensin system (RAS). This condition manifests prenatally with oligohydramninos and postnatally with persistent anuria, severe refractory hypotension, and defects in skull ossification. CASE PRESENTATION: In this report, we describe a case of a female patient who, despite receiving multi vasopressor treatment, experienced persistent hypotension, ultimately resulting in early death at five days of age. While there was a history of parental consanguinity, no reported family history of renal disease existed. Blood samples from the parents and the remaining DNA sample of the patient underwent Whole Genome Sequencing (WGS). The genetic analysis revealed a rare homozygous loss of function variant (NM_000685.5; c.415C > T; p.Arg139*) in the Angiotensin II Receptor Type 1 (AGTR1) gene. CONCLUSION: This case highlights the consequence of loss-of-function variants in AGTR1 gene leading to RTD, which is characterized by high mortality rate at birth or during the neonatal period. Furthermore, we provide a comprehensive review of previously reported variants in the AGTR1 gene, which is the least encountered genetic cause of RTD, along with their associated clinical features.

The link between glycemic control measures and eye microvascular complications in a clinical cohort of type 2 diabetes with microRNA-223-3p signature.

BACKGROUND: Type 2 diabetes (T2D) is a critical healthcare challenge and priority in Qatar which is listed amongst the top 10 countries in the world, with its prevalence presently at 17% double the global average. MicroRNAs (miRNAs) are implicated in the pathogenesis of (T2D) and long-term microvascular complications including diabetic retinopathy (DR). METHODS: In this study, a T2D cohort that accurately matches the characteristics of the general population was employed to find microRNA (miRNA) signatures that are correlated with glycemic and ß cell function measurements. Targeted miRNA profiling was performed in (471) T2D individuals with or without DR and (491) (non-diabetic) healthy controls from the Qatar Biobank. Discovery analysis identified 20 differentially expressed miRNAs in T2D compared to controls, of which miR-223-3p was significantly upregulated (fold change:5.16, p = 3.6e-02) and positively correlated with glucose and hemoglobin A1c (HbA1c) levels (p-value = 9.88e-04 and 1.64e-05, respectively), but did not show any significant associations with insulin or C-peptide. Accordingly, we performed functional validation using a miR-223-3p mimic (overexpression) under control and hyperglycemia-induced conditions in a zebrafish model. RESULTS: Over-expression of miR-223-3p alone was associated with significantly higher glucose (42.7 mg/dL, n = 75 vs 38.7 mg/dL, n = 75, p = 0.02) and degenerated retinal vasculature, and altered retinal morphology involving changes in the ganglion cell layer and inner and outer nuclear layers. Assessment of retinal angiogenesis revealed significant upregulation in the expression of vascular endothelial growth factor and its receptors, including kinase insert domain receptor. Further, the pancreatic markers, pancreatic and duodenal homeobox 1, and the insulin gene expressions were upregulated in the miR-223-3p group. CONCLUSION: Our zebrafish model validates a novel correlation between miR-223-3p and DR development. Targeting miR-223-3p in T2D patients may serve as a promising therapeutic strategy to control DR in at-risk individuals.

Assessing the genetic burden of familial hypercholesterolemia in a large middle eastern biobank.

BACKGROUND: The genetic architecture underlying Familial Hypercholesterolemia (FH) in Middle Eastern Arabs is yet to be fully described, and approaches to assess this from population-wide biobanks are important for public health planning and personalized medicine. METHODS: We evaluate the pilot phase cohort (n = 6,140 adults) of the Qatar Biobank (QBB) for FH using the Dutch Lipid Clinic Network (DLCN) criteria, followed by an in-depth characterization of all genetic alleles in known dominant (LDLR, APOB, and PCSK9) and recessive (LDLRAP1, ABCG5, ABCG8, and LIPA) FH-causing genes derived from whole-genome sequencing (WGS). We also investigate the utility of a globally established 12-SNP polygenic risk score to predict FH individuals in this cohort with Arab ancestry. RESULTS: Using DLCN criteria, we identify eight (0.1%) 'definite', 41 (0.7%) 'probable' and 334 (5.4%) 'possible' FH individuals, estimating a prevalence of 'definite or probable' FH in the Qatari cohort of ~ 1:125. We identify ten previously known pathogenic single-nucleotide variants (SNVs) and 14 putatively novel SNVs, as well as one novel copy number variant in PCSK9. Further, despite the modest sample size, we identify one homozygote for a known pathogenic variant (ABCG8, p. Gly574Arg, global MAF = 4.49E-05) associated with Sitosterolemia 2. Finally, calculation of polygenic risk scores found that individuals with 'definite or probable' FH have a significantly higher LDL-C SNP score than 'unlikely' individuals (p = 0.0003), demonstrating its utility in Arab populations. CONCLUSION: We design and implement a standardized approach to phenotyping a population biobank for FH risk followed by systematically identifying known variants and assessing putative novel variants contributing to FH burden in Qatar. Our results motivate similar studies in population-level biobanks - especially those with globally under-represented ancestries - and highlight the importance of genetic screening programs for early detection and management of individuals with high FH risk in health systems.

Diagnosis and treatment of type 1 diabetes at the dawn of the personalized medicine era.

Type 1 diabetes affects millions of people globally and requires careful management to avoid serious long-term complications, including heart and kidney disease, stroke, and loss of sight. The type 1 diabetes patient cohort is highly heterogeneous, with individuals presenting with disease at different stages and severities, arising from distinct etiologies, and overlaying varied genetic backgrounds. At present, the "one-size-fits-all" treatment for type 1 diabetes is exogenic insulin substitution therapy, but this approach fails to achieve optimal blood glucose control in many individuals. With advances in our understanding of early-stage diabetes development, diabetes stratification, and the role of genetics, type 1 diabetes is a promising candidate for a personalized medicine approach, which aims to apply "the right therapy at the right time, to the right patient". In the case of type 1 diabetes, great efforts are now being focused on risk stratification for diabetes development to enable pre-clinical detection, and the application of treatments such as gene therapy, to prevent pancreatic destruction in a sub-set of patients. Alongside this, breakthroughs in stem cell therapies hold great promise for the regeneration of pancreatic tissues in some individuals. Here we review the recent initiatives in the field of personalized medicine for type 1 diabetes, including the latest discoveries in stem cell and gene therapy for the disease, and current obstacles that must be overcome before the dream of personalized medicine for all type 1 diabetes patients can be realized.

Genome sequencing unveils mutational landscape of the familial Mediterranean fever: Potential implications of IL33/ST2 signalling.

Familial Mediterranean fever (FMF) is the most common auto-inflammatory disease. It is transmitted as autosomal recessive trait with mutations in MEditerranean FeVer (MEFV) gene. Despite a typical clinical expression, many patients have either a single or no mutation in MEFV. The current work is aimed to revisit the genetic landscape of FMF disease using high-coverage whole genome sequencing. In atypical patients (carrying a single or no mutation in MEFV), we revealed many rare variants in genes associated with auto-inflammatory disorders, and more interestingly, we discovered a novel variant ( a 2.1-Kb deletion) in exon 11 of IL1RL1 gene, present only in patients. To validate and screen this patient-specific variant, a tandem of allele-specific PCR and quantitative real-time PCR was performed in 184 FMF patients and 218 healthy controls and we demonstrated that the novel deletion was absent in controls and was present in more than 19% of patients. This study sheds more light on the mutational landscape of FMF. Our discovery of a disease-specific variant in IL1RL1 gene may constitute a novel genetic marker for FMF. This finding suggesting a potential role of the IL33/ST2 signalling in the disease pathogenicity highlights a new paradigm in FMF pathophysiology.

Rare variants at KCNJ2 are associated with LDL-cholesterol levels in a cross-population study.

Leveraging whole genome sequencing data of 1751 individuals from the UK and 2587 Qatari subjects, we suggest here an association of rare variants mapping to the sour taste-associated gene KCNJ2 with reduced low-density lipoprotein cholesterol (LDL-C, P = 2.10 × 10-12) and with a 22% decreased dietary trans-fat intake. This study identifies a novel candidate rare locus for LDL-C, adding insights into the genetic architecture of a complex trait implicated in cardiovascular disease.

Burden of Mendelian disorders in a large Middle Eastern biobank.

BACKGROUND: Genome sequencing of large biobanks from under-represented ancestries provides a valuable resource for the interrogation of Mendelian disease burden at world population level, complementing small-scale familial studies. METHODS: Here, we interrogate 6045 whole genomes from Qatar-a Middle Eastern population with high consanguinity and understudied mutational burden-enrolled at the national Biobank and phenotyped for 58 clinically-relevant quantitative traits. We examine a curated set of 2648 Mendelian genes from 20 panels, annotating known and novel pathogenic variants and assessing their penetrance and impact on the measured traits. RESULTS: We find that 62.5% of participants are carriers of at least 1 known pathogenic variant relating to recessive conditions, with homozygosity observed in 1 in 150 subjects (0.6%) for which Peninsular Arabs are particularly enriched versus other ancestries (5.8-fold). On average, 52.3 loss-of-function variants were found per genome, 6.5 of which affect a known Mendelian gene. Several variants annotated in ClinVar/HGMD as pathogenic appeared at intermediate frequencies in this cohort (1-3%), highlighting Arab founder effect, while others have exceedingly high frequencies (> 5%) prompting reconsideration as benign. Furthermore, cumulative gene burden analysis revealed 56 genes having gene carrier frequency > 1/50, including 5 ACMG Tier 3 panel genes which would be candidates for adding to newborn screening in the country. Additionally, leveraging 58 biobank traits, we systematically assess the impact of novel/rare variants on phenotypes and discover 39 candidate large-effect variants associating with extreme quantitative traits. Furthermore, through rare variant burden testing, we discover 13 genes with high mutational load, including 5 with impact on traits relevant to disease conditions, including metabolic disorder and type 2 diabetes, consistent with the high prevalence of these conditions in the region. CONCLUSIONS: This study on the first phase of the growing Qatar Genome Program cohort provides a comprehensive resource from a Middle Eastern population to understand the global mutational burden in Mendelian genes and their impact on traits in seemingly healthy individuals in high consanguinity settings.

Network-based identification and prioritization of key transcriptional factors of diabetic kidney disease.

Diabetic nephropathy (DN) is one of the most established microvascular complications of diabetes and a key cause of end-stage renal disease. It is well established that gene susceptibility to DN plays a critical role in disease pathophysiology. Therefore, many genetic studies have been performed to categorize candidate genes in prominent diabetic cohorts, aiming to investigate DN pathogenesis and etiology. In this study, we performed a meta-analysis on the expression profiles of GSE1009, GSE30122, GSE96804, GSE99340, GSE104948, GSE104954, and GSE111154 to identify critical transcriptional factors associated with DN progression. The analysis was conducted for all individual datasets for each kidney tissue (glomerulus, tubules, and kidney cortex). We identified distinct clusters of susceptibility genes that were dysregulated in a renal compartment-specific pattern. Further, we recognized a small but a closely connected set of these susceptibility genes enriched for podocyte differentiation, several of which were characterized as genes encoding critical transcriptional factors (TFs) involved in DN development and podocyte function. To validate the role of identified TFs in DN progression, we functionally validated the three main TFs (DACH1, LMX1B, and WT1) identified through differential gene expression and network analysis using the hyperglycemic zebrafish model. We report that hyperglycemia-induced altered gene expression of the key TF genes leads to morphological abnormalities in zebrafish glomeruli, pronephric tubules, proximal and distal ducts. This study demonstrated that altered expression of these TF genes could be associated with hyperglycemia-induced nephropathy and, thus, aids in understanding the molecular drivers, essential genes, and pathways that trigger DN initiation and development.

A Complex Intrachromosomal Rearrangement Disrupting IRF6 in a Family with Popliteal Pterygium and Van der Woude Syndromes.

Clefts of the lip and/or palate (CL/P) are considered the most common form of congenital anomalies occurring either in isolation or in association with other clinical features. Van der woude syndrome (VWS) is associated with about 2% of all CL/P cases and is further characterized by having lower lip pits. Popliteal pterygium syndrome (PPS) is a more severe form of VWS, normally characterized by orofacial clefts, lower lip pits, skin webbing, skeletal anomalies and syndactyly of toes and fingers. Both syndromes are inherited in an autosomal dominant manner, usually caused by heterozygous mutations in the Interferon Regulatory Factor 6 (IRF6) gene. Here we report the case of a two-generation family where the index presented with popliteal pterygium syndrome while both the father and sister had clinical features of van der woude syndrome, but without any point mutations detected by re-sequencing of known gene panels or microarray testing. Using whole genome sequencing (WGS) followed by local de novo assembly, we discover and validate a copy-neutral, 429 kb complex intra-chromosomal rearrangement in the long arm of chromosome 1, disrupting the IRF6 gene. This variant is copy-neutral, novel against publicly available databases, and segregates in the family in an autosomal dominant pattern. This finding suggests that missing heritability in rare diseases may be due to complex genomic rearrangements that can be resolved by WGS and de novo assembly, helping deliver answers to patients where no genetic etiology was identified by other means.

Genomic architecture of autism spectrum disorder in Qatar: The BARAKA-Qatar Study.

BACKGROUND: Autism spectrum disorder (ASD) is a neurodevelopmental condition characterized by impaired social and communication skills, restricted interests, and repetitive behaviors. The prevalence of ASD among children in Qatar was recently estimated to be 1.1%, though the genetic architecture underlying ASD both in Qatar and the greater Middle East has been largely unexplored. Here, we describe the first genomic data release from the BARAKA-Qatar Study-a nationwide program building a broadly consented biorepository of individuals with ASD and their families available for sample and data sharing and multi-omics research. METHODS: In this first release, we present a comprehensive analysis of whole-genome sequencing (WGS) data of the first 100 families (372 individuals), investigating the genetic architecture, including single-nucleotide variants (SNVs), copy number variants (CNVs), tandem repeat expansions (TREs), as well as mitochondrial DNA variants (mtDNA) segregating with ASD in local families. RESULTS: Overall, we identify potentially pathogenic variants in known genes or regions in 27 out of 100 families (27%), of which 11 variants (40.7%) were classified as pathogenic or likely-pathogenic based on American College of Medical Genetics (ACMG) guidelines. Dominant variants, including de novo and inherited, contributed to 15 (55.6%) of these families, consisting of SNVs/indels (66.7%), CNVs (13.3%), TREs (13.3%), and mtDNA variants (6.7%). Moreover, homozygous variants were found in 7 families (25.9%), with a sixfold increase in homozygous burden in consanguineous versus non-consanguineous families (13.6% and 1.8%, respectively). Furthermore, 28 novel ASD candidate genes were identified in 20 families, 23 of which had recurrent hits in MSSNG and SSC cohorts. CONCLUSIONS: This study illustrates the value of ASD studies in under-represented populations and the importance of WGS as a comprehensive tool for establishing a molecular diagnosis for families with ASD. Moreover, it uncovers a significant role for recessive variation in ASD architecture in consanguineous settings and provides a unique resource of Middle Eastern genomes for future research to the global ASD community.

The third international hackathon for applying insights into large-scale genomic composition to use cases in a wide range of organisms.

In October 2021, 59 scientists from 14 countries and 13 U.S. states collaborated virtually in the Third Annual Baylor College of Medicine & DNANexus Structural Variation hackathon. The goal of the hackathon was to advance research on structural variants (SVs) by prototyping and iterating on open-source software. This led to nine hackathon projects focused on diverse genomics research interests, including various SV discovery and genotyping methods, SV sequence reconstruction, and clinically relevant structural variation, including SARS-CoV-2 variants. Repositories for the projects that participated in the hackathon are available at https://github.com/collaborativebioinformatics.

Ethnic-specific association of amylase gene copy number with adiposity traits in a large Middle Eastern biobank.

Studies assessing the impact of amylase genes copy number (CN) on adiposity report conflicting findings in different global populations, likely reflecting the impact of ancestral and ethnic-specific environment and lifestyle on selection at the amylase loci. Here, we leverage population size and detailed adiposity measures from a large population biobank to resolve confounding effects and determine the relationship between salivary (AMY1) and pancreatic (AMY2A) amylase genes CN and adiposity in 2935 Qatari individuals who underwent whole-genome sequencing (WGS) as part of the Qatar Genome Programme. We observe a negative association between AMY1 CNs and trunk fat percentage in the Qatari population (P = 7.50 × 10-3) and show that Qataris of Arab descent have significantly lower CN at AMY1 (P = 1.32 × 10-10) as well as less favorable adiposity and metabolic profiles (P < 1.34 × 10-8) than Qataris with Persian ancestry. Indeed, lower AMY1 CN was associated with increased total and trunk fat percentages in Arabs (P < 4.60 × 10-3) but not in Persians. Notably, overweight and obese Persians reported a significant trend towards dietary restraint following weight gain compared to Arabs (P = 4.29 × 10-5), with AMY1 CN showing negative association with dietary self-restraint (P = 3.22 × 10-3). This study reports an association between amylase gene CN and adiposity traits in a large Middle Eastern population. Importantly, we leverage rich biobank data to demonstrate that the strength of this association varies with ethnicity, and may be influenced by population-specific behaviors that also contribute to adiposity traits.

Thousands of Qatari genomes inform human migration history and improve imputation of Arab haplotypes.

Arab populations are largely understudied, notably their genetic structure and history. Here we present an in-depth analysis of 6,218 whole genomes from Qatar, revealing extensive diversity as well as genetic ancestries representing the main founding Arab genealogical lineages of Qahtanite (Peninsular Arabs) and Adnanite (General Arabs and West Eurasian Arabs). We find that Peninsular Arabs are the closest relatives of ancient hunter-gatherers and Neolithic farmers from the Levant, and that founder Arab populations experienced multiple splitting events 12-20 kya, consistent with the aridification of Arabia and farming in the Levant, giving rise to settler and nomadic communities. In terms of recent genetic flow, we show that these ancestries contributed significantly to European, South Asian as well as South American populations, likely as a result of Islamic expansion over the past 1400 years. Notably, we characterize a large cohort of men with the ChrY J1a2b haplogroup (n = 1,491), identifying 29 unique sub-haplogroups. Finally, we leverage genotype novelty to build a reference panel of 12,432 haplotypes, demonstrating improved genotype imputation for both rare and common alleles in Arabs and the wider Middle East.

Homozygous duplication identified by whole genome sequencing causes LRBA deficiency.

In more than one-third of primary immunodeficiency (PID) patients, extensive genetic analysis including whole-exome sequencing (WES) fails to identify the genetic defect. Whole-genome sequencing (WGS) is able to detect variants missed by other genomics platforms, enabling the molecular diagnosis of otherwise unresolved cases. Here, we report two siblings, offspring of consanguineous parents, who experienced similar severe events encompassing early onset of colitis, lymphoproliferation, and hypogammaglobulinemia, typical of lipopolysaccharide-responsive and beige-like anchor (LRBA) or cytotoxic T lymphocyte antigen 4 (CTLA4) deficiencies. Gene-panel sequencing, comparative genomic hybridization (CGH) array, and WES failed to reveal a genetic aberration in relevant genes. WGS of these patients detected a 12.3 kb homozygous tandem duplication that was absent in control cohorts and is predicted to disrupt the reading frame of the LRBA gene. The variant was validated by PCR and Sanger sequencing, demonstrating the presence of the junction between the reference and the tandem-duplicated sequence. Droplet digital PCR (ddPCR) further confirmed the copy number in the unaffected parents (CN = 3, heterozygous) and affected siblings (CN = 4, homozygous), confirming the expected segregation pattern. In cases of suspected inherited immunodeficiency, WGS may reveal a mutation when other methods such as microarray and WES analysis failed to detect an aberration.

An international virtual hackathon to build tools for the analysis of structural variants within species ranging from coronaviruses to vertebrates.

In October 2020, 62 scientists from nine nations worked together remotely in the Second Baylor College of Medicine & DNAnexus hackathon, focusing on different related topics on Structural Variation, Pan-genomes, and SARS-CoV-2 related research.   The overarching focus was to assess the current status of the field and identify the remaining challenges. Furthermore, how to combine the strengths of the different interests to drive research and method development forward. Over the four days, eight groups each designed and developed new open-source methods to improve the identification and analysis of variations among species, including humans and SARS-CoV-2. These included improvements in SV calling, genotyping, annotations and filtering. Together with advancements in benchmarking existing methods. Furthermore, groups focused on the diversity of SARS-CoV-2. Daily discussion summary and methods are available publicly at  https://github.com/collaborativebioinformatics provides valuable insights for both participants and the research community.