NeMu: a comprehensive pipeline for accurate reconstruction of neutral mutation spectra from evolutionary data

The recognized importance of mutational spectra in molecular evolution is yet to be fully exploited beyond human cancer studies and model organisms. The wealth of intraspecific polymorphism data in the GenBank repository, covering a broad spectrum of genes and species, presents an untapped opportunity for detailed mutational spectrum analysis. Existing methods fall short by ignoring intermediate substitutions on the inner branches of phylogenetic trees and lacking the capability for cross-species mutational comparisons. To address these challenges, we present the NeMu pipeline, available at https://nemu-pipeline.com, a tool grounded in phylogenetic principles designed to provide comprehensive and scalable analysis of mutational spectra. Utilizing extensive sequence data from numerous available genome projects, NeMu rapidly and accurately reconstructs the neutral mutational spectrum. This tool, facilitating the reconstruction of gene- and species-specific mutational spectra, contributes to a deeper understanding of evolutionary mechanisms across the broad spectrum of known species.

Reanalysis of mtDNA mutations of human primordial germ cells (PGCs) reveals NUMT contamination and suggests that selection in PGCs may be positive

The resilience of the mitochondrial genome (mtDNA) to a high mutational pressure depends, in part, on negative purifying selection in the germline. A paradigm in the field has been that such selection, at least in part, takes place in primordial germ cells (PGCs). Specifically, Floros et al. (Nature Cell Biology 20: 144-51) reported an increase in the synonymity of mtDNA mutations (a sign of purifying selection) between early-stage and late-stage PGCs. We re-analyzed Floros' et al. data and determined that their mutational dataset was significantly contaminated with single nucleotide variants (SNVs) derived from a nuclear sequence of mtDNA origin (NUMT) located on chromosome 5. Contamination was caused by co-amplification of the NUMT sequence by cross-specific PCR primers. Importantly, when we removed NUMT-derived SNVs, the evidence of purifying selection was abolished. In addition to bulk PGCs, Floros et al. reported the analysis of single-cell late-stage PGCs, which were amplified with different sets of PCR primers that cannot amplify the NUMT sequence. Accordingly, there were no NUMT-derived SNVs among single PGC mutations. Interestingly, single PGC mutations show adecreaseof synonymity with increased intracellular mutant fraction. More specifically, nonsynonymous mutations show faster intracellular genetic drift towards higher mutant fraction than synonymous ones. This pattern is incompatible with predominantly negative selection. This suggests that germline selection of mtDNA mutations is a complex phenomenon and that the part of this process that takes place in PGCs may be predominantly positive. However counterintuitive, positive germline selection of detrimental mtDNA mutations has been reported previously andpotentially may be evolutionarily advantageous.

Secondary structure of the human mitochondrial genome affects formation of deletions

BACKGROUND

Aging in postmitotic tissues is associated with clonal expansion of somatic mitochondrial deletions, the origin of which is not well understood. Such deletions are often flanked by direct nucleotide repeats, but this alone does not fully explain their distribution. Here, we hypothesized that the close proximity of direct repeats on single-stranded mitochondrial DNA (mtDNA) might play a role in the formation of deletions.

RESULTS

By analyzing human mtDNA deletions in the major arc of mtDNA, which is single-stranded during replication and is characterized by a high number of deletions, we found a non-uniform distribution with a "hot spot" where one deletion breakpoint occurred within the region of 6-9 kb and another within 13-16 kb of the mtDNA. This distribution was not explained by the presence of direct repeats, suggesting that other factors, such as the spatial proximity of these two regions, can be the cause. In silico analyses revealed that the single-stranded major arc may be organized as a large-scale hairpin-like loop with a center close to 11 kb and contacting regions between 6-9 kb and 13-16 kb, which would explain the high deletion activity in this contact zone. The direct repeats located within the contact zone, such as the well-known common repeat with a first arm at 8470-8482 bp (base pair) and a second arm at 13,447-13,459 bp, are three times more likely to cause deletions compared to direct repeats located outside of the contact zone. A comparison of age- and disease-associated deletions demonstrated that the contact zone plays a crucial role in explaining the age-associated deletions, emphasizing its importance in the rate of healthy aging.

CONCLUSIONS

Overall, we provide topological insights into the mechanism of age-associated deletion formation in human mtDNA, which could be used to predict somatic deletion burden and maximum lifespan in different human haplogroups and mammalian species.

Live Birth of a Healthy Child in a Couple with Identical mtDNA Carrying a Pathogenic c.471_477delTTTAAAAinsG Variant in the MOCS2 Gene

Molybdenum cofactor deficiency type B (MOCODB; #252160) is an autosomal recessive metabolic disorder that has only been described in 37 affected patients. In this report, we describe the presence of an in-frame homozygous variant (c.471_477delTTTAAAAinsG) in the MOCS2 gene in an affected child, diagnosed with Ohtahara syndrome according to the clinical manifestations. The analysis of the three-dimensional structure of the protein and the amino acid substitutions suggested the pathogenicity of this mutation. To prevent transmitting this mutation to the next generation, we used preimplantation genetic testing for the monogenic disorders (PGT-M) protocol to select MOCS2 gene mutant-free embryos for transfer in an in vitro fertilization (IVF) program. As a result, a healthy child was born. Interestingly, both parents of the proband shared an identical mitochondrial (mt) DNA control region, assuming their close relationship and thus suggesting that both copies of the nuclear rare variant c.471_477delTTTAAAAinsG may have been transmitted from the same female ancestor. Our estimation of the a priori probability of meeting individuals with the same mtDNA haplotype confirms the assumption of a possible distant maternal relationship among the proband’s direct relatives.

A mitochondria-specific mutational signature of aging: increased rate of A > G substitutions on the heavy strand

The mutational spectrum of the mitochondrial DNA (mtDNA) does not resemble any of the known mutational signatures of the nuclear genome and variation in mtDNA mutational spectra between different organisms is still incomprehensible. Since mitochondria are responsible for aerobic respiration, it is expected that mtDNA mutational spectrum is affected by oxidative damage. Assuming that oxidative damage increases with age, we analyse mtDNA mutagenesis of different species in regards to their generation length. Analysing, (i) dozens of thousands of somatic mtDNA mutations in samples of different ages (ii) 70053 polymorphic synonymous mtDNA substitutions reconstructed in 424 mammalian species with different generation lengths and (iii) synonymous nucleotide content of 650 complete mitochondrial genomes of mammalian species we observed that the frequency of AH > GH substitutions (H: heavy strand notation) is twice bigger in species with high versus low generation length making their mtDNA more AH poor and GH rich. Considering that AH > GH substitutions are also sensitive to the time spent single-stranded (TSSS) during asynchronous mtDNA replication we demonstrated that AH > GH substitution rate is a function of both species-specific generation length and position-specific TSSS. We propose that AH > GH is a mitochondria-specific signature of oxidative damage associated with both aging and TSSS.

Dynamics of the most common pathogenic mtDNA variant m.3243A > G demonstrate frequency-dependency in blood and positive selection in the germline

The A-to-G point mutation at position 3243 in the human mitochondrial genome (m.3243A > G) is the most common pathogenic mtDNA variant responsible for disease in humans. It is widely accepted that m.3243A > G levels decrease in blood with age, and an age correction representing ~ 2% annual decline is often applied to account for this change in mutation level. Here we report that recent data indicate that the dynamics of m.3243A > G are more complex and depend on the mutation level in blood in a bi-phasic way. Consequently, the traditional 2% correction, which is adequate 'on average', creates opposite predictive biases at high and low mutation levels. Unbiased age correction is needed to circumvent these drawbacks of the standard model. We propose to eliminate both biases by using an approach where age correction depends on mutation level in a biphasic way to account for the dynamics of m.3243A > G in blood. The utility of this approach was further tested in estimating germline selection of m.3243A > G. The biphasic approach permitted us to uncover patterns consistent with the possibility of positive selection for m.3243A > G. Germline selection of m.3243A > G shows an 'arching' profile by which selection is positive at intermediate mutant fractions and declines at high and low mutant fractions. We conclude that use of this biphasic approach will greatly improve the accuracy of modelling changes in mtDNA mutation frequencies in the germline and in somatic cells during aging.

Mitochondrial Pseudogenes Suggest Repeated Inter-Species Hybridization among Direct Human Ancestors

The hypothesis that the evolution of humans involves hybridization between diverged species has been actively debated in recent years. We present the following novel evidence in support of this hypothesis: the analysis of nuclear pseudogenes of mtDNA ("NUMTs"). NUMTs are considered "mtDNA fossils" as they preserve sequences of ancient mtDNA and thus carry unique information about ancestral populations. Our comparison of a NUMT sequence shared by humans, chimpanzees, and gorillas with their mtDNAs implies that, around the time of divergence between humans and chimpanzees, our evolutionary history involved the interbreeding of individuals whose mtDNA had diverged as much as ~4.5 Myr prior. This large divergence suggests a distant interspecies hybridization. Additionally, analysis of two other NUMTs suggests that such events occur repeatedly. Our findings suggest a complex pattern of speciation in primate/human ancestors and provide one potential explanation for the mosaic nature of fossil morphology found at the emergence of the hominin lineage. A preliminary version of this manuscript was uploaded to the preprint server BioRxiv in 2017 (10.1101/134502).

Increased prevalence of clonal hematopoiesis of indeterminate potential amongst people living with HIV

People living with human immunodeficiency virus (PLWH) have significantly increased risk for cardiovascular disease in part due to inflammation and immune dysregulation. Clonal hematopoiesis of indeterminate potential (CHIP), the age-related acquisition and expansion of hematopoietic stem cells due to leukemogenic driver mutations, increases risk for both hematologic malignancy and coronary artery disease (CAD). Since increased inflammation is hypothesized to be both a cause and consequence of CHIP, we hypothesized that PLWH have a greater prevalence of CHIP. We searched for CHIP in multi-ethnic cases from the Swiss HIV Cohort Study (SHCS, n = 600) and controls from the Atherosclerosis Risk in the Communities study (ARIC, n = 8111) from blood DNA-derived exome sequences. We observed that HIV is associated with a twofold increase in CHIP prevalence, both in the whole study population and in a subset of 230 cases and 1002 matched controls selected by propensity matching to control for demographic imbalances (SHCS 7%, ARIC 3%, p = 0.005). We also observed that ASXL1 is the most commonly mutated CHIP-associated gene in PLWH. Our results suggest that CHIP may contribute to the excess cardiovascular risk observed in PLWH.

Increased CHIP Prevalence Amongst People Living with HIV

People living with human immunodeficiency virus (PLWH) have significantly increased risk for cardiovascular disease in part due to inflammation and immune dysregulation. Clonal hematopoiesis of indeterminate potential (CHIP), the age-related acquisition and expansion of hematopoietic stem cells due to leukemogenic driver mutations, increases risk for both hematologic malignancy and coronary artery disease (CAD). Since increased inflammation is hypothesized to be both a cause and consequence of CHIP, we hypothesized that PLWH have a greater prevalence of CHIP. We searched for CHIP in multi-ethnic cases from the Swiss HIV Cohort Study (SHCS, n=600) and controls from the Atherosclerosis Risk in the Communities study (ARIC, n=8,111) from blood DNA-derived exome sequences. We observed that HIV is associated with increased CHIP prevalence, both in the whole study population and in a subset of 230 cases and 1002 matched controls selected by propensity matching to control for demographic imbalances (SHCS 7%, ARIC 3%, p=0.005). Additionally, unlike in ARIC, ASXL1 was the most commonly implicated mutated CHIP gene. We propose that CHIP may be one mechanism through which PLWH are at increased risk for CAD. Larger prospective studies should evaluate the hypothesis that CHIP contributes to the excess cardiovascular risk in PLWH.

Molecular design of hypothalamus development

A wealth of specialized neuroendocrine command systems intercalated within the hypothalamus control the most fundamental physiological needs in vertebrates1,2. Nevertheless, we lack a developmental blueprint that integrates the molecular determinants of neuronal and glial diversity along temporal and spatial scales of hypothalamus development3. Here we combine single-cell RNA sequencing of 51,199 mouse cells of ectodermal origin, gene regulatory network (GRN) screens in conjunction with genome-wide association study-based disease phenotyping, and genetic lineage reconstruction to show that nine glial and thirty-three neuronal subtypes are generated by mid-gestation under the control of distinct GRNs. Combinatorial molecular codes that arise from neurotransmitters, neuropeptides and transcription factors are minimally required to decode the taxonomical hierarchy of hypothalamic neurons. The differentiation of γ-aminobutyric acid (GABA) and dopamine neurons, but not glutamate neurons, relies on quasi-stable intermediate states, with a pool of GABA progenitors giving rise to dopamine cells4. We found an unexpected abundance of chemotropic proliferation and guidance cues that are commonly implicated in dorsal (cortical) patterning5 in the hypothalamus. In particular, loss of SLIT-ROBO signalling impaired both the production and positioning of periventricular dopamine neurons. Overall, we identify molecular principles that shape the developmental architecture of the hypothalamus and show how neuronal heterogeneity is transformed into a multimodal neural unit to provide virtually infinite adaptive potential throughout life.

Chromatin three-dimensional interactions mediate genetic effects on gene expression

Studying the genetic basis of gene expression and chromatin organization is key to characterizing the effect of genetic variability on the function and structure of the human genome. Here we unravel how genetic variation perturbs gene regulation using a dataset combining activity of regulatory elements, gene expression, and genetic variants across 317 individuals and two cell types. We show that variability in regulatory activity is structured at the intra- and interchromosomal levels within 12,583 cis-regulatory domains and 30 trans-regulatory hubs that highly reflect the local (that is, topologically associating domains) and global (that is, open and closed chromatin compartments) nuclear chromatin organization. These structures delimit cell type-specific regulatory networks that control gene expression and coexpression and mediate the genetic effects of cis- and trans-acting regulatory variants on genes.

From Normal to Obesity and Back: The Associations between Mitochondrial DNA Copy Number, Gender, and Body Mass Index

Mitochondrial DNA (mtDNA) encodes core subunits of oxidative phosphorylation complexes and, as a result of intricate regulatory crosstalk between nuclear and mitochondrial genomes, the total number of mtDNA copies fits the requirements of each cell type. Deviations from the physiological number of mtDNA copies are expected to be deleterious and might cause some inherited diseases and normal ageing. We studied 46 obese patients with type 2 diabetes (T2DM) one year after a laparoscopic sleeve gastrectomy (LSG) and Roux-en-Y gastric bypass (RYGB). The results were compared with normal-weight patients without T2DM (control group 1) (body mass index (BMI) = 22.5 ± 3.01 kg/m²) and patients with obesity without T2DM (control group 2) (BMI = 36 ± 3.45 kg/m²). We detected an increase of mtDNA copy number in the cells of the buffy coat obtained from peripheral blood, sampled one year after bariatric surgery. We also found that average mtDNA copy number as well as its dynamics (before and after the surgery) are gender-specific. To the best of our knowledge, this is the first evidence for the restoration of mtDNA copy number in obese patients after LSG and RYGB.

Slightly deleterious genomic variants and transcriptome perturbations in Down syndrome embryonic selection

The majority of aneuploid fetuses are spontaneously miscarried. Nevertheless, some aneuploid individuals survive despite the strong genetic insult. Here, we investigate if the survival probability of aneuploid fetuses is affected by the genome-wide burden of slightly deleterious variants. We analyzed two cohorts of live-born Down syndrome individuals (388 genotyped samples and 16 fibroblast transcriptomes) and observed a deficit of slightly deleterious variants on Chromosome 21 and decreased transcriptome-wide variation in the expression level of highly constrained genes. We interpret these results as signatures of embryonic selection, and propose a genetic handicap model whereby an individual bearing an extremely severe deleterious variant (such as aneuploidy) could escape embryonic lethality if the genome-wide burden of slightly deleterious variants is sufficiently low. This approach can be used to study the composition and effect of the numerous slightly deleterious variants in humans and model organisms.

The Immune Signaling Adaptor LAT Contributes to the Neuroanatomical Phenotype of 16p11.2 BP2-BP3 CNVs

Copy-number changes in 16p11.2 contribute significantly to neuropsychiatric traits. Besides the 600 kb BP4-BP5 CNV found in 0.5%-1% of individuals with autism spectrum disorders and schizophrenia and whose rearrangement causes reciprocal defects in head size and body weight, a second distal 220 kb BP2-BP3 CNV is likewise a potent driver of neuropsychiatric, anatomical, and metabolic pathologies. These two CNVs are engaged in complex reciprocal chromatin looping, intimating a functional relationship between genes in these regions that might be relevant to pathomechanism. We assessed the drivers of the distal 16p11.2 duplication by overexpressing each of the nine encompassed genes in zebrafish. Only overexpression of LAT induced a reduction of brain proliferating cells and concomitant microcephaly. Consistently, suppression of the zebrafish ortholog induced an increase of proliferation and macrocephaly. These phenotypes were not unique to zebrafish; Lat knockout mice show brain volumetric changes. Consistent with the hypothesis that LAT dosage is relevant to the CNV pathology, we observed similar effects upon overexpression of CD247 and ZAP70, encoding members of the LAT signalosome. We also evaluated whether LAT was interacting with KCTD13, MVP, and MAPK3, major driver and modifiers of the proximal 16p11.2 600 kb BP4-BP5 syndromes, respectively. Co-injected embryos exhibited an increased microcephaly, suggesting the presence of genetic interaction. Correspondingly, carriers of 1.7 Mb BP1-BP5 rearrangements that encompass both the BP2-BP3 and BP4-BP5 loci showed more severe phenotypes. Taken together, our results suggest that LAT, besides its well-recognized function in T cell development, is a major contributor of the 16p11.2 220 kb BP2-BP3 CNV-associated neurodevelopmental phenotypes.

The genomic landscape of human cellular circadian variation points to a novel role for the signalosome

The importance of natural gene expression variation for human behavior is undisputed, but its impact on circadian physiology remains mostly unexplored. Using umbilical cord fibroblasts, we have determined by genome-wide association how common genetic variation impacts upon cellular circadian function. Gene set enrichment points to differences in protein catabolism as one major source of clock variation in humans. The two most significant alleles regulated expression of COPS7B, a subunit of the COP9 signalosome. We further show that the signalosome complex is imported into the nucleus in timed fashion to stabilize the essential circadian protein BMAL1, a novel mechanism to oppose its proteasome-mediated degradation. Thus, circadian clock properties depend in part upon a genetically-encoded competition between stabilizing and destabilizing forces, and genetic alterations in these mechanisms provide one explanation for human chronotype.

Data on the time of integration of the human mitochondrial pseudogenes (NUMTs) into the nuclear genome

The data and methods presented in this article are supplementing the research article "Integration of mtDNA pseudogenes into the nuclear genome coincides with speciation of the human genus. A hypothesis", DOI: 10.1016/j.mito.2016.12.001 (Gunbin et al., 2017) [1]. Mitochondrial DNA is known to get inserted into nuclear DNA to form NUMTs, i.e. nuclear DNA pseudogenes of the mtDNA. We present here the sequences of selected NUMTs, in which time of integration can be determined with sufficient precision. We report their chromosomal positions , their position within the great ape mtDNA phylogeny, and their times of integration into the nuclear genome. The methods used to generate the data and to control their quality are also presented. The dataset is made publicly available to enable critical or extended analyzes.

Human genome meeting 2016 : Houston, TX, USA. 28 February - 2 March 2016

O1 The metabolomics approach to autism: identification of biomarkers for early detection of autism spectrum disorder A. K. Srivastava, Y. Wang, R. Huang, C. Skinner, T. Thompson, L. Pollard, T. Wood, F. Luo, R. Stevenson O2 Phenome-wide association study for smoking- and drinking-associated genes in 26,394 American women with African, Asian, European, and Hispanic descents R. Polimanti, J. Gelernter O3 Effects of prenatal environment, genotype and DNA methylation on birth weight and subsequent postnatal outcomes: findings from GUSTO, an Asian birth cohort X. Lin, I. Y. Lim, Y. Wu, A. L. Teh, L. Chen, I. M. Aris, S. E. Soh, M. T. Tint, J. L. MacIsaac, F. Yap, K. Kwek, S. M. Saw, M. S. Kobor, M. J. Meaney, K. M. Godfrey, Y. S. Chong, J. D. Holbrook, Y. S. Lee, P. D. Gluckman, N. Karnani, GUSTO study group O4 High-throughput identification of specific qt interval modulating enhancers at the SCN5A locus A. Kapoor, D. Lee, A. Chakravarti O5 Identification of extracellular matrix components inducing cancer cell migration in the supernatant of cultivated mesenchymal stem cells C. Maercker, F. Graf, M. Boutros O6 Single cell allele specific expression (ASE) IN T21 and common trisomies: a novel approach to understand DOWN syndrome and other aneuploidies G. Stamoulis, F. Santoni, P. Makrythanasis, A. Letourneau, M. Guipponi, N. Panousis, M. Garieri, P. Ribaux, E. Falconnet, C. Borel, S. E. Antonarakis O7 Role of microRNA in LCL to IPSC reprogramming S. Kumar, J. Curran, J. Blangero O8 Multiple enhancer variants disrupt gene regulatory network in Hirschsprung disease S. Chatterjee, A. Kapoor, J. Akiyama, D. Auer, C. Berrios, L. Pennacchio, A. Chakravarti O9 Metabolomic profiling for the diagnosis of neurometabolic disorders T. R. Donti, G. Cappuccio, M. Miller, P. Atwal, A. Kennedy, A. Cardon, C. Bacino, L. Emrick, J. Hertecant, F. Baumer, B. Porter, M. Bainbridge, P. Bonnen, B. Graham, R. Sutton, Q. Sun, S. Elsea O10 A novel causal methylation network approach to Alzheimer’s disease Z. Hu, P. Wang, Y. Zhu, J. Zhao, M. Xiong, David A Bennett O11 A microRNA signature identifies subtypes of triple-negative breast cancer and reveals MIR-342-3P as regulator of a lactate metabolic pathway A. Hidalgo-Miranda, S. Romero-Cordoba, S. Rodriguez-Cuevas, R. Rebollar-Vega, E. Tagliabue, M. Iorio, E. D’Ippolito, S. Baroni O12 Transcriptome analysis identifies genes, enhancer RNAs and repetitive elements that are recurrently deregulated across multiple cancer types B. Kaczkowski, Y. Tanaka, H. Kawaji, A. Sandelin, R. Andersson, M. Itoh, T. Lassmann, the FANTOM5 consortium, Y. Hayashizaki, P. Carninci, A. R. R. Forrest O13 Elevated mutation and widespread loss of constraint at regulatory and architectural binding sites across 11 tumour types C. A. Semple O14 Exome sequencing provides evidence of pathogenicity for genes implicated in colorectal cancer E. A. Rosenthal, B. Shirts, L. Amendola, C. Gallego, M. Horike-Pyne, A. Burt, P. Robertson, P. Beyers, C. Nefcy, D. Veenstra, F. Hisama, R. Bennett, M. Dorschner, D. Nickerson, J. Smith, K. Patterson, D. Crosslin, R. Nassir, N. Zubair, T. Harrison, U. Peters, G. Jarvik, NHLBI GO Exome Sequencing Project O15 The tandem duplicator phenotype as a distinct genomic configuration in cancer F. Menghi, K. Inaki, X. Woo, P. Kumar, K. Grzeda, A. Malhotra, H. Kim, D. Ucar, P. Shreckengast, K. Karuturi, J. Keck, J. Chuang, E. T. Liu O16 Modeling genetic interactions associated with molecular subtypes of breast cancer B. Ji, A. Tyler, G. Ananda, G. Carter O17 Recurrent somatic mutation in the MYC associated factor X in brain tumors H. Nikbakht, M. Montagne, M. Zeinieh, A. Harutyunyan, M. Mcconechy, N. Jabado, P. Lavigne, J. Majewski O18 Predictive biomarkers to metastatic pancreatic cancer treatment J. B. Goldstein, M. Overman, G. Varadhachary, R. Shroff, R. Wolff, M. Javle, A. Futreal, D. Fogelman O19 DDIT4 gene expression as a prognostic marker in several malignant tumors L. Bravo, W. Fajardo, H. Gomez, C. Castaneda, C. Rolfo, J. A. Pinto O20 Spatial organization of the genome and genomic alterations in human cancers K. C. Akdemir, L. Chin, A. Futreal, ICGC PCAWG Structural Alterations Group O21 Landscape of targeted therapies in solid tumors S. Patterson, C. Statz, S. Mockus O22 Genomic analysis reveals novel drivers and progression pathways in skin basal cell carcinoma S. N. Nikolaev, X. I. Bonilla, L. Parmentier, B. King, F. Bezrukov, G. Kaya, V. Zoete, V. Seplyarskiy, H. Sharpe, T. McKee, A. Letourneau, P. Ribaux, K. Popadin, N. Basset-Seguin, R. Ben Chaabene, F. Santoni, M. Andrianova, M. Guipponi, M. Garieri, C. Verdan, K. Grosdemange, O. Sumara, M. Eilers, I. Aifantis, O. Michielin, F. de Sauvage, S. Antonarakis O23 Identification of differential biomarkers of hepatocellular carcinoma and cholangiocarcinoma via transcriptome microarray meta-analysis S. Likhitrattanapisal O24 Clinical validity and actionability of multigene tests for hereditary cancers in a large multi-center study S. Lincoln, A. Kurian, A. Desmond, S. Yang, Y. Kobayashi, J. Ford, L. Ellisen O25 Correlation with tumor ploidy status is essential for correct determination of genome-wide copy number changes by SNP array T. L. Peters, K. R. Alvarez, E. F. Hollingsworth, D. H. Lopez-Terrada O26 Nanochannel based next-generation mapping for interrogation of clinically relevant structural variation A. Hastie, Z. Dzakula, A. W. Pang, E. T. Lam, T. Anantharaman, M. Saghbini, H. Cao, BioNano Genomics O27 Mutation spectrum in a pulmonary arterial hypertension (PAH) cohort and identification of associated truncating mutations in TBX4 C. Gonzaga-Jauregui, L. Ma, A. King, E. Berman Rosenzweig, U. Krishnan, J. G. Reid, J. D. Overton, F. Dewey, W. K. Chung O28 NORTH CAROLINA macular dystrophy (MCDR1): mutations found affecting PRDM13 K. Small, A. DeLuca, F. Cremers, R. A. Lewis, V. Puech, B. Bakall, R. Silva-Garcia, K. Rohrschneider, M. Leys, F. S. Shaya, E. Stone O29 PhenoDB and genematcher, solving unsolved whole exome sequencing data N. L. Sobreira, F. Schiettecatte, H. Ling, E. Pugh, D. Witmer, K. Hetrick, P. Zhang, K. Doheny, D. Valle, A. Hamosh O30 Baylor-Johns Hopkins Center for Mendelian genomics: a four year review S. N. Jhangiani, Z. Coban Akdemir, M. N. Bainbridge, W. Charng, W. Wiszniewski, T. Gambin, E. Karaca, Y. Bayram, M. K. Eldomery, J. Posey, H. Doddapaneni, J. Hu, V. R. Sutton, D. M. Muzny, E. A. Boerwinkle, D. Valle, J. R. Lupski, R. A. Gibbs O31 Using read overlap assembly to accurately identify structural genetic differences in an ashkenazi jewish trio S. Shekar, W. Salerno, A. English, A. Mangubat, J. Bruestle O32 Legal interoperability: a sine qua non for international data sharing A. Thorogood, B. M. Knoppers, Global Alliance for Genomics and Health - Regulatory and Ethics Working Group O33 High throughput screening platform of competent sineups: that can enhance translation activities of therapeutic target H. Takahashi, K. R. Nitta, A. Kozhuharova, A. M. Suzuki, H. Sharma, D. Cotella, C. Santoro, S. Zucchelli, S. Gustincich, P. Carninci O34 The undiagnosed diseases network international (UDNI): clinical and laboratory research to meet patient needs J. J. Mulvihill, G. Baynam, W. Gahl, S. C. Groft, K. Kosaki, P. Lasko, B. Melegh, D. Taruscio O36 Performance of computational algorithms in pathogenicity predictions for activating variants in oncogenes versus loss of function mutations in tumor suppressor genes R. Ghosh, S. Plon O37 Identification and electronic health record incorporation of clinically actionable pharmacogenomic variants using prospective targeted sequencing S. Scherer, X. Qin, R. Sanghvi, K. Walker, T. Chiang, D. Muzny, L. Wang, J. Black, E. Boerwinkle, R. Weinshilboum, R. Gibbs O38 Melanoma reprogramming state correlates with response to CTLA-4 blockade in metastatic melanoma T. Karpinets, T. Calderone, K. Wani, X. Yu, C. Creasy, C. Haymaker, M. Forget, V. Nanda, J. Roszik, J. Wargo, L. Haydu, X. Song, A. Lazar, J. Gershenwald, M. Davies, C. Bernatchez, J. Zhang, A. Futreal, S. Woodman O39 Data-driven refinement of complex disease classification from integration of heterogeneous functional genomics data in GeneWeaver E. J. Chesler, T. Reynolds, J. A. Bubier, C. Phillips, M. A. Langston, E. J. Baker O40 A general statistic framework for genome-based disease risk prediction M. Xiong, L. Ma, N. Lin, C. Amos O41 Integrative large-scale causal network analysis of imaging and genomic data and its application in schizophrenia studies N. Lin, P. Wang, Y. Zhu, J. Zhao, V. Calhoun, M. Xiong O42 Big data and NGS data analysis: the cloud to the rescue O. Dobretsberger, M. Egger, F. Leimgruber O43 Cpipe: a convergent clinical exome pipeline specialised for targeted sequencing S. Sadedin, A. Oshlack, Melbourne Genomics Health Alliance O44 A Bayesian classification of biomedical images using feature extraction from deep neural networks implemented on lung cancer data V. A. A. Antonio, N. Ono, Clark Kendrick C. Go O45 MAV-SEQ: an interactive platform for the Management, Analysis, and Visualization of sequence data Z. Ahmed, M. Bolisetty, S. Zeeshan, E. Anguiano, D. Ucar O47 Allele specific enhancer in EPAS1 intronic regions may contribute to high altitude adaptation of Tibetans C. Zeng, J. Shao O48 Nanochannel based next-generation mapping for structural variation detection and comparison in trios and populations H. Cao, A. Hastie, A. W. Pang, E. T. Lam, T. Liang, K. Pham, M. Saghbini, Z. Dzakula O49 Archaic introgression in indigenous populations of Malaysia revealed by whole genome sequencing Y. Chee-Wei, L. Dongsheng, W. Lai-Ping, D. Lian, R. O. Twee Hee, Y. Yunus, F. Aghakhanian, S. S. Mokhtar, C. V. Lok-Yung, J. Bhak, M. Phipps, X. Shuhua, T. Yik-Ying, V. Kumar, H. Boon-Peng O50 Breast and ovarian cancer prevention: is it time for population-based mutation screening of high risk genes? I. Campbell, M.-A. Young, P. James, Lifepool O53 Comprehensive coverage from low DNA input using novel NGS library preparation methods for WGS and WGBS C. Schumacher, S. Sandhu, T. Harkins, V. Makarov O54 Methods for large scale construction of robust PCR-free libraries for sequencing on Illumina HiSeqX platform H. DoddapaneniR. Glenn, Z. Momin, B. Dilrukshi, H. Chao, Q. Meng, B. Gudenkauf, R. Kshitij, J. Jayaseelan, C. Nessner, S. Lee, K. Blankenberg, L. Lewis, J. Hu, Y. Han, H. Dinh, S. Jireh, K. Walker, E. Boerwinkle, D. Muzny, R. Gibbs O55 Rapid capture methods for clinical sequencing J. Hu, K. Walker, C. Buhay, X. Liu, Q. Wang, R. Sanghvi, H. Doddapaneni, Y. Ding, N. Veeraraghavan, Y. Yang, E. Boerwinkle, A. L. Beaudet, C. M. Eng, D. M. Muzny, R. A. Gibbs O56 A diploid personal human genome model for better genomes from diverse sequence data K. C. C. Worley, Y. Liu, D. S. T. Hughes, S. C. Murali, R. A. Harris, A. C. English, X. Qin, O. A. Hampton, P. Larsen, C. Beck, Y. Han, M. Wang, H. Doddapaneni, C. L. Kovar, W. J. Salerno, A. Yoder, S. Richards, J. Rogers, J. R. Lupski, D. M. Muzny, R. A. Gibbs O57 Development of PacBio long range capture for detection of pathogenic structural variants Q. Meng, M. Bainbridge, M. Wang, H. Doddapaneni, Y. Han, D. Muzny, R. Gibbs O58 Rhesus macaques exhibit more non-synonymous variation but greater impact of purifying selection than humans R. A. Harris, M. Raveenedran, C. Xue, M. Dahdouli, L. Cox, G. Fan, B. Ferguson, J. Hovarth, Z. Johnson, S. Kanthaswamy, M. Kubisch, M. Platt, D. Smith, E. Vallender, R. Wiseman, X. Liu, J. Below, D. Muzny, R. Gibbs, F. Yu, J. Rogers O59 Assessing RNA structure disruption induced by single-nucleotide variation J. Lin, Y. Zhang, Z. Ouyang P1 A meta-analysis of genome-wide association studies of mitochondrial dna copy number A. Moore, Z. Wang, J. Hofmann, M. Purdue, R. Stolzenberg-Solomon, S. Weinstein, D. Albanes, C.-S. Liu, W.-L. Cheng, T.-T. Lin, Q. Lan, N. Rothman, S. Berndt P2 Missense polymorphic genetic combinations underlying down syndrome susceptibility E. S. Chen P4 The evaluation of alteration of ELAM-1 expression in the endometriosis patients H. Bahrami, A. Khoshzaban, S. Heidari Keshal P5 Obesity and the incidence of apolipoprotein E polymorphisms in an assorted population from Saudi Arabia population K. K. R. Alharbi P6 Genome-associated personalized antithrombotical therapy for patients with high risk of thrombosis and bleeding M. Zhalbinova, A. Akilzhanova, S. Rakhimova, M. Bekbosynova, S. Myrzakhmetova P7 Frequency of Xmn1 polymorphism among sickle cell carrier cases in UAE population M. Matar P8 Differentiating inflammatory bowel diseases by using genomic data: dimension of the problem and network organization N. Mili, R. Molinari, Y. Ma, S. Guerrier P9 Vulnerability of genetic variants to the risk of autism among Saudi children N. Elhawary, M. Tayeb, N. Bogari, N. Qotb P10 Chromatin profiles from ex vivo purified dopaminergic neurons establish a promising model to support studies of neurological function and dysfunction S. A. McClymont, P. W. Hook, L. A. Goff, A. McCallion P11 Utilization of a sensitized chemical mutagenesis screen to identify genetic modifiers of retinal dysplasia in homozygous Nr2e3^rd7 mice Y. Kong, J. R. Charette, W. L. Hicks, J. K. Naggert, L. Zhao, P. M. Nishina P12 Ion torrent next generation sequencing of recessive polycystic kidney disease in Saudi patients B. M. Edrees, M. Athar, F. A. Al-Allaf, M. M. Taher, W. Khan, A. Bouazzaoui, N. A. Harbi, R. Safar, H. Al-Edressi, A. Anazi, N. Altayeb, M. A. Ahmed, K. Alansary, Z. Abduljaleel P13 Digital expression profiling of Purkinje neurons and dendrites in different subcellular compartments A. Kratz, P. Beguin, S. Poulain, M. Kaneko, C. Takahiko, A. Matsunaga, S. Kato, A. M. Suzuki, N. Bertin, T. Lassmann, R. Vigot, P. Carninci, C. Plessy, T. Launey P14 The evolution of imperfection and imperfection of evolution: the functional and functionless fractions of the human genome D. Graur P16 Species-independent identification of known and novel recurrent genomic entities in multiple cancer patients J. Friis-Nielsen, J. M. Izarzugaza, S. Brunak P18 Discovery of active gene modules which are densely conserved across multiple cancer types reveal their prognostic power and mutually exclusive mutation patterns B. S. Soibam P19 Whole exome sequencing of dysplastic leukoplakia tissue indicates sequential accumulation of somatic mutations from oral precancer to cancer D. Das, N. Biswas, S. Das, S. Sarkar, A. Maitra, C. Panda, P. Majumder P21 Epigenetic mechanisms of carcinogensis by hereditary breast cancer genes J. J. Gruber, N. Jaeger, M. Snyder P22 RNA direct: a novel RNA enrichment strategy applied to transcripts associated with solid tumors K. Patel, S. Bowman, T. Davis, D. Kraushaar, A. Emerman, S. Russello, N. Henig, C. Hendrickson P23 RNA sequencing identifies gene mutations for neuroblastoma K. Zhang P24 Participation of SFRP1 in the modulation of TMPRSS2-ERG fusion gene in prostate cancer cell lines M. Rodriguez-Dorantes, C. D. Cruz-Hernandez, C. D. P. Garcia-Tobilla, S. Solorzano-Rosales P25 Targeted Methylation Sequencing of Prostate Cancer N. Jäger, J. Chen, R. Haile, M. Hitchins, J. D. Brooks, M. Snyder P26 Mutant TPMT alleles in children with acute lymphoblastic leukemia from México City and Yucatán, Mexico S. Jiménez-Morales, M. Ramírez, J. Nuñez, V. Bekker, Y. Leal, E. Jiménez, A. Medina, A. Hidalgo, J. Mejía P28 Genetic modifiers of Alström syndrome J. Naggert, G. B. Collin, K. DeMauro, R. Hanusek, P. M. Nishina P31 Association of genomic variants with the occurrence of angiotensin-converting-enzyme inhibitor (ACEI)-induced coughing among Filipinos E. M. Cutiongco De La Paz, R. Sy, J. Nevado, P. Reganit, L. Santos, J. D. Magno, F. E. Punzalan , D. Ona , E. Llanes, R. L. Santos-Cortes , R. Tiongco, J. Aherrera, L. Abrahan, P. Pagauitan-Alan; Philippine Cardiogenomics Study Group P32 The use of “humanized” mouse models to validate disease association of a de novo GARS variant and to test a novel gene therapy strategy for Charcot-Marie-Tooth disease type 2D K. H. Morelli, J. S. Domire, N. Pyne, S. Harper, R. Burgess P34 Molecular regulation of chondrogenic human induced pluripotent stem cells M. A. Gari, A. Dallol, H. Alsehli, A. Gari, M. Gari, A. Abuzenadah P35 Molecular profiling of hematologic malignancies: implementation of a variant assessment algorithm for next generation sequencing data analysis and clinical reporting M. Thomas, M. Sukhai, S. Garg, M. Misyura, T. Zhang, A. Schuh, T. Stockley, S. Kamel-Reid P36 Accessing genomic evidence for clinical variants at NCBI S. Sherry, C. Xiao, D. Slotta, K. Rodarmer, M. Feolo, M. Kimelman, G. Godynskiy, C. O’Sullivan, E. Yaschenko P37 NGS-SWIFT: a cloud-based variant analysis framework using control-accessed sequencing data from DBGAP/SRA C. Xiao, E. Yaschenko, S. Sherry P38 Computational assessment of drug induced hepatotoxicity through gene expression profiling C. Rangel-Escareño, H. Rueda-Zarate P40 Flowr: robust and efficient pipelines using a simple language-agnostic approach;ultraseq; fast modular pipeline for somatic variation calling using flowr S. Seth, S. Amin, X. Song, X. Mao, H. Sun, R. G. Verhaak, A. Futreal, J. Zhang P41 Applying “Big data” technologies to the rapid analysis of heterogenous large cohort data S. J. Whiite, T. Chiang, A. English, J. Farek, Z. Kahn, W. Salerno, N. Veeraraghavan, E. Boerwinkle, R. Gibbs P42 FANTOM5 web resource for the large-scale genome-wide transcription start site activity profiles of wide-range of mammalian cells T. Kasukawa, M. Lizio, J. Harshbarger, S. Hisashi, J. Severin, A. Imad, S. Sahin, T. C. Freeman, K. Baillie, A. Sandelin, P. Carninci, A. R. R. Forrest, H. Kawaji, The FANTOM Consortium P43 Rapid and scalable typing of structural variants for disease cohorts W. Salerno, A. English, S. N. Shekar, A. Mangubat, J. Bruestle, E. Boerwinkle, R. A. Gibbs P44 Polymorphism of glutathione S-transferases and sulphotransferases genes in an Arab population A. H. Salem, M. Ali, A. Ibrahim, M. Ibrahim P46 Genetic divergence of CYP3A5*3 pharmacogenomic marker for native and admixed Mexican populations J. C. Fernandez-Lopez, V. Bonifaz-Peña, C. Rangel-Escareño, A. Hidalgo-Miranda, A. V. Contreras P47 Whole exome sequence meta-analysis of 13 white blood cell, red blood cell, and platelet traits L. Polfus, CHARGE and NHLBI Exome Sequence Project Working Groups P48 Association of adipoq gene with type 2 diabetes and related phenotypes in african american men and women: The jackson heart study S. Davis, R. Xu, S. Gebeab, P Riestra, A Gaye, R. Khan, J. Wilson, A. Bidulescu P49 Common variants in casr gene are associated with serum calcium levels in koreans S. H. Jung, N. Vinayagamoorthy, S. H. Yim, Y. J. Chung P50 Inference of multiple-wave population admixture by modeling decay of linkage disequilibrium with multiple exponential functions Y. Zhou, S. Xu P51 A Bayesian framework for generalized linear mixed models in genome-wide association studies X. Wang, V. Philip, G. Carter P52 Targeted sequencing approach for the identification of the genetic causes of hereditary hearing impairment A. A. Abuzenadah, M. Gari, R. Turki, A. Dallol P53 Identification of enhancer sequences by ATAC-seq open chromatin profiling A. Uyar, A. Kaygun, S. Zaman, E. Marquez, J. George, D. Ucar P54 Direct enrichment for the rapid preparation of targeted NGS libraries C. L. Hendrickson, A. Emerman, D. Kraushaar, S. Bowman, N. Henig, T. Davis, S. Russello, K. Patel P56 Performance of the Agilent D5000 and High Sensitivity D5000 ScreenTape assays for the Agilent 4200 Tapestation System R. Nitsche, L. Prieto-Lafuente P57 ClinVar: a multi-source archive for variant interpretation M. Landrum, J. Lee, W. Rubinstein, D. Maglott P59 Association of functional variants and protein physical interactions of human MUTY homolog linked with familial adenomatous polyposis and colorectal cancer syndrome Z. Abduljaleel, W. Khan, F. A. Al-Allaf, M. Athar , M. M. Taher, N. Shahzad P60 Modification of the microbiom constitution in the gut using chicken IgY antibodies resulted in a reduction of acute graft-versus-host disease after experimental bone marrow transplantation A. Bouazzaoui, E. Huber, A. Dan, F. A. Al-Allaf, W. Herr, G. Sprotte, J. Köstler, A. Hiergeist, A. Gessner, R. Andreesen, E. Holler P61 Compound heterozygous mutation in the LDLR gene in Saudi patients suffering severe hypercholesterolemia F. Al-Allaf, A. Alashwal, Z. Abduljaleel, M. Taher, A. Bouazzaoui, H. Abalkhail, A. Al-Allaf, R. Bamardadh, M. Athar

Genomic analysis identifies new drivers and progression pathways in skin basal cell carcinoma

Basal cell carcinoma (BCC) of the skin is the most common malignant neoplasm in humans. BCC is primarily driven by the Sonic Hedgehog (Hh) pathway. However, its phenotypic variation remains unexplained. Our genetic profiling of 293 BCCs found the highest mutation rate in cancer (65 mutations/Mb). Eighty-five percent of the BCCs harbored mutations in Hh pathway genes (PTCH1, 73% or SMO, 20% (P = 6.6 × 10(-8)) and SUFU, 8%) and in TP53 (61%). However, 85% of the BCCs also harbored additional driver mutations in other cancer-related genes. We observed recurrent mutations in MYCN (30%), PPP6C (15%), STK19 (10%), LATS1 (8%), ERBB2 (4%), PIK3CA (2%), and NRAS, KRAS or HRAS (2%), and loss-of-function and deleterious missense mutations were present in PTPN14 (23%), RB1 (8%) and FBXW7 (5%). Consistent with the mutational profiles, N-Myc and Hippo-YAP pathway target genes were upregulated. Functional analysis of the mutations in MYCN, PTPN14 and LATS1 suggested their potential relevance in BCC tumorigenesis.

DNA-Methylation Patterns in Trisomy 21 Using Cells from Monozygotic Twins

DNA methylation is essential in mammalian development. We have hypothesized that methylation differences induced by trisomy 21 (T21) contribute to the phenotypic characteristics and heterogeneity in Down syndrome (DS). In order to determine the methylation differences in T21 without interference of the interindividual genomic variation, we have used fetal skin fibroblasts from monozygotic (MZ) twins discordant for T21. We also used skin fibroblasts from MZ twins concordant for T21, normal MZ twins without T21, and unrelated normal and T21 individuals. Reduced Representation Bisulfite Sequencing (RRBS) revealed 35 differentially methylated promoter regions (DMRs) (Absolute methylation differences = 25%, FDR < 0.001) in MZ twins discordant for T21 that have also been observed in comparison between unrelated normal and T21 individuals. The identified DMRs are enriched for genes involved in embryonic organ morphogenesis (FDR = 1.60 e -03) and include genes of the HOXB and HOXD clusters. These DMRs are maintained in iPS cells generated from this twin pair and are correlated with the gene expression changes. We have also observed an increase in DNA methylation level in the T21 methylome compared to the normal euploid methylome. This observation is concordant with the up regulation of DNA methyltransferase enzymes (DNMT3B and DNMT3L) and down regulation of DNA demethylation enzymes (TET2 and TET3) observed in the iPSC of the T21 versus normal twin. Altogether, the results of this study highlight the epigenetic effects of the extra chromosome 21 in T21 on loci outside of this chromosome that are relevant to DS associated phenotypes.

When man got his mtDNA deletions?

Somatic mtDNA mutations and deletions in particular are known to clonally expand within cells, eventually reaching detrimental intracellular concentrations. The possibility that clonal expansion is a slow process taking a lifetime had prompted an idea that founder mutations of mutant clones that cause mitochondrial dysfunction in the aged tissue might have originated early in life. If, conversely, expansion was fast, founder mutations should predominantly originate later in life. This distinction is important: indeed, from which mutations should we protect ourselves – those of early development/childhood or those happening at old age? Recently, high-resolution data describing the distribution of mtDNA deletions have been obtained using a novel, highly efficient method (Taylor et al., 2014). These data have been interpreted as supporting predominantly early origin of founder mutations. Re-analysis of the data implies that the data actually better fit mostly late origin of founders, although more research is clearly needed to resolve the controversy.

Genetic and epigenetic regulation of human lincRNA gene expression

Large intergenic noncoding RNAs (lincRNAs) are still poorly functionally characterized. We analyzed the genetic and epigenetic regulation of human lincRNA expression in the GenCord collection by using three cell types from 195 unrelated European individuals. We detected a considerable number of cis expression quantitative trait loci (cis-eQTLs) and demonstrated that the genetic regulation of lincRNA expression is independent of the regulation of neighboring protein-coding genes. lincRNAs have relatively more cis-eQTLs than do equally expressed protein-coding genes with the same exon number. lincRNA cis-eQTLs are located closer to transcription start sites (TSSs) and their effect sizes are higher than cis-eQTLs found for protein-coding genes, suggesting that lincRNA expression levels are less constrained than that of protein-coding genes. Additionally, lincRNA cis-eQTLs can influence the expression level of nearby protein-coding genes and thus could be considered as QTLs for enhancer activity. Enrichment of expressed lincRNA promoters in enhancer marks provides an additional argument for the involvement of lincRNAs in the regulation of transcription in cis. By investigating the epigenetic regulation of lincRNAs, we observed both positive and negative correlations between DNA methylation and gene expression (expression quantitative trait methylation [eQTMs]), as expected, and found that the landscapes of passive and active roles of DNA methylation in gene regulation are similar to protein-coding genes. However, lincRNA eQTMs are located closer to TSSs than are protein-coding gene eQTMs. These similarities and differences in genetic and epigenetic regulation between lincRNAs and protein-coding genes contribute to the elucidation of potential functions of lincRNAs.

The complex SNP and CNV genetic architecture of the increased risk of congenital heart defects in Down syndrome

Congenital heart defect (CHD) occurs in 40% of Down syndrome (DS) cases. While carrying three copies of chromosome 21 increases the risk for CHD, trisomy 21 itself is not sufficient to cause CHD. Thus, additional genetic variation and/or environmental factors could contribute to the CHD risk. Here we report genomic variations that in concert with trisomy 21, determine the risk for CHD in DS. This case-control GWAS includes 187 DS with CHD (AVSD = 69, ASD = 53, VSD = 65) as cases, and 151 DS without CHD as controls. Chromosome 21–specific association studies revealed rs2832616 and rs1943950 as CHD risk alleles (adjusted genotypic P-values <0.05). These signals were confirmed in a replication cohort of 92 DS-CHD cases and 80 DS-without CHD (nominal P-value 0.0022). Furthermore, CNV analyses using a customized chromosome 21 aCGH of 135K probes in 55 DS-AVSD and 53 DS-without CHD revealed three CNV regions associated with AVSD risk (FDR ≤ 0.05). Two of these regions that are located within the previously identified CHD region on chromosome 21 were further confirmed in a replication study of 49 DS-AVSD and 45 DS- without CHD (FDR ≤ 0.05). One of these CNVs maps near the RIPK4 gene, and the second includes the ZBTB21 (previously ZNF295) gene, highlighting the potential role of these genes in the pathogenesis of CHD in DS. We propose that the genetic architecture of the CHD risk of DS is complex and includes trisomy 21, and SNP and CNV variations in chromosome 21. In addition, a yet-unidentified genetic variation in the rest of the genome may contribute to this complex genetic architecture.

Accumulation of slightly deleterious mutations in mitochondrial protein-coding genes of large versus small mammals

After the effective size of a population, N(e), declines, some slightly deleterious amino acid replacements which were initially suppressed by purifying selection become effectively neutral and can reach fixation. Here we investigate this phenomenon for a set of all 13 mitochondrial protein-coding genes from 110 mammalian species. By using body mass as a proxy for N(e), we show that large mammals (i.e., those with low N(e)) as compared with small ones (in our sample these are, on average, 369.5 kg and 275 g, respectively) have a 43% higher rate of accumulation of nonsynonymous nucleotide substitutions relative to synonymous substitutions, and an 8-40% higher rate of accumulation of radical amino acid substitutions relative to conservative substitutions, depending on the type of amino acid classification. These higher rates result in a 6% greater amino acid dissimilarity between modern species and their most recent reconstructed ancestors in large versus small mammals. Because nonsynonymous substitutions are likely to be more harmful than synonymous substitutions, and radical amino acid substitutions are likely to be more harmful than conservative ones, our results suggest that large mammals experience less efficient purifying selection than small mammals. Furthermore, because in the course of mammalian evolution body size tends to increase and, consequently, N(e) tends to decline, evolution of mammals toward large body size may involve accumulation of slightly deleterious mutations in mitochondrial protein-coding genes, which may contribute to decline or extinction of large mammals.

Purifying selection in mitochondria, free-living and obligate intracellular proteobacteria

BACKGROUND

The effectiveness of elimination of slightly deleterious mutations depends mainly on drift and recombination frequency. Here we analyze the influence of these two factors on the strength of the purifying selection in mitochondrial and proteobacterial orthologous genes taking into account the differences in the organism lifestyles.

RESULTS

(I) We found that the probability of fixation of nonsynonymous substitutions (Kn/Ks) in mitochondria is significantly lower compared to obligate intracellular bacteria and even marginally significantly lower compared to free-living bacteria. The comparison of bacteria of different lifestyles demonstrates more effective elimination of slightly deleterious mutations in (II) free-living bacteria as compared to obligate intracellular species and in (III) obligate intracellular parasites as compared to obligate intracellular symbionts. (IV) Finally, we observed that the level of the purifying selection (i.e. 1-Kn/Ks) increases with the density of mobile elements in bacterial genomes.

CONCLUSION

This study shows that the comparison of patterns of molecular evolution of orthologous genes between ecologically different groups of organisms allow to elucidate the genetic consequences of their various lifestyles. Comparing the strength of the purifying selection among proteobacteria with different lifestyles we obtained results, which are in concordance with theoretical expectations: (II) low effective population size and level of recombination in obligate intracellular proteobacteria lead to less effective elimination of mutations compared to free-living relatives; (III) rare horizontal transmissions, i.e. effectively zero recombination level in symbiotic obligate intracellular bacteria leads to less effective purifying selection than in parasitic obligate intracellular bacteria; (IV) the increased frequency of recombination in bacterial genomes with high mobile element density leads to a more effective elimination of slightly deleterious mutations. At the same time, (I) more effective purifying selection in relatively small populations of nonrecombining mitochondria as compared to large populations of recombining proteobacteria was unexpected. We hypothesize that additional features such as the high number of protein-protein interactions or female germ-cell atresia increase evolutionary constraints and maintain the effective purifying selection in mitochondria, but more work is needed to definitely establish these additional features.