Estimating variance components in population scale family trees.

The rapid digitization of genealogical and medical records enables the assembly of extremely large pedigree records spanning millions of individuals and trillions of pairs of relatives. Such pedigrees provide the opportunity to investigate the sociological and epidemiological history of human populations in scales much larger than previously possible. Linear mixed models (LMMs) are routinely used to analyze extremely large animal and plant pedigrees for the purposes of selective breeding. However, LMMs have not been previously applied to analyze population-scale human family trees. Here, we present Sparse Cholesky factorIzation LMM (Sci-LMM), a modeling framework for studying population-scale family trees that combines techniques from the animal and plant breeding literature and from human genetics literature. The proposed framework can construct a matrix of relationships between trillions of pairs of individuals and fit the corresponding LMM in several hours. We demonstrate the capabilities of Sci-LMM via simulation studies and by estimating the heritability of longevity and of reproductive fitness (quantified via number of children) in a large pedigree spanning millions of individuals and over five centuries of human history. Sci-LMM provides a unified framework for investigating the epidemiological history of human populations via genealogical records.

Genetic linkage analysis of a large family identifies FIGN as a candidate modulator of reduced penetrance in heritable pulmonary arterial hypertension.

BACKGROUND: Mapping the genetic component of molecular mechanisms responsible for the reduced penetrance (RP) of rare disorders constitutes one of the most challenging problems in human genetics. Heritable pulmonary arterial hypertension (PAH) is one such disorder characterised by rare mutations mostly occurring in the bone morphogenetic protein receptor type 2 (BMPR2) gene and a wide heterogeneity of penetrance modifier mechanisms. Here, we analyse 32 genotyped individuals from a large Iberian family of 65 members, including 22 carriers of the pathogenic BMPR2 mutation c.1472G>A (p.Arg491Gln), 8 of them diagnosed with PAH by right-heart catheterisation, leading to an RP rate of 36.4%. METHODS: We performed a linkage analysis on the genotyping data to search for genetic modifiers of penetrance. Using functional genomics data, we characterised the candidate region identified by linkage analysis. We also predicted the haplotype segregation within the family. RESULTS: We identified a candidate chromosome region in 2q24.3, 38 Mb upstream from BMPR2, with significant linkage (LOD=4.09) under a PAH susceptibility model. This region contains common variants associated with vascular aetiology and shows functional evidence that the putative genetic modifier is located in the upstream distal promoter of the fidgetin (FIGN) gene. CONCLUSION: Our results suggest that the genetic modifier acts through FIGN transcriptional regulation, whose expression variability would contribute to modulating heritable PAH. This finding may help to advance our understanding of RP in PAH across families sharing the p.Arg491Gln pathogenic mutation in BMPR2.

Multikernel linear mixed models for complex phenotype prediction.

Linear mixed models (LMMs) and their extensions have recently become the method of choice in phenotype prediction for complex traits. However, LMM use to date has typically been limited by assuming simple genetic architectures. Here, we present multikernel linear mixed model (MKLMM), a predictive modeling framework that extends the standard LMM using multiple-kernel machine learning approaches. MKLMM can model genetic interactions and is particularly suitable for modeling complex local interactions between nearby variants. We additionally present MKLMM-Adapt, which automatically infers interaction types across multiple genomic regions. In an analysis of eight case-control data sets from the Wellcome Trust Case Control Consortium and more than a hundred mouse phenotypes, MKLMM-Adapt consistently outperforms competing methods in phenotype prediction. MKLMM is as computationally efficient as standard LMMs and does not require storage of genotypes, thus achieving state-of-the-art predictive power without compromising computational feasibility or genomic privacy.

Accurate liability estimation improves power in ascertained case-control studies.

Linear mixed models (LMMs) have emerged as the method of choice for confounded genome-wide association studies. However, the performance of LMMs in nonrandomly ascertained case-control studies deteriorates with increasing sample size. We propose a framework called LEAP (liability estimator as a phenotype; https://github.com/omerwe/LEAP) that tests for association with estimated latent values corresponding to severity of phenotype, and we demonstrate that this can lead to a substantial power increase.

Phospholipase A2-activating protein is associated with a novel form of leukoencephalopathy.

Leukoencephalopathies are a group of white matter disorders related to abnormal formation, maintenance, and turnover of myelin in the central nervous system. These disorders of the brain are categorized according to neuroradiological and pathophysiological criteria. Herein, we have identified a unique form of leukoencephalopathy in seven patients presenting at ages 2 to 4 months with progressive microcephaly, spastic quadriparesis, and global developmental delay. Clinical, metabolic, and imaging characterization of seven patients followed by homozygosity mapping and linkage analysis were performed. Next generation sequencing, bioinformatics, and segregation analyses followed, to determine a loss of function sequence variation in the phospholipase A2-activating protein encoding gene (PLAA). Expression and functional studies of the encoded protein were performed and included measurement of prostaglandin E2 and cytosolic phospholipase A2 activity in membrane fractions of fibroblasts derived from patients and healthy controls. Plaa-null mice were generated and prostaglandin E2 levels were measured in different tissues. The novel phenotype of our patients segregated with a homozygous loss-of-function sequence variant, causing the substitution of leucine at position 752 to phenylalanine, in PLAA, which causes disruption of the protein's ability to induce prostaglandin E2 and cytosolic phospholipase A2 synthesis in patients' fibroblasts. Plaa-null mice were perinatal lethal with reduced brain levels of prostaglandin E2 The non-functional phospholipase A2-activating protein and the associated neurological phenotype, reported herein for the first time, join other complex phospholipid defects that cause leukoencephalopathies in humans, emphasizing the importance of this axis in white matter development and maintenance.

Induced pluripotent stem cells as a model for telomeric abnormalities in ICF type I syndrome.

Human telomeric regions are packaged as constitutive heterochromatin, characterized by extensive subtelomeric DNA methylation and specific histone modifications. ICF (immunodeficiency, centromeric instability, facial anomalies) type I patients carry mutations in DNA methyltransferase 3B (DNMT3B) that methylates de novo repetitive sequences during early embryonic development. ICF type I patient fibroblasts display hypomethylated subtelomeres, abnormally short telomeres and premature senescence. In order to study the molecular mechanism by which the failure to de novo methylate subtelomeres results in accelerated telomere shortening, we generated induced pluripotent stem cells (iPSCs) from 3 ICF type I patients. Telomeres were elongated in ICF-iPSCs during reprogramming, and the senescence phenotype was abolished despite sustained subtelomeric hypomethylation and high TERRA levels. Fibroblast-like cells (FLs) isolated from differentiated ICF-iPSCs maintained abnormally high TERRA levels, and telomeres in these cells shortened at an accelerated rate, leading to early senescence, thus recapitulating the telomeric phenotype of the parental fibroblasts. These findings demonstrate that the abnormal telomere phenotype associated with subtelomeric hypomethylation is overridden in cells expressing telomerase, therefore excluding telomerase inhibition by TERRA as a central mechanism responsible for telomere shortening in ICF syndrome. The data in the current study lend support to the use of ICF-iPSCs for modeling of phenotypic and molecular defects in ICF syndrome and for unraveling the mechanism whereby subtelomeric hypomethylation is linked to accelerated telomeric loss in this syndrome.

Seborrhea-like dermatitis with psoriasiform elements caused by a mutation in ZNF750, encoding a putative C2H2 zinc finger protein.

We describe an Israeli Jewish Moroccan family presenting with autosomal dominant seborrhea-like dermatosis with psoriasiform elements, including enhanced keratinocyte proliferation, parakeratosis, follicular plugging, Pityrosporum ovale overgrowth and dermal CD4 lymphocyte infiltrate. We mapped the disease gene to a 0.5-cM region overlapping the PSORS2 locus (17q25) and identified a frameshift mutation in ZNF750, which encodes a putative C2H2 zinc finger protein. ZNF750 is normally expressed in keratinocytes but not in fibroblasts and is barely detectable in CD4 lymphocytes.

A system for exact and approximate genetic linkage analysis of SNP data in large pedigrees.

MOTIVATION: The use of dense single nucleotide polymorphism (SNP) data in genetic linkage analysis of large pedigrees is impeded by significant technical, methodological and computational challenges. Here we describe Superlink-Online SNP, a new powerful online system that streamlines the linkage analysis of SNP data. It features a fully integrated flexible processing workflow comprising both well-known and novel data analysis tools, including SNP clustering, erroneous data filtering, exact and approximate LOD calculations and maximum-likelihood haplotyping. The system draws its power from thousands of CPUs, performing data analysis tasks orders of magnitude faster than a single computer. By providing an intuitive interface to sophisticated state-of-the-art analysis tools coupled with high computing capacity, Superlink-Online SNP helps geneticists unleash the potential of SNP data for detecting disease genes. RESULTS: Computations performed by Superlink-Online SNP are automatically parallelized using novel paradigms, and executed on unlimited number of private or public CPUs. One novel service is large-scale approximate Markov Chain-Monte Carlo (MCMC) analysis. The accuracy of the results is reliably estimated by running the same computation on multiple CPUs and evaluating the Gelman-Rubin Score to set aside unreliable results. Another service within the workflow is a novel parallelized exact algorithm for inferring maximum-likelihood haplotyping. The reported system enables genetic analyses that were previously infeasible. We demonstrate the system capabilities through a study of a large complex pedigree affected with metabolic syndrome. AVAILABILITY: Superlink-Online SNP is freely available for researchers at http://cbl-hap.cs.technion.ac.il/superlink-snp. The system source code can also be downloaded from the system website. CONTACT: omerw@cs.technion.ac.il SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A founder mutation in Vps37A causes autosomal recessive complex hereditary spastic paraparesis.

BACKGROUND: Members of two seemingly unrelated kindreds of Arab Moslem origin presented with pronounced early onset spastic paraparesis of upper and lower limbs, mild intellectual disability, kyphosis, pectus carinatum and hypertrichosis. METHODS: The authors performed neurological and developmental examinations on the affected individuals. The authors conducted whole genome linkage and haplotype analyses, followed by sequencing of candidate genes; RNA and protein expression studies; and finally proof of principle investigations on knockdown morpholino oligonucleotide injected zebrafish. RESULTS: The authors characterise a novel form of autosomal recessive complex hereditary spastic paraparesis (CHSP). MRI studies of brain and spinal cord were normal. Within a single significantly linked locus the authors ultimately identified a homozygous missense mutation c.1146A>T (p.K382N) in the vacuolar protein sorting 37A (Vps37A) gene, fully penetrant and segregating with the disease in both families. Mobility was significantly reduced in Vps37A knockdown morpholino oligonucleotide injected zebrafish, supporting the causal relationship between mutations in this gene and the phenotype described in the patients of this study. CONCLUSIONS: The authors provide evidence for the involvement of Vps37A, a member of the endosomal sorting complex required for transport (ESCRT) system, in upper motor neuron disease. The ESCRT system has been shown to play a central role in intracellular trafficking, in the maturation of multivesicular bodies and the sorting of ubiquitinated membrane proteins into internal luminal vesicles. Further investigation of mechanisms by which dysfunction of this gene causes CHSP will contribute to the understanding of intracellular trafficking of vesicles by the ESCRT machinery and its relevance to CHSP.

Spinal muscular atrophy with pontocerebellar hypoplasia is caused by a mutation in the VRK1 gene.

The spinal muscular atrophies (SMAs) are a genetically and clinically heterogeneous group of disorders characterized by degeneration and loss of anterior horn cells in the spinal cord, leading to muscle weakness and atrophy. Spinal muscular atrophy with pontocerebellar hypoplasia (SMA-PCH, also known as pontocerebellar hypoplasia type 1 [PCH1]) is one of the rare infantile SMA variants that include additional clinical manifestations, and its genetic basis is unknown. We used a homozygosity mapping and positional cloning approach in a consanguineous family of Ashkenazi Jewish origin and identified a nonsense mutation in the vaccinia-related kinase 1 gene (VRK1) as a cause of SMA-PCH. VRK1, one of three members of the mammalian VRK family, is a serine/threonine kinase that phosphorylates p53 and CREB and is essential for nuclear envelope formation. Its identification as a gene involved in SMA-PCH implies new roles for the VRK proteins in neuronal development and maintenance and suggests the VRK genes as candidates for related phenotypes.

Integration of SNP genotyping confidence scores in IBD inference.

MOTIVATION: High-throughput single nucleotide polymorphism (SNP) arrays have become the standard platform for linkage and association analyses. The high SNP density of these platforms allows high-resolution identification of ancestral recombination events even for distant relatives many generations apart. However, such inference is sensitive to marker mistyping and current error detection methods rely on the genotyping of additional close relatives. Genotyping algorithms provide a confidence score for each marker call that is currently not integrated in existing methods. There is a need for a model that incorporates this prior information within the standard identical by descent (IBD) and association analyses. RESULTS: We propose a novel model that incorporates marker confidence scores within IBD methods based on the Lander-Green Hidden Markov Model. The novel parameter of this model is the joint distribution of confidence scores and error status per array. We estimate this probability distribution by applying a modified expectation-maximization (EM) procedure on data from nuclear families genotyped with Affymetrix 250K SNP arrays. The converged tables from two different genotyping algorithms are shown for a wide range of error rates. We demonstrate the efficacy of our method in refining the detection of IBD signals using nuclear pedigrees and distant relatives. AVAILABILITY: Plinke, a new version of Plink with an extended pairwise IBD inference model allowing per marker error probabilities is freely available at: http://bioinfo.bgu.ac.il/bsu/software/plinke. CONTACT: obirk@bgu.ac.il; markusb@bgu.ac.il SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Gray platelet syndrome: natural history of a large patient cohort and locus assignment to chromosome 3p.

Gray platelet syndrome (GPS) is an inherited bleeding disorder characterized by macrothrombocytopenia and absence of platelet α-granules resulting in typical gray platelets on peripheral smears. GPS is associated with a bleeding tendency, myelofibrosis, and splenomegaly. Reports on GPS are limited to case presentations. The causative gene and underlying pathophysiology are largely unknown. We present the results of molecular genetic analysis of 116 individuals including 25 GPS patients from 14 independent families as well as novel clinical data on the natural history of the disease. The mode of inheritance was autosomal recessive (AR) in 11 and indeterminate in 3 families. Using genome-wide linkage analysis, we mapped the AR-GPS gene to a 9.4-Mb interval on 3p21.1-3p22.1, containing 197 protein-coding genes. Sequencing of 1423 (69%) of the 2075 exons in the interval did not identify the GPS gene. Long-term follow-up data demonstrated the progressive nature of the thrombocytopenia and myelofibrosis of GPS resulting in fatal hemorrhages in some patients. We identified high serum vitamin B(12) as a consistent, novel finding in GPS. Chromosome 3p21.1-3p22.1 has not been previously linked to a platelet disorder; identification of the GPS gene will likely lead to the discovery of novel components of platelet organelle biogenesis. This study is registered at www.clinicaltrials.gov as NCT00069680 and NCT00369421.

Genetic linkage analysis in the presence of germline mosaicism.

Germline mosaicism is a genetic condition in which some germ cells of an individual contain a mutation. This condition violates the assumptions underlying classic genetic analysis and may lead to failure of such analysis. In this work we extend the statistical model used for genetic linkage analysis in order to incorporate germline mosaicism. We develop a likelihood ratio test for detecting whether a genetic trait has been introduced into a pedigree by germline mosaicism. We analyze the statistical properties of this test and evaluate its performance via computer simulations. We demonstrate that genetic linkage analysis has high power to identify linkage in the presence of germline mosaicism when our extended model is used. We further use this extended model to provide solid statistical evidence that the MDN syndrome studied by Genzer-Nir et al. has been introduced by germline mosaicism.

Alopecia, neurological defects, and endocrinopathy syndrome caused by decreased expression of RBM28, a nucleolar protein associated with ribosome biogenesis.

Single-gene disorders offer unique opportunities to shed light upon fundamental physiological processes in humans. We investigated an autosomal-recessive phenotype characterized by alopecia, progressive neurological defects, and endocrinopathy (ANE syndrome). By using homozygosity mapping and candidate-gene analysis, we identified a loss-of-function mutation in RBM28, encoding a nucleolar protein. RBM28 yeast ortholog, Nop4p, was previously found to regulate ribosome biogenesis. Accordingly, electron microscopy revealed marked ribosome depletion and structural abnormalities of the rough endoplasmic reticulum in patient cells, ascribing ANE syndrome to the restricted group of inherited disorders associated with ribosomal dysfunction.

Estimating genome-wide IBD sharing from SNP data via an efficient hidden Markov model of LD with application to gene mapping.

MOTIVATION: Association analysis is the method of choice for studying complex multifactorial diseases. The premise of this method is that affected persons contain some common genomic regions with similar SNP alleles and such areas will be found in this analysis. An important disadvantage of GWA studies is that it does not distinguish between genomic areas that are inherited from a common ancestor [identical by descent (IBD)] and areas that are identical merely by state [identical by state (IBS)]. Clearly, areas that can be marked with higher probability as IBD and have the same correlation with the disease status of identical areas that are more probably only IBS, are better candidates to be causative, and yet this distinction is not encoded in standard association analysis. RESULTS: We develop a factorial hidden Markov model-based algorithm for computing genome-wide IBD sharing. The algorithm accepts as input SNP data of measured individuals and estimates the probability of IBD at each locus for every pair of individuals. For two g-degree relatives, when g > or = 8, the computation yields a precision of IBD tagging of over 50% higher than previous methods for 95% recall. Our algorithm uses a first-order Markovian model for the linkage disequilibrium process and employs a reduction of the state space of the inheritance vector from being exponential in g to quadratic. The higher accuracy along with the reduced time complexity marks our method as a feasible means for IBD mapping in practical scenarios. AVAILABILITY: A software implementation, called IBDMAP, is freely available at http://bioinfo.cs.technion.ac.il/IBDmap.

Speeding up HMM algorithms for genetic linkage analysis via chain reductions of the state space.

UNLABELLED: We develop an hidden Markov model (HMM)-based algorithm for computing exact parametric and non-parametric linkage scores in larger pedigrees than was possible before. The algorithm is applicable whenever there are chains of persons in the pedigree with no genetic measurements and with unknown affection status. The algorithm is based on shrinking the state space of the HMM considerably using such chains. In a two g-degree cousins pedigree the reduction drops the state space from being exponential in g to being linear in g. For a Finnish family in which two affected children suffer from a rare cold-inducing sweating syndrome, we were able to reduce the state space by more than five orders of magnitude from 2(50) to 2(32). In another pedigree of state-space size of 2(27), used for a study of pituitary adenoma, the state space reduced by a factor of 8.5 and consequently exact linkage scores can now be computed, rather than approximated. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Polyadenylation and degradation of human mitochondrial RNA: the prokaryotic past leaves its mark.

RNA polyadenylation serves a purpose in bacteria and organelles opposite from the role it plays in nuclear systems. The majority of nucleus-encoded transcripts are characterized by stable poly(A) tails at their mature 3' ends, which are essential for stabilization and translation initiation. In contrast, in bacteria, chloroplasts, and plant mitochondria, polyadenylation is a transient feature which promotes RNA degradation. Surprisingly, in spite of their prokaryotic origin, human mitochondrial transcripts possess stable 3'-end poly(A) tails, akin to nucleus-encoded mRNAs. Here we asked whether human mitochondria retain truncated and transiently polyadenylated transcripts in addition to stable 3'-end poly(A) tails, which would be consistent with the preservation of the largely ubiquitous polyadenylation-dependent RNA degradation mechanisms of bacteria and organelles. To this end, using both molecular and bioinformatic methods, we sought and revealed numerous examples of such molecules, dispersed throughout the mitochondrial genome. The broad distribution but low abundance of these polyadenylated truncated transcripts strongly suggests that polyadenylation-dependent RNA degradation occurs in human mitochondria. The coexistence of this system with stable 3'-end polyadenylation, despite their seemingly opposite effects, is so far unprecedented in bacteria and other organelles.

Polyadenylation of ribosomal RNA in human cells.

The addition of poly(A)-tails to RNA is a process common to almost all organisms. In eukaryotes, stable poly(A)-tails, important for mRNA stability and translation initiation, are added to the 3' ends of most nuclear-encoded mRNAs, but not to rRNAs. Contrarily, in prokaryotes and organelles, polyadenylation stimulates RNA degradation. Recently, polyadenylation of nuclear-encoded transcripts in yeast was reported to promote RNA degradation, demonstrating that polyadenylation can play a double-edged role for RNA of nuclear origin. Here we asked whether in human cells ribosomal RNA can undergo polyadenylation. Using both molecular and bioinformatic approaches, we detected non-abundant polyadenylated transcripts of the 18S and 28S rRNAs. Interestingly, many of the post-transcriptionally added tails were composed of heteropolymeric poly(A)-rich sequences containing the other nucleotides in addition to adenosine. These polyadenylated RNA fragments are most likely degradation intermediates, as primer extension (PE) analysis revealed the presence of distal fragmented molecules, some of which matched the polyadenylation sites of the proximal cleavage products revealed by oligo(dT) and circled RT-PCR. These results suggest the presence of a mechanism to degrade ribosomal RNAs in human cells, that possibly initiates with endonucleolytic cleavages and involves the addition of poly(A) or poly(A)-rich tails to truncated transcripts, similar to that which operates in prokaryotes and organelles.

Quantitative analysis of population-scale family trees with millions of relatives.

Family trees have vast applications in fields as diverse as genetics, anthropology, and economics. However, the collection of extended family trees is tedious and usually relies on resources with limited geographical scope and complex data usage restrictions. We collected 86 million profiles from publicly available online data shared by genealogy enthusiasts. After extensive cleaning and validation, we obtained population-scale family trees, including a single pedigree of 13 million individuals. We leveraged the data to partition the genetic architecture of human longevity and to provide insights into the geographical dispersion of families. We also report a simple digital procedure to overlay other data sets with our resource.