Dupuytren´s disease, a fibromatosis of the connective tissue in the palm, is a common complex disease with a strong genetic component. Up to date nine genetic loci have been found to be associated with the disease. Six of these loci contain genes that code for Wnt signalling proteins. In spite of this striking first insight into the genetic factors in Dupuytren´s disease, much of the inherited risk in Dupuytren´s disease still needs to be discovered. The already identified loci jointly explain ~1% of the heritability in this disease. To further elucidate the genetic basis of Dupuytren´s disease, we performed a genome-wide meta-analysis combining three genome-wide association study (GWAS) data sets, comprising 1,580 cases and 4,480 controls. We corroborated all nine previously identified loci, six of these with genome-wide significance (p-value < 5x10-8). In addition, we identified 14 new suggestive loci (p-value < 10-5). Intriguingly, several of these new loci contain genes associated with Wnt signalling and therefore represent excellent candidates for replication. Next, we compared whole-transcriptome data between patient- and control-derived tissue samples and found the Wnt/β-catenin pathway to be the top deregulated pathway in patient samples. We then conducted network and pathway analyses in order to identify protein networks that are enriched for genes highlighted in the GWAS meta-analysis and expression data sets. We found further evidence that the Wnt signalling pathways in conjunction with other pathways may play a critical role in Dupuytren´s disease.

So far, the annotation of translation initiation sites (TISs) has been based mostly upon bioinformatics rather than experimental evidence. We adapted ribosomal footprinting to puromycin-treated cells to generate a transcriptome-wide map of TISs in a human monocytic cell line. A neural network was trained on the ribosomal footprints observed at previously annotated AUG translation initiation codons (TICs), and used for the ab initio prediction of TISs in 5062 transcripts with sufficient sequence coverage. Functional interpretation suggested 2994 novel upstream open reading frames (uORFs) in the 5' UTR, 1406 uORFs overlapping with the coding sequence, and 546 N-terminal protein extensions. The TIS detection method was validated on the basis of previously published alternative TISs and uORFs. Among primates, TICs in newly annotated TISs were significantly more conserved than control codons, both for AUGs and near-cognate codons. The transcriptome-wide map of novel candidate TISs derived as part of the study will shed further light on the way in which human proteome diversity is influenced by alternative translation initiation and regulation.

Studies on the association of fatty acid-binding protein 2 (FABP2) A54T and promoter polymorphism, and type 2 diabetes mellitus, insulin, and triglyceride levels are controversial. The aim of this study was to investigate the interfering effect of FABP2 A54T and promoter polymorphism on the postprandial response to a mixed meal and an oral glucose load. Seven hundred men from the Metabolic Intervention Cohort Kiel underwent a standard glucose tolerance test and a standardized mixed meal test and were genotyped by use of the Taqman method. When calculated independently from promoter variability, postprandial triglyceride levels were significantly higher and postprandial insulin sensitivity (homeostasis model assessment index) was lower in homozygous carriers of FABP2 T54T compared with carriers of the FABP2 exon wild-type allele (FABP2 A54A and A54T). This confirms previous findings. The effect of the exon T54T genotype on triglyceride levels and insulin sensitivity, however, was dependent on promoter variability. We found a significant increase in postprandial triglyceride levels and a decrease in insulin sensitivity due to T54T only in the presence of the homozygous B genotype at the promoter polymorphism. Similar results were obtained after oral glucose tolerance test. Reporter gene assays indicated a higher responsiveness to peroxisome proliferator-activating receptor-gamma (PPAR-gamma)/retinoid X receptor (RXR) of FABP2 promoter B vs promoter A. Synergism between a higher inducibility of FABP2 expression and a higher activity of T54 variant may explain higher postprandial triglycerides in case of combined genotype (promoter B + T54). This interference and different linkage disequilibrium between FABP2 exon and promoter polymorphisms may explain the different results obtained in different cohorts.

To analyze family-based whole-genome sequence (WGS) data for complex traits, we developed a rare variant (RV) non-parametric linkage (NPL) analysis method, which has advantages over association methods. The RV-NPL differs from the NPL in that RVs are analyzed, and allele sharing among affected relative-pairs is estimated only for minor alleles. Analyzing families can increase power because causal variants with familial aggregation usually have larger effect sizes than those underlying sporadic diseases. Differing from association analysis, for NPL only affected individuals are analyzed, which can increase power, since unaffected family members can be susceptibility variant carriers. RV-NPL is robust to population substructure and admixture, inclusion of nonpathogenic variants, as well as allelic and locus heterogeneity and can readily be applied outside of coding regions. In contrast to analyzing common variants using NPL, where loci localize to large genomic regions (e.g., >50 Mb), mapped regions are well defined for RV-NPL. Using simulation studies, we demonstrate that RV-NPL is substantially more powerful than applying traditional NPL methods to analyze RVs. The RV-NPL was applied to analyze 107 late-onset Alzheimer disease (LOAD) pedigrees of Caribbean Hispanic and European ancestry with WGS data, and statistically significant linkage (LOD ≥ 3.8) was found with RVs in PSMF1 and PTPN21 which have been shown to be involved in LOAD etiology. Additionally, nominally significant linkage was observed with RVs in ABCA7, ACE, EPHA1, and SORL1, genes that were previously reported to be associated with LOAD. RV-NPL is an ideal method to elucidate the genetic etiology of complex familial diseases.

Clinically identified genetic variants in ion channels can be benign or cause disease by increasing or decreasing the protein function. As a consequence, therapeutic decision-making is challenging without molecular testing of each variant. Our biophysical knowledge of ion-channel structures and function is just emerging, and it is currently not well understood which amino acid residues cause disease when mutated. We sought to systematically identify biological properties associated with variant pathogenicity across all major voltage and ligand-gated ion-channel families. We collected and curated 3049 pathogenic variants from hundreds of neurodevelopmental and other disorders and 12 546 population variants for 30 ion channel or channel subunits for which a high-quality protein structure was available. Using a wide range of bioinformatics approaches, we computed 163 structural features and tested them for pathogenic variant enrichment. We developed a novel 3D spatial distance scoring approach that enables comparisons of pathogenic and population variant distribution across protein structures. We discovered and independently replicated that several pore residue properties and proximity to the pore axis were most significantly enriched for pathogenic variants compared to population variants. Using our 3D scoring approach, we showed that the strongest pathogenic variant enrichment was observed for pore-lining residues and alpha-helix residues within 5Å distance from the pore axis centre and not involved in gating. Within the subset of residues located at the pore, the hydrophobicity of the pore was the feature most strongly associated with variant pathogenicity. We also found an association between the identified properties and both clinical phenotypes and functional in vitro assays for voltage-gated sodium channels (SCN1A, SCN2A, SCN8A) and N-methyl-D-aspartate receptor (GRIN1, GRIN2A, GRIN2B) encoding genes. In an independent expert-curated dataset of 1422 neurodevelopmental disorder pathogenic patient variants and 679 electrophysiological experiments, we show that pore axis distance is associated with seizure age of onset and cognitive performance as well as differential gain versus loss-of-channel function. In summary, we identified biological properties associated with ion-channel malfunction and show that these are correlated with in vitro functional readouts and clinical phenotypes in patients with neurodevelopmental disorders. Our results suggest that clinical decision support algorithms that predict variant pathogenicity and function are feasible in the future.

Analyses of few gene-sets in epilepsy showed a potential to unravel key disease associations. We set out to investigate the burden of ultra-rare variants (URVs) in a comprehensive range of biologically informed gene-sets presumed to be implicated in epileptogenesis.

Pathogenic copy-number variants (CNVs) can cause a heterogeneous spectrum of rare and severe disorders. However, most CNVs are benign and are part of natural variation in human genomes. CNV pathogenicity classification, genotype-phenotype analyses, and therapeutic target identification are challenging and time-consuming tasks that require the integration and analysis of information from multiple scattered sources by experts.

High frequency electrosurgery has a key role in the broadening application of liver surgery. Its molecular signature, i.e. the metabolites evolving from electrocauterization which may inhibit hepatic wound healing, have not been systematically studied.

The analysis of structural variants, in particular of copy-number variations (CNVs), has proven valuable in unraveling the genetic basis of human diseases. Hence, a large number of algorithms have been developed for the detection of CNVs in SNP array signal intensity data. Using the European and African HapMap trio data, we undertook a comparative evaluation of six commonly used CNV detection software tools, namely Affymetrix Power Tools (APT), QuantiSNP, PennCNV, GLAD, R-gada and VEGA, and assessed their level of pair-wise prediction concordance. The tool-specific CNV prediction accuracy was assessed in silico by way of intra-familial validation. Software tools differed greatly in terms of the number and length of the CNVs predicted as well as the number of markers included in a CNV. All software tools predicted substantially more deletions than duplications. Intra-familial validation revealed consistently low levels of prediction accuracy as measured by the proportion of validated CNVs (34-60%). Moreover, up to 20% of apparent family-based validations were found to be due to chance alone. Software using Hidden Markov models (HMM) showed a trend to predict fewer CNVs than segmentation-based algorithms albeit with greater validity. PennCNV yielded the highest prediction accuracy (60.9%). Finally, the pairwise concordance of CNV prediction was found to vary widely with the software tools involved. We recommend HMM-based software, in particular PennCNV, rather than segmentation-based algorithms when validity is the primary concern of CNV detection. QuantiSNP may be used as an additional tool to detect sets of CNVs not detectable by the other tools. Our study also reemphasizes the need for laboratory-based validation, such as qPCR, of CNVs predicted in silico.

Leprosy, a chronic infectious disease caused by Mycobacterium leprae (M. leprae), was very common in Europe till the 16th century. Here, we perform an ancient DNA study on medieval skeletons from Denmark that show lesions specific for lepromatous leprosy (LL). First, we test the remains for M. leprae DNA to confirm the infection status of the individuals and to assess the bacterial diversity. We assemble 10 complete M. leprae genomes that all differ from each other. Second, we evaluate whether the human leukocyte antigen allele DRB1*15:01, a strong LL susceptibility factor in modern populations, also predisposed medieval Europeans to the disease. The comparison of genotype data from 69 M. leprae DNA-positive LL cases with those from contemporary and medieval controls reveals a statistically significant association in both instances. In addition, we observe that DRB1*15:01 co-occurs with DQB1*06:02 on a haplotype that is a strong risk factor for inflammatory diseases today.

The definition of haplotype blocks of single-nucleotide polymorphisms (SNPs) has been proposed so that the haplotypes can be used as markers in association studies and to efficiently describe human genetic variation. The International Haplotype Map (HapMap) project to construct a comprehensive catalog of haplotypic variation in humans is underway. However, a number of factors have already been shown to influence the definition of blocks, including the population studied and the sample SNP density. Here, we examine the effect that marker selection has on the definition of blocks and the pattern of haplotypes by using comparable but complementary SNP sets and a number of block definition methods in various genomic regions and populations that were provided by the Encyclopedia of DNA Elements (ENCODE) project. We find that the chosen SNP set has a profound effect on the block-covered sequence and block borders, even at high marker densities. Our results question the very concept of discrete haplotype blocks and the possibility of generalizing block findings from the HapMap project. We comparatively apply the block-free tagging-SNP approach and discuss both the haplotype approach and the tagging-SNP approach as means to efficiently catalog genetic variation.

Polymorphisms in the interleukin-1 (IL1) gene have been suggested to influence transcription of IL1A (interleukin-1alpha) and IL1B (interleukin-1beta) and thereby the pathophysiology of periodontitis. This case-control association study on 415 northern European Caucasian patients with aggressive periodontitis (AgP) and 874 healthy controls was conducted to examine 10 single-nucleotide polymorphisms (SNPs) in the genes of the IL1 cluster for association with IL1A, IL1B, CKAP2L (cytoskeleton-associated protein 2-like), and IL1RN (IL-1 receptor antagonist). The results do not support an association between variants in the IL1 gene cluster and AgP. This case-control study had at least 95% power to detect genuine associations with variants carrying relative risks of at least 1.5 for heterozygous carriers and 2.25 for homozygous carriers. Previous reports of an association between IL1 promoter SNPs and periodontitis might reflect subpopulation effects and have to be interpreted with care.

Whole-genome and whole-exome sequencing of individual patients allow the study of rare and potentially causative genetic variation. In this study, we sequenced DNA of a trio comprising a boy with very-early-onset inflammatory bowel disease (veoIBD) and his unaffected parents. We identified a rare, X-linked missense variant in the NAPDH oxidase NOX1 gene (c.C721T, p.R241C) in heterozygous state in the mother and in hemizygous state in the patient. We discovered that, in addition, the patient was homozygous for a common missense variant in the CYBA gene (c.T214C, p.Y72H). CYBA encodes the p22phox protein, a cofactor for NOX1. Functional assays revealed reduced cellular ROS generation and antibacterial capacity of NOX1 and p22phox variants in intestinal epithelial cells. Moreover, the identified NADPH oxidase complex variants affected NOD2-mediated immune responses, and p22phox was identified as a novel NOD2 interactor. In conclusion, we detected missense variants in a veoIBD patient that disrupt the host response to bacterial challenges and reduce protective innate immune signaling via NOD2. We assume that the patient's individual genetic makeup favored disturbed intestinal mucosal barrier function.

Guided by the practice of classical epidemiology, research into the genetic basis of complex disease has usually taken for granted the dictum that causative mutations are invariably over-represented among clinically affected as compared to unaffected individuals. However, we show that this supposition is not true and that a mutation contributing to the etiology of a complex disease can, under certain circumstances, be depleted among patients. Populations with defined disease prevalence were repeatedly simulated under a Wright-Fisher model, assuming various types of population history and genotype-phenotype relationship. For each simulation, the resulting mutation-specific population frequencies and odds ratios (ORs) were evaluated. In addition, the relationship between mutation frequency and OR was studied using real data from the NIH GWAS catalogue of reported phenotype associations of single-nucleotide polymorphisms (SNPs). While rare diseases (prevalence <1%) were found to be consistently caused by rare mutations with ORs>1, up to 20% of mutations causing a pandemic disease (prevalence 10-20%) had ORs<1, and their population frequency ranged from 0% to 100%. Moreover, simulation-based ORs exhibited a wide distribution, irrespective of mutation frequency. In conclusion, a substantial proportion of mutations causing common complex diseases may appear 'protective' in genetic epidemiological studies and hence would normally tend to be excluded, albeit erroneously, from further study. This apparently paradoxical result is explicable in terms of mutual confounding of the respective genotype-phenotype relationships due to a negative correlation between causal mutations induced by their common gene genealogy. As would be predicted by our findings, a significant negative correlation became apparent in published genome-wide association studies between the OR of genetic variants associated with a particular disease and the prevalence of that disease.

Epileptic encephalopathies are a phenotypically and genetically heterogeneous group of severe epilepsies accompanied by intellectual disability and other neurodevelopmental features. Using next-generation sequencing, we identified four different de novo mutations in KCNA2, encoding the potassium channel KV1.2, in six isolated patients with epileptic encephalopathy (one mutation recurred three times independently). Four individuals presented with febrile and multiple afebrile, often focal seizure types, multifocal epileptiform discharges strongly activated by sleep, mild to moderate intellectual disability, delayed speech development and sometimes ataxia. Functional studies of the two mutations associated with this phenotype showed almost complete loss of function with a dominant-negative effect. Two further individuals presented with a different and more severe epileptic encephalopathy phenotype. They carried mutations inducing a drastic gain-of-function effect leading to permanently open channels. These results establish KCNA2 as a new gene involved in human neurodevelopmental disorders through two different mechanisms, predicting either hyperexcitability or electrical silencing of KV1.2-expressing neurons.

Psoriasis (PS) and Crohn disease (CD) have been shown to be epidemiologically, pathologically, and therapeutically connected, but little is known about their shared genetic causes. We performed meta-analyses of five published genome-wide association studies on PS (2,529 cases and 4,955 controls) and CD (2,142 cases and 5,505 controls), followed up 20 loci that showed strongest evidence for shared disease association and, furthermore, tested cross-disease associations for previously reported PS and CD risk alleles in additional 6,115 PS cases, 4,073 CD cases, and 10,100 controls. We identified seven susceptibility loci outside the human leukocyte antigen region (9p24 near JAK2, 10q22 at ZMIZ1, 11q13 near PRDX5, 16p13 near SOCS1, 17q21 at STAT3, 19p13 near FUT2, and 22q11 at YDJC) shared between PS and CD with genome-wide significance (p < 5 × 10(-8)) and confirmed four already established PS and CD risk loci (IL23R, IL12B, REL, and TYK2). Three of the shared loci are also genome-wide significantly associated with PS alone (10q22 at ZMIZ1, p(rs1250544) = 3.53 × 10(-8), 11q13 near PRDX5, p(rs694739) = 3.71 × 10(-09), 22q11 at YDJC, p(rs181359) = 8.02 × 10(-10)). In addition, we identified one susceptibility locus for CD (16p13 near SOCS1, p(rs4780355) = 4.99 × 10(-8)). Refinement of association signals identified shared genome-wide significant associations for exonic SNPs at 10q22 (ZMIZ1) and in silico expression quantitative trait locus analyses revealed that the associations at ZMIZ1 and near SOCS1 have a potential functional effect on gene expression. Our results show the usefulness of joint analyses of clinically distinct immune-mediated diseases and enlarge the map of shared genetic risk loci.

Numerous studies of human populations in Europe and Asia have revealed a concordance between their extant genetic structure and the prevailing regional pattern of geography and language. For native South Americans, however, such evidence has been lacking so far. Therefore, we examined the relationship between Y-chromosomal genotype on the one hand, and male geographic origin and linguistic affiliation on the other, in the largest study of South American natives to date in terms of sampled individuals and populations. A total of 1,011 individuals, representing 50 tribal populations from 81 settlements, were genotyped for up to 17 short tandem repeat (STR) markers and 16 single nucleotide polymorphisms (Y-SNPs), the latter resolving phylogenetic lineages Q and C. Virtually no structure became apparent for the extant Y-chromosomal genetic variation of South American males that could sensibly be related to their inter-tribal geographic and linguistic relationships. This continent-wide decoupling is consistent with a rapid peopling of the continent followed by long periods of isolation in small groups. Furthermore, for the first time, we identified a distinct geographical cluster of Y-SNP lineages C-M217 (C3*) in South America. Such haplotypes are virtually absent from North and Central America, but occur at high frequency in Asia. Together with the locally confined Y-STR autocorrelation observed in our study as a whole, the available data therefore suggest a late introduction of C3* into South America no more than 6,000 years ago, perhaps via coastal or trans-Pacific routes. Extensive simulations revealed that the observed lack of haplogroup C3* among extant North and Central American natives is only compatible with low levels of migration between the ancestor populations of C3* carriers and non-carriers. In summary, our data highlight the fact that a pronounced correlation between genetic and geographic/cultural structure can only be expected under very specific conditions, most of which are likely not to have been met by the ancestors of native South Americans.

Acrocallosal syndrome (ACLS) is an autosomal recessive neurodevelopmental disorder caused by KIF7 defects and belongs to the heterogeneous group of ciliopathies related to Joubert syndrome (JBTS). While ACLS is characterized by macrocephaly, prominent forehead, depressed nasal bridge, and hypertelorism, facial dysmorphism has not been emphasized in JBTS cohorts with molecular diagnosis. To evaluate the specificity and etiology of ACLS craniofacial features, we performed whole exome or targeted Sanger sequencing in patients with the aforementioned overlapping craniofacial appearance but variable additional ciliopathy features followed by functional studies. We found (likely) pathogenic variants of KIF7 in 5 out of 9 families, including the original ACLS patients, and delineated 1000 to 4000-year-old Swiss founder alleles. Three of the remaining families had (likely) pathogenic variants in the JBTS gene C5orf42, and one patient had a novel de novo frameshift variant in SHH known to cause autosomal dominant holoprosencephaly. In accordance with the patients' craniofacial anomalies, we showed facial midline widening after silencing of C5orf42 in chicken embryos. We further supported the link between KIF7, SHH, and C5orf42 by demonstrating abnormal primary cilia and diminished response to a SHH agonist in fibroblasts of C5orf42-mutated patients, as well as axonal pathfinding errors in C5orf42-silenced chicken embryos similar to those observed after perturbation of Shh signaling. Our findings, therefore, suggest that beside the neurodevelopmental features, macrocephaly and facial widening are likely more general signs of disturbed SHH signaling. Nevertheless, long-term follow-up revealed that C5orf42-mutated patients showed catch-up development and fainting of facial features contrary to KIF7-mutated patients.

Genotype imputation for single nucleotide polymorphisms (SNPs) has been shown to be a powerful means to include genetic markers in exploratory genetic association studies without having to genotype them, and is becoming a standard procedure. A number of different software programs are available. In our experience, user-friendliness is often the deciding factor in the choice of software to solve a particular task. We therefore evaluated the usability of three publicly available imputation programs: BEAGLE, IMPUTE and MACH. We found all three programs to perform well with HapMap reference data, with little effort needed for data preparation and subsequent association analysis. Each of them has different strengths and weaknesses, however, and none is optimal for all situations.

AIDS has changed from a mostly male-specific health problem to one that predominantly affects females. Although sex differences in HIV-1 susceptibility are beyond doubt, the extent to which sex affects the onset and progression of AIDS has remained elusive. Here, we provide evidence for an influence of X chromosomal variation on the course of retroviral infection, both in HIV-1-infected patients and in the rhesus macaque model of AIDS. A two-stage, microsatellite-based GWAS of SIV-infected monkeys revealed MHC class I markers and a hitherto-unknown X chromosomal locus as being associated with a nominal score measuring progression to AIDS (Fisher's exact p < 10(-6)). The X chromosomal association was subsequently confirmed in HIV-1-infected patients with published SNP genotype data. SNP rs5968255, located at human Xq21.1 in a conserved sequence element near the RPS6KA6 and CYLC1 genes, was identified as a significant genetic determinant of disease progression in females (ANOVA p = 8.8 x 10(-5)), but not in males (p = 0.19). Heterozygous female carriers of the C allele showed significantly slower CD4 cell decline and a lower viral load at set point than TT homozygous females and than males. Inspection of HapMap revealed that the CT genotype is significantly more frequent among Asians than among Europeans or Africans. Our results suggest that, in addition to the individual innate and adaptive immunity status, sex-linked genetic variation impacts upon the rate of progression to AIDS. Elucidating the mechanisms underlying this sex-specific effect will promote the development of antiretroviral therapies with high efficacy in both sexes.