@article {847, title = {Chromosome-scale, haplotype-resolved assembly of human genomes.}, journal = {Nat Biotechnol}, volume = {39}, year = {2021}, month = {2021 03}, pages = {309-312}, abstract = {

Haplotype-resolved or phased genome assembly provides a complete picture of genomes and their complex genetic variations. However, current algorithms for phased assembly either do not generate chromosome-scale phasing or require pedigree information, which limits their application. We present a method named diploid assembly (DipAsm) that uses long, accurate reads and long-range conformation data for single individuals to generate a chromosome-scale phased assembly within 1 day. Applied to four public human genomes, PGP1, HG002, NA12878 and HG00733, DipAsm produced haplotype-resolved assemblies with minimum contig length needed to cover 50\% of the known genome (NG50) up to 25 Mb and phased ~99.5\% of heterozygous sites at 98-99\% accuracy, outperforming other approaches in terms of both contiguity and phasing completeness. We demonstrate the importance of chromosome-scale phased assemblies for the discovery of structural variants (SVs), including thousands of new transposon insertions, and of highly polymorphic and medically important regions such as the human leukocyte antigen (HLA) and killer cell immunoglobulin-like receptor (KIR) regions. DipAsm will facilitate high-quality precision medicine and studies of individual haplotype variation and population diversity.

}, keywords = {Algorithms, Chromosomes, Human, Genome, Human, Haplotypes, Heterozygote, Humans, Polymorphism, Single Nucleotide}, issn = {1546-1696}, doi = {10.1038/s41587-020-0711-0}, author = {Garg, Shilpa and Fungtammasan, Arkarachai and Carroll, Andrew and Chou, Mike and Schmitt, Anthony and Zhou, Xiang and Mac, Stephen and Peluso, Paul and Hatas, Emily and Ghurye, Jay and Maguire, Jared and Mahmoud, Medhat and Cheng, Haoyu and Heller, David and Zook, Justin M and Moemke, Tobias and Marschall, Tobias and Sedlazeck, Fritz J and Aach, John and Chin, Chen-Shan and Church, George M and Li, Heng} } @article {898, title = {Haplotype-resolved diverse human genomes and integrated analysis of structural variation.}, journal = {Science}, volume = {372}, year = {2021}, month = {2021 04 02}, abstract = {

Long-read and strand-specific sequencing technologies together facilitate the de novo assembly of high-quality haplotype-resolved human genomes without parent-child trio data. We present 64 assembled haplotypes from 32 diverse human genomes. These highly contiguous haplotype assemblies (average minimum contig length needed to cover 50\% of the genome: 26 million base pairs) integrate all forms of genetic variation, even across complex loci. We identified 107,590 structural variants (SVs), of which 68\% were not discovered with short-read sequencing, and 278 SV hotspots (spanning megabases of gene-rich sequence). We characterized 130 of the most active mobile element source elements and found that 63\% of all SVs arise through homology-mediated mechanisms. This resource enables reliable graph-based genotyping from short reads of up to 50,340 SVs, resulting in the identification of 1526 expression quantitative trait loci as well as SV candidates for adaptive selection within the human population.

}, keywords = {Female, Genetic Variation, Genome, Human, Genotype, Haplotypes, High-Throughput Nucleotide Sequencing, Humans, INDEL Mutation, Interspersed Repetitive Sequences, Male, Population Groups, Quantitative Trait Loci, Retroelements, Sequence Analysis, DNA, Sequence Inversion, Whole Genome Sequencing}, issn = {1095-9203}, doi = {10.1126/science.abf7117}, author = {Ebert, Peter and Audano, Peter A and Zhu, Qihui and Rodriguez-Martin, Bernardo and Porubsky, David and Bonder, Marc Jan and Sulovari, Arvis and Ebler, Jana and Zhou, Weichen and Serra Mari, Rebecca and Yilmaz, Feyza and Zhao, Xuefang and Hsieh, PingHsun and Lee, Joyce and Kumar, Sushant and Lin, Jiadong and Rausch, Tobias and Chen, Yu and Ren, Jingwen and Santamarina, Martin and H{\"o}ps, Wolfram and Ashraf, Hufsah and Chuang, Nelson T and Yang, Xiaofei and Munson, Katherine M and Lewis, Alexandra P and Fairley, Susan and Tallon, Luke J and Clarke, Wayne E and Basile, Anna O and Byrska-Bishop, Marta and Corvelo, Andr{\'e} and Evani, Uday S and Lu, Tsung-Yu and Chaisson, Mark J P and Chen, Junjie and Li, Chong and Brand, Harrison and Wenger, Aaron M and Ghareghani, Maryam and Harvey, William T and Raeder, Benjamin and Hasenfeld, Patrick and Regier, Allison A and Abel, Haley J and Hall, Ira M and Flicek, Paul and Stegle, Oliver and Gerstein, Mark B and Tubio, Jose M C and Mu, Zepeng and Li, Yang I and Shi, Xinghua and Hastie, Alex R and Ye, Kai and Chong, Zechen and Sanders, Ashley D and Zody, Michael C and Talkowski, Michael E and Mills, Ryan E and Devine, Scott E and Lee, Charles and Korbel, Jan O and Marschall, Tobias and Eichler, Evan E} } @article {950, title = {8q24 genetic variation and comprehensive haplotypes altering familial risk of prostate cancer.}, journal = {Nat Commun}, volume = {11}, year = {2020}, month = {2020 03 23}, pages = {1523}, abstract = {

The 8q24 genomic locus is tied to the origin of numerous cancers. We investigate its contribution to hereditary prostate cancer (HPC) in independent study populations of the Nashville Familial Prostate Cancer Study and International Consortium for Prostate Cancer Genetics (combined: 2,836 HPC cases, 2,206 controls of European ancestry). Here we report 433 variants concordantly associated with HPC in both study populations, accounting for 9\% of heritability and modifying age of diagnosis as well as aggressiveness; 183 reach genome-wide significance. The variants comprehensively distinguish independent risk-altering haplotypes overlapping the 648 kb locus (three protective, and four risk (peak odds ratios: 1.5, 4, 5, and 22)). Sequence of the near-Mendelian haplotype reveals eleven causal mutation candidates. We introduce a linkage disequilibrium-based algorithm discerning eight independent sentinel variants, carrying considerable risk prediction ability (AUC = 0.625) for a single locus. These findings elucidate 8q24 locus structure and correlates for clinical prediction of prostate cancer risk.

}, keywords = {Aged, Case-Control Studies, Chromosomes, Human, Pair 8, European Continental Ancestry Group, Genetic Loci, Genetic Predisposition to Disease, Genome-Wide Association Study, Haplotypes, Heterozygote, Humans, Linkage Disequilibrium, Male, Middle Aged, Mutation, Prostatic Neoplasms, Protective Factors, Risk Factors}, issn = {2041-1723}, doi = {10.1038/s41467-020-15122-1}, author = {Dupont, William D and Breyer, Joan P and Plummer, W Dale and Chang, Sam S and Cookson, Michael S and Smith, Joseph A and Blue, Elizabeth E and Bamshad, Michael J and Smith, Jeffrey R} } @article {756, title = {Congenital myasthenic syndrome with mild intellectual disability caused by a recurrent SLC25A1 variant.}, journal = {Eur J Hum Genet}, volume = {28}, year = {2020}, month = {2020 03}, pages = {373-377}, abstract = {

Congenital myasthenic syndromes (CMS) are a clinically and genetically heterogeneous group of disorders caused by mutations which lead to impaired neuromuscular transmission. SLC25A1 encodes a mitochondrial citrate carrier, associated mainly with the severe neurometabolic disease combined D-2- and L-2-hydroxyglutaric aciduria (D/L-2-HGA). We previously reported a single family with a homozygous missense variant in SLC25A1 with a phenotype restricted to relatively mild CMS with intellectual disability, but to date no additional cases of this CMS subtype had been reported. Here, we performed whole exome sequencing (WES) in three additional and unrelated families presenting with CMS and mild intellectual disability to identify the underlying causative gene. The WES analysis revealed the presence of a homozygous c.740G>A; p.(Arg247Gln) missense SLC25A1 variant, the same SLC25A1 variant as identified in the original family with this phenotype. Electron microscopy of muscle from two cases revealed enlarged and accumulated mitochondria. Haplotype analysis performed in two unrelated families suggested that this variant is a result of recurrent mutation and not a founder effect. This suggests that p.(Arg247Gln) is associated with a relatively mild CMS phenotype with subtle mitochondrial abnormalities, while other variants in this gene cause more severe neurometabolic disease. In conclusion, the p.(Arg247Gln) SLC25A1 variant should be considered in patients presenting with a presynaptic CMS phenotype, particularly with accompanying intellectual disability.

}, keywords = {Adult, Female, Haplotypes, Homozygote, Humans, Intellectual Disability, Male, Mitochondrial Proteins, Muscle, Skeletal, Mutation, Missense, Myasthenic Syndromes, Congenital, Organic Anion Transporters}, issn = {1476-5438}, doi = {10.1038/s41431-019-0506-2}, author = {Balaraju, Sunitha and Topf, Ana and McMacken, Grace and Kumar, Veeramani Preethish and Pechmann, Astrid and Roper, Helen and Vengalil, Seena and Polavarapu, Kiran and Nashi, Saraswati and MAHjan, Niranjan Prakash and Barbosa, Ines A and Deshpande, Charu and Taylor, Robert W and Cossins, Judith and Beeson, David and Laurie, Steven and Kirschner, Janbernd and Horvath, Rita and McFarland, Robert and Nalini, Atchayaram and Lochm{\"u}ller, Hanns} } @article {774, title = {TBX6 missense variants expand the mutational spectrum in a non-Mendelian inheritance disease.}, journal = {Hum Mutat}, volume = {41}, year = {2020}, month = {2020 01}, pages = {182-195}, abstract = {

Congenital scoliosis (CS) is a birth defect with variable clinical and anatomical manifestations due to spinal malformation. The genetic etiology underlying about 10\% of CS cases in the Chinese population is compound inheritance by which the gene dosage is reduced below that of haploinsufficiency. In this genetic model, the trait manifests as a result of the combined effect of a rare variant and common pathogenic variant allele at a locus. From exome sequencing (ES) data of 523 patients in Asia and two patients in Texas, we identified six TBX6 gene-disruptive variants from 11 unrelated CS patients via ES and in vitro functional testing. The in trans mild hypomorphic allele was identified in 10 of the 11 subjects; as anticipated these 10 shared a similar spinal deformity of hemivertebrae. The remaining case has a homozygous variant in TBX6 (c.418C>T) and presents a more severe spinal deformity phenotype. We found decreased transcriptional activity and abnormal cellular localization as the molecular mechanisms for TBX6 missense loss-of-function alleles. Expanding the mutational spectrum of TBX6 pathogenic alleles enabled an increased molecular diagnostic detection rate, provided further evidence for the gene dosage-dependent genetic model underlying CS, and refined clinical classification.

}, keywords = {Alleles, Cell Line, Female, Gene Expression, Genes, Reporter, Genetic Association Studies, Genetic Predisposition to Disease, Genotype, Haplotypes, Humans, Inheritance Patterns, Male, Models, Molecular, Molecular Diagnostic Techniques, Mutation, Missense, Phenotype, Protein Conformation, Radiography, Sequence Analysis, DNA, Spine, Structure-Activity Relationship, T-Box Domain Proteins, Whole Exome Sequencing}, issn = {1098-1004}, doi = {10.1002/humu.23907}, author = {Chen, Weisheng and Lin, Jiachen and Wang, Lianlei and Li, Xiaoxin and Zhao, Sen and Liu, Jiaqi and Akdemir, Zeynep C and Zhao, Yanxue and Du, Renqian and Ye, Yongyu and Song, Xiaofei and Zhang, Yuanqiang and Yan, Zihui and Yang, Xinzhuang and Lin, Mao and Shen, Jianxiong and Wang, Shengru and Gao, Na and Yang, Ying and Liu, Ying and Li, Wenli and Liu, Jia and Zhang, Na and Yang, Xu and Xu, Yuan and Zhang, Jianguo and Delgado, Mauricio R and Posey, Jennifer E and Qiu, Guixing and Rios, Jonathan J and Liu, Pengfei and Wise, Carol A and Zhang, Feng and Wu, Zhihong and Lupski, James R and Wu, Nan} } @article {526, title = {The coexistence of copy number variations (CNVs) and single nucleotide polymorphisms (SNPs) at a locus can result in distorted calculations of the significance in associating SNPs to disease.}, journal = {Hum Genet}, volume = {137}, year = {2018}, month = {2018 Jul}, pages = {553-567}, abstract = {

With the recent advance in genome-wide association studies (GWAS), disease-associated single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) have been extensively reported. Accordingly, the issue of incorrect identification of recombination events that can induce the distortion of multi-allelic or hemizygous variants has received more attention. However, the potential distorted calculation bias or significance of a detected association in a GWAS due to the coexistence of CNVs and SNPs in the same genomic region may remain under-recognized. Here we performed the association study within a congenital scoliosis (CS) cohort whose genetic etiology was recently elucidated as a compound inheritance model, including mostly one rare variant deletion CNV null allele and one common variant non-coding hypomorphic haplotype of the TBX6 gene. We demonstrated that the existence of a deletion in TBX6 led to an overestimation of the contribution of the SNPs on the hypomorphic allele. Furthermore, we generalized a model to explain the calculation bias, or distorted significance calculation for an association study, that can be {\textquoteright}induced{\textquoteright} by CNVs at a locus. Meanwhile, overlapping between the disease-associated SNPs from published GWAS and common CNVs (overlap 10\%) and pathogenic/likely pathogenic CNVs (overlap 99.69\%) was significantly higher than the random distribution (p < 1 {\texttimes} 10 and p = 0.034, respectively), indicating that such co-existence of CNV and SNV alleles might generally influence data interpretation and potential outcomes of a GWAS. We also verified and assessed the influence of colocalizing CNVs to the detection sensitivity of disease-associated SNP variant alleles in another adolescent idiopathic scoliosis (AIS) genome-wide association study. We proposed that detecting co-existent CNVs when evaluating the association signals between SNPs and disease traits could improve genetic model analyses and better integrate GWAS with robust Mendelian principles.

}, keywords = {Adolescent, Congenital Abnormalities, DNA Copy Number Variations, Female, Genetic Predisposition to Disease, Genome, Human, Genome-Wide Association Study, Genomics, Genotype, Haplotypes, Humans, Phenotype, Polymorphism, Single Nucleotide, Scoliosis}, issn = {1432-1203}, doi = {10.1007/s00439-018-1910-3}, author = {Liu, Jiaqi and Zhou, Yangzhong and Liu, Sen and Song, Xiaofei and Yang, Xin-Zhuang and Fan, Yanhui and Chen, Weisheng and Akdemir, Zeynep Coban and Yan, Zihui and Zuo, Yuzhi and Du, Renqian and Liu, Zhenlei and Yuan, Bo and Zhao, Sen and Liu, Gang and Chen, Yixin and Zhao, Yanxue and Lin, Mao and Zhu, Qiankun and Niu, Yuchen and Liu, Pengfei and Ikegawa, Shiro and Song, You-Qiang and Posey, Jennifer E and Qiu, Guixing and Zhang, Feng and Wu, Zhihong and Lupski, James R and Wu, Nan} } @article {628, title = {Dissecting the Causal Mechanism of X-Linked Dystonia-Parkinsonism by Integrating Genome and Transcriptome Assembly.}, journal = {Cell}, volume = {172}, year = {2018}, month = {2018 02 22}, pages = {897-909.e21}, abstract = {

X-linked Dystonia-Parkinsonism (XDP) is a Mendelian~neurodegenerative disease that is endemic to the Philippines and is associated with a founder haplotype. We integrated multiple genome and transcriptome assembly technologies to narrow the causal mutation to the TAF1 locus, which included a SINE-VNTR-Alu (SVA) retrotransposition into intron 32 of the gene. Transcriptome analyses identified decreased expression of the canonical cTAF1 transcript among XDP probands, and de novo assembly~across multiple pluripotent stem-cell-derived neuronal lineages discovered aberrant TAF1 transcription that involved alternative splicing and intron retention (IR) in proximity to the SVA that was anti-correlated with overall TAF1 expression. CRISPR/Cas9 excision of the SVA rescued this XDP-specific~transcriptional signature and normalized TAF1 expression in probands. These data suggest an SVA-mediated aberrant transcriptional mechanism associated with XDP and may provide a roadmap for layered technologies and integrated assembly-based analyses for other unsolved Mendelian disorders.

}, keywords = {Alternative Splicing, Alu Elements, Base Sequence, Cohort Studies, CRISPR-Cas Systems, Dystonic Disorders, Family, Female, Genetic Diseases, X-Linked, Genetic Loci, Genome, Human, Haplotypes, High-Throughput Nucleotide Sequencing, Histone Acetyltransferases, Humans, Induced Pluripotent Stem Cells, Introns, Male, Minisatellite Repeats, Models, Genetic, Nerve Degeneration, Neural Stem Cells, Neurons, RNA, Messenger, Short Interspersed Nucleotide Elements, TATA-Binding Protein Associated Factors, Transcription Factor TFIID, Transcriptome}, issn = {1097-4172}, doi = {10.1016/j.cell.2018.02.011}, author = {Aneichyk, Tatsiana and Hendriks, William T and Yadav, Rachita and Shin, David and Gao, Dadi and Vaine, Christine A and Collins, Ryan L and Domingo, Aloysius and Currall, Benjamin and Stortchevoi, Alexei and Multhaupt-Buell, Trisha and Penney, Ellen B and Cruz, Lilian and Dhakal, Jyotsna and Brand, Harrison and Hanscom, Carrie and Antolik, Caroline and Dy, Marisela and Ragavendran, Ashok and Underwood, Jason and Cantsilieris, Stuart and Munson, Katherine M and Eichler, Evan E and Acu{\~n}a, Patrick and Go, Criscely and Jamora, R Dominic G and Rosales, Raymond L and Church, Deanna M and Williams, Stephen R and Garcia, Sarah and Klein, Christine and M{\"u}ller, Ulrich and Wilhelmsen, Kirk C and Timmers, H T Marc and Sapir, Yechiam and Wainger, Brian J and Henderson, Daniel and Ito, Naoto and Weisenfeld, Neil and Jaffe, David and Sharma, Nutan and Breakefield, Xandra O and Ozelius, Laurie J and Bragg, D Cristopher and Talkowski, Michael E} } @article {490, title = {Variants in EXOSC9 Disrupt the RNA Exosome and Result in Cerebellar Atrophy with Spinal Motor Neuronopathy.}, journal = {Am J Hum Genet}, volume = {102}, year = {2018}, month = {2018 05 03}, pages = {858-873}, abstract = {

The exosome is a conserved multi-protein complex that is essential for correct RNA processing. Recessive variants in exosome components EXOSC3, EXOSC8, and RBM7 cause various constellations of pontocerebellar hypoplasia (PCH), spinal muscular atrophy (SMA), and central nervous system demyelination. Here, we report on four unrelated affected individuals with recessive variants in EXOSC9 and the effect of the variants on the function of the RNA exosome in~vitro in affected individuals{\textquoteright} fibroblasts and skeletal muscle and in~vivo in zebrafish. The clinical presentation was severe, early-onset, progressive SMA-like motor neuronopathy, cerebellar atrophy, and in one affected individual, congenital fractures of the long bones. Three affected individuals of different ethnicity carried the homozygous c.41T>C (p.Leu14Pro) variant, whereas one affected individual was compound heterozygous for c.41T>C (p.Leu14Pro) and c.481C>T (p.Arg161). We detected reduced EXOSC9 in fibroblasts and skeletal muscle and observed a reduction of the whole multi-subunit exosome complex on blue-native polyacrylamide gel electrophoresis. RNA sequencing of fibroblasts and skeletal muscle detected significant >2-fold changes in genes involved in neuronal development and cerebellar and motor neuron degeneration, demonstrating the widespread effect of the variants. Morpholino oligonucleotide knockdown and CRISPR/Cas9-mediated mutagenesis of exosc9 in zebrafish recapitulated aspects of the human phenotype, as they have in other zebrafish models of exosomal disease. Specifically, portions of the cerebellum and hindbrain were absent, and motor neurons failed to develop and migrate properly. In summary, we show that variants in EXOSC9 result in a neurological syndrome combining cerebellar atrophy and spinal motoneuronopathy, thus expanding the list of human exosomopathies.

}, keywords = {Amino Acid Sequence, Animals, Atrophy, Base Sequence, Cerebellum, Child, Preschool, Exosome Multienzyme Ribonuclease Complex, Exosomes, Female, Fibroblasts, Gene Knockdown Techniques, Genetic Variation, Haplotypes, Humans, Infant, Male, Motor Neurons, Muscle, Skeletal, Pedigree, RNA-Binding Proteins, Spinal Cord, Zebrafish}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2018.03.011}, author = {Burns, David T and Donkervoort, Sandra and M{\"u}ller, Juliane S and Knierim, Ellen and Bharucha-Goebel, Diana and Faqeih, Eissa Ali and Bell, Stephanie K and AlFaifi, Abdullah Y and Monies, Dorota and Millan, Francisca and Retterer, Kyle and Dyack, Sarah and MacKay, Sara and Morales-Gonzalez, Susanne and Giunta, Michele and Munro, Benjamin and Hudson, Gavin and Scavina, Mena and Baker, Laura and Massini, Tara C and Lek, Monkol and Hu, Ying and Ezzo, Daniel and Alkuraya, Fowzan S and Kang, Peter B and Griffin, Helen and Foley, A Reghan and Schuelke, Markus and Horvath, Rita and B{\"o}nnemann, Carsten G} } @article {83, title = {Mitral regurgitation as a phenotypic manifestation of nonphotosensitive trichothiodystrophy due to a splice variant in MPLKIP.}, journal = {BMC Med Genet}, volume = {17}, year = {2016}, month = {2016 Feb 16}, pages = {13}, abstract = {

BACKGROUND: 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.

METHODS: 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.

RESULTS: 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.

CONCLUSION: 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.

}, keywords = {Adaptor Proteins, Signal Transducing, Adolescent, Adult, Alleles, Asian Continental Ancestry Group, Child, Cloning, Molecular, Exome, Female, Genetic Linkage, Haplotypes, HEK293 Cells, Homozygote, Humans, Introns, Male, Mitral Valve Insufficiency, Pakistan, Pedigree, Phenotype, Polymorphism, Single Nucleotide, RNA Splicing, Sequence Analysis, DNA, Trichothiodystrophy Syndromes, Young Adult}, issn = {1471-2350}, doi = {10.1186/s12881-016-0275-5}, author = {Shah, Khadim and Ali, Raja Hussain and Ansar, Muhammad and Lee, Kwanghyuk and Chishti, Muhammad Salman and Abbe, Izoduwa and Li, Biao and Smith, Joshua D and Nickerson, Deborah A and Shendure, Jay and Coucke, Paul J and Steyaert, Wouter and Bamshad, Michael J and Santos-Cortez, Regie Lyn P and Leal, Suzanne M and Ahmad, Wasim} } @article {43, title = {PADRE: Pedigree-Aware Distant-Relationship Estimation.}, journal = {Am J Hum Genet}, volume = {99}, year = {2016}, month = {2016 Jul 07}, pages = {154-62}, abstract = {

Accurate estimation of shared ancestry is an important component of many genetic studies; current prediction tools accurately estimate pairwise genetic relationships up to the ninth degree. Pedigree-aware distant-relationship estimation (PADRE) combines relationship likelihoods generated by estimation of recent shared ancestry (ERSA) with likelihoods from family networks reconstructed by pedigree reconstruction and identification of a maximum unrelated set (PRIMUS), improving the power to detect distant relationships between pedigrees. Using PADRE, we estimated relationships from simulated pedigrees and three extended pedigrees, correctly predicting 20\% more fourth- through ninth-degree simulated relationships than when using ERSA alone. By leveraging pedigree information, PADRE can even identify genealogical relationships between individuals who are genetically unrelated. For example, although 95\% of 13(th)-degree relatives are genetically unrelated, in simulations, PADRE correctly predicted 50\% of 13(th)-degree relationships to within one degree of relatedness. The improvement in prediction accuracy was consistent between simulated and actual pedigrees. We also applied PADRE to the HapMap3 CEU samples and report new cryptic relationships and validation of previously described relationships between families. PADRE greatly expands the range of relationships that can be estimated by using genetic data in pedigrees.

}, keywords = {Algorithms, Female, Haplotypes, Humans, Male, Models, Genetic, Pedigree, Reproducibility of Results}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2016.05.020}, author = {Staples, Jeffrey and Witherspoon, David J and Jorde, Lynn B and Nickerson, Deborah A and Below, Jennifer E and Huff, Chad D} } @article {123, title = {PRIMUS: improving pedigree reconstruction using mitochondrial and Y haplotypes.}, journal = {Bioinformatics}, volume = {32}, year = {2016}, month = {2016 Feb 15}, pages = {596-8}, abstract = {

UNLABELLED: PRIMUS is a pedigree reconstruction algorithm that uses estimates of genome-wide identity by descent to reconstruct pedigrees consistent with observed genetic data. However, when genetic data for individuals within a pedigree are missing, often multiple pedigrees can be reconstructed that fit the data. We report a major expansion of PRIMUS that uses mitochondrial (mtDNA) and non-recombining Y chromosome (NRY) haplotypes to eliminate many pedigree structures that are inconsistent with the genetic data. We demonstrate that discordances in mtDNA and NRY haplotypes substantially reduce the number of potential pedigrees, and often lead to the identification of the correct pedigree.

AVAILABILITY AND IMPLEMENTATION: We have implemented PRIMUS updates in PERL and it is available at primus.gs.washington.edu.

}, keywords = {Algorithms, Chromosomes, Human, Y, Computer Simulation, DNA, Mitochondrial, Genetics, Population, Haplotypes, Humans, Linkage Disequilibrium, Pedigree, Software}, issn = {1367-4811}, doi = {10.1093/bioinformatics/btv618}, author = {Staples, Jeffrey and Ekunwe, Lynette and Lange, Ethan and Wilson, James G and Nickerson, Deborah A and Below, Jennifer E} } @article {178, title = {Collapsed haplotype pattern method for linkage analysis of next-generation sequence data.}, journal = {Eur J Hum Genet}, volume = {23}, year = {2015}, month = {2015 Dec}, pages = {1739-43}, abstract = {

Recent advances in next-generation sequencing (NGS) make it possible to directly sequence genomes and exomes of individuals with Mendelian diseases and screen sequence data for causal variants. With the reduction in cost of NGS, DNA samples from entire families can be sequenced and linkage analysis can be performed directly using NGS data. Inspired by {\textquoteright}burden{\textquoteright} tests, which are used for complex trait rare variant association studies, we developed the collapsed haplotype pattern (CHP) method for linkage analysis. Using data from several deafness genes we demonstrate that the CHP method is substantially more powerful than analyzing individual variants. Unlike applying NGS data filtering approaches, the CHP method provides statistical evidence of a gene{\textquoteright}s involvement in disease etiology and is also less likely to exclude causal variants in the presence of phenocopies and/or reduced penetrance. The CHP method was implemented in the SEQLinkage software package, which can perform linkage analysis on NGS data or can generate data compatible with many linkage analysis programs, reviving them for use in NGS era.

}, keywords = {Chromosome Mapping, Genetic Linkage, Haplotypes, Hearing Loss, High-Throughput Nucleotide Sequencing, Humans, Sequence Analysis, DNA, Software}, issn = {1476-5438}, doi = {10.1038/ejhg.2015.64}, author = {Wang, Gao T and Zhang, Di and Li, Biao and Dai, Hang and Leal, Suzanne M} } @article {107, title = {Comparative Genomic Analyses of the Human NPHP1 Locus Reveal Complex Genomic Architecture and Its Regional Evolution in Primates.}, journal = {PLoS Genet}, volume = {11}, year = {2015}, month = {2015 Dec}, pages = {e1005686}, abstract = {

Many loci in the human genome harbor complex genomic structures that can result in susceptibility to genomic rearrangements leading to various genomic disorders. Nephronophthisis 1 (NPHP1, MIM$\#$ 256100) is an autosomal recessive disorder that can be caused by defects of NPHP1; the gene maps within the human 2q13 region where low copy repeats (LCRs) are abundant. Loss of function of NPHP1 is responsible for approximately 85\% of the NPHP1 cases-about 80\% of such individuals carry a large recurrent homozygous NPHP1 deletion that occurs via nonallelic homologous recombination (NAHR) between two flanking directly oriented ~45 kb LCRs. Published data revealed a non-pathogenic inversion polymorphism involving the NPHP1 gene flanked by two inverted ~358 kb LCRs. Using optical mapping and array-comparative genomic hybridization, we identified three potential novel structural variant (SV) haplotypes at the NPHP1 locus that may protect a haploid genome from the NPHP1 deletion. Inter-species comparative genomic analyses among primate genomes revealed massive genomic changes during evolution. The aggregated data suggest that dynamic genomic rearrangements occurred historically within the NPHP1 locus and generated SV haplotypes observed in the human population today, which may confer differential susceptibility to genomic instability and the NPHP1 deletion within a personal genome. Our study documents diverse SV haplotypes at a complex LCR-laden human genomic region. Comparative analyses provide a model for how this complex region arose during primate evolution, and studies among humans suggest that intra-species polymorphism may potentially modulate an individual{\textquoteright}s susceptibility to acquiring disease-associated alleles.

}, keywords = {Adaptor Proteins, Signal Transducing, Alleles, Animals, Comparative Genomic Hybridization, Evolution, Molecular, Gene Dosage, Gene Rearrangement, Genome, Human, Genomic Structural Variation, Haplotypes, Humans, Kidney Diseases, Cystic, Membrane Proteins, Primates}, issn = {1553-7404}, doi = {10.1371/journal.pgen.1005686}, author = {Yuan, Bo and Liu, Pengfei and Gupta, Aditya and Beck, Christine R and Tejomurtula, Anusha and Campbell, Ian M and Gambin, Tomasz and Simmons, Alexandra D and Withers, Marjorie A and Harris, R Alan and Rogers, Jeffrey and Schwartz, David C and Lupski, James R} } @article {152, title = {Rare A2ML1 variants confer susceptibility to otitis media.}, journal = {Nat Genet}, volume = {47}, year = {2015}, month = {2015 Aug}, pages = {917-20}, abstract = {

A duplication variant within the middle ear-specific gene A2ML1 cosegregates with otitis media in an indigenous Filipino pedigree (LOD score = 7.5 at reduced penetrance) and lies within a founder haplotype that is also shared by 3 otitis-prone European-American and Hispanic-American children but is absent in non-otitis-prone children and >62,000 next-generation sequences. We identified seven additional A2ML1 variants in six otitis-prone children. Collectively, our studies support a role for A2ML1 in the pathophysiology of otitis media.

}, keywords = {alpha-Macroglobulins, Animals, Base Sequence, Child, Cochlea, Exome, Family Health, Female, Gene Duplication, Gene Frequency, Genetic Predisposition to Disease, Genotype, Haplotypes, Humans, Male, Mice, Inbred C57BL, Models, Molecular, Otitis Media, Pedigree, Principal Component Analysis, Protein Conformation, Sequence Analysis, DNA}, issn = {1546-1718}, doi = {10.1038/ng.3347}, author = {Santos-Cortez, Regie Lyn P and Chiong, Charlotte M and Reyes-Quintos, Ma Rina T and Tantoco, Ma Leah C and Wang, Xin and Acharya, Anushree and Abbe, Izoduwa and Giese, Arnaud P and Smith, Joshua D and Allen, E Kaitlynn and Li, Biao and Cutiongco-de la Paz, Eva Maria and Garcia, Marieflor Cristy and Llanes, Erasmo Gonzalo D V and Labra, Patrick John and Gloria-Cruz, Teresa Luisa I and Chan, Abner L and Wang, Gao T and Daly, Kathleen A and Shendure, Jay and Bamshad, Michael J and Nickerson, Deborah A and Patel, Janak A and Riazuddin, Saima and Sale, Michele M and Chonmaitree, Tasnee and Ahmed, Zubair M and Abes, Generoso T and Leal, Suzanne M} } @article {279, title = {Rare-variant extensions of the transmission disequilibrium test: application to autism exome sequence data.}, journal = {Am J Hum Genet}, volume = {94}, year = {2014}, month = {2014 Jan 02}, pages = {33-46}, abstract = {

Many population-based rare-variant (RV) association tests, which aggregate variants across a region, have been developed to analyze sequence data. A drawback of analyzing population-based data is that it is difficult to adequately control for population substructure and admixture, and spurious associations can occur. For RVs, this problem can be substantial, because the spectrum of rare variation can differ greatly between populations. A solution is to analyze parent-child trio data, by using the transmission disequilibrium test (TDT), which is robust to population substructure and admixture. We extended the TDT to test for RV associations using four commonly used methods. We demonstrate that for all RV-TDT methods, using proper analysis strategies, type I error is well-controlled even when there are high levels of population substructure or admixture. For trio data, unlike for population-based data, RV allele-counting association methods will lead to inflated type I errors. However type I errors can be properly controlled by obtaining p values empirically through haplotype permutation. The power of the RV-TDT methods was evaluated and compared to the analysis of case-control data with a number of genetic and disease models. The RV-TDT was also used to analyze exome data from 199 Simons Simplex Collection autism trios and an association was observed with variants in ABCA7. Given the problem of adequately controlling for population substructure and admixture in RV association studies and the growing number of sequence-based trio studies, the RV-TDT is extremely beneficial to elucidate the involvement of RVs in the etiology of complex traits.

}, keywords = {Alleles, Autistic Disorder, Computer Simulation, Exome, Gene Frequency, Genetic Association Studies, Genetic Predisposition to Disease, Genetic Variation, Haplotypes, Humans, Linkage Disequilibrium, Models, Genetic, Phenotype, Sequence Analysis, DNA}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2013.11.021}, author = {He, Zongxiao and O{\textquoteright}Roak, Brian J and Smith, Joshua D and Wang, Gao and Hooker, Stanley and Santos-Cortez, Regie Lyn P and Li, Biao and Kan, Mengyuan and Krumm, Nik and Nickerson, Deborah A and Shendure, Jay and Eichler, Evan E and Leal, Suzanne M} } @article {280, title = {Effect of catechol-o-methyltransferase-gene (COMT) variants on experimental and acute postoperative pain in 1,000 women undergoing surgery for breast cancer.}, journal = {Anesthesiology}, volume = {119}, year = {2013}, month = {2013 Dec}, pages = {1422-33}, abstract = {

BACKGROUND: Catechol-O-methyltransferase (COMT) metabolizes catecholamines in different tissues. Polymorphisms in COMT gene can attenuate COMT activity and increase sensitivity to pain. Human studies exploring the effect of COMT polymorphisms on pain sensitivity have mostly included small, heterogeneous samples and have ignored several important single nucleotide polymorphisms (SNPs). This study examines the effect of COMT polymorphisms on experimental and postoperative pain phenotypes in a large ethnically homogeneous female patient cohort.

METHODS: Intensity of cold (+2-4{\textdegree}C) and heat (+48{\textdegree}C) pain and tolerance to cold pain were assessed in 1,000 patients scheduled for breast cancer surgery. Acute postoperative pain and oxycodone requirements were recorded. Twenty-two COMT SNPs were genotyped and their association with six pain phenotypes analyzed with linear regression.

RESULTS: There was no association between any of the tested pain phenotypes and SNP rs4680. The strongest association signals were seen between rs165774 and heat pain intensity as well as rs887200 and cold pain intensity. In both cases, minor allele carriers reported less pain. Neither of these results remained significant after strict multiple testing corrections. When analyzed further, the effect of rs887200 was, however, shown to be significant and consistent throughout the cold pressure test. No evidence of association between the SNPs and postoperative oxycodone consumption was found.

CONCLUSIONS: SNPs rs887200 and rs165774 located in the untranslated regions of the gene had the strongest effects on pain sensitivity. Their effect on pain is described here for the first time. These results should be confirmed in further studies and the potential functional mechanisms of the variants studied.

}, keywords = {Acute Pain, Adult, Aged, Analgesics, Opioid, Anesthesia, Breast Neoplasms, Catechol O-Methyltransferase, Cold Temperature, Databases, Genetic, Female, Follow-Up Studies, Genetic Variation, Genome-Wide Association Study, Genotype, Haplotypes, Hot Temperature, Humans, Middle Aged, Oxycodone, Pain Measurement, Pain, Postoperative, Polymorphism, Single Nucleotide}, issn = {1528-1175}, doi = {10.1097/ALN.0000000000000013}, author = {Kambur, Oleg and Kaunisto, Mari A and Tikkanen, Emmi and Leal, Suzanne M and Ripatti, Samuli and Kalso, Eija A} } @article {299, title = {Mutations in KARS, encoding lysyl-tRNA synthetase, cause autosomal-recessive nonsyndromic hearing impairment DFNB89.}, journal = {Am J Hum Genet}, volume = {93}, year = {2013}, month = {2013 Jul 11}, pages = {132-40}, abstract = {

Previously, DFNB89, a locus associated with autosomal-recessive nonsyndromic hearing impairment (ARNSHI), was mapped to chromosomal region 16q21-q23.2 in three unrelated, consanguineous Pakistani families. Through whole-exome sequencing of a hearing-impaired individual from each family, missense mutations were identified at highly conserved residues of lysyl-tRNA synthetase (KARS): the c.1129G>A (p.Asp377Asn) variant was found in one family, and the c.517T>C (p.Tyr173His) variant was found in the other two families. Both variants were predicted to be damaging by multiple bioinformatics tools. The two variants both segregated with the nonsyndromic-hearing-impairment phenotype within the three families, and neither mutation was identified in ethnically matched controls or within variant databases. Individuals homozygous for KARS mutations had symmetric, severe hearing impairment across all frequencies but did not show evidence of auditory or limb neuropathy. It has been demonstrated that KARS is expressed in hair cells of zebrafish, chickens, and mice. Moreover, KARS has strong localization to the spiral ligament region of the cochlea, as well as to Deiters{\textquoteright} cells, the sulcus epithelium, the basilar membrane, and the surface of the spiral limbus. It is hypothesized that KARS variants affect aminoacylation in inner-ear cells by interfering with binding activity to tRNA or p38 and with tetramer formation. The identification of rare KARS variants in ARNSHI-affected families defines a gene that is associated with ARNSHI.

}, keywords = {Animals, Case-Control Studies, Chickens, Cochlea, Computational Biology, Consanguinity, Female, Genetic Linkage, Hair Cells, Auditory, Inner, Haplotypes, Hearing Loss, Homozygote, Humans, Lysine-tRNA Ligase, Male, Mice, Mutation, Missense, Pedigree, Transfer RNA Aminoacylation, Zebrafish}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2013.05.018}, author = {Santos-Cortez, Regie Lyn P and Lee, Kwanghyuk and Azeem, Zahid and Antonellis, Patrick J and Pollock, Lana M and Khan, Saadullah and Andrade-Elizondo, Paula B and Chiu, Ilene and Adams, Mark D and Basit, Sulman and Smith, Joshua D and Nickerson, Deborah A and McDermott, Brian M and Ahmad, Wasim and Leal, Suzanne M} } @article {317, title = {Digenic inheritance of an SMCHD1 mutation and an FSHD-permissive D4Z4 allele causes facioscapulohumeral muscular dystrophy type 2.}, journal = {Nat Genet}, volume = {44}, year = {2012}, month = {2012 Dec}, pages = {1370-4}, abstract = {

Facioscapulohumeral dystrophy (FSHD) is characterized by chromatin relaxation of the D4Z4 macrosatellite array on chromosome 4 and expression of the D4Z4-encoded DUX4 gene in skeletal muscle. The more common form, autosomal dominant FSHD1, is caused by contraction of the D4Z4 array, whereas the genetic determinants and inheritance of D4Z4 array contraction-independent FSHD2 are unclear. Here, we show that mutations in SMCHD1 (encoding structural maintenance of chromosomes flexible hinge domain containing 1) on chromosome 18 reduce SMCHD1 protein levels and segregate with genome-wide D4Z4 CpG hypomethylation in human kindreds. FSHD2 occurs in individuals who inherited both the SMCHD1 mutation and a normal-sized D4Z4 array on a chromosome 4 haplotype permissive for DUX4 expression. Reducing SMCHD1 levels in skeletal muscle results in D4Z4 contraction-independent DUX4 expression. Our study identifies SMCHD1 as an epigenetic modifier of the D4Z4 metastable epiallele and as a causal genetic determinant of FSHD2 and possibly other human diseases subject to epigenetic regulation.

}, keywords = {Adult, Aged, Chromosomal Proteins, Non-Histone, Chromosomes, Human, Pair 18, CpG Islands, DNA Methylation, Epigenesis, Genetic, Female, Haplotypes, Heredity, Homeodomain Proteins, Humans, Male, Middle Aged, Muscular Dystrophy, Facioscapulohumeral, Mutation}, issn = {1546-1718}, doi = {10.1038/ng.2454}, author = {Lemmers, Richard J L F and Tawil, Rabi and Petek, Lisa M and Balog, Judit and Block, Gregory J and Santen, Gijs W E and Amell, Amanda M and van der Vliet, Patrick J and Almomani, Rowida and Straasheijm, Kirsten R and Krom, Yvonne D and Klooster, Rinse and Sun, Yu and den Dunnen, Johan T and Helmer, Quinta and Donlin-Smith, Colleen M and Padberg, George W and van Engelen, Baziel G M and de Greef, Jessica C and Aartsma-Rus, Annemieke M and Frants, Rune R and de Visser, Marianne and Desnuelle, Claude and Sacconi, Sabrina and Filippova, Galina N and Bakker, Bert and Bamshad, Michael J and Tapscott, Stephen J and Miller, Daniel G and van der Maarel, Silv{\`e}re M} }