Deborah A. Nickerson

Joubert syndrome (JBTS) is a recessive neurodevelopmental ciliopathy characterized by a pathognomonic hindbrain malformation. All known JBTS genes encode proteins involved in the structure or function of primary cilia, ubiquitous antenna-like organelles essential for cellular signal transduction. Here, we used the recently identified JBTS-associated protein armadillo repeat motif–containing 9 (ARMC9) in tandem-affinity purification and yeast 2-hybrid screens to identify a ciliary module whose dysfunction underlies JBTS. In addition to the known JBTS-associated proteins CEP104 and CSPP1, we identified coiled-coil domain containing 66 (CCDC66) and TOG array regulator of axonemal microtubules 1 (TOGARAM1) as ARMC9 interaction partners. We found that TOGARAM1 variants cause JBTS and disrupt TOGARAM1 interaction with ARMC9. Using a combination of protein interaction analyses, characterization of patient-derived fibroblasts, and analysis of CRISPR/Cas9-engineered zebrafish and hTERT-RPE1 cells, we demonstrated that dysfunction of ARMC9 or TOGARAM1 resulted in short cilia with decreased axonemal acetylation and polyglutamylation, but relatively intact transition zone function. Aberrant serum-induced ciliary resorption and cold-induced depolymerization in ARMC9 and TOGARAM1 patient cell lines suggest a role for this new JBTS-associated protein module in ciliary stability.

Baseline lung function, quantified as forced expiratory volume in the first second of exhalation (FEV ₁ ), is a standard diagnostic criterion used by clinicians to identify and classify lung diseases. Using whole genome sequencing data from the National Heart, Lung, and Blood Institute TOPMed project, we identified a novel genetic association with FEV ₁ on chromosome 12 in 867 African American children with asthma (p = 1.26 × 10 ⁻⁸ , β = 0.302). Conditional analysis within 1 Mb of the tag signal (rs73429450) yielded one major and two other weaker independent signals within this peak. We explored statistical and functional evidence for all variants in linkage disequilibrium with the three independent signals and yielded 9 variants as the most likely candidates responsible for the association with FEV ₁ . Hi-C data and eQTL analysis demonstrated that these variants physically interacted with KITLG (aka SCF) and their minor alleles were associated with increased expression of KITLG gene in nasal epithelial cells. Gene-by-air-pollution interaction analysis found that the candidate variant rs58475486 interacted with past-year SO ₂ exposure (p = 0.003, β = 0.32). This study identified a novel protective genetic association with FEV ₁ , possibly mediated through KITLG , in African American children with asthma.

Nonphotosensitive trichothiodystrophy (TTDN) is a rare autosomal recessive disorder of neuroectodermal origin. The condition is marked by hair abnormalities, intellectual impairment, nail dystrophies and susceptibility to infections but with no UV sensitivity. We identified three consanguineous Pakistani families with varied TTDN features and used homozygosity mapping, linkage analysis, and Sanger and exome sequencing in order to identify pathogenic variants. Haplotype analysis was performed and haplotype age estimated. A splicing assay was used to validate the effect of the MPLKIP splice variant on expression. Affected individuals from all families exhibit several TTDN features along with a heart-specific feature, i.e. mitral regurgitation. Exome sequencing in the probands from families ED168 and ED241 identified a homozygous splice mutation c.339 + 1G > A within MPLKIP. The same splice variant co-segregates with TTDN in a third family ED210. The MPLKIP splice variant was not found in public databases, e.g. the Exome Aggregation Consortium, and in unrelated Pakistani controls. Functional analysis of the splice variant confirmed intron retention, which leads to protein truncation and loss of a phosphorylation site. Haplotype analysis identified a 585.1-kb haplotype which includes the MPLKIP variant, supporting the existence of a founder haplotype that is estimated to be 25,900 years old. This study extends the allelic and phenotypic spectra of MPLKIP-related TTDN, to include a splice variant that causes cardiomyopathy as part of the TTDN phenotype.

To identify polymorphisms associated with variability of iron overload severity in HFE‐associated hemochromatosis, we performed exome sequencing of DNA from 35 male HFE C282Y homozygotes with either markedly increased iron stores (n = 22; cases) or with normal or mildly increased iron stores (n = 13; controls). The 35 participants, residents of the United States, Canada, and Australia, reported no or light alcohol consumption. Sequencing data included 82,068 single‐nucleotide variants, and 10,337 genes were tested for a difference between cases and controls. A variant in the GNPAT gene showed the most significant association with severe iron overload (P = 3 × 10⁻⁶; P = 0.033 by the likelihood ratio test after correction for multiple comparisons). Sixteen of twenty‐two participants with severe iron overload had glyceronephosphate O‐acyltransferase (GNPAT) polymorphism p.D519G (rs11558492; 15 heterozygotes, one homozygote). No control participant had this polymorphism. To examine functional consequences of GNPAT deficiency, we performed small interfering RNA–based knockdown of GNPAT in the human liver‐derived cell line, HepG2/C3A. This knockdown resulted in a >17‐fold decrease in expression of the messenger RNA encoding the iron‐regulatory hormone, hepcidin. Conclusion: GNPAT p.D519G is associated with a high‐iron phenotype in HFE C282Y homozygotes and may participate in hepcidin regulation. (Hepatology 2015;62:429–439

Efforts to define the genetic architecture underlying variable statin response have met with limited success, possibly because previous studies were limited to effect based on a single dose. We lever...

The emergence of massively parallel DNA sequencing platforms has made resequencing an affordable approach to study genetic variation. However, the cost of whole genome resequencing remains too high to apply to large numbers of human samples. Genomic partitioning methods allow enrichment for regions of interest at a scale that is matched to the throughput of the new sequencing platforms. We review general categories of methods for genomic partitioning including multiplex PCR, capture-by-circularization, and capture-by-hybridization. Parameters that are relevant to the performance of any given method include multiplexity, specificity, uniformity, input requirements, scalability, and cost. The successful development of genomic partitioning strategies will be key to taking full advantage of massively parallel sequencing, at least until resequencing of complete mammalian genomes becomes widely affordable.

Common genetic polymorphism may explain a portion of the heritable risk for common diseases, so considerable effort has been devoted to finding and typing common single-nucleotide polymorphisms (SNPs) in the human genome. Many SNPs show correlated genotypes, or linkage disequilibrium (LD), suggesting that only a subset of all SNPs (known as tagging SNPs, or tagSNPs) need to be genotyped for disease association studies. Based on the genetic differences that exist among human populations, most tagSNP sets are defined in a single population and applied only in populations that are closely related. To improve the efficiency of multi-population analyses, we have developed an algorithm called MultiPop-TagSelect that finds a near-minimal union of population-specific tagSNP sets across an arbitrary number of populations. We present this approach as an extension of LD-select, a tagSNP selection method that uses a greedy algorithm to group SNPs into bins based on their pairwise association patterns, although the MultiPop-TagSelect algorithm could be used with any SNP tagging approach that allows choices between nearly equivalent SNPs. We evaluate the algorithm by considering tagSNP selection in candidate-gene resequencing data and lower density whole-chromosome data. Our analysis reveals that an exhaustive search is often intractable, while the developed algorithm can quickly and reliably find near-optimal solutions even for difficult tagSNP selection problems. Using populations of African, Asian, and European ancestry, we also show that an optimal multi-population set of tagSNPs can be substantially smaller (up to 44%) than a typical set obtained through independent or sequential selection.

Pinpointing genetic associations in the human genome relies heavily on the accuracy of the underlying genotype data. Null alleles can generate significant inaccuracies in genotype data and can negatively affect the statistical power of a study. Existing quality control (QC) tests, including tests of Hardy–Weinberg equilibrium, are not sensitive enough to detect the presence of even moderately frequent null alleles in the data. We show that direct analysis of raw data from a quantitative genotyping platform can detect up to 75% of null alleles, even at frequencies below the sensitivity of more traditional methods. Detecting unexpected null alleles not only has benefits in QC of genotype data but may also be valuable in detecting rare, functional null alleles that would otherwise be missed.