SPEAQeasy: a scalable pipeline for expression analysis and quantification for R/bioconductor-powered RNA-seq analyses.

BACKGROUND: RNA sequencing (RNA-seq) is a common and widespread biological assay, and an increasing amount of data is generated with it. In practice, there are a large number of individual steps a researcher must perform before raw RNA-seq reads yield directly valuable information, such as differential gene expression data. Existing software tools are typically specialized, only performing one step-such as alignment of reads to a reference genome-of a larger workflow. The demand for a more comprehensive and reproducible workflow has led to the production of a number of publicly available RNA-seq pipelines. However, we have found that most require computational expertise to set up or share among several users, are not actively maintained, or lack features we have found to be important in our own analyses. RESULTS: In response to these concerns, we have developed a Scalable Pipeline for Expression Analysis and Quantification (SPEAQeasy), which is easy to install and share, and provides a bridge towards R/Bioconductor downstream analysis solutions. SPEAQeasy is portable across computational frameworks (SGE, SLURM, local, docker integration) and different configuration files are provided ( http://research.libd.org/SPEAQeasy/ ). CONCLUSIONS: SPEAQeasy is user-friendly and lowers the computational-domain entry barrier for biologists and clinicians to RNA-seq data processing as the main input file is a table with sample names and their corresponding FASTQ files. The goal is to provide a flexible pipeline that is immediately usable by researchers, regardless of their technical background or computing environment.

Label-free quantitative proteomics in serum reveals candidate biomarkers associated with low bone mineral density in Mexican postmenopausal women.

Postmenopausal osteoporosis is a public health problem leading to an increased risk of fractures, negatively impacting women's health. The absence of sensitive and specific biomarkers for early detection of osteoporosis represents a substantial challenge for improving patient management. Herein, we aimed to identify potential candidate proteins associated with low bone mineral density (BMD) in postmenopausal women from the Mexican population. Serum samples from postmenopausal women (40 with normal BMD, 40 with osteopenia (OS), and 20 with osteoporosis (OP)) were analyzed by label-free LC-MS/MS quantitative proteomics. Proteome profiling revealed significant differences between the OS and OP groups compared to individuals with normal BMD. A quantitative comparison of proteins between groups indicated 454 differentially expressed proteins (DEPs). Compared to normal BMD, 14 and 214 DEPs were found in OS and OP groups, respectively, while 226 DEPs were identified between OS and OP groups. The protein-protein interaction and enrichment analysis of DEPs were closely linked to the bone mineral content, skeletal morphology, and immune response activation. Based on their role in bone metabolism, a panel of 12 candidate biomarkers was selected, of which 1 DEP (RYR1) was found upregulated in the OS and OP groups, 8 DEPs (APOA1, SHBG, FETB, MASP1, PTK2B, KNG1, GSN, and B2M) were upregulated in OP and 3 DEPs (APOA2, RYR3, and HBD) were downregulated in OS or OP. The proteomic analysis described here may help discover new and potentially non-invasive biomarkers for the early diagnosis of osteoporosis in postmenopausal women.

F-actin nanostructures rearrangements and regulation are essential for SARS-CoV-2 particle production in host pulmonary cells.

Our study focused on deciphering the role of F-actin and related regulatory factors during SARS-CoV-2 particle production and transmission in human pulmonary cells. Quantitative high-resolution microscopies revealed that the late phases of SARS-CoV-2 infection induce a strong rearrangement of F-actin nanostructures dependent on the viral M, E, and N structural proteins. Intracellular vesicles containing viral components are labeled with Rab7 and Lamp1 and are surrounded by F-actin ring-shaped structures, suggesting their role in viral trafficking toward the cell membrane for virus release. Furthermore, filopodia-like nanostructures were loaded with viruses, potentially facilitating their egress and transmission between lung cells. Gene expression analysis revealed the involvement of alpha-actinins under the regulation of the protein kinase N (PKN). The use of a PKN inhibitor efficiently reduces virus particle production, restoring endoplasmic reticulum and F-actin cellular shape. Our results highlight an important role of F-actin rearrangements during the productive phases of SARS-CoV-2 particles.

Correction: Whole genome variation in 27 Mexican indigenous populations, demographic and biomedical insights.

[This corrects the article DOI: 10.1371/journal.pone.0249773.].

Whole genome variation in 27 Mexican indigenous populations, demographic and biomedical insights.

There has been limited study of Native American whole genome diversity to date, which impairs effective implementation of personalized medicine and a detailed description of its demographic history. Here we report high coverage whole genome sequencing of 76 unrelated individuals, from 27 indigenous groups across Mexico, with more than 97% average Native American ancestry. On average, each individual has 3.26 million Single Nucleotide Variants and short indels, that together comprise a catalog of 9,737,152 variants, 44,118 of which are novel. We report 497 common Single Nucleotide Variants (with allele frequency > 5%) mapped to drug responses and 316,577 in enhancer or promoter elements; interestingly we found some of these enhancer variants in PPARG, a nuclear receptor involved in highly prevalent health problems in Mexican population, such as obesity, diabetes, and insulin resistance. By detecting signals of positive selection we report 24 enriched key pathways under selection, most of them related to immune mechanisms. No missense variants in ACE2, the receptor responsible for the entry of the SARS CoV-2 virus, were found in any individual. Population genomics and phylogenetic analyses demonstrated stratification in a Northern-Central-Southern axis, with major substructure in the Central region. The Seri, a northern group with the most genetic divergence in our study, showed a distinctive genomic context with the most novel variants, and the most population specific genotypes. Genome-wide analysis showed that the average haplotype blocks are longer in Native Mexicans than in other world populations. With this dataset we describe previously undetected population level variation in Native Mexicans, helping to reduce the gap in genomic data representation of such groups.