Search results

WOS: 000373197500020 ; PubMed ID: 27016271 ; BACKGROUND AND OBJECTIVE: Developmental language disorder (DLD) is a highly prevalent neurodevelopmental disorder associated with negative outcomes in different domains; the etiology of DLD is unknown. To investigate the genetic underpinnings of DLD, we performed genome-wide association and whole exome sequencing studies in a geographically isolated population with a substantially elevated prevalence of the disorder (ie, the AZ sample). METHODS: DNA samples were collected from 359 individuals for the genome-wide association study and from 12 severely affected individuals for whole exome sequencing. Multifaceted phenotypes, representing major domains of expressive language functioning, were derived from collected speech samples. RESULTS: Gene-based analyses revealed a significant association between SETBP1 and complexity of linguistic output (P = 5.47 x 10(-7)). The analysis of exome variants revealed coding sequence variants in 14 genes, most of which play a role in neural development. Targeted enrichment analysis implicated myocyte enhancer factor-2 (MEF2)-regulated genes in DLD in the AZ population. The main findings were successfully replicated in an independent cohort of children at risk for related disorders (n = 37). CONCLUSIONS: MEF2-regulated pathways were identified as potential candidate pathways in the etiology of DLD. Several genes (including the candidate SETBP1 and other MEF2-related genes) seem to jointly influence certain, but not all, facets of the DLD phenotype. Even when genetic and environmental diversity is reduced, DLD is best conceptualized as etiologically complex. Future research should establish whether the signals detected in the AZ population can be replicated in other samples and languages and provide further characterization of the identified pathway. ; National Institute of Health [R01 DC007665, P50 HD052120]; NIH Centers for Mendelian Genomics [5U54HG006504]; National Science Foundation Integrative Graduate Education and Research Traineeship grant [114399]; Government of the Russian Federation [14.Z50.31.0027]; National Institutes of Health (NIH) ; Supported by National Institute of Health grants R01 DC007665 (Dr Grigorenko, Principal Investigator) and P50 HD052120 (Richard Wagner, Principal Investigator), NIH Centers for Mendelian Genomics (5U54HG006504), National Science Foundation Integrative Graduate Education and Research Traineeship grant 114399 (Dr Magnuson, Principal Investigator), and grant 14.Z50.31.0027 from the Government of the Russian Federation (Dr Grigorenko, Principal Investigator). Funded by the National Institutes of Health (NIH).

WOS: 000373197500020 ; PubMed ID: 27016271 ; BACKGROUND AND OBJECTIVE: Developmental language disorder (DLD) is a highly prevalent neurodevelopmental disorder associated with negative outcomes in different domains; the etiology of DLD is unknown. To investigate the genetic underpinnings of DLD, we performed genome-wide association and whole exome sequencing studies in a geographically isolated population with a substantially elevated prevalence of the disorder (ie, the AZ sample). METHODS: DNA samples were collected from 359 individuals for the genome-wide association study and from 12 severely affected individuals for whole exome sequencing. Multifaceted phenotypes, representing major domains of expressive language functioning, were derived from collected speech samples. RESULTS: Gene-based analyses revealed a significant association between SETBP1 and complexity of linguistic output (P = 5.47 x 10(-7)). The analysis of exome variants revealed coding sequence variants in 14 genes, most of which play a role in neural development. Targeted enrichment analysis implicated myocyte enhancer factor-2 (MEF2)-regulated genes in DLD in the AZ population. The main findings were successfully replicated in an independent cohort of children at risk for related disorders (n = 37). CONCLUSIONS: MEF2-regulated pathways were identified as potential candidate pathways in the etiology of DLD. Several genes (including the candidate SETBP1 and other MEF2-related genes) seem to jointly influence certain, but not all, facets of the DLD phenotype. Even when genetic and environmental diversity is reduced, DLD is best conceptualized as etiologically complex. Future research should establish whether the signals detected in the AZ population can be replicated in other samples and languages and provide further characterization of the identified pathway. ; National Institute of Health [R01 DC007665, P50 HD052120]; NIH Centers for Mendelian Genomics [5U54HG006504]; National Science Foundation Integrative Graduate Education and Research Traineeship grant [114399]; Government of the Russian Federation [14.Z50.31.0027]; National Institutes of Health (NIH) ; Supported by National Institute of Health grants R01 DC007665 (Dr Grigorenko, Principal Investigator) and P50 HD052120 (Richard Wagner, Principal Investigator), NIH Centers for Mendelian Genomics (5U54HG006504), National Science Foundation Integrative Graduate Education and Research Traineeship grant 114399 (Dr Magnuson, Principal Investigator), and grant 14.Z50.31.0027 from the Government of the Russian Federation (Dr Grigorenko, Principal Investigator). Funded by the National Institutes of Health (NIH).

Leukocyte telomere length (LTL) is a heritable biomarker of genomic aging. In this study, we perform a genome-wide meta-analysis of LTL by pooling densely genotyped and imputed association results across large-scale European-descent studies including up to 78,592 individuals. We identify 49 genomic regions at a false dicovery rate (FDR) 350,000 UK Biobank participants suggest that genetically shorter telomere length increases the risk of hypothyroidism and decreases the risk of thyroid cancer, lymphoma, and a range of proliferative conditions. Our results replicate previously reported associations with increased risk of coronary artery disease and lower risk for multiple cancer types. Our findings substantially expand current knowledge on genes that regulate LTL and their impact on human health and disease. ; The ENGAGE Project was funded under the European Union Framework 7 – Health Theme (HEALTH-F4-2007- 201413). The InterAct project received funding from the European Union (Integrated Project LSHM-CT-2006-037197 in the Framework Programme 6 of the European Community). The EPIC-CVD study was supported by core funding from the UK Medical Research Council (MR/L003120/1), the British Heart Foundation (RG/13/13/30194; RG/18/13/33946), the European Commission Framework Programme 7 (HEALTH-F2-2012-279233), and the National Institute for Health Research [Cambridge Biomedical Research Centre at the Cambridge University Hospitals NHS Foundation Trust]. C.P.N is funded by the BHF. V.C., C.P.N. and N.J.S. are supported by the NIHR Leicester Cardiovascular Biomedical Research Centre and N.J.S. holds an NIHR Senior Investigator award. Chen Li is support by a 4-year Wellcome Trust PhD Studentship; CL, LAL, NJW are funded by the Medical Research Council (MC_UU_12015/1). NJW is an NIHR Senior Investigator. JD is funded by the National Institute for Health Research [Senior Investigator Award]. Cohort specific and further acknowledgements are given in the Supplemental Data.

As part of the Genes, Environment and Development Initiative, the Minnesota Center for Twin and Family Research (MCTFR) undertook a genome-wide association study, which we describe here. A total of 8,405 research participants, clustered in four-member families, have been successfully genotyped on 527,829 single nucleotide polymorphism (SNP) markers using Illumina's Human660W-Quad array. Quality control screening of samples and markers as well as SNP imputation procedures are described. We also describe methods for ancestry control and how the familial clustering of the MCTFR sample can be accounted for in the analysis using a Rapid Feasible Generalized Least Squares algorithm. The rich longitudinal MCTFR assessments provide numerous opportunities for collaboration.

BACKGROUND: We examined the associations between germline variants and breast cancer mortality using a large meta-analysis of women of European ancestry. METHODS: Meta-analyses included summary estimates based on Cox models of twelve datasets using ~10.4 million variants for 96,661 women with breast cancer and 7697 events (breast cancer-specific deaths). Oestrogen receptor (ER)-specific analyses were based on 64,171 ER-positive (4116) and 16,172 ER-negative (2125) patients. We evaluated the probability of a signal to be a true positive using the Bayesian false discovery probability (BFDP). RESULTS: We did not find any variant associated with breast cancer-specific mortality at P < 5 × 10-8. For ER-positive disease, the most significantly associated variant was chr7:rs4717568 (BFDP = 7%, P = 1.28 × 10-7, hazard ratio [HR] = 0.88, 95% confidence interval [CI] = 0.84-0.92); the closest gene is AUTS2. For ER-negative disease, the most significant variant was chr7:rs67918676 (BFDP = 11%, P = 1.38 × 10-7, HR = 1.27, 95% CI = 1.16-1.39); located within a long intergenic non-coding RNA gene (AC004009.3), close to the HOXA gene cluster. CONCLUSIONS: We uncovered germline variants on chromosome 7 at BFDP < 15% close to genes for which there is biological evidence related to breast cancer outcome. However, the paucity of variants associated with mortality at genome-wide significance underpins the challenge in providing genetic-based individualised prognostic information for breast cancer patients. ; BCAC is funded by Cancer Research UK [C1287/A16563 and C1287/A10118], the European Union's Horizon 2020 Research and Innovation Programme (Grant numbers 634935 and 633784 for BRIDGES and B-CAST, respectively), and by the European Community's Seventh Framework Programme under grant agreement number 223175 (Grant number HEALTH-F2-2009-223175) (COGS).

Breastfeeding has been an important survival trait during human history, though it has long been recognized that individuals differ in their exact breastfeeding behavior. Here our aims were, first, to explore to what extent genetic and environmental influences contributed to the individual differences in breastfeeding behavior; second, to detect possible genetic variants related to breastfeeding; and lastly, to test if the genetic variants associated with breastfeeding have been previously found to be related with breast size. Data were collected from a large community-based cohort of Australian twins, with 3,364 women participating in the twin modelling analyses and 1,521 of them included in the genome-wide association study (GWAS). Monozygotic (MZ) twin correlations (rMZ= 0.52, 95% CI 0.46–0.57) were larger than dizygotic (DZ) twin correlations (rDZ= 0.35, 95% CI 0.25–0.43) and the best-fitting model was the one composed by additive genetics and unique environmental factors, explaining 53% and 47% of the variance in breastfeeding behavior, respectively. No breastfeeding-related genetic variants reached genome-wide significance. The polygenic risk score analyses showed no significant results, suggesting breast size does not influence breastfeeding. This study confers a replication of a previous one exploring the sources of variance of breastfeeding and, to our knowledge, is the first one to conduct a GWAS on breastfeeding and look at the overlap with variants for breast size.

Altres ajuts: The Stroke Genetics Network (SiGN) study was funded by a cooperative agreement grant from the US National Institute of Neurological Disorders and Stroke (NINDS), NIH (U01 NS069208 and R01 NS100178). SAHLSIS was supported by the Swedish Heart and Lung Foundation (HLF-20160316), the Swedish Research Council (K2014-64X-14605-12-5), the Swedish Stroke Association, the Swedish state (under the "Avtal om Läkarutbildning och Medicinsk Forskning, ALF") (ALFGBG-720081). Australian Stroke Genetics Collaboration study was supported by the National Health and Medical Research Council, Australia. Stroke Pharmacogenomics and Genetics group was supported by Invictus plus network, Generation project, and Miguel Servet programme from Instituto de Salud Carlos III, GODs project and Epigenesis project from Marató de TV3 Foundation and Agaur from Generalitat de Catalunya Government. Arne Lindgren was supported by the Swedish Heart and Lung Foundation, Region Skåne, Skåne University Hospital, the Freemasons Lodge of Instruction EOS in Lund, Lund University, the Foundation of Färs & Frosta-one of Spar-banken Skåne's ownership Foundations, and the Swedish Stroke Association. Martin Söderholm was supported by grants from the Swedish Stroke Association, the Foundation of Färs & Frosta-one of Sparbanken Skåne's ownership Foundations, and the Swedish government (under the "Avtal om Läkarutbildning och Medicinsk Forskning, ALF"). Annie Pedersen was supported by grants from the Swedish government (under the "Avtal om Läkarutbildning och Medi-cinsk Forskning, ALF") and the Gothenburg Foundation for Neurological Research. Natalia Rost was in part supported by NIH-NINDS (R01NS086905 and R01NS082285). Daniel Strbian was supported by the Finnish Subsidiary Governmental Fund (VTR). The authors thank NINDS for funding the genotyping of patients included in the SiGN study (U01 NS069208 and R01 NS100178) and Sólveig Grétarsdóttir for genotyping a subsample of the SAHLSIS cohort. ; ObjectiveTo discover common genetic variants ...

1 página-. Resumen del póster presentado en el Seminario celebrado en Tullm, Austria, entre el 8 y el 11 de abril de 2018. ; In current study, GWAS has been performed for the first time for detecting high-resolution QTL for resistance to fumonisin accumulation in maize kernel. A subset of 270 inbred lines from a maize diversity panel (composed of 302 inbred lines) that represents much of the diversity available in public breeding sector around the world was evaluated under inoculation with Fusarium verticillioides. Significant associations between polymorphisms at SNPs in bins 1.07, 1.09, 2.08, 3.02, 3.04, 3.05, 3.06, 3.08, 3.09, 4.02, 4.05, 5.02, 6.07, 7.05, 8.07, 9.03 and fumonisin accumulation (determined by ELISA tests) in maize kernel were detected. At least, four associations should be behind novel QTL because QTL were not found in that bin or in adjacent bins in previous studies (Robertson-Hoyt et al. 2006; Maschietto et al. 2017). These novel QTLs for fumonisin content would be at bins 3.02, 5.02, 7.05, and 8.03. Genes likely involved in ROS production, ROS-scavenging and ROS-detoxification, in signaling and/or regulation of plant hormone-mediated responses and lipid metabolism were proposed as candidate genes behind significant associations. ; This research was funded by the "Plan Estatal de Ciencia y Tecnologia de España" (project AGL2015-67313-C2-1-R, co-financed with European Union funds under the FEDER program) and by Autonomous Government of Galicia, Spain (projects IN607A/013 and ED431F 2016/014). ; Peer reviewed

Background: Osteoporosis is diagnosed by the measurement of bone mineral density, which is a highly heritable and multifactorial trait. We aimed to identify genetic loci that are associated with bone mineral density. Methods: In this genome-wide association study, we identified the most promising of 314 075 single nucleotide polymorphisms (SNPs) in 2094 women in a UK study. We then tested these SNPs for replication in 6463 people from three other cohorts in western Europe. We also investigated allelic expression in lymphoblast cell lines. We tested the association between the replicated SNPs and osteoporotic fractures with data from two studies. Findings: We identified genome-wide evidence for an association between bone mineral density and two SNPs (p<5×10-8). The SNPs were rs4355801, on chromosome 8, near to the TNFRSF11B (osteoprotegerin) gene, and rs3736228, on chromosome 11 in the LRP5 (lipoprotein-receptor-related protein) gene. A non-synonymous SNP in the LRP5 gene was associated with decreased bone mineral density (rs3736228, p=6·3×10-12for lumbar spine and p=1·9×10-4for femoral neck) and an increased risk of both osteoporotic fractures (odds ratio [OR] 1·3, 95% CI 1·09-1·52, p=0·002) and osteoporosis (OR 1·3, 1·08-1·63, p=0·008). Three SNPs near the TNFRSF11B gene were associated with decreased bone mineral density (top SNP, rs4355801: p=7·6×10-10for lumbar spine and p=3·3×10-8for femoral neck) and increased risk of osteoporosis (OR 1·2, 95% CI 1·01-1·42, p=0·038). For carriers of the risk allele at rs4355801, expression of TNFRSF11B in lymphoblast cell lines was halved (p=3·0×10-6). 1883 (22%) of 8557 people were at least heterozygous for these risk alleles, and these alleles had a cumulative association with bone mineral density (trend p=2·3×10-17). The presence of both risk alleles increased the risk of osteoporotic fractures (OR 1·3, 1·08-1·63, p=0·006) and this effect was independent of bone mineral density. Interpretation: Two gene variants of key biological proteins increase the risk of osteoporosis and osteoporotic fracture. The combined effect of these risk alleles on fractures is similar to that of most well-replicated environmental risk factors, and they are present in more than one in five white people, suggesting a potential role in screening. Funding: Wellcome Trust, European Commission, NWO Investments, Arthritis Research Campaign, Chronic Disease Research Foundation, Canadian Institutes of Health Research, European Society for Clinical and Economic Aspects of Osteoporosis, Genome Canada, Genome Quebéc, Canada Research Chairs, National Health and Medical Research Council of Australia, and European Union.

Abstract
Soybean is increasingly becoming one of the main oilseeds crops in Kazakhstan, with a major growing area in the south-eastern part of the country. The comparative assessment of small-size (108) and medium-size (276 accessions) collections of soybean originated from five regions of the world was made in field trials in Almaty region (south-east of Kazakhstan). The full collection was genotyped using the Illumina iSelect 6K SNP (single nucleotide polymorphism) array that allowed observation of 4651 polymorphic SNPs. The genome-wide association study (GWAS) was conducted using the following six important agronomic traits: flowering time, seed maturation time, plant height, number of fertile nodes, number of seeds per node, and thousand seed weight. In total, GWAS using TASSEL package identified 56 stable QTNs (quantitative trait nucleotdies) for all the studied traits, including 42 QTNs detected using the medium-size collection. It was shown that the medium-size population, in comparison with the small-size population, was preferable for the detection of significant associations. A comparative analysis of the data with previously published results suggested that 36 of the stable 56 QTNs were presumably novel associations. The revealed QTNs can be used in the molecular breeding projects for soybean yield improvement.

Mexican Americans (MAs) and American Indians (AIs) constitute conspicuously understudied groups with respect to risk for post-traumatic stress disorder (PTSD), especially in light of findings showing racial/ethnic differences in trauma exposure and risk for PTSD. The purpose of this study was to examine genetic influences on PTSD in two minority cohorts. A genome-wide association study (GWAS) with sum PTSD symptoms for trauma-exposed subjects was run in each cohort. Six highly correlated variants in olfactory receptor family 11 subfamily L member 1 (OR11L1) were suggestively associated with PTSD in the MA cohort. These associations remained suggestively significant after permutation testing. A signal in a nearby olfactory receptor on chromosome 1, olfactory receptor family 2 subfamily L member 13 (OR2L13), tagged by rs151319968, was nominally associated with PTSD in the AI sample. Although no variants were significantly associated after correction for multiple testing in a meta-analysis of the two cohorts, pathway analysis of the top hits showed an enrichment cluster of terms related to sensory transduction, olfactory receptor activity, G-protein coupled receptors, and membrane. As previous studies have proposed a role for olfaction in PTSD, our results indicate this influence may be partially driven by genetic variation in the olfactory system.

Ye-Jin Lee,1 SeungHo Choi,2 Sung-Youn Kwon,2 Yunhwan Lee,2 Jung Kyu Lee,3 Eun Young Heo,3 Hee Soon Chung,3,4 Deog Kyeom Kim3,4 1Division of Pulmonary, Allergy and Critical Care Medicine, Department of Internal Medicine, Hallym University Kangdong Sacred Heart Hospital, Seoul, Korea; 2Department of Internal Medicine, Healthcare Research Institute, Healthcare System Gangnam Center, Seoul National University Hospital, Seoul 135-984 Korea; 3Division of Pulmonary and Critical Care Medicine, Seoul Metropolitan Government-Seoul National University Boramae Medical Center, Seoul, Korea; 4Department of Internal Medicine, Seoul National University College of Medicine, Seoul, KoreaCorrespondence: Deog Kyeom KimDivision of Pulmonary and Critical Care Medicine, Seoul Metropolitan Government-Seoul National University Boramae Medical Center, Boramae-Gil 41, Dongjak-gu, Seoul 156-707, Republic of KoreaTel +82-2-870-3207Fax +82-2-831-0714Email kimdkmd@gmail.comBackground: Identifying the genetic basis of airflow limitation is one of the most interesting issues for understanding chronic obstructive pulmonary disease (COPD) pathophysiology. Several studies have shown that some genetic variants associated with COPD have been identified in genome-wide association study (GWAS), especially in patients with moderate to severe COPD; genetic susceptibility for airflow limitation in the early COPD phase has not been widely studied.Objective: We investigated the genetic variants in early COPD.Methods: The present study analyzed Gene-environment interaction and phenotype (GENIE) cohort that included participants who received health screening examination. The association between single nucleotide polymorphism (SNP) and susceptibility to early COPD (FEV1 predicted ≥ 50% and FEV1/FVC < 0.7) was tested.Results: A total of 130 patients with early COPD and 3478 controls (1700 ever smokers and 1778 never smokers) were recruited. When compared with the total controls, certain SNPs (rs2818103, rs875033, rs9354627, rs34552148) on chromosome 6 were included at the top of our list (p= 5.6 × 10– 7 ∼ 9.6 × 10– 6) although they did not reach genome-wide significance. When compared with the never smoker controls, two SNPs (rs2857210, rs2621419) of the HLA-DQB2 gene class were persistently associated with susceptibility to early COPD.Conclusion: Certain SNPs located on chromosome 6 or the HLA-DQB2 gene were the top-scoring SNPs for the association with susceptibility to early COPD in the Korean GENIE cohort.Keywords: early chronic obstructive pulmonary disease, genome-wide association study, single nucleotide polymorphism, SNP, HLA-DQ gene

AbstractHuman height and body mass index are influenced by a large number of genes, each with small effects, along with environment. To identify common genetic variants associated with these traits, we performed genome-wide association studies in 11,536 individuals composed of Australian twins, family members, and unrelated individuals at ∼550,000 genotyped SNPs. We identified a single genome-wide significant variant for height (Pvalue = 1.06 × 10–9) located inHHIP, a well-replicated height-associated gene. Suggestive levels of association were found for other known genes associated with height (Pvalues < 1 × 10–6):ADAMTSL3,EFEMP1,GPR126, andHMGA2; and BMI (Pvalues < 1 × 10–4):FTOandMC4R. Together, these variants explain less than 2% of total phenotypic variation for height and 0.5% for BMI.

Genome-wide association studies (GWAS) of attention-deficit/hyperactivity disorder (ADHD) offer the benefit of a hypothesis-free approach to measuring the quantitative effect of genetic variants on affection status. Generally the findings of GWAS relying on ADHD status have been non-significant, but the one study using quantitative measures of symptoms found SLC9A9 and SLC6A1 were associated with inattention and hyperactivity–impulsivity. Accordingly, we performed a GWAS using quantitative measures of each ADHD subtype measured with the Strengths and Weaknesses of ADHD and Normal Behaviour (SWAN) scale in two community-based samples. This scale captures the full range of attention and kinetic behavior; from high levels of attention and appropriate activity to the inattention and hyperactivity–impulsivity associated with ADHD within two community-based samples. Our discovery sample comprised 1,851 participants (mean age = 22.8 years [4.8]; 50.6% female), while our replication sample comprised 155 participants (mean age = 26.3 years [3.1]; 68.4% females). Age, sex, age × sex, and age2 were included as covariates and the results from each sample were combined using meta-analysis, then analyzed with a gene-based test to estimate the combined effect of markers within genes. We compare our results with markers that have previously been found to have a strong association with ADHD symptoms. Neither the GWAS nor subsequent meta-analyses yielded genome-wide significant results; the strongest effect was observed at rs2110267 (4.62 × 10−7) for symptoms of hyperactivity–impulsivity. The strongest effect in the gene-based test was for GPR139 on symptoms of inattention (6.40 × 10−5). Replication of this study with larger samples will add to our understanding of the genetic etiology of ADHD.

Background: Systemic lupus erythematosus (SLE) is a common systemic autoimmune disease with a complex genetic inheritance. Genome-wide association studies (GWAS) have significantly increased the number of significant loci associated with SLE risk. To date, however, established loci account for less than 30% of the disease heritability and additional risk variants have yet to be identified. Here we performed a GWAS followed by a meta-analysis to identify new genome-wide significant loci for SLE. Methods: We genotyped a cohort of 907 patients with SLE (cases) and 1524 healthy controls from Spain and performed imputation using the 1000 Genomes reference data. We tested for association using logistic regression with correction for the principal components of variation. Meta-analysis of the association results was subsequently performed on 7,110,321 variants using genetic data from a large cohort of 4036 patients with SLE and 6959 controls of Northern European ancestry. Genetic association was also tested at the pathway level after removing the effect of known risk loci using PASCAL software. Results: We identified five new loci associated with SLE at the genome-wide level of significance (p < 5 × 10? 8): GRB2, SMYD3, ST8SIA4, LAT2 and ARHGAP27. Pathway analysis revealed several biological processes significantly associated with SLE risk: B cell receptor signaling (p = 5.28 × 10? 6), CTLA4 co-stimulation during T cell activation (p = 3.06 × 10? 5), interleukin-4 signaling (p = 3.97 × 10? 5) and cell surface interactions at the vascular wall (p = 4.63 × 10? 5). Conclusions: Our results identify five novel loci for SLE susceptibility, and biologic pathways associated via multiple low-effect-size loci. ; This study was funded by the Spanish Ministry of Economy and Competitiveness (grant numbers: PSE-010000-2006-6 and IPT-010000-2010-36). This work has been also sponsored by the grant SEV-2011-00067 of Severo Ochoa Program, awarded by the Spanish Government. This work was supported by an EFSD/Lilly research fellowship. Josep M. Mercader was supported by Sara Borrell Fellowship from the Instituto Carlos III. Sílvia Bonàs was awarded an FI-DGR Fellowship from FI-DGR 2013 from Agència de Gestió d'Ajuts Universitaris i de Recerca (AGAUR, Generalitat de Catalunya). The study sponsors had no role in the design of the study, collection, analysis or interpretation of the data.