Publications

Risk stratification in pulmonary arterial hypertension (PAH) is essential to provide more aggressive treatment for patients at higher risk. Nevertheless, recently introduced simplified prognostic tools neglect the genetic background. Additionally, pulmonary veno-oclusive disease (PVOD) has never been considered in risk assessment strategies.

Genetic recombination is a common evolutionary mechanism that produces molecular diversity. However, its consequences on protein folding stability have not attracted the same attention as in the case of point mutations. Here, we studied the effects of homologous recombination on the computationally predicted protein folding stability for several protein families, finding less detrimental effects than we previously expected. Although recombination can affect multiple protein sites, we found that the fraction of recombined proteins that are eliminated by negative selection because of insufficient stability is not significantly larger than the corresponding fraction of proteins produced by mutation events. Indeed, although recombination disrupts epistatic interactions, the mean stability of recombinant proteins is not lower than that of their parents. On the other hand, the difference of stability between recombined proteins is amplified with respect to the parents, promoting phenotypic diversity. As a result, at least one third of recombined proteins present stability between those of their parents, and a substantial fraction have higher or lower stability than those of both parents. As expected, we found that parents with similar sequences tend to produce recombined proteins with stability close to that of the parents. Finally, the simulation of protein evolution along the ancestral recombination graph with empirical substitution models commonly used in phylogenetics, which ignore constraints on protein folding stability, showed that recombination favors the decrease of folding stability, supporting the convenience of adopting structurally constrained models when possible for inferences of protein evolutionary histories with recombination.

Despite the increased recognition of (T-BOX transcription factor 4)-associated pulmonary arterial hypertension (PAH), genotype-phenotype associations are lacking and may provide important insights. To compile and functionally characterize all variants reported to date and undertake a comprehensive genotype-phenotype analysis. We assembled a multicenter cohort of 137 patients harboring monoallelic variants and assessed the pathogenicity of missense variation ( = 42) using a novel luciferase reporter assay containing T-BOX binding motifs. We sought genotype-phenotype correlations and undertook a comparative analysis with patients with PAH with (Bone Morphogenetic Protein Receptor type 2) causal variants ( = 162) or no identified variants in PAH-associated genes ( = 741) genotyped via the National Institute for Health Research BioResource-Rare Diseases. Functional assessment of missense variants led to the novel finding of gain-of-function effects associated with older age at diagnosis of lung disease compared with loss-of-function effects ( = 0.038). Variants located in the T-BOX and nuclear localization domains were associated with earlier presentation ( = 0.005) and increased incidence of interstitial lung disease ( = 0.003). Event-free survival (death or transplantation) was shorter in the T-BOX group ( = 0.022), although age had a significant effect in the hazard model ( = 0.0461). Carriers of variants were diagnosed at a younger age ( < 0.001) and had worse baseline lung function (FEV, FVC) ( = 0.009) than the and no identified causal variant groups. We demonstrated that syndrome is not strictly the result of haploinsufficiency but can also be caused by gain of function. The pleiotropic effects of in lung disease may be in part explained by the differential effect of pathogenic mutations located in critical protein domains.

Cancer still poses a global threat, since a lot of tumors remain untreatable despite all the available chemotherapeutic drugs, whose side effects, it must also be noted, still raise concerns. The antitumoral properties of marine seaweeds make them a potential source of new, less toxic, and more active antitumoral agents. Furthermore, these natural extracts can be combined with nanotechnology to increase their efficacy and improve targeting. In this work, a (CT) aqueous extract was employed for the green synthesis of gold nanoparticles (Au@CT). The complete characterization of Au@CT was performed by UV-Vis spectroscopy, Fourier transform infrared (FTIR) spectroscopy, Zeta potential, electron microscopy, X-ray powder diffraction (XRD), high-performance steric exclusion chromatography (HPSEC), and by the determination of their antioxidant capacity. The antiproliferative activity of Au@CT was then tested in hepatic (HEPG-2) and pancreatic (BxPC-3) cell lines. Their potential capacity as enhancers of gemcitabine, a drug frequently used to treat both types of tumors, was also tested. The activity of Au@CT was compared to the activity of the CT extract alone. A synergistic effect with gemcitabine was proven for HEPG-2. Our results showed that gold nanoparticles synthesized from seaweed extracts with antitumoral activity could be a good gemcitabine enhancer.

The genomic profiling of circulating tumor cells (CTCs) in the bloodstream should provide clinically relevant information on therapeutic efficacy and help predict cancer survival. Here, we contrasted the genomic profiles of CTC pools recovered from metastatic colorectal cancer (mCRC) patients using different enrichment strategies (CellSearch, Parsortix, and FACS). Mutations inferred in the CTC pools differed depending on the enrichment strategy and, in all cases, represented a subset of the mutations detected in the matched primary tumor samples. However, the CTC pools from Parsortix, and in part, CellSearch, showed diversity estimates, mutational signatures, and drug-suitability scores remarkably close to those found in matching primary tumor samples. In addition, FACS CTC pools were enriched in apparent sequencing artifacts, leading to much higher genomic diversity estimates. Our results highlight the utility of CTCs to assess the genomic heterogeneity of individual tumors and help clinicians prioritize drugs in mCRC.

We present SIEVE, a statistical method for the joint inference of somatic variants and cell phylogeny under the finite-sites assumption from single-cell DNA sequencing. SIEVE leverages raw read counts for all nucleotides and corrects the acquisition bias of branch lengths. In our simulations, SIEVE outperforms other methods in phylogenetic reconstruction and variant calling accuracy, especially in the inference of homozygous variants. Applying SIEVE to three datasets, one for triple-negative breast (TNBC), and two for colorectal cancer (CRC), we find that double mutant genotypes are rare in CRC but unexpectedly frequent in the TNBC samples.

Recurrence of tumor cells following local and systemic therapy is a significant hurdle in cancer. Most patients with metastatic colorectal cancer (mCRC) will relapse, despite resection of the metastatic lesions. A better understanding of the evolutionary history of recurrent lesions is required to identify the spatial and temporal patterns of metastatic progression and expose the genetic and evolutionary determinants of therapeutic resistance. With this goal in mind, here we leveraged a unique single-cell whole-genome sequencing dataset from recurrent hepatic lesions of an mCRC patient. Our phylogenetic analysis confirms that the treatment induced a severe demographic bottleneck in the liver metastasis but also that a previously diverged lineage survived this surgery, possibly after migration to a different site in the liver. This lineage evolved very slowly for two years under adjuvant drug therapy and diversified again in a very short period. We identified several non-silent mutations specific to this lineage and inferred a substantial contribution of chemotherapy to the overall, genome-wide mutational burden. All in all, our study suggests that mCRC subclones can migrate locally and evade resection, keep evolving despite rounds of chemotherapy, and re-expand explosively.

Due to their rarity and heterogeneity and despite the introduction of molecular features in the current WHO classification, clinical criteria such as those from the European Organization for Research and Treatment of Cancer (EORTC) and the Radiation Therapy Oncology Group (RTOG) are still being used to make treatment decisions in low-grade gliomas (LGG). Patients with diffuse low-grade glioma treated at our institution between 2002 and 2018 were analyzed, retrieving and assessing the degree of consistency between the EORTC and RTOG criteria, as well as the isocitrate dehydrogenase 1 and 2 (IDH) gene mutational status. Likewise, multivariate analyses were performed to ascertain the superiority of any of the factors over the others. One hundred and two patients were included. The degree of consistency between the RTOG and EORTC criteria was 71.6% (K = 0.426; p = 0.0001). Notably, 51.7% of those assigned to low risk by the EORTC were classified as high risk according to the RTOG classification. In multivariate analysis, only complete resection, age > 40 years, size and IDH mutation status were independently correlated with OS. When the RTOG and EORTC scores were entered into the model, only the EORTC model was independently associated with mortality. The degree of consistency between the EORT and RTOG criteria is low. Therefore, there is a need to integrate clinical-molecular scores to improve treatment decisions in LGG.

Advances in phylogenomics and high-throughput sequencing have allowed the reconstruction of deep phylogenetic relationships in the evolution of eukaryotes. Yet, the root of the eukaryotic tree of life remains elusive. The most popular hypothesis in textbooks and reviews is a root between Unikonta (Opisthokonta + Amoebozoa) and Bikonta (all other eukaryotes), which emerged from analyses of a single-gene fusion. Subsequent, highly cited studies based on concatenation of genes supported this hypothesis with some variations or proposed a root within Excavata. However, concatenation of genes does not consider phylogenetically-informative events like gene duplications and losses. A recent study using gene tree parsimony (GTP) suggested the root lies between Opisthokonta and all other eukaryotes, but only including 59 taxa and 20 genes. Here we use GTP with a duplication-loss model in a gene-rich and taxon-rich dataset (i.e., 2,786 gene families from two sets of 155 and 158 diverse eukaryotic lineages) to assess the root, and we iterate each analysis 100 times to quantify tree space uncertainty. We also contrasted our results and discarded alternative hypotheses from the literature using GTP and the likelihood-based method SpeciesRax. Our estimates suggest a root between Fungi or Opisthokonta and all other eukaryotes; but based on further analysis of genome size, we propose that the root between Opisthokonta and all other eukaryotes is the most likely.

The selection of the best-fitting substitution model of molecular evolution is a traditional step for phylogenetic inferences, including ancestral sequence reconstruction (ASR). However, a few recent studies suggested that applying this procedure does not affect the accuracy of phylogenetic tree reconstruction. Here, we revisited this debate topic by analyzing the influence of selection among substitution models of protein evolution, with focus on exchangeability matrices, on the accuracy of ASR using simulated and real data. We found that the selected best-fitting substitution model produces the most accurate ancestral sequences, especially if the data present large genetic diversity. Indeed, ancestral sequences reconstructed under substitution models with similar exchangeability matrices were similar, suggesting that if the selected best-fitting model cannot be used for the reconstruction, applying a model similar to the selected one is preferred. We conclude that selecting among substitution models of protein evolution is recommended for reconstructing accurate ancestral sequences.

Bardet-Biedl syndrome (BBS) is an autosomal recessive ciliopathy characterized by extensive inter- and intra-familial variability, in which oligogenic interactions have been also reported. Our main goal is to elucidate the role of mutational load in the clinical variability of BBS. A cohort of 99 patients from 77 different families with biallelic pathogenic variants in a BBS-associated gene was retrospectively recruited. Human Phenotype Ontology terms were used in the annotation of clinical symptoms. The mutational load in 39 BBS-related genes was studied in index cases using different molecular and next-generation sequencing (NGS) approaches. Candidate allele combinations were analysed using the in silico tools ORVAL and DiGePred. After clinical annotation, 76 out of the 99 cases a priori fulfilled established criteria for diagnosis of BBS or BBS-like. BBS1 alleles, found in 42% of families, were the most represented in our cohort. An increased mutational load was excluded in 41% of the index cases (22/54). Oligogenic inheritance was suspected in 52% of the screened families (23/45), being 40 tested by means of NGS data and 5 only by traditional methods. Together, ORVAL and DiGePred platforms predicted an oligogenic effect in 44% of the triallelic families (10/23). Intrafamilial variable severity could be clinically confirmed in six of the families. Our findings show that the presence of more than two alleles in BBS-associated genes correlated in six families with a more severe phenotype and associated with specific findings, highlighting the role of the mutational load in the management of BBS cases.

Single-nucleotide variants (SNVs) are the most common variations in the human genome. Recently developed methods for SNV detection from single-cell DNA sequencing data, such as SCIΦ and scVILP, leverage the evolutionary history of the cells to overcome the technical errors associated with single-cell sequencing protocols. Despite being accurate, these methods are not scalable to the extensive genomic breadth of single-cell whole-genome (scWGS) and whole-exome sequencing (scWES) data.

Primary cilia are non-motile organelles associated with the cell cycle, which can be found in most vertebrate cell types. Cilia formation occurs through a process called ciliogenesis, which involves several mechanisms including planar cell polarity (PCP) and the Hedgehog (Hh) signaling pathway. Some gene complexes, such as BBSome or CPLANE (ciliogenesis and planar polarity effector), have been linked to ciliogenesis. CPLANE complex is composed of , and , which bind to and for cilia formation. Defects in these genes have been linked to a malfunction of intraflagellar transport and defects in the planar cell polarity, as well as defective activation of the Hedgehog signalling pathway. These faults lead to defective cilium formation, resulting in ciliopathies, including orofacial-digital syndrome (OFDS) and Bardet-Biedl syndrome (BBS). Considering the close relationship, between the CPLANE complex and cilium formation, it can be expected that defects in the genes that encode subunits of the CPLANE complex may be related to other ciliopathies.

The integral utilization of sustainable resources with versatile, efficient and cleaner processes is encouraged. Hydrothermal treatment with subcritical water is a chemical free, tunable and rapid technology providing enhanced yield compared to conventional extraction and was explored for the benign by design extraction and depolymerization of carrageenan from Chondrus crispus. Up to 90% of the seaweed was solubilized operating under nonisothermal regime during heating up to 200 °C and 75.5% crude carrageenan yield was attained at 140 °C. Crude carrageenan could not be precipitated by ethanol from the extracts produced at 180 °C and higher temperatures, but ultrafiltration (100 kDa) of the extract obtained at 160 °C provided comparable recovery yields and similar rheological features to those of the ethanol precipitated product. Operation at 140 °C was preferred based on the higher recovery yield of the biopolymer and the whole extract was suitable for the green synthesis of polycrystalline decahedral quasi-spherical gold nanoparticles with a mean size distribution of 8.4 nm and Z potential value of -40.2 mV. Alternatively, the crude carrageen fraction was used for the formulation of printable biopolymer based gels with suitable mechanical properties, including a relevant gel strength enhancement (about 10-fold) when compared with conventional procedures.

TOPIIA topoisomerases are required for the regulation of DNA topology by DNA cleavage and re-ligation and are important targets of antibiotic and anticancer agents. Humans possess two TOPIIA paralogue genes (TOP2A and TOP2B) with high sequence and structural similarity but distinct cellular functions. Despite their functional and clinical relevance, the evolutionary history of TOPIIA is still poorly understood. Here we show that TOPIIA is highly conserved in Metazoa. We also found that TOPIIA paralogues from jawed and jawless vertebrates had different origins related with tetraploidization events. After duplication, TOP2B evolved under a stronger purifying selection than TOP2A, perhaps promoted by the more specialized role of TOP2B in postmitotic cells. We also detected genetic signatures of positive selection in the highly variable C-terminal domain (CTD), possibly associated with adaptation to cellular interactions. By comparing TOPIIA from modern and archaic humans, we found two amino acid substitutions in the TOP2A CTD, suggesting that TOP2A may have contributed to the evolution of present-day humans, as proposed for other cell cycle-related genes. Finally, we identified six residues conferring resistance to chemotherapy differing between TOP2A and TOP2B. These six residues could be targets for the development of TOP2A-specific inhibitors that would avoid the side effects caused by inhibiting TOP2B. Altogether, our findings clarify the origin, diversification and selection pressures governing the evolution of animal TOPIIA.

Human mitochondria can be genetically distinct within the same individual, a phenomenon known as heteroplasmy. In cancer, this phenomenon seems exacerbated, and most mitochondrial mutations seem to be heteroplasmic. How this genetic variation is arranged within and among normal and tumor cells is not well understood. To address this question, here we sequenced single-cell mitochondrial genomes from multiple normal and tumoral locations in four colorectal cancer patients. Our results suggest that single cells, both normal and tumoral, can carry various mitochondrial haplotypes. Remarkably, this intra-cell heteroplasmy can arise before tumor development and be maintained afterward in specific tumoral cell subpopulations. At least in the colorectal patients studied here, the somatic mutations in the single-cells do not seem to have a prominent role in tumorigenesis.

Long-Term Non-Progressors (LTNPs) are untreated Human Immunodeficiency virus type 1 (HIV-1) infected individuals able to control disease progression for prolonged periods. However, the LTNPs status is temporary, as viral load increases followed by decreases in CD4 + T-cell counts. Control of HIV-1 infection in LTNPs viremic controllers, have been associated with effective immunodominant HIV-1 Gag-CD8 + T-cell responses restricted by protective HLA-B alleles. Individuals carrying HLA-B*14:02 control HIV-1 infection is related to an immunodominant Env-CD8 + T-cell response. Limited data are available on the contribution of HLA-B*14:02 CD8 + T -cells in LTNPs.

The last glacial period (LGP) promoted a loss of genetic diversity in Paleolithic populations of modern humans from diverse regions of the world by range contractions and habitat fragmentation. However, this period also provided some currently submersed lands, such as the Sunda shelf in Southeast Asia (SEA), that could have favored the expansion of our species. Concerning the latter, still little is known about the influence of the lowering sea level on the genetic diversity of current SEA populations. Here, we applied approximate Bayesian computation, based on extensive spatially explicit computer simulations, to evaluate the fitting of mtDNA data from diverse SEA populations with alternative evolutionary scenarios that consider and ignore the LGP and migration through long-distance dispersal (LDD). We found that both the LGP and migration through LDD should be taken into consideration to explain the currently observed genetic diversity in these populations and supported a rapid expansion of first populations throughout SEA. We also found that temporarily available lands caused by the low sea level of the LGP provided additional resources and migration corridors that favored genetic diversity. We conclude that migration through LDD and temporarily available lands during the LGP should be considered to properly understand and model the first expansions of modern humans.

This study reports on the green and cost-efficient synthesis of gold nanoparticles from three different red algae extracts. The nanoparticles synthesized were fully characterized by UV-Vis spectroscopy, HRTEM, and Z-potential. Relevant components occurring in the extracts, such as polysaccharides or phenolic content, were assessed by analytical techniques such as spectrophotometric assays and liquid chromatography. Finally, the antioxidant, antitumoral, and anti-inflammatory potential of both the extracts and the gold nanoparticles synthesized were analyzed in order to determine a possible synergistic effect on the nanoparticles. The results obtained confirmed the obtainment of gold nanoparticles with significant potential as immunotherapeutic agents. The therapeutic potential of these nanoparticles could be higher than that of inert gold nanoparticles loaded with bioactive molecules since the former would allow for higher accumulation into the targeted tissue.

Clonally transmissible cancers are tumour lineages that are transmitted between individuals via the transfer of living cancer cells. In marine bivalves, leukaemia-like transmissible cancers, called hemic neoplasia (HN), have demonstrated the ability to infect individuals from different species. We performed whole-genome sequencing in eight warty venus clams that were diagnosed with HN, from two sampling points located more than 1000 nautical miles away in the Atlantic Ocean and the Mediterranean Sea Coasts of Spain. Mitochondrial genome sequencing analysis from neoplastic animals revealed the coexistence of haplotypes from two different clam species. Phylogenies estimated from mitochondrial and nuclear markers confirmed this leukaemia originated in striped venus clams and later transmitted to clams of the species warty venus, in which it survives as a contagious cancer. The analysis of mitochondrial and nuclear gene sequences supports all studied tumours belong to a single neoplastic lineage that spreads in the Seas of Southern Europe.

Intra-host evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been reported in cases with persistent coronavirus disease 2019 (COVID-19). In this study, we describe a severely immunosuppressed individual with HIV-1/SARS-CoV-2 coinfection with a long-term course of SARS-CoV-2 infection. A 28-year-old man was diagnosed with HIV-1 infection (CD4+ count: 3 cells/µL nd 563000 HIV-1 RNA copies/mL) and simultaneous pneumonia, disseminated infection and SARS-CoV-2 infection. SARS-CoV-2 real-time reverse transcription polymerase chain reaction positivity from nasopharyngeal samples was prolonged for 15 weeks. SARS-CoV-2 was identified as variant Alpha (PANGO lineage B.1.1.7) with mutation S:E484K. Spike-specific T-cell response was similar to HIV-negative controls although enriched in IL-2, and showed disproportionately increased immunological exhaustion marker levels. Despite persistent SARS-CoV-2 infection, adaptive intra-host SARS-CoV-2 evolution, was not identified. Spike-specific T-cell response protected against a severe COVID-19 outcome and the increased immunological exhaustion marker levels might have favoured SARS-CoV-2 persistence.

We introduce CellPhy, a maximum likelihood framework for inferring phylogenetic trees from somatic single-cell single-nucleotide variants. CellPhy leverages a finite-site Markov genotype model with 16 diploid states and considers amplification error and allelic dropout. We implement CellPhy into RAxML-NG, a widely used phylogenetic inference package that provides statistical confidence measurements and scales well on large datasets with hundreds or thousands of cells. Comprehensive simulations suggest that CellPhy is more robust to single-cell genomics errors and outperforms state-of-the-art methods under realistic scenarios, both in accuracy and speed. CellPhy is freely available at https://github.com/amkozlov/cellphy .

Single-cell sequencing has gained popularity in recent years. Despite its numerous applications, single-cell DNA sequencing data is highly error-prone due to technical biases arising from uneven sequencing coverage, allelic dropout, and amplification error. With these artifacts, the identification of somatic genomic variants becomes a challenging task, and over the years, several methods have been developed explicitly for this type of data. Single-cell variant callers implement distinct strategies, make different use of the data, and typically result in many discordant calls when applied to real data. Here, we review current approaches for single-cell variant calling, emphasizing single nucleotide variants. We highlight their potential benefits and shortcomings to help users choose a suitable tool for their data at hand.

is a ubiquitous gene associated with Alström syndrome (ALMS). The main symptoms of ALMS affect multiple organs and tissues, generating at last, multi-organic fibrosis in the lungs, kidneys and liver. TGF-β is one of the main pathways implicated in fibrosis, controlling the cell cycle, apoptosis, cell migration, cell adhesion and epithelial-mesenchymal transition (EMT). Nevertheless, the role of gene in fibrosis generation and other implicated processes such as cell migration or cell adhesion the TGF- β pathway has not been elucidated yet. Initially, we evaluated how depletion of affects different processes like apoptosis, cell cycle and mitochondrial activity in HeLa cells. Then, we performed proteomic profiling with TGF-β stimuli in HeLa -/- cells and validated the results by examining different EMT biomarkers using qPCR. The expression of these EMT biomarkers were also studied in hTERT-BJ-5ta -/-. Finally, we evaluated the SMAD3 and SMAD2 phosphorylation and cell migration capacity in both models. Depletion of generated apoptosis resistance to thapsigargin (THAP) and C2-Ceramide (C2-C), and G2/M cell cycle arrest in HeLa cells. For mitochondrial activity, results did not show significant differences between and . Proteomic results showed inhibition of downstream pathways regulated by TGF-β. The protein-coding genes (PCG) were associated with processes like focal adhesion or cell-substrate adherens junction in HeLa. showed an opposite pattern to what would be expected when activating the EMT in HeLa and BJ-5ta. Finally, in BJ-5ta model a reduced activation of SMAD3 but not SMAD2 were also observed. In HeLa model no alterations in the canonical TGF-β pathway were observed but both cell lines showed a reduction in migration capacity. has a role in controlling the cell cycle and the apoptosis processes. Moreover, the depletion of affects the signal transduction through the TGF-β and other processes like the cell migration and adhesion capacity.

Octocorals possess sclerites, small elements comprised of calcium carbonate (CaCO) that are important diagnostic characters in octocoral taxonomy. Among octocorals, sea pens comprise a unique order (Pennatulacea) that live in a wide range of depths. Habitat depth is considered to be important in the diversification of octocoral species, but a lack of information on sea pens has limited studies on their adaptation and evolution across depth. Here, we aimed to reveal trends of adaptation and evolution of sclerite shapes in sea pens with regards to habitat depth phylogenetic analyses and ancestral reconstruction analyses. Colony form of sea pens is suggested to have undergone convergent evolution and the loss of axis has occurred independently across the evolution of sea pens. Divergences of sea pen taxa and of sclerite forms are suggested to depend on habitat depths. In addition, their sclerite forms may be related to evolutionary history of the sclerite and the surrounding chemical environment as well as water temperature. Three-flanged sclerites may possess the tolerance towards the environment of the deep sea, while plate sclerites are suggested to be adapted towards shallower waters, and have evolved independently multiple times. The common ancestor form of sea pens was predicted to be deep-sea and similar to family Pseudumbellulidae in form, possessing sclerites intermediate in form to those of alcyonaceans and modern sea pens such as spindles, rods with spines, and three-flanged sclerites with serrated edges sclerites, as well as having an axis and bilateral traits.

The reconstruction of genetic material of ancestral organisms constitutes a powerful application of evolutionary biology. A fundamental step in this inference is the ancestral sequence reconstruction (ASR), which can be performed with diverse methodologies implemented in computer frameworks. However, most of these methodologies ignore evolutionary properties frequently observed in microbes, such as genetic recombination and complex selection processes, that can bias the traditional ASR. From a practical perspective, here I review methodologies for the reconstruction of ancestral DNA and protein sequences, with particular focus on microbes, and including biases, recommendations, and software implementations. I conclude that microbial ASR is a complex analysis that should be carefully performed and that there is a need for methods to infer more realistic ancestral microbial sequences.

A detailed understanding of how and when severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission occurs is crucial for designing effective prevention measures. Other than contact tracing, genome sequencing provides information to help infer who infected whom. However, the effectiveness of the genomic approach in this context depends on both (high enough) mutation and (low enough) transmission rates. Today, the level of resolution that we can obtain when describing SARS-CoV-2 outbreaks using just genomic information alone remains unclear. In order to answer this question, we sequenced forty-nine SARS-CoV-2 patient samples from ten local clusters in NW Spain for which partial epidemiological information was available and inferred transmission history using genomic variants. Importantly, we obtained high-quality genomic data, sequencing each sample twice and using unique barcodes to exclude cross-sample contamination. Phylogenetic and cluster analyses showed that consensus genomes were generally sufficient to discriminate among independent transmission clusters. However, levels of intrahost variation were low, which prevented in most cases the unambiguous identification of direct transmission events. After filtering out recurrent variants across clusters, the genomic data were generally compatible with the epidemiological information but did not support specific transmission events over possible alternatives. We estimated the effective transmission bottleneck size to be one to two viral particles for sample pairs whose donor-recipient relationship was likely. Our analyses suggest that intrahost genomic variation in SARS-CoV-2 might be generally limited and that homoplasy and recurrent errors complicate identifying shared intrahost variants. Reliable reconstruction of direct SARS-CoV-2 transmission based solely on genomic data seems hindered by a slow mutation rate, potential convergent events, and technical artifacts. Detailed contact tracing seems essential in most cases to study SARS-CoV-2 transmission at high resolution.

The evolution of structural proteins is generally constrained by the folding stability. However, little is known about the particular capacity of viral proteins to accommodate mutations that can potentially affect the protein stability and, in general, the evolution of the protein stability over time. As an illustrative model case, here, we investigated the evolution of the stability of the human immunodeficiency virus (HIV-1) protease (PR), which is a common HIV-1 drug target, under diverse evolutionary scenarios that include (1) intra-host virus evolution in a cohort of seventy-five patients sampled over time, (2) intra-host virus evolution sampled before and after specific PR-based treatments, and (3) inter-host evolution considering extant and ancestral (reconstructed) PR sequences from diverse HIV-1 subtypes. We also investigated the specific influence of currently known HIV-1 PR resistance mutations on the PR folding stability. We found that the HIV-1 PR stability fluctuated over time within a constant and wide range in any studied evolutionary scenario, accommodating multiple mutations that partially affected the stability while maintaining activity. We did not identify relationships between change of PR stability and diverse clinical parameters such as viral load, CD4 T-cell counts, and a surrogate of time from infection. Counterintuitively, we predicted that nearly half of the studied HIV-1 PR resistance mutations do not significantly decrease stability, which, together with compensatory mutations, would allow the protein to adapt without requiring dramatic stability changes. We conclude that the HIV-1 PR presents a wide structural plasticity to acquire molecular adaptations without affecting the overall evolution of stability.

Pulmonary Arterial Hypertension (PAH) is a rare disease caused by the obliteration of the pulmonary arterioles, increasing pulmonary vascular resistance and eventually causing right heart failure. Endothelin-1 (EDN1) is a vasoconstrictor peptide whose levels are indicators of disease progression and its pathway is one of the most common targeted by current treatments. We sequenced the EDN1 untranslated regions of a small subset of patients with PAH, predicted the effect , and used a luciferase assay with the different genotypes to analyze its influence on gene expression. Finally, we used siRNAs against the major transcription factors (TFs) predicted for these regions [peroxisome proliferator-activated receptor γ (PPARγ), Krüppel-Like Factor 4 (KLF4), and vitamin D receptor (VDR)] to assess EDN1 expression in cell culture and validate the binding sites. First, we detected a single nucleotide polymorphism (SNP) in the 5' untranslated region (UTR; rs397751713) and another in the 3'regulatory region (rs2859338) that altered luciferase activity depending on their genotype. We determined that KLF4/PPARγ could bind to the rs397751713 and VDR to rs2859338. By using siRNAs and luciferase assays, we determined that PPARγ binds differentially to rs397751713. PPARγ and VDR Knock-Down (KD) increased the mRNA levels and EDN1 production in porcine aortic endothelial cells (PAECs), while PPARγ and KLF4 KD increased the EDN1 production in HeLa. In conclusion, common variants in EDN1 regulatory regions could alter EDN1 levels. We were able to validate that PPARγ binds in rs397751713 and is a key regulator of EDN1. In addition, KLF4 and VDR regulate EDN1 production in a cell-dependent manner, but VDR does not bind directly to the regions we studied.

The evolutionary processes of mutation and recombination, upon which selection operates, are fundamental to understand the observed molecular diversity. Unlike nucleotide sequences, the estimation of the recombination rate in protein sequences has been little explored, neither implemented in evolutionary frameworks, despite protein sequencing methods are largely used.

A coalescent model of a sample of size n is derived from a birth-death process that originates at a random time in the past from a single founder individual. Over time, the descendants of the founder evolve into a population of large (infinite) size from which a sample of size n is taken. The parameters and time of the birth-death process are scaled in N, the size of the present-day population, while letting N→∞, similarly to how the standard Kingman coalescent process arises from the Wright-Fisher model. The model is named the Limit Birth-Death (LBD) coalescent model. Simulations from the LBD coalescent model with sample size n are computationally slow compared to standard coalescent models. Therefore, we suggest different approximations to the LBD coalescent model assuming the population size is a deterministic function of time rather than a stochastic process. Furthermore, we introduce a hybrid LBD coalescent model, that combines the exactness of the LBD coalescent model model with the speed of the approximations.

Diverse phylogenetic methods require a substitution model of evolution that should mimic, as accurately as possible, the real substitution process. At the protein level, empirical substitution models have traditionally been based on a large number of different proteins from particular taxonomic levels. However, these models assume that all of the proteins of a taxonomic level evolve under the same substitution patterns. We believe that this assumption is highly unrealistic and should be relaxed by considering protein-specific substitution models that account for protein-specific selection processes. In order to test this hypothesis, we inferred and evaluated four new empirical substitution models for the protease and integrase of HIV and other viruses. We found that these models more accurately fit, compared with any of the currently available empirical substitution models, the evolutionary process of these proteins. We conclude that evolutionary inferences from protein sequences are more accurate if they are based on protein-specific substitution models rather than taxonomic-specific (generalist) substitution models. We also present four new empirical substitution models of protein evolution that could be useful for phylogenetic inferences of viral protease and integrase.

Extensive transmission of SARS-CoV-2 during the COVID-19 pandemic allowed the generation of thousands of mutations within its genome. While several of these become rare, others largely increase in prevalence, potentially jeopardizing the sensitivity of PCR-based diagnostics. Taking advantage of SARS-CoV-2 genomic knowledge, we designed a one-step probe-based multiplex RT-qPCR (OmniSARS2) to simultaneously detect short fragments of the SARS-CoV-2 genome in ORF1ab, E gene and S gene. Comparative genomics of the most common SARS-CoV-2 lineages, other human betacoronavirus and alphacoronavirus, was the basis for this design, targeting both highly conserved regions across SARS-CoV-2 lineages and variable or absent in other viruses. The highest analytical sensitivity of this method for SARS-CoV-2 detection was 94.2 copies/mL at 95% detection probability (~1 copy per total reaction volume) for the S gene assay, matching the most sensitive available methods. In vitro specificity tests, performed using reference strains, showed no cross-reactivity with other human coronavirus or common pathogens. The method was compared with commercially available methods and detected the virus in clinical samples encompassing different SARS-CoV-2 lineages, including B.1, B.1.1, B.1.177 or B.1.1.7 and rarer lineages. OmniSARS2 revealed a sensitive and specific viral detection method that is less likely to be affected by lineage evolution oligonucleotide-sample mismatch, of relevance to ensure the accuracy of COVID-19 molecular diagnostic methods.

The COVID-19 pandemic is promoting online education. In this concern, not only online teaching of wet laboratory practices is being challenging, but also online teaching of practices in bioinformatics. Here, I discuss the problematic of online teaching practices in bioinformatics and provide fundamental guidelines, especially oriented for beginner educators, that can help in the time of the pandemic.

The evolution of protein-coding genes is usually driven by selective processes, which favor some evolutionary trajectories over others, optimizing the subsequent protein stability and activity. The analysis of selection in this type of genetic data is broadly performed with the metric nonsynonymous/synonymous substitution rate ratio (dN/dS). However, most of the well-established methodologies to estimate this metric make crucial assumptions, such as lack of recombination or invariable codon frequencies along genes, which can bias the estimation. Here, we review the most relevant biases in the dN/dS estimation and provide a detailed guide to estimate this metric using state-of-the-art procedures that account for such biases, along with illustrative practical examples and recommendations. We also discuss the traditional interpretation of the estimated dN/dS emphasizing the importance of considering complementary biological information such as the role of the observed substitutions on the stability and function of proteins. This review is oriented to help evolutionary biologists that aim to accurately estimate selection in protein-coding sequences.

A new species of Etelis is described based on 16 specimens collected from the Red Sea and Western Australia, with confirmed genetic records throughout the Indo-West Pacific. It is similar to and was often misidentified as Etelis carbunculus Cuvier, with both species sharing the diagnostic character of low number of developed gill rakers. Nonetheless, the two species are genetically divergent and differ morphologically in adult body length; proportions of eye, snout, cheek and caudal fin; shape of head, opercular spine and sagittal otolith; and coloration of the tip of the upper caudal fin. Etelis boweni has a wide Indo-west Pacific distribution that largely overlaps with E. carbunculus, and the two species are often caught on the same fishing line.

DNA topoisomerase III beta (TOP3B) is unique by operating on both DNA and RNA substrates to regulate gene expression and genomic stability. Mutations in human TOP3B are linked to neurodevelopmental and cognitive disorders, highlighting its relevance for human health. Despite the emerging importance of TOP3B, its precise cellular functions and evolutionary history remain poorly understood. Here, we show that TOP3B is conserved across main metazoan groups and evolved under strong purifying selection. Subdomain IV was identified as the most conserved TOP3B region, in agreement with its role in providing the structural foundation of the protein. On the contrary, subdomain II is the less conserved, possibly because it is the most structurally flexible region of all TOP3B regions. Interestingly, TOP3B residue at position 472, previously associated with schizophrenia, is highly variable across animals, suggesting a more specific role in humans and related species. Finally, we show that all TOP3B CXXC zinc finger motifs previously identified at the protein C-terminal region are retained across metazoans. We also found that the two major methylation sites known to regulate TOP3B activity are located in the most conserved region of the C-terminal arginine-glycine-glycine (RGG) box, suggesting that a similar regulatory mechanism may operate throughout animals. Overall, our results provide a better understanding of the evolution and functional roles of TOP3B.

Carrageenan was extracted from Mastocarpus stellatus using hot water extraction under atmospheric and pressurized conditions. The influence of heating temperature during a non-isothermal heating profile up to temperatures in the range 70-190 °C was studied to evaluate the extraction yields and properties of the carrageenan fraction. Under the selected conditions (130 °C), extracted carrageenan (CMs) was used for the green synthesis of gold nanoparticles (AuNPs). After the optimization of the reaction conditions, the synthesized gold nanoparticles (Au@CMs) were characterized by UV-Vis spectroscopy, Z potential measurements, electron microscopy, and X-ray diffraction analysis, which confirmed the formation of spherical, polycrystalline, and negatively charged nanoparticles with a mean diameter of 14.3 ± 2.1 nm. The study conducted by scanning transmission electron microscopy, energy dispersive X-ray analysis and mapping confirmed the presence of carrageenan stabilizing AuNPs. Finally, Fourier transformed infrared spectroscopy was performed to analyze the functional groups of CMs involved in the reduction and stabilization of AuNPs. The selective cytotoxicity and the antioxidant activity of the Au@CMs were evaluated in different cell lines and compared to the CMs. Au@CMs showed an improved antioxidant capacity in cells under oxidative stress and the induction of apoptosis in a monocytic cell line, while no antitumor effect was observed in a lung endothelial cell line.

Pulmonary Arterial Hypertension (PAH) is a severe complication of Connective Tissue Disease (CTD), with remarkable morbidity and mortality. However, the molecular and genetic basis of CTD-PAH remains incompletely understood. This study aimed to screen for genetic defects in a cohort of patients with CTD-PAH, using a PAH-specific panel of 35 genes. During recruitment, 79 patients were studied, including 59 Systemic Sclerosis patients (SSc) and 69 females. Disease-associated variants were observed in nine patients: 4 pathogenic/likely pathogenic variants in 4 different genes ( and ) and 5 Variants of Unknown Significance (VUS) in 4 genes ( and ). One patient with mixed CTD had a frameshift pathogenic variant in . Two patients with SSc-PAH carried variants in . A patient diagnosed with Systemic Lupus Erythematous (SLE) presented a pathogenic nonsense variant in . Another patient with SSc-PAH presented a pathogenic variant in . Four patients with SSc-PAH carried a VUS in and , respectively. These findings suggest that genetic factors may contribute to Pulmonary Vascular Disease (PVD) in CTD patients.

Over the last years, there has been an increasing trend towards the use of environmentally friendly processes to synthesize nanomaterials. In the case of nanomedicine, the use of bionanofactories with associated biological properties, such as seaweed, has emerged as a promising field of work due to the possibility they open for both the preservation of those properties in the nanomaterials synthesized and/or the reduction of their toxicity. In the present study, gold (Au@SP) and silver (Ag@SP) nanoparticles were synthesized using an aqueous extract of Saccorhiza polyschides (SP). Several techniques showed that the nanoparticles formed were spherical and stable, with mean diameters of 14 ± 2 nm for Au@SP and 15 ± 3 nm for Ag@SP. The composition of the biomolecules in the extract and the nanoparticles were also analyzed. The analyses performed indicate that the extract acts as a protective medium, with the particles embedded in it preventing aggregation and coalescence. Au@SP and Ag@SP showed superior immunostimulant and antiproliferative activity on immune and tumor cells, respectively, to that of the SP extract. Moreover, the nanoparticles were able to modulate the release of reactive oxygen species depending on the concentration. Hence, both nanoparticles have a significant therapeutic potential for the treatment of cancer or in immunostimulant therapy.

Recombination between the X and Y human sex chromosomes is limited to the two pseudoautosomal regions (PARs) that present quite distinct evolutionary origins. Despite the crucial importance for male meiosis, genetic diversity patterns and evolutionary dynamics of these regions are poorly understood. In the present study, we analyzed and compared the genetic diversity of the PAR regions using publicly available genomic sequences encompassing both PAR1 and PAR2. Comparisons were performed through allele diversities, linkage disequilibrium status and recombination frequencies within and between X and Y chromosomes. In agreement with previous studies, we confirmed the role of PAR1 as a male-specific recombination hotspot, but also observed similar characteristic patterns of diversity in both regions although male recombination occurs at PAR2 to a much lower extent (at least one recombination event at PAR1 and in ≈1% in normal male meioses at PAR2). Furthermore, we demonstrate that both PARs harbor significantly different allele frequencies between X and Y chromosomes, which could support that recombination is not sufficient to homogenize the pseudoautosomal gene pool or is counterbalanced by other evolutionary forces. Nevertheless, the observed patterns of diversity are not entirely explainable by sexually antagonistic selection. A better understanding of such processes requires new data from intergenerational transmission studies of PARs, which would be decisive on the elucidation of PARs evolution and their role in male-driven heterosomal aneuploidies.

In 1981, the Journal of Molecular Evolution (JME) published an article entitled "Evolutionary trees from DNA sequences: A maximum likelihood approach" by Joseph (Joe) Felsenstein (J Mol Evol 17:368-376, 1981). This groundbreaking work laid the foundation for the emerging field of statistical phylogenetics, providing a tractable way of finding maximum likelihood (ML) estimates of evolutionary trees from DNA sequence data. This paper is the second most cited (more than 9000 citations) in JME after Kimura's (J Mol Evol 16:111-120, 1980) seminal paper on a model of nucleotide substitution (with nearly 20,000 citations). On the occasion of the 50th anniversary of JME, we elaborate on the significance of Felsenstein's ML approach to estimating phylogenetic trees.

The COVID-19 pandemic has sparked an urgent need to uncover the underlying biology of this devastating disease. Though RNA viruses mutate more rapidly than DNA viruses, there are a relatively small number of single nucleotide polymorphisms (SNPs) that differentiate the main SARS-CoV-2 lineages that have spread throughout the world. In this study, we investigated 129 RNA-seq data sets and 6928 consensus genomes to contrast the intra-host and inter-host diversity of SARS-CoV-2. Our analyses yielded three major observations. First, the mutational profile of SARS-CoV-2 highlights intra-host single nucleotide variant (iSNV) and SNP similarity, albeit with differences in C > U changes. Second, iSNV and SNP patterns in SARS-CoV-2 are more similar to MERS-CoV than SARS-CoV-1. Third, a significant fraction of insertions and deletions contribute to the genetic diversity of SARS-CoV-2. Altogether, our findings provide insight into SARS-CoV-2 genomic diversity, inform the design of detection tests, and highlight the potential of iSNVs for tracking the transmission of SARS-CoV-2.

The history of modern humans in the Iberian Peninsula includes a variety of population arrivals sometimes presenting admixture with resident populations. Genetic data from current Iberian populations revealed an overall east-west genetic gradient that some authors interpreted as a direct consequence of the Reconquista, where Catholic Kingdoms expanded their territories toward the south while displacing Muslims. However, this interpretation has not been formally evaluated. Here, we present a qualitative analysis of the causes of the current genetic gradient observed in the Iberian Peninsula using extensive spatially explicit computer simulations based on a variety of evolutionary scenarios. Our results indicate that the Neolithic range expansion clearly produces the orientation of the observed genetic gradient. Concerning the Reconquista (including political borders among Catholic Kingdoms and regions with different languages), if modeled upon a previous Neolithic expansion, it effectively favored the orientation of the observed genetic gradient and shows local isolation of certain regions (i.e., Basques and Galicia). Despite additional evolutionary scenarios could be evaluated to more accurately decipher the causes of the Iberian genetic gradient, here we show that this gradient has a more complex explanation than that previously hypothesized.

Alström syndrome (ALMS) is an ultrarare disease with an estimated prevalence lower than 1 in 1,000,000. It is associated with disease-causing mutations in the Alström syndrome 1 () gene, which codifies for a structural protein of the basal body and centrosomes. The symptomatology involves nystagmus, type 2 diabetes mellitus (T2D), obesity, dilated cardiomyopathy (DCM), neurodegenerative disorders and multiorgan fibrosis. We refined the clinical and genetic diagnosis data of 12 patients from 11 families, all of them from Spain. We also studied the allelic frequency of the different variants present in this cohort and performed a haplotype analysis for the most prevalent allele. The genetic analysis revealed 2 novel homozygous variants located in the exon 8, p.(Glu929Ter) and p.(His1808GlufsTer20) in 2 unrelated patients. These 2 novel variants were classified as pathogenic after an experiment (computer analysis). On the other hand, 2 alleles were detected at a high frequency in our cohort: p.(Tyr1714Ter) (25%) and p.(Ser3872TyrfsTer19) (16.7%). The segregation analysis showed that the pathogenic variant p.(Tyr1714Ter) in 3 families is linked to a rare missense polymorphism, p.(Asn1787Asp). In conclusion, 2 novel pathological mutations have been discovered in homozygosis, as well as a probable founder effect in 3 unrelated families.

Macroalgae are sources of bioactive compounds that are interesting from both a chemical and a medical point of view. Although their use in biomedicine has increased significantly in recent years, tests conducted to date have been mostly related to species from temperate latitudes, with the potential application of Antarctic biodiversity being minor. The wide variety of algae species present on Antarctic coastal areas can be a source of new antioxidants. Bearing this in mind, the brown macroalgae (DA) and the red (IC) were selected for the preparation of aqueous extracts with the aim of analyzing their antioxidant activity. This analysis was performed by determining reducing power, total phenolic content, and 2,2-diphenyl-1-picrylhydrazyl free radical scavenging activity. Furthermore, both extracts were employed to synthesize gold and silver nanoparticles. The nanomaterials were fully characterized by means of UV-Visible spectroscopy, transmission electron microscopy, Z potential measurements, and Fourier transform infrared spectroscopy, which confirmed the formation of stable, spherical nanoparticles with mean diameters of 13.7 ± 3.1 and 17.5 ± 3.7 nm for Ag@DA and Ag@IC and 12.6 ± 1.9 and 12.3 ± 1.6 nm for Au@DA and Au@IC. Antioxidant assays were performed after the synthesis of the nanomaterials to evaluate their possible synergistic effect with the extracts. The results suggest that polysaccharides and proteins may play a key role in the process of reduction and stabilization. Finally, for the sake of comparison, the results obtained for the Antarctic macroalgae and have also been considered in the present work.

Research on gold nanoparticles (AuNPs) occupies a prominent place in the field of biomedicine nowadays, being their putative toxicity and bioactivity areas of major concern. The green synthesis of metallic nanoparticles using extracts from marine organisms allows the avoidance of hazardous production steps while maintaining features of interest, thus enabling the exploitation of their promising bioactivity.

Dementia due to Alzheimer's disease (AD) is a complex neurodegenerative disorder, which much of heritability remains unexplained. At the clinical level, one of the most common physiological alterations is the slowing of oscillatory brain activity, measurable by electroencephalography (EEG). Relative power (RP) at the conventional frequency bands (i.e., delta, theta, alpha, beta-1, and beta-2) can be considered as AD endophenotypes.

In this study, we aimed to evaluate the role of ALMS1 in the morphology of primary cilia and regulation of cellular signaling using a knockdown model of the hTERT-RPE1 cell line. ALMS1 depletion resulted in the formation of longer cilia, which often displayed altered morphology as evidenced by extensive twisting and bending of the axoneme. Transforming growth factor beta/bone morphogenetic protein (TGF-β/BMP) signaling, which is regulated by primary cilia, was similarly affected by ALMS1 depletion as judged by reduced levels of TGFβ-1-mediated activation of SMAD2/3. These results provide novel information on the role of ALMS1 in the function of primary cilia and processing of cellular signaling, which when aberrantly regulated may underlie Alström syndrome.

Different methodological approaches are available to assess DNA methylation biomarkers. In this study, we evaluated two sodium bisulfite conversion-dependent methods, namely pyrosequencing and methylation-specific qPCR (MS-qPCR), with the aim of measuring the closeness of agreement of methylation values between these two methods and its effect when setting a cut-off. Methylation of tumor suppressor gene was evaluated in 80 lung cancer patients from which cytological lymph node samples were obtained. Cluster analyses were used to establish methylated and unmethylated groups for each method. Agreement and concordance between pyrosequencing and MS-qPCR was evaluated with Pearson's correlation, Bland-Altman, Cohen's kappa index and ROC curve analyses. Based on these analyses, cut-offs were derived for MS-qPCR. An acceptable correlation (Pearson's R2 = 0.738) was found between pyrosequencing (PYRmean) and MS-qPCR (NMP; normalized methylation percentage), providing similar clinical results when categorizing data as binary using cluster analysis. Compared to pyrosequencing, MS-qPCR tended to underestimate methylation for values between 0 and 15%, while for methylation >30% overestimation was observed. The estimated cut-off for MS-qPCR data based on cluster analysis, kappa-index agreement and ROC curve analysis were much lower than that derived from pyrosequencing. In conclusion, our results indicate that independently of the approach used for estimating the cut-off, the methylation percentage obtained through MS-qPCR is lower than that calculated for pyrosequencing. These differences in data and therefore in the cut-off should be examined when using methylation biomarkers in the clinical practice.

The Mediterranean mussel Mytilus galloprovincialis is an ecologically and economically relevant edible marine bivalve, highly invasive and resilient to biotic and abiotic stressors causing recurrent massive mortalities in other bivalves. Although these traits have been recently linked with the maintenance of a high genetic variation within natural populations, the factors underlying the evolutionary success of this species remain unclear.

Brain waves, measured by electroencephalography (EEG), are a powerful tool in the investigation of neurophysiological traits and a noninvasive and cost-effective alternative in the diagnostic of some neurological diseases. In order to identify novel Quantitative Trait Loci (QTLs) for brain wave relative power (RP), we collected resting state EEG data in five frequency bands (δ, θ, α, β1, and β2) and genome-wide data in a cohort of 105 patients with late onset Alzheimer's disease (LOAD), 41 individuals with mild cognitive impairment and 45 controls from Iberia, correcting for disease status. One novel association was found with an interesting candidate for a role in brain wave biology, (C-type lectin domain family 16), with a variant at this locus passing the adjusted genome-wide significance threshold after Bonferroni correction. This finding reinforces the importance of immune regulation in brain function. Additionally, at a significance cutoff value of 5 × 10, 18 independent association signals were detected. These signals comprise brain expression Quantitative Loci (eQTLs) in caudate basal ganglia, spinal cord, anterior cingulate cortex and hypothalamus, as well as chromatin interactions in adult and fetal cortex, neural progenitor cells and hippocampus. Moreover, in the set of genes showing signals of association with brain wave RP in our dataset, there is an overrepresentation of loci previously associated with neurological traits and pathologies, evidencing the pleiotropy of the genetic variation modulating brain function.

Colorectal cancer (CRC) development is generally accepted as a sequential process, with genetic mutations determining phenotypic tumor progression. However, matching genetic profiles with histological transition requires the analyses of temporal samples from the same patient at key stages of progression.

Pulmonary Arterial Hypertension (PAH) is a rare and fatal disease where knowledge about its genetic basis continues to increase. In this study, we used targeted panel sequencing in a cohort of 624 adult and pediatric patients from the Spanish PAH registry. We identified 11 rare variants in the ATP-binding Cassette subfamily C member 8 (ABCC8) gene, most of them with splicing alteration predictions. One patient also carried another variant in SMAD1 gene (c.27delinsGTAAAG). We performed an ABCC8 in vitro biochemical analyses using hybrid minigenes to confirm the correct mRNA processing of 3 missense variants (c.211C > T p.His71Tyr, c.298G > A p.Glu100Lys and c.1429G > A p.Val477Met) and the skipping of exon 27 in the novel splicing variant c.3394G > A. Finally, we used structural protein information to further assess the pathogenicity of the variants. The results showed 11 novel changes in ABCC8 and 1 in SMAD1 present in PAH patients. After in silico and in vitro biochemical analyses, we classified 2 as pathogenic (c.3288_3289del and c.3394G > A), 6 as likely pathogenic (c.211C > T, c.1429G > A, c.1643C > T, c.2422C > A, c.2694 + 1G > A, c.3976G > A and SMAD1 c.27delinsGTAAAG) and 3 as Variants of Uncertain Significance (c.298G > A, c.2176G > A and c.3238G > A). In all, we show that coupling in silico tools with in vitro biochemical studies can improve the classification of genetic variants.

Alström Syndrome (ALMS) is an ultra-rare multisystem genetic disorder caused by autosomal recessive variants in the ALMS1 gene, which is located on chromosome 2p13. ALMS is a multisystem, progressive disease characterised by visual disturbance, hearing impairment, cardiomyopathy, childhood obesity, extreme insulin resistance, accelerated non-alcoholic fatty liver disease (NAFLD), renal dysfunction, respiratory disease, endocrine and urologic disorders. Clinical symptoms first appear in infancy with great variability in age of onset and severity. ALMS has an estimated incidence of 1 case per 1,000,000 live births and ethnically or geographically isolated populations have a higher-than-average frequency. The rarity and complexity of the syndrome and the lack of expertise can lead to delayed diagnosis, misdiagnosis and inadequate care. Multidisciplinary and multiprofessional teams of experts are essential for the management of patients with ALMS, as early diagnosis and intervention can slow the progression of multi-organ dysfunctions and improve patient quality of life.These guidelines are intended to define standard of care for patients suspected or diagnosed with ALMS of any age. All information contained in this document has originated from a systematic review of the literature and the experiences of the authors in their care of patients with ALMS. The Appraisal of Guidelines for Research & Evaluation (AGREE II) system was adopted for the development of the guidelines and for defining the related levels of evidence and strengths of recommendations.These guidelines are addressed to: a) specialist centres, other hospital-based medical teams and staffs involved with the care of ALMS patients, b) family physicians and other primary caregivers and c) patients and their families.

Pulmonary arterial hypertension is a very infrequent disease, with a variable etiology and clinical expressivity, making sometimes the clinical diagnosis a challenge. Current classification based on clinical features does not reflect the underlying molecular profiling of these groups. The advance in massive parallel sequencing in PAH has allowed for the describing of several new causative and susceptibility genes related to PAH, improving overall patient diagnosis. In order to address the molecular diagnosis of patients with PAH we designed, validated, and routinely applied a custom panel including 21 genes. Three hundred patients from the National Spanish PAH Registry (REHAP) were included in the analysis. A custom script was developed to annotate and filter the variants. Variant classification was performed according to the ACMG guidelines. Pathogenic and likely pathogenic variants have been found in 15% of the patients with 12% of variants of unknown significance (VUS). We have found variants in patients with connective tissue disease (CTD) and congenital heart disease (CHD). In addition, in a small proportion of patients (1.75%), we observed a possible digenic mode of inheritance. These results stand out the importance of the genetic testing of patients with associated forms of PAH (i.e., CHD and CTD) additionally to the classical IPAH and HPAH forms. Molecular confirmation of the clinical presumptive diagnosis is required in cases with a high clinical overlapping to carry out proper management and follow up of the individuals with the disease.

Proteins are commonly used as molecular targets against pathogens such as viruses and bacteria. However, pathogens can evolve rapidly permitting their populations to increase in protein diversity over time and thus escape to the activity of a molecular therapy. Subsequently, in order to design more durable and robust therapies as well as to understand viral evolution in a host and subsequent transmission, it is central to understand the evolution of pathogen proteins. This understanding can enable the detection of protein regions that can be potential targets for therapies and predict the emergence of molecular resistance against therapies. In this direction, two articles published recently in the Journal of Molecular Evolution investigated the evolution of proteomes of diverse flaviviruses, including Zika virus, Dengue virus and West Nile virus. Here I discuss the importance of considering the evolution of viral proteins, with the use of as realistic as possible models and methods that mimic protein evolution, to improve the design of antiviral therapies.

The COVID-19 pandemic has sparked an urgent need to uncover the underlying biology of this devastating disease. Though RNA viruses mutate more rapidly than DNA viruses, there are a relatively small number of single nucleotide polymorphisms (SNPs) that differentiate the main SARS-CoV-2 clades that have spread throughout the world. In this study, we investigated over 7,000 SARS-CoV-2 datasets to unveil both intrahost and interhost diversity. Our intrahost and interhost diversity analyses yielded three major observations. First, the mutational profile of SARS-CoV-2 highlights iSNV and SNP similarity, albeit with high variability in C>T changes. Second, iSNV and SNP patterns in SARS-CoV-2 are more similar to MERS-CoV than SARS-CoV-1. Third, a significant fraction of small indels fuel the genetic diversity of SARS-CoV-2. Altogether, our findings provide insight into SARS-CoV-2 genomic diversity, inform the design of detection tests, and highlight the potential of iSNVs for tracking the transmission of SARS-CoV-2.

Cavalli-Sforza and coauthors originally explored the genetic variation of modern humans throughout the world and observed an overall east-west genetic gradient in Asia. However, the specific environmental and population genetics processes causing this gradient were not formally investigated and promoted discussion in recent studies. Here we studied the influence of diverse environmental and population genetics processes on Asian genetic gradients and identified which could have produced the observed gradient. To do so, we performed extensive spatially-explicit computer simulations of genetic data under the following scenarios: (a) variable levels of admixture between Paleolithic and Neolithic populations, (b) migration through long-distance dispersal (LDD), (c) Paleolithic range contraction induced by the last glacial maximum (LGM), and (d) Neolithic range expansions from one or two geographic origins (the Fertile Crescent and the Yangzi and Yellow River Basins). Next, we estimated genetic gradients from the simulated data and we found that they were sensible to the analysed processes, especially to the range contraction induced by LGM and to the number of Neolithic expansions. Some scenarios were compatible with the observed east-west genetic gradient, such as the Paleolithic expansion with a range contraction induced by the LGM or two Neolithic range expansions from both the east and the west. In general, LDD increased the variance of genetic gradients among simulations. We interpreted the obtained gradients as a consequence of both allele surfing caused by range expansions and isolation by distance along the vast east-west geographic axis of this continent.

The development of High-Throughput Sequencing (HTS) technologies is having a major impact on the genomic analysis of viral populations. Current HTS platforms can capture nucleic acid variation across millions of genes for both selected amplicons and full viral genomes. HTS has already facilitated the discovery of new viruses, hinted new taxonomic classifications and provided a deeper and broader understanding of their diversity, population and genetic structure. Hence, HTS has already replaced standard Sanger sequencing in basic and applied research fields, but the next step is its implementation as a routine technology for the analysis of viruses in clinical settings. The most likely application of this implementation will be the analysis of viral genomics, because the huge population sizes, high mutation rates and very fast replacement of viral populations have demonstrated the limited information obtained with Sanger technology. In this review, we describe new technologies and provide guidelines for the high-throughput sequencing and genetic and evolutionary analyses of viral populations and metaviromes, including software applications. With the development of new HTS technologies, new and refurbished molecular and bioinformatic tools are also constantly being developed to process and integrate HTS data. These allow assembling viral genomes and inferring viral population diversity and dynamics. Finally, we also present several applications of these approaches to the analysis of viral clinical samples including transmission clusters and outbreak characterization.

Our capacity to study individual cells has enabled a new level of resolution for understanding complex biological systems such as multicellular organisms or microbial communities. Not surprisingly, several methods have been developed in recent years with a formidable potential to investigate the somatic evolution of single cells in both healthy and pathological tissues. However, single-cell sequencing data can be quite noisy due to different technical biases, so inferences resulting from these new methods need to be carefully contrasted. Here, I introduce CellCoal, a software tool for the coalescent simulation of single-cell sequencing genotypes. CellCoal simulates the history of single-cell samples obtained from somatic cell populations with different demographic histories and produces single-nucleotide variants under a variety of mutation models, sequencing read counts, and genotype likelihoods, considering allelic imbalance, allelic dropout, amplification, and sequencing errors, typical of this type of data. CellCoal is a flexible tool that can be used to understand the implications of different somatic evolutionary processes at the single-cell level, and to benchmark dedicated bioinformatic tools for the analysis of single-cell sequencing data. CellCoal is available at https://github.com/dapogon/cellcoal.

A collection of the editors of Journal of Molecular Evolution have gotten together to pose a set of key challenges and future directions for the field of molecular evolution. Topics include challenges and new directions in prebiotic chemistry and the RNA world, reconstruction of early cellular genomes and proteins, macromolecular and functional evolution, evolutionary cell biology, genome evolution, molecular evolutionary ecology, viral phylodynamics, theoretical population genomics, somatic cell molecular evolution, and directed evolution. While our effort is not meant to be exhaustive, it reflects research questions and problems in the field of molecular evolution that are exciting to our editors.

Despite the efforts made to reconstruct the history of modern humans, there are still poorly explored regions that are key for understanding the phylogeography of our species. One of them is the Philippines, which is crucial to unravel the colonization of Southeast Asia and Oceania but where little is known about when and how the first humans arrived. In order to shed light into this settlement, we collected samples from 157 individuals of the Philippines with the four grandparents belonging to the same region and mitochondrial variants older than 20,000 years. Next, we analyzed the hypervariable I mtDNA region by approximate Bayesian computation based on extensive spatially explicit computer simulations to select among several migration routes towards the Philippines and to estimate population genetic parameters of this colonization. We found that the colonization of the Philippines occurred more than 60,000 years ago, with long-distance dispersal and from both north and south migration routes. Our results also suggest an environmental scenario especially optimal for humans, with large carrying capacity and population growth, in comparison to other regions of Asia. In all, our study suggests a rapid expansion of modern humans towards the Philippines that could be associated with the establisment of maritime technologies and favorable environmental conditions.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The nonsynonymous/synonymous substitution rate ratio (dN/dS) is a commonly used parameter to quantify molecular adaptation in protein-coding data. It is known that the estimation of dN/dS can be biased if some evolutionary processes are ignored. In this concern, common ML methods to estimate dN/dS assume invariable codon frequencies among sites, despite this characteristic is rare in nature, and it could bias the estimation of this parameter.

ModelTest-NG is a reimplementation from scratch of jModelTest and ProtTest, two popular tools for selecting the best-fit nucleotide and amino acid substitution models, respectively. ModelTest-NG is one to two orders of magnitude faster than jModelTest and ProtTest but equally accurate and introduces several new features, such as ascertainment bias correction, mixture, and free-rate models, or the automatic processing of single partitions. ModelTest-NG is available under a GNU GPL3 license at https://github.com/ddarriba/modeltest , last accessed September 2, 2019.

The knowledge of hereditary predisposition has changed our understanding of Pulmonary Arterial Hypertension. Genetic testing has been widely extended and the application of Pulmonary Arterial Hypertension specific gene panels has allowed its inclusion in the diagnostic workup and increase the diagnostic ratio compared to the traditional sequencing techniques. This is particularly important in the differential diagnosis between Pulmonary Arterial Hypertension and Pulmonary Venoocclusive Disease.

SPLATCHE3 simulates genetic data under a variety of spatially explicit evolutionary scenarios, extending previous versions of the framework. The new capabilities include long-distance migration, spatially and temporally heterogeneous short-scale migrations, alternative hybridization models, simulation of serial samples of genetic data and a large variety of DNA mutation models. These implementations have been applied independently to various studies, but grouped together in the current version.

How and when tumoral clones start spreading to surrounding and distant tissues is currently unclear. Here we leveraged a model-based evolutionary framework to investigate the demographic and biogeographic history of a colorectal cancer. Our analyses strongly support an early monoclonal metastatic colonization, followed by a rapid population expansion at both primary and secondary sites. Moreover, we infer a hematogenous metastatic spread under positive selection, plus the return of some tumoral cells from the liver back to the colon lymph nodes. This study illustrates how sophisticated techniques typical of organismal evolution can provide a detailed, quantitative picture of the complex tumoral dynamics over time and space.

The molecular clock hypothesis, which states that substitutions accumulate in protein sequences at a constant rate, plays a fundamental role in molecular evolution but it is violated when selective or mutational processes vary with time. Such violations of the molecular clock have been widely investigated for protein sequences, but not yet for protein structures. Here, we introduce a novel statistical test (Significant Clock Violations) and perform a large scale assessment of the molecular clock in the evolution of both protein sequences and structures in three large superfamilies. After validating our method with computer simulations, we find that clock violations are generally consistent in sequence and structure evolution, but they tend to be larger and more significant in structure evolution. Moreover, changes of function assessed through Gene Ontology and InterPro terms are associated with large and significant clock violations in structure evolution. We found that almost one third of significant clock violations are significant in structure evolution but not in sequence evolution, highlighting the advantage to use structure information for assessing accelerated evolution and gathering hints of positive selection. Clock violations between closely related pairs are frequently significant in sequence evolution, consistent with the observed time dependence of the substitution rate attributed to segregation of neutral and slightly deleterious polymorphisms, but not in structure evolution, suggesting that these substitutions do not affect protein structure although they may affect stability. These results are consistent with the view that natural selection, both negative and positive, constrains more strongly protein structures than protein sequences. Our code for computing clock violations is freely available at https://github.com/ugobas/Molecular_clock.

Ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota are widespread in marine and terrestrial habitats, playing a major role in the global nitrogen cycle. However, their evolutionary history remains unexplored, which limits our understanding of their adaptation mechanisms. Here, our comprehensive phylogenomic tree of Thaumarchaeota supports three sequential events: origin of AOA from terrestrial non-AOA ancestors, colonization of the shallow ocean, and expansion to the deep ocean. Careful molecular dating suggests that these events coincided with the Great Oxygenation Event around 2300 million years ago (Mya), and oxygenation of the shallow and deep ocean around 800 and 635-560 Mya, respectively. The first transition was likely enabled by the gain of an aerobic pathway for energy production by ammonia oxidation and biosynthetic pathways for cobalamin and biotin that act as cofactors in aerobic metabolism. The first transition was also accompanied by the loss of dissimilatory nitrate and sulfate reduction, loss of oxygen-sensitive pyruvate oxidoreductase, which reduces pyruvate to acetyl-CoA, and loss of the Wood-Ljungdahl pathway for anaerobic carbon fixation. The second transition involved gain of a K transporter and of the biosynthetic pathway for ectoine, which may function as an osmoprotectant. The third transition was accompanied by the loss of the uvr system for repairing ultraviolet light-induced DNA lesions. We conclude that oxygen availability drove the terrestrial origin of AOA and their expansion to the photic and dark oceans, and that the stressors encountered during these events were partially overcome by gene acquisitions from Euryarchaeota and Bacteria, among other sources.

The evaluation of mediastinal lymph nodes is critical for the correct staging of patients with lung cancer (LC). Endobronchial ultrasound-guided transbronchial needle aspiration (EBUS-TBNA) is a minimally invasive technique for mediastinal staging, though unfortunately lymph node micrometastasis is often missed by cytological analysis. The aim of this study was to evaluate the predictive capacity of methylation biomarkers and provide a classification rule for predicting malignancy in false negative EBUS-TBNA samples. The study included 112 patients with a new or suspected diagnosis of LC that were referred to EBUS-TBNA. Methylation of , and was quantified by nested methylation-specific qPCR in 218 EBUS-TBNA lymph node samples. Cross-validated linear regression models were evaluated to predict malignancy. According to EBUS-TBNA and final diagnosis, 90 samples were true positives for malignancy, 110 were true negatives, and 18 were false negatives. , , and were the methylation markers that better predicted malignancy. The model including sex, age, short axis diameter and standard uptake value of adenopathy, and showed 82.7% cross-validated sensitivity and 82.4% specificity for the detection of malignant lymphadenopathies among negative cytology samples. Our results suggest that the predictive model approach proposed can complement EBUS-TBNA for mediastinal staging.

The multiple genetic approaches available for molecular diagnosis of human diseases have made possible to identify an increasing number of pathogenic genetic changes, particularly with the advent of next generation sequencing (NGS) technologies. However, the main challenge lies in the interpretation of their functional impact, which has resulted in the widespread use of animal models. We describe here the functional modelling of seven BBS loci variants, most of them novel, in zebrafish embryos to validate their in silico prediction of pathogenicity. We show that target knockdown (KD) of known BBS (BBS1, BB5 or BBS6) loci leads to developmental defects commonly associated with ciliopathies, as previously described. These KD pleiotropic phenotypes were rescued by co-injecting human wild type (WT) loci sequence but not with the equivalent mutated mRNAs, providing evidence of the pathogenic effect of these BBS changes. Furthermore, direct assessment of cilia located in Kupffer's vesicle (KV) showed a reduction of ciliary length associated with all the studied variants, thus confirming a deleterious effect. Taken together, our results seem to prove the pathogenicity of the already classified and unclassified new BBS variants, as well as highlight the usefulness of zebrafish as an animal model for in vivo assays in human ciliopathies.

The location of the high mountains of southern Europe has been crucial in the phylogeography of most European species, but how extrinsic (topography of sky islands) and intrinsic features (dispersal dynamics) have interacted to shape the genetic structure in alpine restricted species is still poorly known. Here we investigated the mechanisms explaining the colonisation of Cantabrian sky islands in an endemic flightless grasshopper. We scrutinised the maternal genetic variability and haplotype structure, and we evaluated the fitting of two migration models to understand the extant genetic structure in these populations: Long-distance dispersal (LDD) and gradual distance dispersal (GDD). We found that GDD fits the real data better than the LDD model, with an onset of the expansion matching postglacial expansions after the retreat of the ice sheets. Our findings suggest a scenario with small carrying capacity, migration rates, and population growth rates, being compatible with a slow dispersal process. The gradual expansion process along the Cantabrian sky islands found here seems to be conditioned by the suitability of habitats and the presence of alpine corridors. Our findings shed light on our understanding about how organisms which have adapted to live in alpine habitats with limited dispersal abilities have faced new and suitable environmental conditions.

This is the first study to report on the biocompatible and immunogenic properties of one-pot synthesised gold and silver nanoparticles (Au@UI and Ag@UI) using the macroalgae Ulva intestinalis (UI). The UI aqueous extract, Au@UI, and Ag@UI were obtained under sterile conditions and fully characterized by UV-vis spectroscopy, TEM, HRTEM, STEM and FTIR spectroscopy. Moreover, for the first time, the composition of carbohydrates in the UI extract has been reported along with the changes observed after nanoparticle synthesis by size exclusion chromatography, in order to investigate their possible role in the biosynthetic process. This study suggested that the polysaccharide fraction of the extract is involved in the formation and stabilization of the nanoparticles. The potential toxicity of the samples was evaluated using different cell lines and the hemocompatibility was tested in mouse erythrocytes. In addition, ROS production, complement activation and cytokine release were evaluated to determine the immunogenicity. The results showed that Au@UI and Ag@UI exhibit good biocompatibility and hemocompatibility, with the exception of Ag@UI nanoparticles at high concentration, which were hemolytic. The samples induced ROS release and complement activation, two key mechanisms in innate immunity. The samples also induced the release of cytokines from Th1 and Th2 profiles, and other cytokines implicated in the activation of the immune system. Au@UI and Ag@UI were biocompatible and preserved the immunostimulant properties of the UI extract. Hence, Au@UI and Ag@UI could be useful as adjuvants in vaccine development and promote a balanced Th1 and Th2 immune response mediated by ROS production, cytokine release and complement activation.

Bivalve molluscs of the genus Mytilus are considered a model organism in ecotoxicology and are known to be well adapted to marine ecosystems affected by multiple anthropogenic factors, including pollution. In order to assess whether pollution interferes with the reproductive success of Mytilus and affects the diversity within and between populations, we sequenced the transcriptomes of 72 individuals from 9 populations of Mytilus galloprovincialis collected along a ca. 130-km north-south transect on the Western coast of the Iberian Peninsula. We found that polluted areas are acting as a barrier to gene flow, potentially because of the detrimental effect of anthropogenic chemicals on larvae carried from more pristine environments. Furthermore, we observed an increase in genetic diversity in populations from polluted site, which could be indicative of higher mutagenicity driven by the environment. We propose that a microevolutionary perspective is essential for a full characterization of human activities on the dispersal of M. galloprovincialis and that it should be incorporated into management, and conservation plans and policies in the context of the effects of pollution on populations.

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.

Slo1 and Slo3 encode close paralogues of the Slo potassium (K) channels family. Despite their evolutionary relatedness, Slo1 and Slo3 channels show marked functional differences and evolutionary dynamics. Whereas Slo1 is a highly conserved and widely expressed channel, Slo3 is a rapidly evolving channel restricted to sperm. However, the molecular mechanisms behind the structural-functional differences of Slo1 and Slo3 channels are unknown. In this study, we explored the functional divergence of Slo1 and Slo3 subfamilies in vertebrates and examined the structure-function relationships of our predictions using experimental data. We found that ∼25% of sites between Slo1 and Slo3 underwent altered functional constraints, affecting some residues with important roles in Slo1 channel gating. Because functional divergence was principally generated by accelerated evolution of Slo3 after gene duplication, we explored selective forces behind Slo3 diversification. We observed that Slo3 subjected was principally subjected to relaxation of purifying selection, but we also identified several sites evolving under positive selection in the cytosolic domain of this channel . Concerning Slo1, this channel presented strong purifying selection. Whether residues evolving under different selection in Slo1 and Slo3 are responsible for functional differences observed between these channels, as well as among Slo3 orthologs, remains to be established.

Marine bio-resources are being widely studied as an invaluable source of compounds with therapeutic applicability. In particular, macroalgae contain an extended variety of bioactive compounds with different structures and promising biological applications. In this work, Ulva lactuca L. (hereafter UL) was utilyzed for the synthesis of gold and silver nanoparticles. Full characterization by UV-Vis spectroscopy, TEM, HRTEM and STEM miscroscopies, Z Potential and FTIR spectroscopy was performed. The first time in the scientific literature, the composition of carbohydrates of UL extract and their changes observed after nanoparticles synthesis were explored in order to investigate their possible role in the biosynthetic process. The reducing power, total phenolic content and DPPH scavenging activity of UL extract, Au@UL and Ag@UL nanoparticles were determined. The effects of UL extract, Au@UL and Ag@UL were tested in vitro on the colon cancer cell lines HT-29 and Caco-2, on normal primary neonatal dermal fibroblast cell line PCS-201-010, as well as on normal colon cell line CCD-112CoN. Lastly, the apoptotic activity and cellular uptake evaluation was determined for Au@UL and Ag@UL.

To assess the predictive value of some clinical and biochemical parameters, and of the +49 A/G polymorphism of the CTLA-4 gene, for long-term remission following the withdrawal of antithyroid drugs before starting antithyroid drug therapy.

SummaryThe aim of this study was to determine the capacity of the meiotic spindle (MS) to predict embryo implantation and live birth rates. For this purpose, we performed a broad systematic literature search. Of all publications retrieved, only those in which the implantation rates were related to some characteristics of the MS were evaluated. Despite the different methodology used in all the chosen studies, presence of the MS in oocytes was found to be positively associated with embryo implantation. Moreover, high retardance values, as well as strict criteria of normality in the MS structure, are significantly related to higher embryo implantation numbers and live birth rates.

Phylogenetic inference from protein data is traditionally based on empirical substitution models of evolution that assume that protein sites evolve independently of each other and under the same substitution process. However, it is well known that the structural properties of a protein site in the native state affect its evolution, in particular the sequence entropy and the substitution rate. Starting from the seminal proposal by Halpern and Bruno, where structural properties are incorporated in the evolutionary model through site-specific amino acid frequencies, several models have been developed to tackle the influence of protein structure on sequence evolution. Here we describe stability-constrained substitution (SCS) models that explicitly consider the stability of the native state against both unfolded and misfolded states. One of them, the mean-field model, provides an independent sites approximation that can be readily incorporated in maximum likelihood methods of phylogenetic inference, including ancestral sequence reconstruction. Next, we describe its validation with simulated and real proteins and its limitations and advantages with respect to empirical models that lack site specificity. We finally provide guidelines and recommendations to analyze protein data accounting for stability constraints, including computer simulations and inferences of protein evolution based on maximum likelihood. Some practical examples are included to illustrate these procedures.

is a widespread and common rocky intertidal cephalopod that mates readily in the laboratory, but for which mating behavior has not been reported previously. Four sets of behavioral experiments were recorded wherein three males, small, medium & large in varying order, were introduced to each of six females, for a total of 24 individual females and 12 individual males utilized in the experiments. Video analysis shows that successful mating occurred in each of the mount, reach and beak-to-beak positions. Mating was observed for all males, regardless of size relative to the female, or order of introduction. Females showed preference for the first male to which they were introduced in experimental pairings rather than any specific male trait, and mating time increased significantly with increasing female size. Five novel microsatellite markers were developed and used to test paternity in the eleven broods resulting from these experimental pairings. We found skewed paternity in each brood, with early male precedence and male size being the best predictors of parentage. Multiple paternity was observed in every experimental cross but was estimated to be comparatively low in the field, suggesting that sperm limitation might be common in this species. We saw no evidence of direct sperm competition in , but larger males produced significantly more offspring. This study contributes to the growing research on cephalopod mating systems and indicates that octopus mating dynamics might be more variable and complex than thought previously.

Bivalve molluscs of the genus Mytilus are considered a model organism in ecotoxicology and are known to be well adapted to marine ecosystems affected by multiple anthropogenic factors, including pollution. In order to assess whether pollution interferes with the reproductive success of Mytilus and affects the diversity within and between populations, we sequenced the transcriptomes of 72 individuals from 9 populations of Mytilus galloprovincialis collected along a ca. 130-km north-south transect on the Western coast of the Iberian Peninsula. We found that polluted areas are acting as a barrier to gene flow, potentially because of the detrimental effect of anthropogenic chemicals on larvae carried from more pristine environments. Furthermore, we observed an increase in genetic diversity in populations from polluted site, which could be indicative of higher mutagenicity driven by the environment. We propose that a microevolutionary perspective is essential for a full characterization of human activities on the dispersal of M. galloprovincialis and that it should be incorporated into management, and conservation plans and policies in the context of the effects of pollution on populations.

Slo1 and Slo3 encode close paralogues of the Slo potassium (K) channels family. Despite their evolutionary relatedness, Slo1 and Slo3 channels show marked functional differences and evolutionary dynamics. Whereas Slo1 is a highly conserved and widely expressed channel, Slo3 is a rapidly evolving channel restricted to sperm. However, the molecular mechanisms behind the structural-functional differences of Slo1 and Slo3 channels are unknown. In this study, we explored the functional divergence of Slo1 and Slo3 subfamilies in vertebrates and examined the structure-function relationships of our predictions using experimental data. We found that ∼25% of sites between Slo1 and Slo3 underwent altered functional constraints, affecting some residues with important roles in Slo1 channel gating. Because functional divergence was principally generated by accelerated evolution of Slo3 after gene duplication, we explored selective forces behind Slo3 diversification. We observed that Slo3 subjected was principally subjected to relaxation of purifying selection, but we also identified several sites evolving under positive selection in the cytosolic domain of this channel . Concerning Slo1, this channel presented strong purifying selection. Whether residues evolving under different selection in Slo1 and Slo3 are responsible for functional differences observed between these channels, as well as among Slo3 orthologs, remains to be established.

Marine bio-resources are being widely studied as an invaluable source of compounds with therapeutic applicability. In particular, macroalgae contain an extended variety of bioactive compounds with different structures and promising biological applications. In this work, Ulva lactuca L. (hereafter UL) was utilyzed for the synthesis of gold and silver nanoparticles. Full characterization by UV-Vis spectroscopy, TEM, HRTEM and STEM miscroscopies, Z Potential and FTIR spectroscopy was performed. The first time in the scientific literature, the composition of carbohydrates of UL extract and their changes observed after nanoparticles synthesis were explored in order to investigate their possible role in the biosynthetic process. The reducing power, total phenolic content and DPPH scavenging activity of UL extract, Au@UL and Ag@UL nanoparticles were determined. The effects of UL extract, Au@UL and Ag@UL were tested in vitro on the colon cancer cell lines HT-29 and Caco-2, on normal primary neonatal dermal fibroblast cell line PCS-201-010, as well as on normal colon cell line CCD-112CoN. Lastly, the apoptotic activity and cellular uptake evaluation was determined for Au@UL and Ag@UL.

To assess the predictive value of some clinical and biochemical parameters, and of the +49 A/G polymorphism of the CTLA-4 gene, for long-term remission following the withdrawal of antithyroid drugs before starting antithyroid drug therapy.

SummaryThe aim of this study was to determine the capacity of the meiotic spindle (MS) to predict embryo implantation and live birth rates. For this purpose, we performed a broad systematic literature search. Of all publications retrieved, only those in which the implantation rates were related to some characteristics of the MS were evaluated. Despite the different methodology used in all the chosen studies, presence of the MS in oocytes was found to be positively associated with embryo implantation. Moreover, high retardance values, as well as strict criteria of normality in the MS structure, are significantly related to higher embryo implantation numbers and live birth rates.

European genetic gradients of modern humans were initially interpreted as a consequence of the demic diffusion of expanding Neolithic farmers. However, recent studies showed that these gradients may also be influenced by other evolutionary processes such as population admixture or range contractions. Genetic gradients were observed in the Americas, although their specific evolutionary causes were not investigated. Here we extended the approach used to study genetic gradients in Europe to analyze the influence of diverse evolutionary scenarios on American genetic gradients. Using extensive computer simulations, we evaluated the impact of (i) admixture between expansion waves of modern humans, (ii) the presence of ice-sheets during the last glacial maximum (LGM) and (iii) long-distance dispersal (LDD) events, on the genetic gradients (detected by principal component analysis) of the entire continent, North America and South America. The specific simulation of North and South America showed that genetic gradients are usually orthogonal to the direction of range expansions-either expansions from Bering or posterior re-expansions to recolonize northern regions after ice sheets melting-and we suggest that they result from allele surfing processes. Conversely, our results on the entire continent show a northwest-southeast gradient obtained with any scenario, which we interpreted as a consequence of isolation by distance along the long length of the continent. These findings suggest that distinct genetic gradients can be detected at different regions of the Americas and that subcontinent regions present gradients more sensible to evolutionary and environmental factors (such as LDD and the LGM) than the whole continent.

The originally published article contained errors in Fig. 1 (Decision tree for the selection of a suitable NGS genomic simulator), whereby the labels 'Variants' and 'No Variants' had been switched. The correct figure is presented in this notice.

Cancer is a disease driven by both somatic mutations that increase survival and proliferation of cell lineages and the evolution of genes associated with cancer risk in populations. Several genes associated with cancer in humans, hereafter cancer genes, show evidence of germline positive selection among species. Taking advantage of a large collection of mammalian genomes, we systematically looked for signatures of germline positive selection in 430 cancer genes available in COSMIC. We identified 40 cancer genes with a robust signal of positive selection in mammals. We found evidence for fewer selective constraints-higher number of non-synonymous substitutions per non-synonymous site to the number of synonymous substitutions per synonymous site (dN/dS)-and higher incidence of positive selection-more positively selected sites-in cancer genes bearing germline and recessive mutations that predispose to cancer. This finding suggests a potential association between relaxed selection, positive selection, and risk of hereditary cancer. On the other hand, we did not find significant differences in terms of tissue or gene type. Human cancer genes under germline positive selection in mammals are significantly enriched in the processes of DNA repair, with high presence of Fanconi anaemia/Breast Cancer A (FA/BRCA) pathway components and T cell proliferation genes. We also show that the inferred positively selected sites in the two genes with the strongest signal of positive selection, i.e., and , are in regions of functional relevance, which could be relevant to cancer susceptibility.

A reconciliation is an annotation of the nodes of a gene tree with evolutionary events-for example, speciation, gene duplication, transfer, loss, etc.-along with a mapping onto a species tree. Many algorithms and software produce or use reconciliations but often using different reconciliation formats, regarding the type of events considered or whether the species tree is dated or not. This complicates the comparison and communication between different programs.

Selecting among alternative scenarios of human evolution is nowadays a common methodology to investigate the history of our species. This strategy is usually based on computer simulations of genetic data under different evolutionary scenarios, followed by a fitting of the simulated data with the real data. A recent trend in the investigation of ancestral evolutionary processes of modern humans is the application of genetic gradients as a measure of fitting, since evolutionary processes such as range expansions, range contractions, and population admixture (among others) can lead to different genetic gradients. In addition, this strategy allows the analysis of the genetic causes of the observed genetic gradients. Here, we review recent findings on the selection among alternative scenarios of human evolution based on simulated genetic gradients, including pros and cons. First, we describe common methodologies to simulate genetic gradients and apply them to select among alternative scenarios of human evolution. Next, we review previous studies on the influence of range expansions, population admixture, last glacial period, and migration with long-distance dispersal on genetic gradients for some regions of the world. Finally, we discuss this analytical approach, including technical limitations, required improvements, and advice. Although here we focus on human evolution, this approach could be extended to study other species.

Next-generation sequencing (NGS), also known as high-throughput sequencing, is changing the field of microbial genomics research. NGS allows for a more comprehensive analysis of the diversity, structure and composition of microbial genes and genomes compared to the traditional automated Sanger capillary sequencing at a lower cost. NGS strategies have expanded the versatility of standard and widely used typing approaches based on nucleotide variation in several hundred DNA sequences and a few gene fragments (MLST, MLVA, rMLST and cgMLST). NGS can now accommodate variation in thousands or millions of sequences from selected amplicons to full genomes (WGS, NGMLST and HiMLST). To extract signals from high-dimensional NGS data and make valid statistical inferences, novel analytic and statistical techniques are needed. In this review, we describe standard and new approaches for microbial sequence typing at gene and genome levels and guidelines for subsequent analysis, including methods and computational frameworks. We also present several applications of these approaches to some disciplines, namely genotyping, phylogenetics and molecular epidemiology.

Mutation and recombination drive the evolution of most pathogens by generating the genetic variants upon which selection operates. Those variants can, for example, confer resistance to host immune systems and drug therapies or lead to epidemic outbreaks. Given their importance, diverse evolutionary studies have investigated the abundance and consequences of mutation and recombination in pathogen populations. However, some controversies persist regarding the contribution of each evolutionary force to the development of particular phenotypic observations (e.g., drug resistance). In this study, we revise the importance of mutation and recombination in the evolution of pathogens at both intra-host and inter-host levels. We also describe state-of-the-art analytical methodologies to detect and quantify these two evolutionary forces, including biases that are often ignored in evolutionary studies. Finally, we present some of our former studies involving pathogenic taxa where mutation and recombination played crucial roles in the recovery of pathogenic fitness, the generation of interspecific genetic diversity, or the design of centralized vaccines. This review also illustrates several common controversies and pitfalls in the analysis and in the evaluation and interpretation of mutation and recombination outcomes.

Pulmonary arterial hypertension associated with systemic lupus erythematosus (PAH-SLE) is a rare disease with a low incidence rate. In this study, PAH related genes and genetic modifiers were characterised molecularly in patients with PAH-SLE.

Advances in sequencing technologies have made it feasible to obtain massive datasets for phylogenomic inference, often consisting of large numbers of loci from multiple species and individuals. The phylogenomic analysis of next-generation sequencing (NGS) data requires a complex computational pipeline where multiple technical and methodological decisions are necessary that can influence the final tree obtained, like those related to coverage, assembly, mapping, variant calling and/or phasing.

Alström syndrome is a rare autosomal recessive inherited disorder caused by mutations in the ALMS1 gene.

Querying cancer genomes at single-cell resolution is expected to provide a powerful framework to understand in detail the dynamics of cancer evolution. However, given the high costs currently associated with single-cell sequencing, together with the inevitable technical noise arising from single-cell genome amplification, cost-effective strategies that maximize the quality of single-cell data are critically needed. Taking advantage of previously published single-cell whole-genome and whole-exome cancer datasets, we studied the impact of sequencing depth and sampling effort towards single-cell variant detection.

Millepora (Cnidaria: Hydrozoa: Milleporidae) spp. are distributed throughout shallow subtropical and tropical marine environments in the Indo-Pacific and Caribbean-Atlantic, and have traditionally been identified using pore characteristics and colony form. Until now, representatives of Millepora spp. on the island of Okinawa-jima, Japan, have been divided into five species; three branching species (Millepora intricata, M. tenera, M. dichotoma), one species with plate-like morphology (M. platyphylla), and one encrusting species (M. exaesa). There have been only a few reports from the Indo-Pacific that have studied the genetic diversity within Millepora spp., although phylogenetic analyses in the Caribbean-Atlantic have proven useful in delimiting closely-related species, while demonstrating that morphologically-based identification systems may have problems. In the present study, we sought to clarify taxonomic confusion of Millepora spp. in the Pacific by using sequence data of the ribosomal internal transcribed spacer (ITS-rDNA) of specimens from Okinawa, Japan and other localities (Johnston Atoll, Great Barrier Reef). Four separate clades were recovered from the ITS-rDNA analyses. Although we examined specimens of all three branching Millepora spp. previously reported from Okinawa-jima Island, in our phylogenetic analyses they were concentrated within a single clade, with only three specimens in other clades. Encrusting Millepora specimens were found within all clades, although it should be noted all species initially start as encrusting forms, and plate-like specimens were found within three clades. Our data also point to the existence of a previously unknown lineage within Millepora characterized by its ability to overgrow live scleractinian corals.

In this study, we analysed the possible influence of the c.419-43delT BMPR2 variant in patients with Graves' disease (GD), in a molecular basis, focusing our efforts on possible alterations in the mRNA processing and synthesis. The molecular assessment of this variant in patients with GD would shed light on the association between the BMPR2 gene and the disease. The variant was detected in 18%, 55% and 10% of patients with pulmonary arterial hypertension, GD and in general population, respectively. Patients with GD fold change showed increased BMPR2 expression when matched against the controls, with a mean of 4.21 ± 1.73 (P = 0.001); BMPR2 was overexpressed in the analysed cell cycle stages. Fold change analysis of variant carriers and non-carriers showed slight overexpression and differences between phases, but none of them were statistically significant. BMPR2 expression was confirmed in the lymphoblastoid cell lines (LCLs) with a molecular weight of 115 kD, and no differences between variant carriers and non-carriers were detected. To conclude, the BMPR2 variant c.419-19delT appears in high frequency in patients with GD, and independently of its presence, BMPR2 is overexpressed in the LCLs from the GD patients tested. This increase could be paired with the described decreased expression of transforming growth factor-β1 in thyroid tissue from patients with GD.

Protein structures strongly influence molecular evolution. In particular, the evolutionary rate of a protein site depends on the number of its native contacts. Stability-constrained models of protein evolution consider this influence of protein structure on evolution by predicting the effect of mutations on the stability of the native state, but they currently neglect how mutations affect the protein structure. These models predict that buried protein sites with more native contacts are more constrained by natural selection and less variable, as observed. Nevertheless, previous work did not consider the stability against compact misfolded conformations, although it is known that the negative design that destabilizes these misfolded conformations influences protein evolution significantly. Here, we show that stability-constrained models that consider misfolding predict that site-specific sequence entropy and substitution rate peak at amphiphilic sites with an intermediate number of contacts, as these sites are less constrained than exposed sites with few contacts whose hydrophobicity must be limited. This result holds both for a mean-field model with independent sites and for a pairwise model that takes as a reference the wild-type sequence, but it contrasts with the observations that indicate that the entropy and the substitution rate decrease monotonically with the number of contacts. Our work suggests that stability-constrained models overestimate the tolerance of amphiphilic sites against mutations, either because of the limits of the free energy function or, more importantly in our opinion, because they do not consider how mutations perturb the native protein structure.

Inherited syndromic retinopathies are a highly heterogeneous group of diseases that involve retinal anomalies and systemic manifestations. They include retinal ciliopathies, other well-defined clinical syndromes presenting with retinal alterations and cases of non-specific multisystemic diseases. The heterogeneity of these conditions makes molecular and clinical characterization of patients challenging in daily clinical practice. We explored the capacity of targeted resequencing and copy-number variation analysis to improve diagnosis of a heterogeneous cohort of 47 patients mainly comprising atypical cases that did not clearly fit a specific clinical diagnosis. Thirty-three likely pathogenic variants were identified in 18 genes (ABCC6, ALMS1, BBS1, BBS2, BBS12, CEP41, CEP290, IFT172, IFT27, MKKS, MYO7A, OTX2, PDZD7, PEX1, RPGRIP1, USH2A, VPS13B, and WDPCP). Molecular findings and additional clinical reassessments made it possible to accurately characterize 14 probands (30% of the total). Notably, clinical refinement of complex phenotypes was achieved in 4 cases, including 2 de novo OTX2-related syndromes, a novel phenotypic association for the ciliary CEP41 gene, and the co-existence of biallelic USH2A variants and a Koolen-de-Vries syndrome-related 17q21.31 microdeletion. We demonstrate that combining next-generation sequencing and CNV analysis is a comprehensive and useful approach to unravel the extensive phenotypic and genotypic complexity of inherited syndromic retinopathies.

In this work, fungi present in the grapevine's phyllosphere collected from the main demarcated wine regions of Portugal were identified, and their phylogenetic relationships were analysed. A total of 46 vine samples (leaves and berries) were collected from different parts of the country, being isolated a total of 117 fungal colonies that were identified to the genus level and sequenced in the following genetic regions: internal transcribed spacer region and 18S rRNA and β-tubulin gene. Next, a phylogenetic tree reconstruction for each genetic region was built. The isolates retrieved from environmental samples belonged to the genera Alternaria (31%), Cladosporium (21%), Penicillium (19%), Aspergillus (7%) and Epicoccum (3%). No genetic signatures of exchange of genetic material were detected, and consequently, the reconstructed phylogenetic trees allowed to distinguish between these different species/genera. In the fungal composition of the Vitis vinifera phyllosphere, several potential pathogens were identified that can be associated with decreases in crop productivity. Knowledge of fungi identification and genetic diversity is pivotal for the development of more adequate crop management strategies. Furthermore, this information will provide guidelines for a more specific and wiser use of fungicides.

In recent years, the management of food waste processing has emerged as a major concern. One such type of food waste, grape pomace, has been shown to be a great source of bioactive compounds which might be used for more environmentally - friendly processes for the synthesis of nanomaterials. In this study, grape pomace of Vitis vinifera has been used for the obtainment of an aqueous extract. Firstly, the reducing activity, total phenolic content and DPPH scavenging activity of the aqueous extract were determined. Then, the aqueous extract was used for the synthesis of gold and silver nanoparticles. The formation of spherical and stable nanoparticles with mean diameters of 35.3±5.2nm for Au@GP and 42.9±6.4nm for Ag@GP was confirmed by UV-vis spectroscopy and transmission electron microscopy. Furthermore, the functional group of biomolecules present in grape pomace extract, Au@GP and Ag@GP, were characterized by Fourier transform infrared spectroscopy prior to and after the synthesis, in order to obtain information about the biomolecules involved in the reducing and stabilization process. This study is the first to deal with the use of Vitis vinifera grape pomace in obtaining gold and silver nanoparticles through an eco-friendly, quick, one-pot synthetic route.

The diversity of DNA sequencing methods and algorithms for genome assemblies presents scientists with a bewildering array of choices. Here, we construct and compare eight candidate assemblies combining overlapping shotgun read data, mate-pair and Chicago libraries and four different genome assemblers to produce a high-quality draft genome of the iconic coral reef Pearlscale Pygmy Angelfish, Centropyge vrolikii (family Pomacanthidae). The best candidate assembly combined all four data types and had a scaffold N50 127.5 times higher than the candidate assembly obtained from shotgun data only. Our best candidate assembly had a scaffold N50 of 8.97 Mb, contig N50 of 189,827, and 97.4% complete for BUSCO v2 (Actinopterygii set) and 95.6% complete for CEGMA matches. These contiguity and accuracy scores are higher than those of any other fish assembly released to date that did not apply linkage map information, including those based on more expensive long-read sequencing data. Our analysis of how different data types improve assembly quality will help others choose the most appropriate de novo genome sequencing strategy based on resources and target applications. Furthermore, the draft genome of the Pearlscale Pygmy angelfish will play an important role in future studies of coral reef fish evolution, diversity and conservation.

The number of amino acids that occupy a given protein site during evolution reflects the selective constraints operating on the site. This evolutionary variability is strongly influenced by the structural properties of the site in the native structure, and it is quantified either through sequence entropy or through substitution rates. However, while the sequence entropy only depends on the equilibrium frequencies of the amino acids, the substitution rate also depends on the exchangeability matrix that describes mutations in the mathematical model of the substitution process. Here we apply two variants of a mathematical model of protein evolution with selection for protein stability, both against unfolding and against misfolding. Exploiting the approximation of independent sites, these models allow computing site-specific substitution processes that satisfy global constraints on folding stability. We find that site-specific substitution rates do not depend only on the selective constraints acting on the site, quantified through its sequence entropy. In fact, polar sites evolve faster than hydrophobic sites even for equal sequence entropy, as a consequence of the fact that polar amino acids are characterized by higher mutational exchangeability than hydrophobic ones. Accordingly, the model predicts that more polar proteins tend to evolve faster. Nevertheless, these results change if we compare proteins that evolve under different mutation biases, such as orthologous proteins in different bacterial genomes. In this case, the substitution rates are faster in genomes that evolve under mutational bias that favor hydrophobic amino acids by preferentially incorporating the nucleotide Thymine that is more frequent in hydrophobic codons. This appearingly contradictory result arises because buried sites occupied by hydrophobic amino acids are characterized by larger selective factors that largely amplify the substitution rate between hydrophobic amino acids, while the selective factors of exposed sites have a weaker effect. Thus, changes in the mutational bias produce deep effects on the biophysical properties of the protein (hydrophobicity) and on its evolutionary properties (sequence entropy and substitution rate) at the same time. The program Prot_evol that implements the two site-specific substitution processes is freely available at https://ub.cbm.uam.es/prot_fold_evol/prot_fold_evol_soft_main.php#Prot_Evol.