Genome sequencing data analysis for rare disease gene discovery.

Rare diseases occur in a smaller proportion of the general population, which is variedly defined as less than 200 000 individuals (US) or in less than 1 in 2000 individuals (Europe). Although rare, they collectively make up to approximately 7000 different disorders, with majority having a genetic origin, and affect roughly 300 million people globally. Most of the patients and their families undergo a long and frustrating diagnostic odyssey. However, advances in the field of genomics have started to facilitate the process of diagnosis, though it is hindered by the difficulty in genome data analysis and interpretation. A major impediment in diagnosis is in the understanding of the diverse approaches, tools and datasets available for variant prioritization, the most important step in the analysis of millions of variants to select a few potential variants. Here we present a review of the latest methodological developments and spectrum of tools available for rare disease genetic variant discovery and recommend appropriate data interpretation methods for variant prioritization. We have categorized the resources based on various steps of the variant interpretation workflow, starting from data processing, variant calling, annotation, filtration and finally prioritization, with a special emphasis on the last two steps. The methods discussed here pertain to elucidating the genetic basis of disease in individual patient cases via trio- or family-based analysis of the genome data. We advocate the use of a combination of tools and datasets and to follow multiple iterative approaches to elucidate the potential causative variant.

Genome-wide association study and trans-ethnic meta-analysis identify novel susceptibility loci for type 2 diabetes mellitus.

BACKGROUND: The genetic basis of type 2 diabetes (T2D) is under-investigated in the Middle East, despite the rapidly growing disease prevalence. We aimed to define the genetic determinants of T2D in Qatar. METHODS: Using whole genome sequencing of 11,436 participants (2765 T2D cases and 8671 controls) from the population-based Qatar Biobank (QBB), we conducted a genome-wide association study (GWAS) of T2D with and without body mass index (BMI) adjustment. RESULTS: We replicated 93 known T2D-associated loci in a BMI-unadjusted model, while 96 known loci were replicated in a BMI-adjusted model. The effect sizes and allele frequencies of replicated SNPs in the Qatari population generally concurred with those from European populations. We identified a locus specific to our cohort located between the APOBEC3H and CBX7 genes in the BMI-unadjusted model. Also, we performed a transethnic meta-analysis of our cohort with a previous GWAS on T2D in multi-ancestry individuals (180,834 T2D cases and 1,159,055 controls). One locus in DYNC2H1 gene reached genome-wide significance in the meta-analysis. Assessing polygenic risk scores derived from European- and multi-ancestries in the Qatari population showed higher predictive performance of the multi-ancestry panel compared to the European panel. CONCLUSION: Our study provides new insights into the genetic architecture of T2D in a Middle Eastern population and identifies genes that may be explored further for their involvement in T2D pathogenesis.

Genome-wide association study and polygenic score assessment of insulin resistance.

Insulin resistance (IR) and beta cell dysfunction are the major drivers of type 2 diabetes (T2D). Genome-Wide Association Studies (GWAS) on IR have been predominantly conducted in European populations, while Middle Eastern populations remain largely underrepresented. We conducted a GWAS on the indices of IR (HOMA2-IR) and beta cell function (HOMA2-%B) in 6,217 non-diabetic individuals from the Qatar Biobank (QBB; Discovery cohort; n = 2170, Replication cohort; n = 4047) with and without body mass index (BMI) adjustment. We also developed polygenic scores (PGS) for HOMA2-IR and compared their performance with a previously derived PGS for HOMA-IR (PGS003470). We replicated 11 loci that have been previously associated with HOMA-IR and 24 loci that have been associated with HOMA-%B, at nominal statistical significance. We also identified a novel locus associated with beta cell function near VEGFC gene, tagged by rs61552983 (P = 4.38 × 10-8). Moreover, our best performing PGS (Q-PGS4; Adj R2 = 0.233 ± 0.014; P = 1.55 x 10-3) performed better than PGS003470 (Adj R2 = 0.194 ± 0.014; P = 5.45 x 10-2) in predicting HOMA2-IR in our dataset. This is the first GWAS on HOMA2 and the first GWAS conducted in the Middle East focusing on IR and beta cell function. Herein, we report a novel locus in VEGFC that is implicated in beta cell dysfunction. Inclusion of under-represented populations in GWAS has potentials to provide important insights into the genetic architecture of IR and beta cell function.

Case Report: Phenotype-Gene Correlation in a Case of Novel Tandem 4q Microduplication With Short Stature, Speech Delay and Microcephaly.

We describe a sporadic case of a pure, tandem, interstitial chromosome 4q duplication, arr[hg19] 4q28.1q32.3 (127,008,069-165,250,477) x3 in a boy born at 36 weeks of gestation. He presented with microcephaly (head circumference <1st percentile), short stature (height <2nd percentile) and poor weight gain (weight <3rd percentile). Hypospadias and horseshoe shaped kidneys were also revealed following a urinary tract ultrasound. Biochemical analysis revealed normal growth hormone and thyroid hormone levels. While gross and fine motor skill development was in line with his age, speech delay was observed. This patient adds to a group of more than 30 cases of pure 4q tandem duplication with common and differing phenotypic presentations. Using a retrospective analysis of previous case studies alongside the current case and bioinformatics analysis of the duplicated region, we deduced the most likely dosage sensitive genes for some of the major phenotypes in the patient. The positive predictive value (PPV) was calculated for each gene and phenotype and was derived by comparing the previously reported patients who have gene duplications and an associated phenotype versus those who had the gene duplications but were unaffected. Thus, the growth retardation phenotype may be associated with NAA15 duplication, speech delay with GRIA2 and microcephaly with PLK4 duplication. Functional studies will help in confirming the observations and elucidating the mechanisms. However, our study highlights the importance of analysing case reports with pure duplications in defining phenotype-gene relationships and in improving our knowledge of the function of precise chromosomal regions.