Transvaginal ultrasound-guided oocyte retrieval in cattle: State-of-the-art and its impact on the in vitro fertilization embryo production outcome.

Transvaginal ultrasound-guided oocyte retrieval (commonly called OPU) and in vitro embryo production (IVP) in cattle has shown significant progress in recent years, in part, as a result of a better understanding of the full potential of these tools by end users. The combination of OPU and IVP (OPU-IVP) has been successfully and widely commercially used worldwide. The main advantages are a greater number of embryos and pregnancies per unit of time, faster genetic progress due to donor quick turn around and more elite sires mating combinations, larger spectrum of female age (calves, prepuberal, heifer, cow) and condition (open, pregnant) from which to retrieve oocytes, a reduced number of sperm (even sexed) required to fertilize the oocytes, among other benefits. OPU-IVP requires significant less donor preparation in comparison to conventional embryo transfer (<50% of usual FSH injections needed) to the extent of no stimulating hormones (FSH) are necessary. Donor synchronization, stimulation, OPU technique, oocyte competence, embryo performance, and its impact on cryopreservation and pregnancy are discussed.

Quantification and Characterization of CTCs and Clusters in Pancreatic Cancer by Means of the Hough Transform Algorithm.

Circulating Tumor Cells (CTCs) are considered a prognostic marker in pancreatic cancer. In this study we present a new approach for counting CTCs and CTC clusters in patients with pancreatic cancer using the IsofluxTM System with the Hough transform algorithm (Hough-IsofluxTM). The Hough-IsofluxTM approach is based on the counting of an array of pixels with a nucleus and cytokeratin expression excluding the CD45 signal. Total CTCs including free and CTC clusters were evaluated in healthy donor samples mixed with pancreatic cancer cells (PCCs) and in samples from patients with pancreatic ductal adenocarcinoma (PDAC). The IsofluxTM System with manual counting was used in a blinded manner by three technicians who used Manual-IsofluxTM as a reference. The accuracy of the Hough-IsofluxTM approach for detecting PCC based on counted events was 91.00% [84.50, 93.50] with a PCC recovery rate of 80.75 ± 16.41%. A high correlation between the Hough-IsofluxTM and Manual-IsofluxTM was observed for both free CTCs and for clusters in experimental PCC (R2 = 0.993 and R2 = 0.902 respectively). However, the correlation rate was better for free CTCs than for clusters in PDAC patient samples (R2 = 0.974 and R2 = 0.790 respectively). In conclusion, the Hough-IsofluxTM approach showed high accuracy for the detection of circulating pancreatic cancer cells. A better correlation rate was observed between Hough-IsofluxTM approach and with the Manual-IsofluxTM for isolated CTCs than for clusters in PDAC patient samples.

Liquid biopsy for KRAS, NRAS and BRAF mutation testing in advanced colorectal cancer patients: the Argentinean experience.

Objective: To analyze the frequency of KRAS, NRAS and BRAF hotspot mutations in circulating tumor DNA (ctDNA) from patients with metastatic colorectal cancer (mCRC). Methods: Observational, descriptive and retrospective study in mCRC patients with available ctDNA-based genotype of KRAS, NRAS and BRAF. Results: The frequencies of plasma mutations for KRAS, NRAS and BRAF were 34% (± 7), 4% (± 3) and 4% (± 3), respectively. Median overall survival of plasma-tested RAS/BRAF-mutated patients was 26.6 months (95% CI: 14.4-not estimable [NE]), while RAS/BRAF wild-type patients did not reach the median survival during follow-up. Median progression-free survival for RAS/BRAF wild-type and RAS/BRAF-mutated patients was 12 (95% CI: 7-NE) and 4 months (95% CI: 4-NE), respectively. Conclusion: Our work supports the utility of KRAS, NRAS and BRAF analysis in liquid biopsy from mCRC patients.

Candidate pharmacological treatments for substance use disorder and suicide identified by gene co-expression network-based drug repositioning.

Patients with substance use disorders (SUD) are at high risk to die by suicide. So far, the neurobiology of the suicide-SUD association has not been elucidated. This study aimed to identify potential pharmacological targets among hub genes from brain gene co-expression networks of individuals with SUD in a suicidal and non-suicidal context. Post-mortem samples from the prefrontal cortex of 79 individuals were analyzed. Individuals were classified into the following groups: suicides with SUD (n = 28), suicides without SUD (n = 23), nonsuicides with SUD (n = 9), nonsuicides without SUD (n = 19). Gene expression profiles were evaluated with the Illumina HumanHT-12 v4 array. Co-expression networks were constructed in WGCNA using the differentially expressed genes found in the comparisons: (a) suicides with and without SUD and (b) nonsuicides with and without SUD. Hub genes were selected for drug-gene interaction testing in the DGIdb database. Among drugs interacting with hub genes in suicides we found MAOA inhibitors and dextromethorphan. In the nonsuicide individuals, we found interactions with eglumegad and antipsychotics (olanzapine, clozapine, loxapine). Modafinil was found to interact with genes in both suicides and nonsuicides. These drugs represent possible candidate treatments for patients with SUD with and without suicidal behavior and their study in each context is encouraged.

Gender imbalance in executive management positions at the Mexican National Institutes of Health.

BACKGROUND: Around the world, there is a significant difference in the proportion of women with access to leadership in healthcare with respect to men. This article studies gender imbalance and wage gap in managerial, executive, and directive job positions at the Mexican National Institutes of Health. METHODS: Cohort data were described using a visual circular representation and modeled using a generalized linear model. Analysis of variance was used to assess model significance, and posterior Fisher's least significant differences were analyzed when appropriate. RESULTS: This study demonstrated that there is a gender imbalance distribution among the hierarchical position at the Mexican National Health Institutes and also exposed that the wage gap exists mainly in the (highest or lowest) ranks in hierarchical order. CONCLUSIONS: Since the majority of the healthcare workforce is female, Mexican women are still underrepresented in executive and directive management positions at national healthcare organizations.

Brain Gene Expression Profiling of Individuals With Dual Diagnosis Who Died by Suicide.

Objective: Dual diagnosis (DD) is the co-occurrence of at least one substance use disorder and one or more mental disorders in a given individual. Despite this comorbidity being highly prevalent and associated with adverse clinical outcomes, its neurobiology remains unclear. Furthermore, patients with DD are at higher risk for suicidal behavior in comparison with single disorder patients. Our objective was to evaluate brain gene expression patterns in individuals with DD who died by suicide. Methods: We compared the gene expression profile in the dorsolateral prefrontal cortex of suicides with DD (n = 10) to the transcriptome of suicides with substance use disorder alone (n = 10), suicides with mood disorders (MD) alone (n = 13), and suicides without mental comorbidities (n = 5). Gene expression profiles were assessed by microarrays. In addition, we performed a brain cell type enrichment to evaluate whether the gene expression profiles could reflect differences in cell type compositions among the groups. Results: When comparing the transcriptome of suicides with DD to suicides with substance use disorder alone and suicides with MD alone, we identified 255 and 172 differentially expressed genes (DEG), respectively. The overlap of DEG between both comparisons (112 genes) highlighted the presence of common disrupted pathways in substance use disorder and MD. When comparing suicides with DD to suicides without mental comorbidities, we identified 330 DEG, mainly enriched in neurogenesis. Cell type enrichment indicated higher levels of glial markers in suicides with DD compared to the other groups. Conclusions: Suicides with DD exhibited a gene expression profile distinct from that of suicides with a single disorder, being substance use disorder or MD, and suicides without mental disorders. Our results suggest alteration in the expression of genes involved in glial specific markers, glutamatergic and GABAergic neurotransmission in suicides with DD compared to suicides with a single disorder and suicides without mental comorbidities. Alterations in the expression of synaptic genes at different levels were found in substance use disorder and MD.

A novel non-parametric method for uncertainty evaluation of correlation-based molecular signatures: its application on PAM50 algorithm.

Motivation: The PAM50 classifier is used to assign patients to the highest correlated breast cancer subtype irrespectively of the obtained value. Nonetheless, all subtype correlations are required to build the risk of recurrence (ROR) score, currently used in therapeutic decisions. Present subtype uncertainty estimations are not accurate, seldom considered or require a population-based approach for this context. Results: Here we present a novel single-subject non-parametric uncertainty estimation based on PAM50's gene label permutations. Simulations results ( n = 5228) showed that only 61% subjects can be reliably 'Assigned' to the PAM50 subtype, whereas 33% should be 'Not Assigned' (NA), leaving the rest to tight 'Ambiguous' correlations between subtypes. The NA subjects exclusion from the analysis improved survival subtype curves discrimination yielding a higher proportion of low and high ROR values. Conversely, all NA subjects showed similar survival behaviour regardless of the original PAM50 assignment. We propose to incorporate our PAM50 uncertainty estimation to support therapeutic decisions. Availability and Implementation: Source code can be found in 'pbcmc' R package at Bioconductor. Contacts: cristobalfresno@gmail.com or efernandez@bdmg.com.ar. Supplementary information: Supplementary data are available at Bioinformatics online.

TarSeqQC: Quality control on targeted sequencing experiments in R.

Targeted sequencing (TS) is growing as a screening methodology used in research and medical genetics to identify genomic alterations causing human diseases. In general, a list of possible genomic variants is derived from mapped reads through a variant calling step. This processing step is usually based on variant coverage, although it may be affected by several factors. Therefore, undercovered relevant clinical variants may not be reported, affecting pathology diagnosis or treatment. Thus, a prior quality control of the experiment is critical to determine variant detection accuracy and to avoid erroneous medical conclusions. There are several quality control tools, but they are focused on issues related to whole-genome sequencing. However, in TS, quality control should assess experiment, gene, and genomic region performances based on achieved coverages. Here, we propose TarSeqQC R package for quality control in TS experiments. The tool is freely available at Bioconductor repository. TarSeqQC was used to analyze two datasets; low-performance primer pools and features were detected, enhancing the quality of experiment results. Read count profiles were also explored, showing TarSeqQC's effectiveness as an exploration tool. Our proposal may be a valuable bioinformatic tool for routinely TS experiments in both research and medical genetics.

The ROCK Inhibitor Fasudil Prevents Chronic Restraint Stress-Induced Depressive-Like Behaviors and Dendritic Spine Loss in Rat Hippocampus.

Background: Dendritic arbor simplification and dendritic spine loss in the hippocampus, a limbic structure implicated in mood disorders, are assumed to contribute to symptoms of depression. These morphological changes imply modifications in dendritic cytoskeleton. Rho GTPases are regulators of actin dynamics through their effector Rho kinase. We have reported that chronic stress promotes depressive-like behaviors in rats along with dendritic spine loss in apical dendrites of hippocampal pyramidal neurons, changes associated with Rho kinase activation. The present study proposes that the Rho kinase inhibitor Fasudil may prevent the stress-induced behavior and dendritic spine loss. Methods: Adult male Sprague-Dawley rats were injected with saline or Fasudil (i.p., 10 mg/kg) starting 4 days prior to and maintained during the restraint stress procedure (2.5 h/d for 14 days). Nonstressed control animals were injected with saline or Fasudil for 18 days. At 24 hours after treatment, forced swimming test, Golgi-staining, and immuno-western blot were performed. Results: Fasudil prevented stress-induced immobility observed in the forced swimming test. On the other hand, Fasudil-treated control animals showed behavioral patterns similar to those of saline-treated controls. Furthermore, we observed that stress induced an increase in the phosphorylation of MYPT1 in the hippocampus, an exclusive target of Rho kinase. This change was accompanied by dendritic spine loss of apical dendrites of pyramidal hippocampal neurons. Interestingly, increased pMYPT1 levels and spine loss were both prevented by Fasudil administration. Conclusion: Our findings suggest that Fasudil may prevent the development of abnormal behavior and spine loss induced by chronic stress by blocking Rho kinase activity.

Effect of spermatozoa motility hyperactivation factors and gamete coincubation duration on in vitro bovine embryo development using flow cytometrically sorted spermatozoa.

The aim of the present study was to evaluate the effects of sperm motility enhancers and different IVF times on cleavage, polyspermy, blastocyst formation, embryo quality and hatching ability. In Experiment 1, sex-sorted X chromosome-bearing Bos taurus spermatozoa were incubated for 30min before 18h fertilisation with hyperactivating factors, namely 10mM caffeine (CA), 5mM theophylline (TH), 10mM caffeine and 5mM theophylline (CA+TH); and untreated spermatozoa (control). In Experiment 2, matured B. taurus oocytes were fertilised using a short (8h) or standard (18h) fertilisation length, comparing two different fertilisation media, namely synthetic oviducal fluid (SOF) fertilisation medium (SOF-FERT) and M199 fertilisation medium (M199-FERT). Cleavage and blastocyst formation rates were significantly higher in the CA+TH group (77% and 27%, respectively) compared with the control group (71% and 21%, respectively). Cleavage rates and blastocyst formation were significantly lower for the shortest fertilisation time (8h) in M199-FERT medium (42% and 12%, respectively). The SOF-FERT medium with an 8h fertilisation time resulted in the highest cleavage rates and blastocyst formation (74% and 29%, respectively). The SOF-FERT medium produced the highest embryo quality (50% Grade 1) and hatching rate (66%). Motility enhancers did not affect polyspermy rates, whereas polyspermy was affected when fertilisation length was extended from 8h (3%) to 18h (9%) and in M199-FERT (14%) compared with SOF-FERT (6%). We conclude that adding the motility enhancers CA and TH to sex sorted spermatozoa and Tyrode's albumin lactate pyruvate (TALP)-Sperm can improve cleavage and embryo development rates without increasing polyspermy. In addition, shortening the oocyte-sperm coincubation time (8h) resulted in similar overall embryo performance rates compared with the prolonged (18h) interval.

Comparison of different fertilisation media for an in vitro maturation?fertilisation?culture system using flow-cytometrically sorted X chromosome-bearing spermatozoa for bovine embryo production.

High demand exists among commercial cattle producers for in vitro-derived bovine embryos fertilised with female sex-sorted spermatozoa from high-value breeding stock. The aim of this study was to evaluate three fertilisation media, namely M199, synthetic oviductal fluid (SOF) and Tyrode's albumin-lactate-pyruvate (TALP), on IVF performance using female sex-sorted spermatozoa. In all, 1143, 1220 and 1041 cumulus-oocyte complexes were fertilised in M199, SOF and TALP, respectively. There were significant differences among fertilisation media (P < 0.05) in cleavage rate (M199 = 57%, SOF = 71% and TALP = 72%), blastocyst formation (M199 = 9%, SOF = 20% and TALP = 19%), proportion of Grade 1 blastocysts (M199 = 15%, SOF = 52% and TALP = 51%), proportion of Grade 3 blastocysts (M199 = 58%, SOF = 21% and TALP = 20%) and hatching rates (M199 = 29%, SOF = 60% and TALP = 65%). The inner cell mass (ICM) and trophectoderm (TE) cells of Day 7 blastocysts were also affected by the fertilisation medium. Embryos derived from SOF and TALP fertilisation media had higher numbers of ICM, TE and total cells than those fertilised in M199. In conclusion, fertilisation media affected cleavage rate, as well as subsequent embryo development, quality and hatching ability. SOF and TALP fertilisation media produced significantly more embryos of higher quality than M199.

RDAVIDWebService: a versatile R interface to DAVID.

SUMMARY: The RDAVIDWebService package provides a class-based interface from R programs/scripts to fully access/control the database for annotation, visualization and integrated discovery, without the need for human interaction on its Web site (http://david.abcc.ncifcrf.gov). The library enhances the database for annotation, visualization and integrated discovery capabilities for Gene Ontology analysis by means of GOstats-based direct acyclic graph conversion methods, in addition to the usual many-genes-to-many-terms visualization. AVAILABILITY AND IMPLEMENTATION: RDAVIDWebService is available as an R package from the Bioconductor project (www.bioconductor.org) and on the authors' Web site (www.bdmg.com.ar) under GPL-2 license, subjected to the terms of use of DAVID (http://david.abcc.ncifcrf.gov/content.jsp?file=WS.html). CONTACT: cfresno@bdmg.com.ar or efernandez@bdmg.com.ar.

Molecular pathology of acute kidney injury in a choline-deficient model and fish oil protective effect.

PURPOSE: The aim of this work was to investigate the potential protective effects of fish oil on the basis of kidney transcriptomic data on a nutritional experimental model. METHODS: Male weanling Wistar rats were divided into four groups and fed choline-deficient (CD) and choline-supplemented (CS) diets with vegetable oil (VO) and menhaden oil (MO): CSVO, CDVO, CSMO and CDMO. Animals were killed after receiving the diets for 6 days. Total RNA was purified from the right kidney and hybridized to Affymetrix GeneChip Rat Gene 1.0 ST Array. Differentially expressed genes were analyzed. RESULTS: All CSVO, CSMO and CDMO rats showed no renal alterations, while all CDVO rats showed renal cortical necrosis. A thorough analysis of the differential expression between groups CSMO and CDMO was carried out. There were no differential genes for p < 0.01. The analysis of the differential expression between groups CSVO and CSMO revealed 32 genes, 11 were over-expressed and 21 were under-expressed in CSMO rats. CONCLUSIONS: This work was part of a large set of experiments and was used in a hypothesis-generating manner. The comprehensive analysis of genetic expression allowed confirming that menhaden oil has a protective effect on this nutritional experimental model and identifying 32 genes that could be responsible for that protection, including Gstp1. These results reveal that gene changes could play a role in renal injury.

Clinical implementation of a cervical cancer screening program via co-testing at a university hospital.

The Human Papillomavirus (HPV) test is a crucial technology for cervical cancer prevention because it enables programs to identify women with high-risk HPV infection who are at risk of developing cervical cancer. Current U.S. Preventive Services Task Force recommendations include cervical cancer screening every three years with cervical cytology alone or every five years with either high-risk HPV testing alone or high-risk HPV testing combined with cytology (co-testing). In Argentina, 7,548 new cervical cancer cases are diagnosed each year with 3,932 deaths attributed to this cause. Our study aims to show the clinical implementation of a cervical cancer screening program by concurrent HPV testing and cervical cytology (co-testing); and to evaluate the possible cervical cancer screening scenarios for Latin America, focusing on their performance and average cost. A cervical cancer screening five year program via co-testing algorithm (Hybrid-2-Capture/cytology) was performed on women aged 30-65 years old at a university hospital. Statistical analysis included a multinomial logistic regression, and two cancer screening classification alternatives were tested (cytology-reflex and HPV-reflex). A total of 2,273 women were included, 91.11% of the participants were double-negative, 2.55% double-positive, 5.90% positive-Hybrid-2-Capture-/negative-cytology, and 0.44% negative-Hybrid-2-Capture/positive-cytology. A thorough follow-up was performed in the positive-Hybrid-2-Capture group. Despite our efforts, 21 (10.93%) were lost, mainly because of changes on their health insurance coverage which excluded them from our screening algorithm. Of the 171 women with positive-Hybrid-2-Capture results and follow-up, 68 (39.77%) cleared the virus infection, 64 (37.43%) showed viral persistence, and 39 (22.81%) were adequately treated after detection via colposcopy/biopsy of histological HSIL (High-Grade Squamous Intraepithelial Lesion). The prevalence of high-risk HPV in this population was 192 women (8.45%), with HSIL histology detection rates of 17 per 1,000 screened women. A multinomial logistic regression analysis was performed over the women with positive-Hybrid-2-Capture considering the follow up (clearance, persistence and HSIL) as dependent variable, and the cytology test results (positive- or negative-cytology and Atypical Squamous Cells of Undetermined Significance, ASC-US) as independent variable. The model supported a direct association between cytology test results and follow up: negative-cytology/clearance, ASC-US/persistence, and positive-cytology/HSIL with the following probabilities of occurrence for these pairs 0.5, 0.647 and 0.647, respectively. Cytology could be considered a prognostic-factor in women with a positive-Hybrid-2-Capture. These findings suggest that the introduction of co-testing could diminish the burden of cervical cancer in low-and middle-income-countries, acting as a tool against inequity in healthcare.

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

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

Rapid Identification of Drug Resistance and Phylogeny in M. tuberculosis, Directly from Sputum Samples.

Tuberculosis (TB) remains one of the most important infectious diseases globally. Establishing a resistance profile from the initial TB diagnosis is a priority. Rapid molecular tests evaluate only the most common genetic variants responsible for resistance to certain drugs, and Whole Genome Sequencing (WGS) needs culture prior to next-generation sequencing (NGS), limiting their clinical value. Targeted sequencing (TS) from clinical samples avoids these drawbacks, providing a signature of genetic markers that can be associated with drug resistance and phylogeny. In this study, a proof-of-concept protocol was developed for detecting genomic variants associated with drug resistance and for the phylogenetic classification of Mycobacterium Tuberculosis (Mtb) in sputum samples. Initially, a set of Mtb reference strains from the WHO were sequenced (WGS and TS). The results from the protocol agreed >95% with WHO reported data and phenotypic drug susceptibility testing (pDST). Lineage genetics results were 100% concordant with those derived from WGS. After that, the TS protocol was applied to sputum samples from TB patients to detect resistance to first- and second-line drugs and derive phylogeny. The accuracy was >90% for all evaluated drugs, except Eto/Pto (77.8%), and 100% were phylogenetically classified. The results indicate that the described protocol, which affords the complete drug resistance profile and phylogeny of Mtb from sputum, could be useful in the clinical area, advancing toward more personalized and more effective treatments in the near future. IMPORTANCE The COVID-19 pandemic negatively affected the progress in accessing essential Tuberculosis (TB) services and reducing the burden of TB disease, resulting in a decreased detection of new cases and increased deaths. Generating molecular diagnostic tests with faster results without losing reliability is considered a priority. Specifically, developing an antimicrobial resistance profile from the initial stages of TB diagnosis is essential to ensure appropriate treatment. Currently available rapid molecular tests evaluate only the most common genetic variants responsible for resistance to certain drugs, limiting their clinical value. In this work, targeted sequencing on sputum samples from TB patients was used to identify Mycobacterium tuberculosis mutations in genes associated with drug resistance and to derive a phylogeny of the infecting strain. This protocol constitutes a proof-of-concept toward the goal of helping clinicians select a timely and appropriate treatment by providing them with actionable information beyond current molecular approaches.

Chromatin topology defines estradiol-primed progesterone receptor and PAX2 binding in endometrial cancer cells.

Estrogen (E2) and Progesterone (Pg), via their specific receptors (ERalpha and PR), are major determinants in the development and progression of endometrial carcinomas, However, their precise mechanism of action and the role of other transcription factors involved are not entirely clear. Using Ishikawa endometrial cancer cells, we report that E2 treatment exposes a set of progestin-dependent PR binding sites which include both E2 and progestin target genes. ChIP-seq results from hormone-treated cells revealed a non-random distribution of PAX2 binding in the vicinity of these estrogen-promoted PR sites. Altered expression of hormone regulated genes in PAX2 knockdown cells suggests a role for PAX2 in fine-tuning ERalpha and PR interplay in transcriptional regulation. Analysis of long-range interactions by Hi-C coupled with ATAC-seq data showed that these regions, that we call 'progestin control regions' (PgCRs), exhibited an open chromatin state even before hormone exposure and were non-randomly associated with regulated genes. Nearly 20% of genes potentially influenced by PgCRs were found to be altered during progression of endometrial cancer. Our findings suggest that endometrial response to progestins in differentiated endometrial tumor cells results in part from binding of PR together with PAX2 to accessible chromatin regions. What maintains these regions open remains to be studied.

A meta-study on transcription factor networks in the suicidal brain.

There is evidence supporting the presence of brain gene expression differences between suicides and non-suicides. Such differences have been implicated in suicide pathophysiology. However, regulatory factors underlying these gene expression differences have not been fully understood. Therefore, the identification of differences in regulatory mechanisms, i.e., transcriptional factors between suicides and non-suicides is crucial for the understanding of suicide neurobiology. In this study, we conducted a transcription factor network meta-study with freely available data from the prefrontal cortex of suicides and non-suicides with different mental disorders, including major depression disorder, bipolar disorder and schizophrenia, as well as healthy controls. Disorder-specific characteristics of suicides and non-suicides transcription factor networks were detected, i.e., the presence of immune response genes in both suicides and non-suicides with major depression disorder networks. Also, we found the presence of ESR1, which has been implicated to give resilience to social stress, in the non-suicides network but not in the suicides with major depression network. Suicides and non-suicides with bipolar disorder shared only three genes in common: FOS, CRY1 and PER2. In addition, we found a higher number of genes involved in immune response in the non-suicides with bipolar disorder compared to the suicides with bipolar disorder network. The suicides and non-suicides with schizophrenia networks exhibited clear differences, including the presence of circadian cycle genes in the suicides with schizophrenia network and their absence in the non-suicides with schizophrenia network. The results of this study provide insight on the regulatory mechanisms underpinning transcriptional changes in the suicidal brain.

Metropolitan age-specific mortality trends at borough and neighborhood level: The case of Mexico City.

Understanding the spatial and temporal patterns of mortality rates in a highly heterogeneous metropolis, is a matter of public policy interest. In this context, there is no, to the best of our knowledge, previous studies that correlate both spatio-temporal and age-specific mortality rates in Mexico City. Spatio-temporal Kriging modeling was used over five age-specific mortality rates (from the years 2000 to 2016 in Mexico City), to gain both spatial (borough and neighborhood) and temporal (year and trimester) data level description. Mortality age-specific patterns have been modeled using multilevel modeling for longitudinal data. Posterior tests were carried out to compare mortality averages between geo-spatial locations. Mortality correlation extends in all study groups for as long as 12 years and as far as 13.27 km. The highest mortality rate takes place in the Cuauhtémoc borough, the commercial, touristic and cultural core downtown of Mexico City. On the contrary, Tlalpan borough is the one with the lowest mortality rates in all the study groups. Post-productive mortality is the first age-specific cause of death, followed by infant, productive, pre-school and scholar groups. The combinations of spatio-temporal Kriging estimation and time-evolution linear mixed-effect models, allowed us to unveil relevant time and location trends that may be useful for public policy planning in Mexico City.

A bioinformatics pipeline for Mycobacterium tuberculosis sequencing that cleans contaminant reads from sputum samples.

Next-Generation Sequencing (NGS) is widely used to investigate genomic variation. In several studies, the genetic variation of Mycobacterium tuberculosis has been analyzed in sputum samples without previous culture, using target enrichment methodologies for NGS. Alignments obtained by different programs generally map the sequences under default parameters, and from these results, it is assumed that only Mycobacterium reads will be obtained. However, variants of interest microorganism in clinical samples can be confused with a vast collection of reads from other bacteria, viruses, and human DNA. Currently, there are no standardized pipelines, and the cleaning success is never verified since there is a lack of rigorous controls to identify and remove reads from other sputum-microorganisms genetically similar to M. tuberculosis. Therefore, we designed a bioinformatic pipeline to process NGS data from sputum samples, including several filters and quality control points to identify and eliminate non-M. tuberculosis reads to obtain a reliable genetic variant report. Our proposal uses the SURPI software as a taxonomic classifier to filter input sequences and perform a mapping that provides the highest percentage of Mycobacterium reads, minimizing the reads from other microorganisms. We then use the filtered sequences to perform variant calling with the GATK software, ensuring the mapping quality, realignment, recalibration, hard-filtering, and post-filter to increase the reliability of the reported variants. Using default mapping parameters, we identified reads of contaminant bacteria, such as Streptococcus, Rhotia, Actinomyces, and Veillonella. Our final mapping strategy allowed a sequence identity of 97.8% between the input reads and the whole M. tuberculosis reference genome H37Rv using a genomic edit distance of three, thus removing 98.8% of the off-target sequences with a Mycobacterium reads loss of 1.7%. Finally, more than 200 unreliable genetic variants were removed during the variant calling, increasing the report's reliability.