Significant interest has emerged in mapping genetic susceptibility for complex traits through whole-genome association studies. These studies rely on the extent of association, i.e., linkage disequilibrium (LD), between single nucleotide polymorphisms (SNPs) across the human genome. LD describes the nonrandom association between SNP pairs and can be used as a metric when designing maximally informative panels of SNPs for association studies in human populations. Using data from the 1.58 million SNPs genotyped by Perlegen, we explored the allele frequency dependence of the LD statistic r(2) both empirically and theoretically. We show that average r(2) values between SNPs unmatched for allele frequency are always limited to much less than 1 (theoretical approximately 0.46 to 0.57 for this dataset). Frequency matching of SNP pairs provides a more sensitive measure for assessing the average decay of LD and generates average r(2) values across nearly the entire informative range (from 0 to 0.89 through 0.95). Additionally, we analyzed the extent of perfect LD (r(2) = 1.0) using frequency-matched SNPs and found significant differences in the extent of LD in genic regions versus intergenic regions. The SNP pairs exhibiting perfect LD showed a significant bias for derived, nonancestral alleles, providing evidence for positive natural selection in the human genome.

Identifying regions of the human genome that have been targets of natural selection will provide important insights into human evolutionary history and may facilitate the identification of complex disease genes. Although the signature that natural selection imparts on DNA sequence variation is difficult to disentangle from the effects of neutral processes such as population demographic history, selective and demographic forces can be distinguished by analyzing multiple loci dispersed throughout the genome. We studied the molecular evolution of 132 genes by comprehensively resequencing them in 24 African-Americans and 23 European-Americans. We developed a rigorous computational approach for taking into account multiple hypothesis tests and demographic history and found that while many apparent selective events can instead be explained by demography, there is also strong evidence for positive or balancing selection at eight genes in the European-American population, but none in the African-American population. Our results suggest that the migration of modern humans out of Africa into new environments was accompanied by genetic adaptations to emergent selective forces. In addition, a region containing four contiguous genes on Chromosome 7 showed striking evidence of a recent selective sweep in European-Americans. More generally, our results have important implications for mapping genes underlying complex human diseases.

Genetic structure in the European American population reflects waves of migration and recent gene flow among different populations. This complex structure can introduce bias in genetic association studies. Using Principal Components Analysis (PCA), we analyze the structure of two independent European American datasets (1,521 individuals-307,315 autosomal SNPs). Individual variation lies across a continuum with some individuals showing high degrees of admixture with non-European populations, as demonstrated through joint analysis with HapMap data. The CEPH Europeans only represent a small fraction of the variation encountered in the larger European American datasets we studied. We interpret the first eigenvector of this data as correlated with ancestry, and we apply an algorithm that we have previously described to select PCA-informative markers (PCAIMs) that can reproduce this structure. Importantly, we develop a novel method that can remove redundancy from the selected SNP panels and show that we can effectively remove correlated markers, thus increasing genotyping savings. Only 150-200 PCAIMs suffice to accurately predict fine structure in European American datasets, as identified by PCA. Simulating association studies, we couple our method with a PCA-based stratification correction tool and demonstrate that a small number of PCAIMs can efficiently remove false correlations with almost no loss in power. The structure informative SNPs that we propose are an important resource for genetic association studies of European Americans. Furthermore, our redundancy removal algorithm can be applied on sets of ancestry informative markers selected with any method in order to select the most uncorrelated SNPs, and significantly decreases genotyping costs.

The short-rib polydactyly syndromes (SRPS) encompass a radiographically and genetically heterogeneous group of skeletal ciliopathies that are characterized by a long narrow chest, short extremities, and variable occurrence of polydactyly. Radiographic abnormalities include undermineralization of the calvarium, shortened and bowed appendicular bones, trident shaped acetabula and polydactyly. In a case of SRPS we identified compound heterozygosity for mutations in IFT52, which encodes a component of the anterograde intraflagellar transport complex. The IFT52 mutant cells synthesized a significantly reduced amount of IFT52 protein, leading to reduced synthesis of IFT74, IFT81, IFT88 and ARL13B, other key anterograde complex members. Ciliogenesis was also disrupted in the mutant cells, with a 60% reduction in the presence of cilia on mutant cells and loss of cilia length regulation for the cells with cilia. These data demonstrate that IFT52 is essential for anterograde complex integrity and for the biosynthesis and maintenance of cilia. The data identify a new locus for SRPS and show that IFT52 mutations result in a ciliopathy with primary effects on the skeleton.

We conducted a genome-wide association analysis of 7 subfractions of low density lipoproteins (LDLs) and 3 subfractions of intermediate density lipoproteins (IDLs) measured by gradient gel electrophoresis, and their response to statin treatment, in 1868 individuals of European ancestry from the Pharmacogenomics and Risk of Cardiovascular Disease study. Our analyses identified four previously-implicated loci (SORT1, APOE, LPA, and CETP) as containing variants that are very strongly associated with lipoprotein subfractions (log(10)Bayes Factor > 15). Subsequent conditional analyses suggest that three of these (APOE, LPA and CETP) likely harbor multiple independently associated SNPs. Further, while different variants typically showed different characteristic patterns of association with combinations of subfractions, the two SNPs in CETP show strikingly similar patterns--both in our original data and in a replication cohort--consistent with a common underlying molecular mechanism. Notably, the CETP variants are very strongly associated with LDL subfractions, despite showing no association with total LDLs in our study, illustrating the potential value of the more detailed phenotypic measurements. In contrast with these strong subfraction associations, genetic association analysis of subfraction response to statins showed much weaker signals (none exceeding log(10)Bayes Factor of 6). However, two SNPs (in APOE and LPA) previously-reported to be associated with LDL statin response do show some modest evidence for association in our data, and the subfraction response proles at the LPA SNP are consistent with the LPA association, with response likely being due primarily to resistance of Lp(a) particles to statin therapy. An additional important feature of our analysis is that, unlike most previous analyses of multiple related phenotypes, we analyzed the subfractions jointly, rather than one at a time. Comparisons of our multivariate analyses with standard univariate analyses demonstrate that multivariate analyses can substantially increase power to detect associations. Software implementing our multivariate analysis methods is available at http://stephenslab.uchicago.edu/software.html.

Isolated hand syndactyly is a common limb malformation with limited known genetic etiology. We used exome sequencing to discover two novel variants, chr11 g.46896373C>G; p.D1403H and chr11 g.46893078G>T; p.Q1564K, in LRP4 in a child with isolated bilateral syndactyly of the third and fourth fingers. Each variant was inherited from a different parent and neither parent was affected. Variants in LRP4 have been previously associated with syndactyly in Cenani-Lenz syndactyly syndrome and Sclerosteosis 2, but have not been reported in individuals with isolated syndactyly. LRP4 inhibits LRP6/LRP5-mediated activation of canonical Wnt signaling and mediates sclerostin-dependent inhibition of bone formation. p.D1403H and p.Q1564K are located within the fourth β-propeller of the extracellular protein domain that has yet to be associated with human disease. Functional analyses of p.D1403H and p.Q1564K show that they significantly decrease LRP4's inhibition of Wnt signaling. These results suggest that variants in the fourth β-propeller of the extracellular protein domain may cause a phenotype distinct from previously characterized LRP4 variants.

We report a severe defect of neuromuscular transmission in a consanguineous patient with a homozygous variant in the laminin α5 subunit gene (LAMA5). The variant c.8046C > T (p.Arg2659Trp) is rare and has a predicted deleterious effect. The affected individual, who also carries a rare homozygous sequence variant in LAMA1, had normal cognitive function, but magnetic resonance brain imaging showed mild volume loss and periventricular T2 prolongation. Repetitive nerve stimulation at 2 Hz showed 50% decrement of compound muscle action potential amplitudes but 250% facilitation immediately after exercise, similar to that seen in Lambert-Eaton myasthenic syndrome. Endplate studies demonstrated a profound reduction of the endplate potential quantal content but normal amplitudes of miniature endplate potentials. Electron microscopy showed endplates with increased postsynaptic folding that were denuded or only partially occupied by small nerve terminals. Expression studies revealed that p.Arg2659Trp caused decreased binding of laminin α5 to SV2A and impaired laminin-521 cell adhesion and cell projection support in primary neuronal cultures. In summary, this report describing severe neuromuscular transmission failure in a patient with a LAMA5 mutation expands the list of phenotypes associated with defects in genes encoding α-laminins.

Acheiropodia, congenital limb truncation, is associated with homozygous deletions in the LMBR1 gene around ZRS, an enhancer regulating SHH during limb development. How these deletions lead to this phenotype is unknown. Using whole-genome sequencing, we fine-mapped the acheiropodia-associated region to 12 kb and show that it does not function as an enhancer. CTCF and RAD21 ChIP-seq together with 4C-seq and DNA FISH identify three CTCF sites within the acheiropodia-deleted region that mediate the interaction between the ZRS and the SHH promoter. This interaction is substituted with other CTCF sites centromeric to the ZRS in the disease state. Mouse knockouts of the orthologous 12 kb sequence have no apparent abnormalities, showcasing the challenges in modelling CTCF alterations in animal models due to inherent motif differences between species. Our results show that alterations in CTCF motifs can lead to a Mendelian condition due to altered enhancer-promoter interactions.

Defects in genes encoding the isoforms of the laminin alpha subunit have been linked to various phenotypic manifestations, including brain malformations, muscular dystrophy, ocular defects, cardiomyopathy, and skin abnormalities. We report here a severe defect of neuromuscular transmission in a consanguineous patient with a homozygous variant in the laminin alpha-5 subunit gene (LAMA5). The variant c.8046C>T (p.Arg2659Trp) is rare and has a predicted deleterious effect. The affected individual, who also carries a rare homozygous sequence variant in LAMA1, had muscle weakness, myopia, and facial tics. Magnetic resonance imaging of brain showed mild volume loss and periventricular T2 prolongation. Repetitive nerve stimulation revealed 50% decrement of compound muscle action potential amplitudes and 250% facilitation immediately after exercise, Endplate studies identified a profound reduction of the endplate potential quantal content and endplates with normal postsynaptic folding that were denuded or partially occupied by small nerve terminals. Expression studies revealed that p.Arg2659Trp caused decreased binding of laminin alpha-5 to SV2A and impaired laminin-521 cell-adhesion and cell projection support in primary neuronal cultures. In summary, this report describing severe neuromuscular transmission failure in a patient with a LAMA5 mutation expands the list of phenotypes associated with defects in genes encoding alpha-laminins.

Large datasets of hundreds to thousands of individuals measuring RNA-seq in observational studies are becoming available. Many popular software packages for analysis of RNA-seq data were constructed to study differences in expression signatures in an experimental design with well-defined conditions (exposures). In contrast, observational studies may have varying levels of confounding transcript-exposure associations; further, exposure measures may vary from discrete (exposed, yes/no) to continuous (levels of exposure), with non-normal distributions of exposure. We compare popular software for gene expression-DESeq2, edgeR and limma-as well as linear regression-based analyses for studying the association of continuous exposures with RNA-seq. We developed a computation pipeline that includes transformation, filtering and generation of empirical null distribution of association P-values, and we apply the pipeline to compute empirical P-values with multiple testing correction. We employ a resampling approach that allows for assessment of false positive detection across methods, power comparison and the computation of quantile empirical P-values. The results suggest that linear regression methods are substantially faster with better control of false detections than other methods, even with the resampling method to compute empirical P-values. We provide the proposed pipeline with fast algorithms in an R package Olivia, and implemented it to study the associations of measures of sleep disordered breathing with RNA-seq in peripheral blood mononuclear cells in participants from the Multi-Ethnic Study of Atherosclerosis.