Genetic data from the extinct giant rat from Tenerife (Canary Islands) points to a recent divergence from mainland relatives.

Evolution of vertebrate endemics in oceanic islands follows a predictable pattern, known as the island rule, according to which gigantism arises in originally small-sized species and dwarfism in large ones. Species of extinct insular giant rodents are known from all over the world. In the Canary Islands, two examples of giant rats, †Canariomys bravoi and †Canariomys tamarani, endemic to Tenerife and Gran Canaria, respectively, disappeared soon after human settlement. The highly derived morphological features of these insular endemic rodents hamper the reconstruction of their evolutionary histories. We have retrieved partial nuclear and mitochondrial data from †C. bravoi and used this information to explore its evolutionary affinities. The resulting dated phylogeny confidently places †C. bravoi within the African grass rat clade (Arvicanthis niloticus). The estimated divergence time, 650 000 years ago (95% higher posterior densities: 373 000-944 000), points toward an island colonization during the Günz-Mindel interglacial stage. †Canariomys bravoi ancestors would have reached the island via passive rafting and then underwent a yearly increase of mean body mass calculated between 0.0015 g and 0.0023 g; this corresponds to fast evolutionary rates (in darwins (d), ranging from 7.09 d to 2.78 d) that are well above those observed for non-insular mammals.

A Common Genetic Origin for Early Farmers from Mediterranean Cardial and Central European LBK Cultures.

The spread of farming out of the Balkans and into the rest of Europe followed two distinct routes: An initial expansion represented by the Impressa and Cardial traditions, which followed the Northern Mediterranean coastline; and another expansion represented by the LBK (Linearbandkeramik) tradition, which followed the Danube River into Central Europe. Although genomic data now exist from samples representing the second migration, such data have yet to be successfully generated from the initial Mediterranean migration. To address this, we generated the complete genome of a 7,400-year-old Cardial individual (CB13) from Cova Bonica in Vallirana (Barcelona), as well as partial nuclear data from five others excavated from different sites in Spain and Portugal. CB13 clusters with all previously sequenced early European farmers and modern-day Sardinians. Furthermore, our analyses suggest that both Cardial and LBK peoples derived from a common ancient population located in or around the Balkan Peninsula. The Iberian Cardial genome also carries a discernible hunter-gatherer genetic signature that likely was not acquired by admixture with local Iberian foragers. Our results indicate that retrieving ancient genomes from similarly warm Mediterranean environments such as the Near East is technically feasible.

Genetic evidence for patrilocal mating behavior among Neandertal groups.

The remains of 12 Neandertal individuals have been found at the El Sidrón site (Asturias, Spain), consisting of six adults, three adolescents, two juveniles, and one infant. Archaeological, paleontological, and geological evidence indicates that these individuals represent all or part of a contemporaneous social group of Neandertals, who died at around the same time and later were buried together as a result of a collapse of an underground karst. We sequenced phylogenetically informative positions of mtDNA hypervariable regions 1 and 2 from each of the remains. Our results show that the 12 individuals stem from three different maternal lineages, accounting for seven, four, and one individual(s), respectively. Using a Y-chromosome assay to confirm the morphological determination of sex for each individual, we found that, although the three adult males carried the same mtDNA lineage, each of the three adult females carried different mtDNA lineages. These findings provide evidence to indicate that Neandertal groups not only were small and characterized by low genetic diversity but also were likely to have practiced patrilocal mating behavior.

The response to fasting and refeeding reveals functional regulation of lipoprotein lipase proteoforms.

Lipoprotein lipase (LPL) is responsible for the intravascular catabolism of triglyceride-rich lipoproteins and plays a central role in whole-body energy balance and lipid homeostasis. As such, LPL is subject to tissue-specific regulation in different physiological conditions, but the mechanisms of this regulation remain incompletely characterized. Previous work revealed that LPL comprises a set of proteoforms with different isoelectric points, but their regulation and functional significance have not been studied thus far. Here we studied the distribution of LPL proteoforms in different rat tissues and their regulation under physiological conditions. First, analysis by two-dimensional electrophoresis and Western blot showed different patterns of LPL proteoforms (i.e., different pI or relative abundance of LPL proteoforms) in different rat tissues under basal conditions, which could be related to the tissue-specific regulation of the enzyme. Next, the comparison of LPL proteoforms from heart and brown adipose tissue between adults and 15-day-old rat pups, two conditions with minimal regulation of LPL in these tissues, yielded virtually the same tissue-specific patterns of LPL proteoforms. In contrast, the pronounced downregulation of LPL activity observed in white adipose tissue during fasting is accompanied by a prominent reconfiguration of the LPL proteoform pattern. Furthermore, refeeding reverts this downregulation of LPL activity and restores the pattern of LPL proteoforms in this tissue. Importantly, this reversible proteoform-specific regulation during fasting and refeeding indicates that LPL proteoforms are functionally diverse. Further investigation of potential differences in the functional properties of LPL proteoforms showed that all proteoforms exhibit lipolytic activity and have similar heparin-binding affinity, although other functional aspects remain to be investigated. Overall, this study demonstrates the ubiquity, differential distribution and specific regulation of LPL proteoforms in rat tissues and underscores the need to consider the existence of LPL proteoforms for a complete understanding of LPL regulation under physiological conditions.

Using population-based data to evaluate the impact of adherence to endocrine therapy on survival in breast cancer through the web-application BreCanSurvPred.

We show how the use and interpretation of population-based cancer survival indicators can help oncologists talk with breast cancer (BC) patients about the relationship between their prognosis and their adherence to endocrine therapy (ET). The study population comprised a population-based cohort of estrogen receptor positive BC patients (N = 1268) diagnosed in Girona and Tarragona (Northeastern Spain) and classified according to HER2 status (+ / -), stage at diagnosis (I/II/III) and five-year cumulative adherence rate (adherent > 80%; non-adherent ≤ 80%). Cox regression analysis was performed to identify significant prognostic factors for overall survival, whereas relative survival (RS) was used to estimate the crude probability of death due to BC (PBC). Stage and adherence to ET were the significant factors for predicting all-cause mortality. Compared to stage I, risk of death increased in stage II (hazard ratio [HR] 2.24, 95% confidence interval [CI]: 1.51-3.30) and stage III (HR 5.11, 95% CI 3.46-7.51), and it decreased with adherence to ET (HR 0.57, 95% CI 0.41-0.59). PBC differences were higher in non-adherent patients compared to adherent ones and increased across stages: stage I: 6.61% (95% CI 0.05-13.20); stage II: 9.77% (95% CI 0.59-19.01), and stage III: 22.31% (95% CI 6.34-38.45). The age-adjusted survival curves derived from this modeling were implemented in the web application BreCanSurvPred ( https://pdocomputation.snpstats.net/BreCanSurvPred ). Web applications like BreCanSurvPred can help oncologists discuss the consequences of non-adherence to prescribed ET with patients.

Mapping co-ancestry connections between the genome of a Medieval individual and modern Europeans.

Historical genetic links among similar populations can be difficult to establish. Identity by descent (IBD) analyses find genomic blocks that represent direct genealogical relationships among individuals. However, this method has rarely been applied to ancient genomes because IBD stretches are progressively fragmented by recombination and thus not recognizable after few tens of generations. To explore such genealogical relationships, we estimated long IBD blocks among modern Europeans, generating networks to uncover the genetic structures. We found that Basques, Sardinians, Icelanders and Orcadians form, each of them, highly intraconnected sub-clusters in a European network, indicating dense genealogical links within small, isolated populations. We also exposed individual genealogical links -such as the connection between one Basque and one Icelandic individual- that cannot be uncovered with other, widely used population genetics methods such as PCA or ADMIXTURE. Moreover, using ancient DNA technology we sequenced a Late Medieval individual (Barcelona, Spain) to high genomic coverage and identified IBD blocks shared between her and modern Europeans. The Medieval IBD blocks are statistically overrepresented only in modern Spaniards, which is the geographically closest population. This approach can be used to produce a fine-scale reflection of shared ancestry across different populations of the world, offering a direct genetic link from the past to the present.

Missing data imputation and synthetic data simulation through modeling graphical probabilistic dependencies between variables (ModGraProDep): An application to breast cancer survival.

BACKGROUND: Two common issues may arise in certain population-based breast cancer (BC) survival studies: I) missing values in a survivals' predictive variable, such as "Stage" at diagnosis, and II) small sample size due to "imbalance class problem" in certain subsets of patients, demanding data modeling/simulation methods. METHODS: We present a procedure, ModGraProDep, based on graphical modeling (GM) of a dataset to overcome these two issues. The performance of the models derived from ModGraProDep is compared with a set of frequently used classification and machine learning algorithms (Missing Data Problem) and with oversampling algorithms (Synthetic Data Simulation). For the Missing Data Problem we assessed two scenarios: missing completely at random (MCAR) and missing not at random (MNAR). Two validated BC datasets provided by the cancer registries of Girona and Tarragona (northeastern Spain) were used. RESULTS: In both MCAR and MNAR scenarios all models showed poorer prediction performance compared to three GM models: the saturated one (GM.SAT) and two with penalty factors on the partial likelihood (GM.K1 and GM.TEST). However, GM.SAT predictions could lead to non-reliable conclusions in BC survival analysis. Simulation of a "synthetic" dataset derived from GM.SAT could be the worst strategy, but the use of the remaining GMs models could be better than oversampling. CONCLUSION: Our results suggest the use of the GM-procedure presented for one-variable imputation/prediction of missing data and for simulating "synthetic" BC survival datasets. The "synthetic" datasets derived from GMs could be also used in clinical applications of cancer survival data such as predictive risk analysis.

Malaria was a weak selective force in ancient Europeans.

Malaria, caused by Plasmodium parasites, is thought to be one of the strongest selective forces that has shaped the genome of modern humans and was endemic in Europe until recent times. Due to its eradication around mid-twentieth century, the potential selective history of malaria in European populations is largely unknown. Here, we screen 224 ancient European genomes from the Upper Palaeolithic to the post-Roman period for 22 malaria-resistant alleles in twelve genes described in the literature. None of the most specific mutations for malaria resistance, like those at G6PD, HBB or Duffy blood group, have been detected among the available samples, while many other malaria-resistant alleles existed well before the advent of agriculture. We detected statistically significant differences between ancient and modern populations for the ATP2B4, FCGR2B and ABO genes and we found evidence of selection at IL-10 and ATP2B4 genes. However it is unclear whether malaria is the causative agent, because these genes are also involved in other immunological challenges. These results suggest that the selective force represented by malaria was relatively weak in Europe, a fact that could be associated to a recent historical introduction of the severe malaria pathogen.

WISCOD: a statistical web-enabled tool for the identification of significant protein coding regions.

Classically, gene prediction programs are based on detecting signals such as boundary sites (splice sites, starts, and stops) and coding regions in the DNA sequence in order to build potential exons and join them into a gene structure. Although nowadays it is possible to improve their performance with additional information from related species or/and cDNA databases, further improvement at any step could help to obtain better predictions. Here, we present WISCOD, a web-enabled tool for the identification of significant protein coding regions, a novel software tool that tackles the exon prediction problem in eukaryotic genomes. WISCOD has the capacity to detect real exons from large lists of potential exons, and it provides an easy way to use global P value called expected probability of being a false exon (EPFE) that is useful for ranking potential exons in a probabilistic framework, without additional computational costs. The advantage of our approach is that it significantly increases the specificity and sensitivity (both between 80% and 90%) in comparison to other ab initio methods (where they are in the range of 70-75%). WISCOD is written in JAVA and R and is available to download and to run in a local mode on Linux and Windows platforms.

Mitochondrial DNA from El Mirador cave (Atapuerca, Spain) reveals the heterogeneity of Chalcolithic populations.

Previous mitochondrial DNA analyses on ancient European remains have suggested that the current distribution of haplogroup H was modeled by the expansion of the Bell Beaker culture (ca 4,500-4,050 years BP) out of Iberia during the Chalcolithic period. However, little is known on the genetic composition of contemporaneous Iberian populations that do not carry the archaeological tool kit defining this culture. Here we have retrieved mitochondrial DNA (mtDNA) sequences from 19 individuals from a Chalcolithic sample from El Mirador cave in Spain, dated to 4,760-4,200 years BP and we have analyzed the haplogroup composition in the context of modern and ancient populations. Regarding extant African, Asian and European populations, El Mirador shows affinities with Near Eastern groups. In different analyses with other ancient samples, El Mirador clusters with Middle and Late Neolithic populations from Germany, belonging to the Rössen, the Salzmünde and the Baalberge archaeological cultures but not with contemporaneous Bell Beakers. Our analyses support the existence of a common genetic signal between Western and Central Europe during the Middle and Late Neolithic and points to a heterogeneous genetic landscape among Chalcolithic groups.

Genetic analysis of high bone mass cases from the BARCOS cohort of Spanish postmenopausal women.

The aims of the study were to establish the prevalence of high bone mass (HBM) in a cohort of Spanish postmenopausal women (BARCOS) and to assess the contribution of LRP5 and DKK1 mutations and of common bone mineral density (BMD) variants to a HBM phenotype. Furthermore, we describe the expression of several osteoblast-specific and Wnt-pathway genes in primary osteoblasts from two HBM cases. A 0.6% of individuals (10/1600) displayed Z-scores in the HBM range (sum Z-score >4). While no mutation in the relevant exons of LRP5 was detected, a rare missense change in DKK1 was found (p.Y74F), which cosegregated with the phenotype in a small pedigree. Fifty-five BMD SNPs from Estrada et al. [NatGenet 44:491-501,2012] were genotyped in the HBM cases to obtain risk scores for each individual. In this small group of samples, Z-scores were found inversely related to risk scores, suggestive of a polygenic etiology. There was a single exception, which may be explained by a rare penetrant genetic variant, counterbalancing the additive effect of the risk alleles. The expression analysis in primary osteoblasts from two HBM cases and five controls suggested that IL6R, DLX3, TWIST1 and PPARG are negatively related to Z-score. One HBM case presented with high levels of RUNX2, while the other displayed very low SOX6. In conclusion, we provide evidence of lack of LRP5 mutations and of a putative HBM-causing mutation in DKK1. Additionally, we present SNP genotyping and expression results that suggest additive effects of several genes for HBM.

Genomic analysis of the blood attributed to Louis XVI (1754-1793), king of France.

A pyrographically decorated gourd, dated to the French Revolution period, has been alleged to contain a handkerchief dipped into the blood of the French king Louis XVI (1754-1793) after his beheading but recent analyses of living males from two Bourbon branches cast doubts on its authenticity. We sequenced the complete genome of the DNA contained in the gourd at low coverage (~2.5×) with coding sequences enriched at a higher ~7.3× coverage. We found that the ancestry of the gourd's genome does not seem compatible with Louis XVI's known ancestry. From a functional perspective, we did not find an excess of alleles contributing to height despite being described as the tallest person in Court. In addition, the eye colour prediction supported brown eyes, while Louis XVI had blue eyes. This is the first draft genome generated from a person who lived in a recent historical period; however, our results suggest that this sample may not correspond to the alleged king.

An integrative computational analysis provides evidence for FBN1-associated network deregulation in trisomy 21.

Although approximately 50% of Down Syndrome (DS) patients have heart abnormalities, they exhibit an overprotection against cardiac abnormalities related with the connective tissue, for example a lower risk of coronary artery disease. A recent study reported a case of a person affected by DS who carried mutations in FBN1, the gene causative for a connective tissue disorder called Marfan Syndrome (MFS). The fact that the person did not have any cardiac alterations suggested compensation effects due to DS. This observation is supported by a previous DS meta-analysis at the molecular level where we have found an overall upregulation of FBN1 (which is usually downregulated in MFS). Additionally, that result was cross-validated with independent expression data from DS heart tissue. The aim of this work is to elucidate the role of FBN1 in DS and to establish a molecular link to MFS and MFS-related syndromes using a computational approach. To reach that, we conducted different analytical approaches over two DS studies (our previous meta-analysis and independent expression data from DS heart tissue) and revealed expression alterations in the FBN1 interaction network, in FBN1 co-expressed genes and FBN1-related pathways. After merging the significant results from different datasets with a Bayesian approach, we prioritized 85 genes that were able to distinguish control from DS cases. We further found evidence for several of these genes (47%), such as FBN1, DCN, and COL1A2, being dysregulated in MFS and MFS-related diseases. Consequently, we further encourage the scientific community to take into account FBN1 and its related network for the study of DS cardiovascular characteristics.

North African populations carry the signature of admixture with Neandertals.

One of the main findings derived from the analysis of the Neandertal genome was the evidence for admixture between Neandertals and non-African modern humans. An alternative scenario is that the ancestral population of non-Africans was closer to Neandertals than to Africans because of ancient population substructure. Thus, the study of North African populations is crucial for testing both hypotheses. We analyzed a total of 780,000 SNPs in 125 individuals representing seven different North African locations and searched for their ancestral/derived state in comparison to different human populations and Neandertals. We found that North African populations have a significant excess of derived alleles shared with Neandertals, when compared to sub-Saharan Africans. This excess is similar to that found in non-African humans, a fact that can be interpreted as a sign of Neandertal admixture. Furthermore, the Neandertal's genetic signal is higher in populations with a local, pre-Neolithic North African ancestry. Therefore, the detected ancient admixture is not due to recent Near Eastern or European migrations. Sub-Saharan populations are the only ones not affected by the admixture event with Neandertals.

Three strategies to stabilise nearly monodispersed silver nanoparticles in aqueous solution.

Silver nanoparticles are extensively used due to their chemical and physical properties and promising applications in areas such as medicine and electronics. Controlled synthesis of silver nanoparticles remains a major challenge due to the difficulty in producing long-term stable particles of the same size and shape in aqueous solution. To address this problem, we examine three strategies to stabilise aqueous solutions of 15 nm citrate-reduced silver nanoparticles using organic polymeric capping, bimetallic core-shell and bimetallic alloying. Our results show that these strategies drastically improve nanoparticle stability by distinct mechanisms. Additionally, we report a new role of polymer functionalisation in preventing further uncontrolled nanoparticle growth. For bimetallic nanoparticles, we attribute the presence of a higher valence metal on the surface of the nanoparticle as one of the key factors for improving their long-term stability. Stable silver-based nanoparticles, free of organic solvents, will have great potential for accelerating further environmental and nanotoxicity studies.PACS: 81.07.-b; 81.16.Be; 82.70.Dd.

COL1A1 haplotypes and hip fracture.

Fragility fractures resulting from low-trauma events such as a fall from standing height are associated with osteoporosis and are very common in older people, especially women. Three single nucleotide polymorphisms (SNPs) at the COL1A1 gene (rs1107946, rs11327935, and rs1800012) have been widely studied and previously associated with bone mineral density (BMD) and fracture. A rare haplotype (T-delT-T) of these three SNPs was found to be greatly overrepresented in fractured individuals compared with nonfractured controls, thus becoming a good candidate for predicting increased fracture risk. The aim of our study was to assess the association of this haplotype with fracture risk in Spanish individuals. We recruited two independent groups of ∼100 patients with hip fracture (a total of 203 individuals) and compared the genotype and haplotype distributions of the three SNPs in the fractured patients with those of 397 control individuals from the BARCOS Spanish cohort. We found no association with risk of fracture at the genotype level for any of the SNPs, and no differences in the SNP frequencies between the two groups. At the haplotype level, we found no association between the T-delT-T haplotype and fracture. However, we observed a small but significant (p = 0.03) association with another rare haplotype, G-insT-T, which was slightly overrepresented in the patient group.

Fragmentation of contaminant and endogenous DNA in ancient samples determined by shotgun sequencing; prospects for human palaeogenomics.

BACKGROUND: Despite the successful retrieval of genomes from past remains, the prospects for human palaeogenomics remain unclear because of the difficulty of distinguishing contaminant from endogenous DNA sequences. Previous sequence data generated on high-throughput sequencing platforms indicate that fragmentation of ancient DNA sequences is a characteristic trait primarily arising due to depurination processes that create abasic sites leading to DNA breaks. METHODOLOGY/PRINCIPALS FINDINGS: To investigate whether this pattern is present in ancient remains from a temperate environment, we have 454-FLX pyrosequenced different samples dated between 5,500 and 49,000 years ago: a bone from an extinct goat (Myotragus balearicus) that was treated with a depurinating agent (bleach), an Iberian lynx bone not subjected to any treatment, a human Neolithic sample from Barcelona (Spain), and a Neandertal sample from the El Sidrón site (Asturias, Spain). The efficiency of retrieval of endogenous sequences is below 1% in all cases. We have used the non-human samples to identify human sequences (0.35 and 1.4%, respectively), that we positively know are contaminants. CONCLUSIONS: We observed that bleach treatment appears to create a depurination-associated fragmentation pattern in resulting contaminant sequences that is indistinguishable from previously described endogenous sequences. Furthermore, the nucleotide composition pattern observed in 5' and 3' ends of contaminant sequences is much more complex than the flat pattern previously described in some Neandertal contaminants. Although much research on samples with known contaminant histories is needed, our results suggest that endogenous and contaminant sequences cannot be distinguished by the fragmentation pattern alone.