Kaviar: an accessible system for testing SNV novelty

Mutations of RagA GTPase in mTORC1 Pathway Are Associated with Autosomal Dominant Cataracts

Cataracts are a significant public health problem with no proven methods for prevention. Discovery of novel disease mechanisms to delineate new therapeutic targets is of importance in cataract prevention and therapy. Herein, we report that mutations in the RagA GTPase (RRAGA), a key regulator of the mechanistic rapamycin complex 1 (mTORC1), are associated with autosomal dominant cataracts. We performed whole exome sequencing in a family with autosomal dominant juvenile-onset cataracts, and identified a novel p.Leu60Arg mutation in RRAGA that co-segregated with the disease, after filtering against the dbSNP database, and at least 123,000 control chromosomes from public and in-house exome databases. In a follow-up direct screening of RRAGA in another 22 families and 142 unrelated patients with congenital or juvenile-onset cataracts, RRAGA was found to be mutated in two unrelated patients (p.Leu60Arg and c.-16G>A respectively). Functional studies in human lens epithelial cells revealed that the RRAGA mutations exerted deleterious effects on mTORC1 signaling, including increased relocation of RRAGA to the lysosomes, up-regulated mTORC1 phosphorylation, down-regulated autophagy, altered cell growth or compromised promoter activity. These data indicate that the RRAGA mutations, associated with autosomal dominant cataracts, play a role in the disease by acting through disruption of mTORC1 signaling.

A group of guanine nucleotide-binding molecules called Rag GTPases are known to play a crucial role in regulation of mTORC1 signaling cascade. In the current study, whole exome sequencing has led to the identification of the RagA GTPase (RRAGA) gene for cataracts and we proceeded to study properties of the RRAGA protein. We captured and sequenced the whole exome for four affected patients from a family with autosomal dominant juvenile-onset posterior cataracts, and found a novel rare mutation in RagA GTPase (RRAGA). To validate this finding, we then sequenced more families and patients, and observed RRAGA mutations in unrelated patients with related phenotypes, suggesting that RRAGA could be mutated in congenital and juvenile-onset cataracts. We further demonstrated supporting evidence that in human lens epithelial cells the RRAGA mutations exerted deleterious effects on relocation of RRAGA to the lysosomes, mTORC1 phosphorylation, autophagy and cell growth. This study gives important new insight into the roles of RRAGA and mTROC1 signaling in the etiology of cataracts.

Genomic architecture of inflammatory bowel disease in five families with multiple affected individuals

Currently, the best clinical predictor for inflammatory bowel disease (IBD) is family history. Over 163 sequence variants have been associated with IBD in genome-wide association studies, but they have weak effects and explain only a fraction of the observed heritability. It is expected that additional variants contribute to the genomic architecture of IBD, possibly including rare variants with effect sizes larger than the identified common variants. Here we applied a family study design and sequenced 38 individuals from five families, under the hypothesis that families with multiple IBD-affected individuals harbor one or more risk variants that (i) are shared among affected family members, (ii) are rare and (iii) have substantial effect on disease development. Our analysis revealed not only novel candidate risk variants but also high polygenic risk scores for common known risk variants in four out of the five families. Functional analysis of our top novel variant in the remaining family, a rare missense mutation in the ubiquitin ligase TRIM11, suggests that it leads to increased nuclear factor of kappa light chain enhancer in B-cells (NF-κB) signaling. We conclude that an accumulation of common weak-effect variants accounts for the high incidence of IBD in most, but not all families we analyzed and that a family study design can identify novel rare variants conferring risk for IBD with potentially large effect size, such as the TRIM11 p.H414Y mutation.

Privacy Risks from Genomic Data-Sharing Beacons

The human genetics community needs robust protocols that enable secure sharing of genomic data from participants in genetic research. Beacons are web servers that answer allele-presence queries—such as “Do you have a genome that has a specific nucleotide (e.g., A) at a specific genomic position (e.g., position 11,272 on chromosome 1)?”—with either “yes” or “no.” Here, we show that individuals in a beacon are susceptible to re-identification even if the only data shared include presence or absence information about alleles in a beacon. Specifically, we propose a likelihood-ratio test of whether a given individual is present in a given genetic beacon. Our test is not dependent on allele frequencies and is the most powerful test for a specified false-positive rate. Through simulations, we showed that in a beacon with 1,000 individuals, re-identification is possible with just 5,000 queries. Relatives can also be identified in the beacon. Re-identification is possible even in the presence of sequencing errors and variant-calling differences. In a beacon constructed with 65 European individuals from the 1000 Genomes Project, we demonstrated that it is possible to detect membership in the beacon with just 250 SNPs. With just 1,000 SNP queries, we were able to detect the presence of an individual genome from the Personal Genome Project in an existing beacon. Our results show that beacons can disclose membership and implied phenotypic information about participants and do not protect privacy a priori. We discuss risk mitigation through policies and standards such as not allowing anonymous pings of genetic beacons and requiring minimum beacon sizes.

Big biomedical data as the key resource for discovery science

Modern biomedical data collection is generating exponentially more data in a multitude of formats. This flood of complex data poses significant opportunities to discover and understand the critical interplay among such diverse domains as genomics, proteomics, metabolomics, and phenomics, including imaging, biometrics, and clinical data. The Big Data for Discovery Science Center is taking an "-ome to home" approach to discover linkages between these disparate data sources by mining existing databases of proteomic and genomic data, brain images, and clinical assessments. In support of this work, the authors developed new technological capabilities that make it easy for researchers to manage, aggregate, manipulate, integrate, and model large amounts of distributed data. Guided by biological domain expertise, the Center's computational resources and software will reveal relationships and patterns, aiding researchers in identifying biomarkers for the most confounding conditions and diseases, such as Parkinson's and Alzheimer's.

Whole-genome sequencing suggests a chemokine gene cluster that modifies age at onset in familial Alzheimer's disease

We have sequenced the complete genomes of 72 individuals affected with early-onset familial Alzheimer's disease caused by an autosomal dominant, highly penetrant mutation in the presenilin-1 (PSEN1) gene, and performed genome-wide association testing to identify variants that modify age at onset (AAO) of Alzheimer's disease. Our analysis identified a haplotype of single-nucleotide polymorphisms (SNPs) on chromosome 17 within a chemokine gene cluster associated with delayed onset of mild-cognitive impairment and dementia. Individuals carrying this haplotype had a mean AAO of mild-cognitive impairment at 51.0 ± 5.2 years compared with 41.1 ± 7.4 years for those without these SNPs. This haplotype thus appears to modify Alzheimer's AAO, conferring a large (~10 years) protective effect. The associated locus harbors several chemokines including eotaxin-1 encoded by CCL11, and the haplotype includes a missense polymorphism in this gene. Validating this association, we found plasma eotaxin-1 levels were correlated with disease AAO in an independent cohort from the University of California San Francisco Memory and Aging Center. In this second cohort, the associated haplotype disrupted the typical age-associated increase of eotaxin-1 levels, suggesting a complex regulatory role for this haplotype in the general population. Altogether, these results suggest eotaxin-1 as a novel modifier of Alzheimer's disease AAO and open potential avenues for therapy.

Heterozygous Loss-of-Function Mutations in DLL4 Cause Adams-Oliver Syndrome

Adams-Oliver syndrome (AOS) is a rare developmental disorder characterized by the presence of aplasia cutis congenita (ACC) of the scalp vertex and terminal limb-reduction defects. Cardiovascular anomalies are also frequently observed. Mutations in five genes have been identified as a cause for AOS prior to this report. Mutations in EOGT and DOCK6 cause autosomal-recessive AOS, whereas mutations in ARHGAP31, RBPJ, and NOTCH1 lead to autosomal-dominant AOS. Because RBPJ, NOTCH1, and EOGT are involved in NOTCH signaling, we hypothesized that mutations in other genes involved in this pathway might also be implicated in AOS pathogenesis. Using a candidate-gene-based approach, we prioritized DLL4, a critical NOTCH ligand, due to its essential role in vascular development in the context of cardiovascular features in AOS-affected individuals. Targeted resequencing of the DLL4 gene with a custom enrichment panel in 89 independent families resulted in the identification of seven mutations. A defect in DLL4 was also detected in two families via whole-exome or genome sequencing. In total, nine heterozygous mutations in DLL4 were identified, including two nonsense and seven missense variants, the latter encompassing four mutations that replace or create cysteine residues, which are most likely critical for maintaining structural integrity of the protein. Affected individuals with DLL4 mutations present with variable clinical expression with no emerging genotype-phenotype correlations. Our findings demonstrate that DLL4 mutations are an additional cause of autosomal-dominant AOS or isolated ACC and provide further evidence for a key role of NOTCH signaling in the etiology of this disorder.

Sequence-Level Analysis of the Major European Huntington Disease Haplotype

Huntington disease (HD) reflects the dominant consequences of a CAG-repeat expansion in HTT. Analysis of common SNP-based haplotypes has revealed that most European HD subjects have distinguishable HTT haplotypes on their normal and disease chromosomes and that ∼50% of the latter share the same major HD haplotype. We reasoned that sequence-level investigation of this founder haplotype could provide significant insights into the history of HD and valuable information for gene-targeting approaches. Consequently, we performed whole-genome sequencing of HD and control subjects from four independent families in whom the major European HD haplotype segregates with the disease. Analysis of the full-sequence-based HTT haplotype indicated that these four families share a common ancestor sufficiently distant to have permitted the accumulation of family-specific variants. Confirmation of new CAG-expansion mutations on this haplotype suggests that unlike most founders of human disease, the common ancestor of HD-affected families with the major haplotype most likely did not have HD. Further, availability of the full sequence data validated the use of SNP imputation to predict the optimal variants for capturing heterozygosity in personalized allele-specific gene-silencing approaches. As few as ten SNPs are capable of revealing heterozygosity in more than 97% of European HD subjects. Extension of allele-specific silencing strategies to the few remaining homozygous individuals is likely to be achievable through additional known SNPs and discovery of private variants by complete sequencing of HTT. These data suggest that the current development of gene-based targeting for HD could be extended to personalized allele-specific approaches in essentially all HD individuals of European ancestry.

State of the Human Proteome in 2014/2015 As Viewed through PeptideAtlas: Enhancing Accuracy and Coverage through the AtlasProphet

The Human PeptideAtlas is a compendium of the highest quality peptide identifications from over 1000 shotgun mass spectrometry proteomics experiments collected from many different laboratories, all reanalyzed through a uniform processing pipeline. The latest 2015-03 build contains substantially more input data than past releases, is mapped to a recent version of our merged reference proteome, and uses improved informatics processing and the development of the AtlasProphet to provide the highest quality results. Within the set of ∼20,000 neXtProt primary entries, 14,070 (70%) are confidently detected in the latest build, 5% are ambiguous, 9% are redundant, leaving the total percentage of proteins for which there are no mapping detections at just 16% (3166), all derived from over 133 million peptide-spectrum matches identifying more than 1 million distinct peptides using AtlasProphet to characterize and classify the protein matches. Improved handling for detection and presentation of single amino-acid variants (SAAVs) reveals the detection of 5326 uniquely mapping SAAVs across 2794 proteins. With such a large amount of data, the control of false positives is a challenge. We present the methodology and results for maintaining rigorous quality along with a discussion of the implications of the remaining sources of errors in the build.

An Efficient Search Algorithm for Finding Genomic-Range Overlaps Based on the Maximum Range Length

Efficient search algorithms for finding genomic-range overlaps are essential for various bioinformatics applications. A majority of fast algorithms for searching the overlaps between a query range (e.g., a genomic variant) and a set of N reference ranges (e.g., exons) has time complexity of O(k + logN), where kdenotes a term related to the length and location of the reference ranges. Here, we present a simple but efficient algorithm that reduces k, based on the maximum reference range length. Specifically, for a given query range and the maximum reference range length, the proposed method divides the reference range set into three subsets: always, potentially, and never overlapping. Therefore, search effort can be reduced by excluding never overlapping subset. We demonstrate that the running time of the proposed algorithm is proportional to potentially overlapping subset size, that is proportional to the maximum reference range length if all the other conditions are the same. Moreover, an implementation of our algorithm was 13.8 to 30.0 percent faster than one of the fastest range search methods available when tested on various genomic-range data sets. The proposed algorithm has been incorporated into a disease-linked variant prioritization pipeline for WGS (http://gnome.tchlab.org) and its implementation is available at http://ml.ssu.ac.kr/gSearch.

Alternating Hemiplegia of Childhood: Retrospective Genetic Study and Genotype-Phenotype Correlations in 187 Subjects from the US AHCF Registry

Mutations in ATP1A3 cause Alternating Hemiplegia of Childhood (AHC) by disrupting function of the neuronal Na+/K+ ATPase. Published studies to date indicate 2 recurrent mutations, D801N and E815K, and a more severe phenotype in the E815K cohort. We performed mutation analysis and retrospective genotype-phenotype correlations in all eligible patients with AHC enrolled in the US AHC Foundation registry from 1997-2012. Clinical data were abstracted from standardized caregivers’ questionnaires and medical records and confirmed by expert clinicians. We identified ATP1A3 mutations by Sanger and whole genome sequencing, and compared phenotypes within and between 4 groups of subjects, those with D801N, E815K, other ATP1A3 or no ATP1A3 mutations. We identified heterozygous ATP1A3 mutations in 154 of 187 (82%) AHC patients. Of 34 unique mutations, 31 (91%) are missense, and 16 (47%) had not been previously reported. Concordant with prior studies, more than 2/3 of all mutations are clustered in exons 17 and 18. Of 143 simplex occurrences, 58 had D801N (40%), 38 had E815K (26%) and 11 had G937R (8%) mutations. Patients with an E815K mutation demonstrate an earlier age of onset, more severe motor impairment and a higher prevalence of status epilepticus. This study further expands the number and spectrum of ATP1A3 mutations associated with AHC and confirms a more deleterious effect of the E815K mutation on selected neurologic outcomes. However, the complexity of the disorder and the extensive phenotypic variability among subgroups merits caution and emphasizes the need for further studies.

Leveraging ancestry to improve causal variant identification in exome sequencing for monogenic disorders

Recent breakthroughs in exome-sequencing technology have made possible the identification of many causal variants of monogenic disorders. Although extremely powerful when closely related individuals (eg, child and parents) are simultaneously sequenced, sequencing of a single case is often unsuccessful due to the large number of variants that need to be followed up for functional validation. Many approaches filter out common variants above a given frequency threshold (eg, 1%), and then prioritize the remaining variants according to their functional, structural and conservation properties. Here we present methods that leverage the genetic structure across different populations to improve filtering performance while accounting for the finite sample size of the reference panels. We show that leveraging genetic structure reduces the number of variants that need to be followed up by 16% in simulations and by up to 38% in empirical data of 20 exomes from individuals with monogenic disorders for which the causal variants are known.

Rare variants in neuronal excitability genes influence risk for bipolar disorder

We sequenced the genomes of 200 individuals from 41 families multiply affected with bipolar disorder (BD) to identify contributions of rare variants to genetic risk. We initially focused on 3,087 candidate genes with known synaptic functions or prior evidence from genome-wide association studies. BD pedigrees had an increased burden of rare variants in genes encoding neuronal ion channels, including subunits of GABAA receptors and voltage-gated calcium channels. Four uncommon coding and regulatory variants also showed significant association, including a missense variant in GABRA6. Targeted sequencing of 26 of these candidate genes in an additional 3,014 cases and 1,717 controls confirmed rare variant associations in ANK3, CACNA1B, CACNA1C, CACNA1D, CACNG2, CAMK2A, and NGF. Variants in promoters and 5' and 3' UTRs contributed more strongly than coding variants to risk for BD, both in pedigrees and in the case-control cohort. The genes and pathways identified in this study regulate diverse aspects of neuronal excitability. We conclude that rare variants in neuronal excitability genes contribute to risk for BD.

Hoyeraal-Hreidarsson syndrome caused by a germline mutation in the TEL patch of the telomere protein TPP1

Germline mutations in telomere biology genes cause dyskeratosis congenita (DC), an inherited bone marrow failure and cancer predisposition syndrome. Hoyeraal-Hreidarsson syndrome (HH) is a clinically severe variant of DC. Using exome sequencing, Kocak et al. identified mutations in ACD (encoding TPP1), a component of the telomeric shelterin complex, in one family affected by HH. Characterization of the mutations revealed that the single-amino-acid deletion affecting the TEL patch surface of the TPP1 protein significantly compromises both telomerase recruitment and processivity.

Germline mutations in telomere biology genes cause dyskeratosis congenita (DC), an inherited bone marrow failure and cancer predisposition syndrome. DC is a clinically heterogeneous disorder diagnosed by the triad of dysplastic nails, abnormal skin pigmentation, and oral leukoplakia; Hoyeraal-Hreidarsson syndrome (HH), a clinically severe variant of DC, also includes cerebellar hypoplasia, immunodeficiency, and intrauterine growth retardation. Approximately 70% of DC cases are associated with a germline mutation in one of nine genes, the products of which are all involved in telomere biology. Using exome sequencing, we identified mutations in Adrenocortical Dysplasia Homolog (ACD) (encoding TPP1), a component of the telomeric shelterin complex, in one family affected by HH. The proband inherited a deletion from his father and a missense mutation from his mother, resulting in extremely short telomeres and a severe clinical phenotype. Characterization of the mutations revealed that the single-amino-acid deletion affecting the TEL patch surface of the TPP1 protein significantly compromises both telomerase recruitment and processivity, while the missense mutation in the TIN2-binding region of TPP1 is not as clearly deleterious to TPP1 function. Our results emphasize the critical roles of the TEL patch in proper stem cell function and demonstrate that TPP1 is the second shelterin component (in addition to TIN2) to be implicated in DC.

Mutations in NOTCH1 cause Adams-Oliver syndrome

Notch signaling determines and reinforces cell fate in bilaterally symmetric multicellular eukaryotes. Despite the involvement of Notch in many key developmental systems, human mutations in Notch signaling components have mainly been described in disorders with vascular and bone effects. Here, we report five heterozygous NOTCH1 variants in unrelated individuals with Adams-Oliver syndrome (AOS), a rare disease with major features of aplasia cutis of the scalp and terminal transverse limb defects. Using whole-genome sequencing in a cohort of 11 families lacking mutations in the four genes with known roles in AOS pathology (ARHGAP31, RBPJ, DOCK6, and EOGT), we found a heterozygous de novo 85 kb deletion spanning the NOTCH1 5' region and three coding variants (c.1285T>C [p.Cys429Arg], c.4487G>A [p.Cys1496Tyr], and c.5965G>A [p.Asp1989Asn]), two of which are de novo, in four unrelated probands. In a fifth family, we identified a heterozygous canonical splice-site variant (c.743-1 G>T) in an affected father and daughter. These variants were not present in 5,077 in-house control genomes or in public databases. In keeping with the prominent developmental role described for Notch1 in mouse vasculature, we observed cardiac and multiple vascular defects in four of the five families. We propose that the limb and scalp defects might also be due to a vasculopathy in NOTCH1-related AOS. Our results suggest that mutations in NOTCH1 are the most common cause of AOS and add to a growing list of human diseases that have a vascular and/or bony component and are caused by alterations in the Notch signaling pathway.

Diffuse angiopathy in Adams-Oliver syndrome associated with truncating DOCK6 mutations

Adams-Oliver syndrome (AOS) is a rare malformation syndrome characterized by the presence of two anomalies: aplasia cutis congenita of the scalp and transverse terminal limb defects. Many affected individuals also have additional malformations, including a variety of intracranial anomalies such as periventricular calcification in keeping with cerebrovascular microbleeds, impaired neuronal migration, epilepsy, and microcephaly. Cardiac malformations can be present, as can vascular dysfunction in the forms of cutis marmorata telangiectasia congenita, pulmonary vein stenoses, and abnormal hepatic microvasculature. Elucidated genetic causes include four genes in different pathways, leading to a model of AOS as a multi-pathway disorder. We identified an infant with mild aplasia cutis congenita and terminal transverse limb defects, developmental delay and a severe, diffuse angiopathy with incomplete microvascularization. Whole-genome sequencing documented two rare truncating variants in DOCK6, a gene associated with a type of autosomal recessive AOS that recurrently features periventricular calcification and impaired neurodevelopment. We highlight an unexpectedly high frequency of likely deleterious mutations in this gene in the general population, relative to the rarity of the disease, and discuss possible explanations for this discrepancy.

Whole-exome sequencing and functional studies identify RPS29 as a novel gene mutated in multicase Diamond-Blackfan anemia families

Diamond-Blackfan anemia (DBA) is a cancer-prone inherited bone marrow failure syndrome. Approximately half of DBA patients have a germ-line mutation in a ribosomal protein gene. We used whole-exome sequencing to identify disease-causing genes in 2 large DBA families. After filtering, 1 nonsynonymous mutation (p.I31F) in the ribosomal protein S29 (RPS29[AUQ1]) gene was present in all 5 DBA-affected individuals and the obligate carrier, and absent from the unaffected noncarrier parent in 1 DBA family. A second DBA family was found to have a different nonsynonymous mutation (p.I50T) in RPS29. Both mutations are amino acid substitutions in exon 2 predicted to be deleterious and resulted in haploinsufficiency of RPS29 expression compared with wild-type RPS29 expression from an unaffected control. The DBA proband with the p.I31F RPS29 mutation had a pre-ribosomal RNA (rRNA) processing defect compared with the healthy control. We demonstrated that both RPS29 mutations failed to rescue the defective erythropoiesis in the rps29(-/-) mutant zebra fish DBA model. RPS29 is a component of the small 40S ribosomal subunit and essential for rRNA processing and ribosome biogenesis. We uncovered a novel DBA causative gene, RPS29, and showed that germ-line mutations in RPS29 can cause a defective erythropoiesis phenotype using a zebra fish model.

A recessive founder mutation in regulator of telomere elongation helicase 1, RTEL1, underlies severe immunodeficiency and features of Hoyeraal Hreidarsson syndrome

Dyskeratosis congenita (DC) is a heterogeneous inherited bone marrow failure and cancer predisposition syndrome in which germline mutations in telomere biology genes account for approximately one-half of known families. Hoyeraal Hreidarsson syndrome (HH) is a clinically severe variant of DC in which patients also have cerebellar hypoplasia and may present with severe immunodeficiency and enteropathy. We discovered a germline autosomal recessive mutation in RTEL1, a helicase with critical telomeric functions, in two unrelated families of Ashkenazi Jewish (AJ) ancestry. The affected individuals in these families are homozygous for the same mutation, R1264H, which affects three isoforms of RTEL1. Each parent was a heterozygous carrier of one mutant allele. Patient-derived cell lines revealed evidence of telomere dysfunction, including significantly decreased telomere length, telomere length heterogeneity, and the presence of extra-chromosomal circular telomeric DNA. In addition, RTEL1 mutant cells exhibited enhanced sensitivity to the interstrand cross-linking agent mitomycin C. The molecular data and the patterns of inheritance are consistent with a hypomorphic mutation in RTEL1 as the underlying basis of the clinical and cellular phenotypes. This study further implicates RTEL1 in the etiology of DC/HH and immunodeficiency, and identifies the first known homozygous autosomal recessive disease-associated mutation in RTEL1.

Germline mutations of regulator of telomere elongation helicase 1, RTEL1, in Dyskeratosis congenita

Dyskeratosis congenita (DC) is an inherited bone marrow failure and cancer predisposition syndrome caused by aberrant telomere biology. The classic triad of dysplastic nails, abnormal skin pigmentation, and oral leukoplakia is diagnostic of DC, but substantial clinical heterogeneity exists; the clinically severe variant Hoyeraal Hreidarsson syndrome (HH) also includes cerebellar hypoplasia, severe immunodeficiency, enteropathy, and intrauterine growth retardation. Germline mutations in telomere biology genes account for approximately one-half of known DC families. Using exome sequencing, we identified mutations in RTEL1, a helicase with critical telomeric functions, in two families with HH. In the first family, two siblings with HH and very short telomeres inherited a premature stop codon from their mother who has short telomeres. The proband from the second family has HH and inherited a premature stop codon in RTEL1 from his father and a missense mutation from his mother, who also has short telomeres. In addition, inheritance of only the missense mutation led to very short telomeres in the proband's brother. Targeted sequencing identified a different RTEL1 missense mutation in one additional DC proband who has bone marrow failure and short telomeres. Both missense mutations affect the helicase domain of RTEL1, and three in silico prediction algorithms suggest that they are likely deleterious. The nonsense mutations both cause truncation of the RTEL1 protein, resulting in loss of the PIP box; this may abrogate an important protein-protein interaction. These findings implicate a new telomere biology gene, RTEL1, in the etiology of DC.

Reduced sensitivity of the ferroportin Q248H mutant to physiological concentrations of hepcidin

Ferroportin Q248H mutation has an allele frequency of 2.2-13.4% in African populations and is associated with a mild tendency to increased serum ferritin in the general population. Some investigators have reported that ferroportin Q248H is degraded after exposure to hepcidin in exactly the same manner as wild-type ferroportin, but supraphysiological concentrations of hepcidin were used. The aim of our study was to determine whether ferroportin Q248H may have reduced sensitivity to physiological concentrations of hepcidin. The sensitivity of ferroportin Q248H to hepcidin was determined in 293T cells transiently expressing ferroportin using immunoblotting and fluorescence analysis. Ferritin concentrations were measured in these cells and also in human primary monocytes derived from humans with different ferroportin genotypes. The effect of Q248H on serum iron measures was examined in patients with sickle cell anemia. Immunoblotting and fluorescence analysis showed decreased sensitivity of ferroportin Q248H to physiological concentrations of hepcidin. Lower ferritin concentrations were observed after incubation with iron and hepcidin in 293T cells expressing ferroportin Q248H and in primary monocytes from ferroportin Q248H subjects. In sickle cell anemia, ferroportin Q248H heterozygotes had lower serum ferritin concentrations than wild-type subjects, consistent with enhanced iron release by macrophage ferroportin Q248H. A clinical benefit of ferroportin Q248H was suggested by lower echocardiographic estimates of pulmonary artery pressure in patients carrying mutant alleles. In conclusion, our results suggest that ferroportin Q248H protein is resistant to physiological concentrations of hepcidin and that this mutation has discernible effects on iron metabolism-related clinical complications of sickle cell anemia. They provide a mechanistic explanation for the effect of ferroportin Q248H on iron status in individuals of African descent and suggest that these changes in iron metabolism may be beneficial under certain disease-specific circumstances. (ClinicalTrials.gov Identifier:NCT00011648).

gSearch: a fast and flexible general search tool for whole-genome sequencing

BACKGROUND

Various processes such as annotation and filtering of variants or comparison of variants in different genomes are required in whole-genome or exome analysis pipelines. However, processing different databases and searching among millions of genomic loci is not trivial.

RESULTS

gSearch compares sequence variants in the Genome Variation Format (GVF) or Variant Call Format (VCF) with a pre-compiled annotation or with variants in other genomes. Its search algorithms are subsequently optimized and implemented in a multi-threaded manner. The proposed method is not a stand-alone annotation tool with its own reference databases. Rather, it is a search utility that readily accepts public or user-prepared reference files in various formats including GVF, Generic Feature Format version 3 (GFF3), Gene Transfer Format (GTF), VCF and Browser Extensible Data (BED) format. Compared to existing tools such as ANNOVAR, gSearch runs more than 10 times faster. For example, it is capable of annotating 52.8 million variants with allele frequencies in 6 min.

AVAILABILITY

gSearch is available at http://ml.ssu.ac.kr/gSearch. It can be used as an independent search tool or can easily be integrated to existing pipelines through various programming environments such as Perl, Ruby and Python.

Novel rare variants in congenital cardiac arrhythmia genes are frequent in drug-induced torsades de pointes

Marked prolongation of the QT interval and polymorphic ventricular tachycardia following medication (drug-induced long QT syndrome, diLQTS) is a severe adverse drug reaction (ADR) that phenocopies congenital long QT syndrome (cLQTS) and is one of the leading causes for drug withdrawal and relabeling. We evaluated the frequency of rare non-synonymous variants in genes contributing to the maintenance of heart rhythm in cases of diLQTS using targeted capture coupled to next-generation sequencing. Eleven of 31 diLQTS subjects (36%) carried a novel missense mutation in genes with known congenital arrhythmia associations or with a known cLQTS mutation. In the 26 Caucasian subjects, 23% carried a highly conserved rare variant predicted to be deleterious to protein function in these genes compared with only 2-4% in public databases (P<0.003). We conclude that the rare variation in genes responsible for congenital arrhythmia syndromes is frequent in diLQTS. Our findings demonstrate that diLQTS is a pharmacogenomic syndrome predisposed by rare genetic variants.