The genetic architecture of autism spectrum disorder involves the interplay of common and rare variants and their impact on hundreds of genes. Using exome sequencing, here we show that analysis of rare coding variation in 3,871 autism cases and 9,937 ancestry-matched or parental controls implicates 22 autosomal genes at a false discovery rate (FDR) < 0.05, plus a set of 107 autosomal genes strongly enriched for those likely to affect risk (FDR < 0.30). These 107 genes, which show unusual evolutionary constraint against mutations, incur de novo loss-of-function mutations in over 5% of autistic subjects. Many of the genes implicated encode proteins for synaptic formation, transcriptional regulation and chromatin-remodelling pathways. These include voltage-gated ion channels regulating the propagation of action potentials, pacemaking and excitability-transcription coupling, as well as histone-modifying enzymes and chromatin remodellers-most prominently those that mediate post-translational lysine methylation/demethylation modifications of histones.

In the Iberian Peninsula, which includes mainly Spain and Portugal, large genomic rearrangements (LGRs) of BRCA1 and BRCA2 have respectively been found in up to 2.33% and 8.4% of families with hereditary breast and/or ovarian cancer (HBOC) that lack point mutations and small indels. In Galicia (Northwest Spain), the spectrum and frequency of BRCA1/BRCA2 point mutations differs from the rest of the Iberian populations. However, to date there are no Galician frequency reports of BRCA1/BRCA2 LGRs. Here we used multiplex ligation-dependent probe amplification (MLPA) to screen 651 Galician index cases (out of the 830 individuals referred for genetic analysis) without point mutations or small indels. We identified three different BRCA1 LGRs in four families. Two of them have been previously classified as pathogenic LGRs: the complete deletion of BRCA1 (identified in two unrelated families) and the deletion of exons 1 to 13. We also identified the duplication of exons 1 and 2 that is a LGR with unknown pathogenicity. Determination of the breakpoints of the BRCA1 LGRs using CNV/SNP arrays and sequencing identified them as NG_005905.2:g.70536_180359del, NG_005905.2:g.90012_97270dup, and NC_000017.10:g.41230935_41399840delinsAluSx1, respectively; previous observations of BRCA1 exon1-24del, exon1-2dup, and exon1-13del LGRs have not characterized them in such detail. All the BRCA1 LGRs arose from unequal homologous recombination events involving Alu elements. We also detected, by sequencing, one BRCA2 LGR, the Portuguese founder mutation c.156_157insAluYa5. The low frequency of BRCA1 LGRs within BRCA1 mutation carriers in Galicia (2.34%, 95% CI: 0.61-7.22) seems to differ from the Spanish population (9.93%, 95% CI: 6.76-14.27, P-value = 0.013) and from the rest of the Iberian population (9.76%, 95% CI: 6.69-13.94, P-value = 0.014).

Determine the prevalence and compare some genetic markers involved in addictive behavior in a group of addicts to derivative of coca (cocaine/crack) or heroin and a control group of non-addicted people matched for gender, age and ethnicity.

Capecitabine is an oral 5-fluorouracil (5-FU) pro-drug commonly used to treat colorectal carcinoma and other tumours. About 35% of patients experience dose-limiting toxicity. The few proven genetic biomarkers of 5-FU toxicity are rare variants and polymorphisms, respectively, at candidate loci dihydropyrimidine dehydrogenase (DPYD) and thymidylate synthase (TYMS).

The genetic background of Argentineans is a mosaic of different continental ancestries. From colonial to present times, the genetic contribution of Europeans and sub-Saharan Africans has superposed to or replaced the indigenous genetic 'stratum'. A sample of 384 individuals representing different Argentinean provinces was collected and genotyped for the first and the second mitochondrial DNA (mtDNA) hypervariable regions, and selectively genotyped for mtDNA SNPs. This data was analyzed together with additional 440 profiles from rural and urban populations plus 304 from Native American Argentineans, all available from the literature. A worldwide database was used for phylogeographic inferences, inter-population comparisons, and admixture analysis. Samples identified as belonging to hg (hg) H2a5 were sequenced for the entire mtDNA genome.

Only a few genetic studies have been carried out to date in Bolivia. However, some of the most important (pre)historical enclaves of South America were located in these territories. Thus, the (sub)-Andean region of Bolivia was part of the Inca Empire, the largest state in Pre-Columbian America. We have genotyped the first hypervariable region (HVS-I) of 720 samples representing the main regions in Bolivia, and these data have been analyzed in the context of other pan-American samples (>19,000 HVS-I mtDNAs). Entire mtDNA genome sequencing was also undertaken on selected Native American lineages. Additionally, a panel of 46 Ancestry Informative Markers (AIMs) was genotyped in a sub-set of samples. The vast majority of the Bolivian mtDNAs (98.4%) were found to belong to the main Native American haplogroups (A: 14.3%, B: 52.6%, C: 21.9%, D: 9.6%), with little indication of sub-Saharan and/or European lineages; however, marked patterns of haplogroup frequencies between main regions exist (e.g. haplogroup B: Andean [71%], Sub-Andean [61%], Llanos [32%]). Analysis of entire genomes unraveled the phylogenetic characteristics of three Native haplogroups: the pan-American haplogroup B2b (originated ∼21.4 thousand years ago [kya]), A2ah (∼5.2 kya), and B2o (∼2.6 kya). The data suggest that B2b could have arisen in North California (an origin even in the north most region of the American continent cannot be disregarded), moved southward following the Pacific coastline and crossed Meso-America. Then, it most likely spread into South America following two routes: the Pacific path towards Peru and Bolivia (arriving here at about ∼15.2 kya), and the Amazonian route of Venezuela and Brazil southwards. In contrast to the mtDNA, Ancestry Informative Markers (AIMs) reveal a higher (although geographically variable) European introgression in Bolivians (25%). Bolivia shows a decreasing autosomal molecular diversity pattern along the longitudinal axis, from the Altiplano to the lowlands. Both autosomes and mtDNA revealed a low impact (1-2%) of a sub-Saharan component in Bolivians.

Between the 15th and 19th centuries ad, the Atlantic slave trade resulted in the forced movement of approximately 13 million people from Africa, mainly to the Americas. Only approximately 11 million survived the passage, and many more died in the early years of captivity. We have studied 481 mitochondrial DNAs (mtDNAs) of recent African ancestry in the Americas and in Eurasia, in an attempt to trace them back to particular regions of Africa. Our results show that mtDNAs in America and Eurasia can, in many cases, be traced to broad geographical regions within Africa, largely in accordance with historical evidence, and raise the possibility that a greater resolution may be possible in the future. However, they also indicate that, at least for the moment, considerable caution is warranted when assessing claims to be able to trace the ancestry of particular lineages to a particular locality within modern-day Africa.

Thyroid cancer is a common malignant disease of the endocrine system with increasing incidence rates over the last few decades. In this study, we sought to analyze the possible association of 45 single nucleotide polymorphisms (SNPs) with thyroid cancer in a population from Rio Grande do Norte, Brazil.

The Iberian Peninsula is a well-delimited geographic region with a rich and complex human history. However, the causes of its genetic structure and past migratory dynamics are not yet fully understood. In order to shed light on them, here we evaluated the gene flow and genetic structure throughout the Iberian Peninsula with spatially explicit modelling applied to a georeferenced genetic dataset composed of genome-wide SNPs from 746 individuals belonging to 17 different regions of the Peninsula. We found contrasting patterns of genetic structure throughout Iberia. In particular, we identified strong patterns of genetic differentiation caused by relevant barriers to gene flow in northern regions and, on the other hand, a large genetic similarity in central and southern regions. In addition, our results showed a preferential north to south migratory dynamics and suggest a sex-biased dispersal in Mediterranean and southern regions. The estimated genetic patterns did not fit with the geographical relief of the Iberian landscape and they rather seem to follow political and linguistic territorial boundaries.

Mosaic loss of chromosome Y (mLOY) is a common ageing-related somatic event and has been previously associated with Alzheimer's disease (AD). However, mLOY estimation from genotype microarray data only reflects the mLOY degree of subjects at the moment of DNA sampling. Therefore, mLOY phenotype associations with AD can be severely age-confounded in the context of genome-wide association studies. Here, we applied Mendelian randomisation to construct an age-independent mLOY polygenic risk score (mloy-PRS) using 114 autosomal variants. The mloy-PRS instrument was associated with an 80% increase in mLOY risk per standard deviation unit (p = 4.22 × 10-20) and was orthogonal with age. We found that a higher genetic risk for mLOY was associated with faster progression to AD in men with mild cognitive impairment (hazard ratio (HR) = 1.23, p = 0.01). Importantly, mloy-PRS had no effect on AD conversion or risk in the female group, suggesting that these associations are caused by the inherent loss of the Y chromosome. Additionally, the blood mLOY phenotype in men was associated with increased cerebrospinal fluid levels of total tau and phosphorylated tau181 in subjects with mild cognitive impairment and dementia. Our results strongly suggest that mLOY is involved in AD pathogenesis.

Autism spectrum disorders (ASD) is a complex neurodevelopmental disorder that may significantly impact on the affected individual's life. Common variation (SNPs) could explain about 50% of ASD heritability. Despite this fact and the large size of the last GWAS meta-analysis, it is believed that hundreds of risk genes in ASD have yet to be discovered. New tools, such as TWAS (Transcriptome Wide Association Studies) which integrate tissue expression and genetic data, are a great approach to identify new ASD susceptibility genes. The main goal of this study is to use UTMOST with the publicly available summary statistics from the largest ASD GWAS meta-analysis as genetic input. In addition, an in silico biological characterization for the novel associated loci was performed. Our results have shown the association of 4 genes at the brain level (CIPC, PINX1, NKX2-2, and PTPRE) and have highlighted the association of NKX2-2, MANBA, ERI1, and MITF at the gastrointestinal level. The gastrointestinal associations are quite relevant given the well-established but unexplored relationship between ASD and gastrointestinal symptoms. Cross-tissue analysis has shown the association of NKX2-2 and BLK. UTMOST-associated genes together with their in silico biological characterization seems to point to different biological mechanisms underlying ASD etiology. Thus, it would not be restricted to brain tissue and it will involve the participation of other body tissues such as the gastrointestinal.

Attention-Deficit Hyperactivity Disorder (ADHD) is a complex neurodevelopmental disorder (NDD) which may significantly impact on the affected individual's life. ADHD is acknowledged to have a high heritability component (70-80%). Recently, a meta-analysis of GWAS (Genome Wide Association Studies) has demonstrated the association of several independent loci. Our main aim here, is to apply PASCAL (pathway scoring algorithm), a new gene-based analysis (GBA) method, to the summary statistics obtained in this meta-analysis. PASCAL will take into account the linkage disequilibrium (LD) across genomic regions in a different way than the most commonly employed GBA methods (MAGMA or VEGAS (Versatile Gene-based Association Study)). In addition to PASCAL analysis a gene network and an enrichment analysis for KEGG and GO terms were carried out. Moreover, GENE2FUNC tool was employed to create gene expression heatmaps and to carry out a (DEG) (Differentially Expressed Gene) analysis using GTEX v7 and BrainSpan data.

Metabolic plasticity is the ability of a biological system to adapt its metabolic phenotype to different environmental stressors. We used a whole-body and tissue-specific phenotypic, functional, proteomic, metabolomic and transcriptomic approach to systematically assess metabolic plasticity in diet-induced obese mice after a combined nutritional and exercise intervention. Although most obesity and overnutrition-related pathological features were successfully reverted, we observed a high degree of metabolic dysfunction in visceral white adipose tissue, characterized by abnormal mitochondrial morphology and functionality. Despite two sequential therapeutic interventions and an apparent global healthy phenotype, obesity triggered a cascade of events in visceral adipose tissue progressing from mitochondrial metabolic and proteostatic alterations to widespread cellular stress, which compromises its biosynthetic and recycling capacity. In humans, weight loss after bariatric surgery showed a transcriptional signature in visceral adipose tissue similar to our mouse model of obesity reversion. Overall, our data indicate that obesity prompts a lasting metabolic fingerprint that leads to a progressive breakdown of metabolic plasticity in visceral adipose tissue.

This study aimed to identify the genetic defect in a multigenerational family presenting an autosomal dominant myopathy with histological features of congenital fiber type disproportion. Linkage analysis and genetic sequencing identified, in all affected members of the family, the c.5807A>G heterozygous mutation in MYH7, which encodes the slow/β-cardiac myosin heavy chain. This mutation causes skeletal but not cardiac involvement. Myosin heavy chain expression pattern was also characterized by immunohistochemistry, western blot and q-PCR in muscle biopsies from two patients aged 25 and 62, respectively. While only congenital fiber type disproportion was observed in the younger patient, older patient's biopsy presented aggregates of slow myosin heavy chains, in fiber sub-sarcolemmal region. These clinico-pathologic findings suggest a novel phenotype within the emerging group of hereditary myosin myopathies, which in this family presents typical characteristics of congenital fiber type disproportion in early stages and later evolves to myosin storage myopathy.

The relationship between variants in SLCO1B1 and SLCO2B1 genes and lipid-lowering response to atorvastatin was investigated.

Colorectal cancer is an important cause of mortality in the developed world. Hereditary forms are due to germ-line mutations in APC, MUTYH, and the mismatch repair genes, but many cases present familial aggregation but an unknown inherited cause. The hypothesis of rare high-penetrance mutations in new genes is a likely explanation for the underlying predisposition in some of these familial cases.

In the last four years, Genome-Wide Association Studies (GWAS) have identified sixteen low-penetrance polymorphisms on fourteen different loci associated with colorectal cancer (CRC). Due to the low risks conferred by known common variants, most of the 35% broad-sense heritability estimated by twin studies remains unexplained. Recently our group performed a case-control study for eight Single Nucleotide Polymorphisms (SNPs) in 4 CRC genes. The present investigation is a follow-up of that study. We have genotyped six SNPs that showed a positive association and carried out a meta-analysis based on eight additional studies comprising in total more than 8000 cases and 6000 controls. The estimated recessive odds ratio for one of the SNPs, rs3219489 (MUTYH Q338H), decreased from 1.52 in the original Swedish study, to 1.18 in the Swedish replication, and to 1.08 in the initial meta-analysis. Since the corresponding summary probability value was 0.06, we decided to retrieve additional information for this polymorphism. The incorporation of six further studies resulted in around 13000 cases and 13000 controls. The newly updated OR was 1.03. The results from the present large, multicenter study illustrate the possibility of decreasing effect sizes with increasing samples sizes. Phenotypic heterogeneity, differential environmental exposures, and population specific linkage disequilibrium patterns may explain the observed difference of genetic effects between Sweden and the other investigated cohorts.

Epigenetics are thought to play a major role in the carcinogenesis of multiple sporadic colorectal cancers (CRC). Previous studies have suggested concordant DNA hypermethylation between tumor pairs. However, only a few methylation markers have been analyzed. This study was aimed at describing the epigenetic signature of multiple CRC using a genome-scale DNA methylation profiling. We analyzed 12 patients with synchronous CRC and 29 age-, sex-, and tumor location-paired patients with solitary tumors from the EPICOLON II cohort. DNA methylation profiling was performed using the Illumina Infinium HM27 DNA methylation assay. The most significant results were validated by Methylight. Tumors samples were also analyzed for the CpG Island Methylator Phenotype (CIMP); KRAS and BRAF mutations and mismatch repair deficiency status. Functional annotation clustering was performed. We identified 102 CpG sites that showed significant DNA hypermethylation in multiple tumors with respect to the solitary counterparts (difference in β value ≥0.1). Methylight assays validated the results for 4 selected genes (p = 0.0002). Eight out of 12(66.6%) multiple tumors were classified as CIMP-high, as compared to 5 out of 29(17.2%) solitary tumors (p = 0.004). Interestingly, 76 out of the 102 (74.5%) hypermethylated CpG sites found in multiple tumors were also seen in CIMP-high tumors. Functional analysis of hypermethylated genes found in multiple tumors showed enrichment of genes involved in different tumorigenic functions. In conclusion, multiple CRC are associated with a distinct methylation phenotype, with a close association between tumor multiplicity and CIMP-high. Our results may be important to unravel the underlying mechanism of tumor multiplicity.

Characterizing the spatial patterns of genetic diversity in human populations has a wide range of applications, from detecting genetic mutations associated with disease to inferring human history. Current approaches, including the widely used principal-component analysis, are not suited for the analysis of linked markers, and local and long-range linkage disequilibrium (LD) can dramatically reduce the accuracy of spatial localization when unaccounted for. To overcome this, we have introduced an approach that performs spatial localization of individuals on the basis of their genetic data and explicitly models LD among markers by using a multivariate normal distribution. By leveraging external reference panels, we derive closed-form solutions to the optimization procedure to achieve a computationally efficient method that can handle large data sets. We validate the method on empirical data from a large sample of European individuals from the POPRES data set, as well as on a large sample of individuals of Spanish ancestry. First, we show that by modeling LD, we achieve accuracy superior to that of existing methods. Importantly, whereas other methods show decreased performance when dense marker panels are used in the inference, our approach improves in accuracy as more markers become available. Second, we show that accurate localization of genetic data can be achieved with only a part of the genome, and this could potentially enable the spatial localization of admixed samples that have a fraction of their genome originating from a given continent. Finally, we demonstrate that our approach is resistant to distortions resulting from long-range LD regions; such distortions can dramatically bias the results when unaccounted for.

We examined the associations between germline variants and breast cancer mortality using a large meta-analysis of women of European ancestry.