VarSome: the human genomic variant search engine.

SUMMARY: VarSome.com is a search engine, aggregator and impact analysis tool for human genetic variation and a community-driven project aiming at sharing global expertise on human variants. AVAILABILITY AND IMPLEMENTATION: VarSome is freely available at http://varsome.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Natural variation in genome architecture among 205 Drosophila melanogaster Genetic Reference Panel lines.

The Drosophila melanogaster Genetic Reference Panel (DGRP) is a community resource of 205 sequenced inbred lines, derived to improve our understanding of the effects of naturally occurring genetic variation on molecular and organismal phenotypes. We used an integrated genotyping strategy to identify 4,853,802 single nucleotide polymorphisms (SNPs) and 1,296,080 non-SNP variants. Our molecular population genomic analyses show higher deletion than insertion mutation rates and stronger purifying selection on deletions. Weaker selection on insertions than deletions is consistent with our observed distribution of genome size determined by flow cytometry, which is skewed toward larger genomes. Insertion/deletion and single nucleotide polymorphisms are positively correlated with each other and with local recombination, suggesting that their nonrandom distributions are due to hitchhiking and background selection. Our cytogenetic analysis identified 16 polymorphic inversions in the DGRP. Common inverted and standard karyotypes are genetically divergent and account for most of the variation in relatedness among the DGRP lines. Intriguingly, variation in genome size and many quantitative traits are significantly associated with inversions. Approximately 50% of the DGRP lines are infected with Wolbachia, and four lines have germline insertions of Wolbachia sequences, but effects of Wolbachia infection on quantitative traits are rarely significant. The DGRP complements ongoing efforts to functionally annotate the Drosophila genome. Indeed, 15% of all D. melanogaster genes segregate for potentially damaged proteins in the DGRP, and genome-wide analyses of quantitative traits identify novel candidate genes. The DGRP lines, sequence data, genotypes, quality scores, phenotypes, and analysis and visualization tools are publicly available.

Genomic variation and its impact on gene expression in Drosophila melanogaster.

Understanding the relationship between genetic and phenotypic variation is one of the great outstanding challenges in biology. To meet this challenge, comprehensive genomic variation maps of human as well as of model organism populations are required. Here, we present a nucleotide resolution catalog of single-nucleotide, multi-nucleotide, and structural variants in 39 Drosophila melanogaster Genetic Reference Panel inbred lines. Using an integrative, local assembly-based approach for variant discovery, we identify more than 3.6 million distinct variants, among which were more than 800,000 unique insertions, deletions (indels), and complex variants (1 to 6,000 bp). While the SNP density is higher near other variants, we find that variants themselves are not mutagenic, nor are regions with high variant density particularly mutation-prone. Rather, our data suggest that the elevated SNP density around variants is mainly due to population-level processes. We also provide insights into the regulatory architecture of gene expression variation in adult flies by mapping cis-expression quantitative trait loci (cis-eQTLs) for more than 2,000 genes. Indels comprise around 10% of all cis-eQTLs and show larger effects than SNP cis-eQTLs. In addition, we identified two-fold more gene associations in males as compared to females and found that most cis-eQTLs are sex-specific, revealing a partial decoupling of the genomic architecture between the sexes as well as the importance of genetic factors in mediating sex-biased gene expression. Finally, we performed RNA-seq-based allelic expression imbalance analyses in the offspring of crosses between sequenced lines, which revealed that the majority of strong cis-eQTLs can be validated in heterozygous individuals.

Automated protein-DNA interaction screening of Drosophila regulatory elements.

Drosophila melanogaster has one of the best characterized metazoan genomes in terms of functionally annotated regulatory elements. To explore how these elements contribute to gene regulation, we need convenient tools to identify the proteins that bind to them. Here we describe the development and validation of a high-throughput yeast one-hybrid platform, which enables screening of DNA elements versus an array of full-length, sequence-verified clones containing over 85% of predicted Drosophila transcription factors. Using six well-characterized regulatory elements, we identified 33 transcription factor-DNA interactions of which 27 were previously unidentified. To simultaneously validate these interactions and locate the binding sites of involved transcription factors, we implemented a powerful microfluidics-based approach that enabled us to retrieve DNA-occupancy data for each transcription factor throughout the respective target DNA elements. Finally, we biologically validated several interactions and identified two new regulators of sine oculis gene expression and hence eye development.

WebPrInSeS: automated full-length clone sequence identification and verification using high-throughput sequencing data.

High-throughput sequencing (HTS) is revolutionizing our ability to obtain cheap, fast and reliable sequence information. Many experimental approaches are expected to benefit from the incorporation of such sequencing features in their pipeline. Consequently, software tools that facilitate such an incorporation should be of great interest. In this context, we developed WebPrInSeS, a web server tool allowing automated full-length clone sequence identification and verification using HTS data. WebPrInSeS encompasses two separate software applications. The first is WebPrInSeS-C which performs automated sequence verification of user-defined open-reading frame (ORF) clone libraries. The second is WebPrInSeS-E, which identifies positive hits in cDNA or ORF-based library screening experiments such as yeast one- or two-hybrid assays. Both tools perform de novo assembly using HTS data from any of the three major sequencing platforms. Thus, WebPrInSeS provides a highly integrated, cost-effective and efficient way to sequence-verify or identify clones of interest. WebPrInSeS is available at http://webprinses.epfl.ch/ and is open to all users.

Non-familial cardiomyopathies in Lebanon: exome sequencing results for five idiopathic cases.

BACKGROUND: Cardiomyopathies affect more than 0.5% of the general population. They are associated with high risk of sudden cardiac death, which can result from either heart failure or electrical abnormalities. Although different mechanisms underlie the various types of cardiomyopathies, a principal pathology is common to all and is usually at the level of the cardiac muscle. With a relatively high incidence rate in most countries, and a subsequent major health burden on both the families and governments, cardiomyopathies are gaining more attention by researchers and pharmaceutical companies as well as health government bodies. In Lebanon, there is no official data about the spectrum of the diseases in terms of their respective prevalence, clinical, or genetic profiles. METHODS: We used exome sequencing to unravel the genetic basis of idiopathic cases of cardiomyopathies in Lebanon, a relatively small country with high rates of consanguineous marriages. RESULTS: Five cases were diagnosed with different forms of cardiomyopathies, and exome sequencing revealed the presence of already documented or novel mutations in known genes in three cases: LMNA for an Emery Dreifuss Muscular Dystrophy case, PKP2 for an arrhythmogenic right ventricle dysplasia case, and MYPN for a dilated cardiomyopathy case. Interestingly two brothers with hypertrophic cardiomyopathy have a novel missense variation in NPR1, the gene encoding the natriuretic peptides receptor type I, not reported previously to be causing cardiomyopathies. CONCLUSION: Our results unravel novel mutations in known genes implicated in cardiomyopathies in Lebanon. Changes in clinical management however, require genetic profiling of a larger cohort of patients.

Integrating clinical and genetic approaches in the diagnosis of 46,XY disorders of sex development.

46,XY differences and/or disorders of sex development (DSD) are clinically and genetically heterogeneous conditions. Although complete androgen insensitivity syndrome has a strong genotype-phenotype correlation, the other types of 46,XY DSD are less well defined, and thus, the precise diagnosis is challenging. This study focused on comparing the relationship between clinical assessment and genetic findings in a cohort of well-phenotyped patients with 46,XY DSD. The study was an analysis of clinical investigations followed by genetic testing performed on 35 patients presenting to a single center. The clinical assessment included external masculinization score (EMS), endocrine profiling and radiological evaluation. Array-comparative genomic hybridization (array-CGH) and sequencing of DSD-related genes were performed. Using an integrated approach, reaching the definitive diagnosis was possible in 12 children. The correlation between clinical and genetic findings was higher in patients with a more severe phenotype (median EMS 2.5 vs 6; P = 0.04). However, in 13 children, at least one variant of uncertain significance was identified, and most times this variant did not correspond to the original clinical diagnosis. In three patients, the genetic studies guided further clinical assessment which resulted in a reclassification of initial clinical diagnosis. Furthermore, we identified eight patients harboring variants in more than one DSD genes, which was not seen in controls (2.5%; P = 0.0003). In summary, taking into account potential challenges in reaching the definitive diagnosis in 46,XY DSD, only integrated approach seems to be the best routine practice.

A recurrent splice-site mutation in EPHA2 causing congenital posterior nuclear cataract.

Intoduction: Inherited cataract, opacification of the lens, is the most common worldwide cause of blindness in children. We aimed to identify the genetic cause of autosomal dominant (AD) posterior nuclear cataract in a four generation British family. METHODS: Whole genome sequence (WGS) was performed on two affected and one unaffected individual of the family and further validated by direct sequencing. Haplotype analysis was performed via genotying. RESULTS: A splice-site mutation c.2826-9G>A in the gene EPHA2, encoding EPH receptor A2 was identified and found to co-segregate with disease. CONCLUSIONS: We have identified a recurrent splice-site mutation c.2826-9G>A in EPHA2 causing isolated posterior nuclear cataract, providing evidence of further phenotypic heterogeneity associated with this variant.

Genetic, molecular and physiological basis of variation in Drosophila gut immunocompetence.

Gut immunocompetence involves immune, stress and regenerative processes. To investigate the determinants underlying inter-individual variation in gut immunocompetence, we perform enteric infection of 140 Drosophila lines with the entomopathogenic bacterium Pseudomonas entomophila and observe extensive variation in survival. Using genome-wide association analysis, we identify several novel immune modulators. Transcriptional profiling further shows that the intestinal molecular state differs between resistant and susceptible lines, already before infection, with one transcriptional module involving genes linked to reactive oxygen species (ROS) metabolism contributing to this difference. This genetic and molecular variation is physiologically manifested in lower ROS activity, lower susceptibility to ROS-inducing agent, faster pathogen clearance and higher stem cell activity in resistant versus susceptible lines. This study provides novel insights into the determinants underlying population-level variability in gut immunocompetence, revealing how relatively minor, but systematic genetic and transcriptional variation can mediate overt physiological differences that determine enteric infection susceptibility.

GETPrime: a gene- or transcript-specific primer database for quantitative real-time PCR.

The vast majority of genes in humans and other organisms undergo alternative splicing, yet the biological function of splice variants is still very poorly understood in large part because of the lack of simple tools that can map the expression profiles and patterns of these variants with high sensitivity. High-throughput quantitative real-time polymerase chain reaction (qPCR) is an ideal technique to accurately quantify nucleic acid sequences including splice variants. However, currently available primer design programs do not distinguish between splice variants and also differ substantially in overall quality, functionality or throughput mode. Here, we present GETPrime, a primer database supported by a novel platform that uniquely combines and automates several features critical for optimal qPCR primer design. These include the consideration of all gene splice variants to enable either gene-specific (covering the majority of splice variants) or transcript-specific (covering one splice variant) expression profiling, primer specificity validation, automated best primer pair selection according to strict criteria and graphical visualization of the latter primer pairs within their genomic context. GETPrime primers have been extensively validated experimentally, demonstrating high transcript specificity in complex samples. Thus, the free-access, user-friendly GETPrime database allows fast primer retrieval and visualization for genes or groups of genes of most common model organisms, and is available at http://updepla1srv1.epfl.ch/getprime/. Database URL: http://deplanckelab.epfl.ch.