Molecular evolution of HoxA13 and the multiple origins of limbless morphologies in amphibians and reptiles

Developmental processes and their results, morphological characters, are inherited through transmission of genes regulating development. While there is ample evidence that cis-regulatory elements tend to be modular, with sequence segments dedicated to different roles, the situation for proteins is less clear, being particularly complex for transcription factors with multiple functions. Some motifs mediating protein-protein interactions may be exclusive to particular developmental roles, but it is also possible that motifs are mostly shared among different processes. Here we focus on HoxA13, a protein essential for limb development. We asked whether the HoxA13 amino acid sequence evolved similarly in three limbless clades: Gymnophiona, Amphisbaenia and Serpentes. We explored variation in ω (dN/dS) using a maximum-likelihood framework and HoxA13sequences from 47 species. Comparisons of evolutionary models provided low ω global values and no evidence that HoxA13 experienced relaxed selection in limbless clades. Branch-site models failed to detect evidence for positive selection acting on any site along branches of Amphisbaena and Gymnophiona, while three sites were identified in Serpentes. Examination of alignments did not reveal consistent sequence differences between limbed and limbless species. We conclude that HoxA13 has no modules exclusive to limb development, which may be explained by its involvement in multiple developmental processes

Recommendations for accurate genotyping of SARS-CoV-2 using amplicon-based sequencing of clinical samples.

Genotyping of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been instrumental in monitoring viral evolution and transmission during the pandemic. The quality of the sequence data obtained from these genotyping efforts depends on several factors, including the quantity/integrity of the input material, the technology, and laboratory-specific implementation. The current lack of guidelines for SARS-CoV-2 genotyping leads to inclusion of error-containing genome sequences in genomic epidemiology studies. We aimed to establish clear and broadly applicable recommendations for reliable virus genotyping. We established and used a sequencing data analysis workflow that reliably identifies and removes technical artefacts; such artefacts can result in miscalls when using alternative pipelines to process clinical samples and synthetic viral genomes with an amplicon-based genotyping approach. We evaluated the impact of experimental factors, including viral load and sequencing depth, on correct sequence determination. We found that at least 1000 viral genomes are necessary to confidently detect variants in the SARS-CoV-2 genome at frequencies of ≥10%. The broad applicability of our recommendations was validated in over 200 clinical samples from six independent laboratories. The genotypes we determined for clinical isolates with sufficient quality cluster by sampling location and period. Our analysis also supports the rise in frequencies of 20A.EU1 and 20A.EU2, two recently reported European strains whose dissemination was facilitated by travel during the summer of 2020. We present much-needed recommendations for the reliable determination of SARS-CoV-2 genome sequences and demonstrate their broad applicability in a large cohort of clinical samples

Tight associations between transcription promoter type and epigenetic variation in histone positioning and modification

Abstract Background Transcription promoters are fundamental genomic cis-elements controlling gene expression. They can be classified into two types by the degree of imprecision of their transcription start sites: peak promoters, which initiate transcription from a narrow genomic region; and broad promoters, which initiate transcription from a wide-ranging region. Eukaryotic transcription initiation is suggested to be associated with the genomic positions and modifications of nucleosomes. For instance, it has been recently shown that histone with H3K9 acetylation (H3K9ac) is more likely to be distributed around broad promoters rather than peak promoters; it can thus be inferred that there is an association between histone H3K9 and promoter architecture. Results Here, we performed a systematic analysis of transcription promoters and gene expression, as well as of epigenetic histone behaviors, including genomic position, stability within the chromatin, and several modifications. We found that, in humans, broad promoters, but not peak promoters, generally had significant associations with nucleosome positioning and modification. Specifically, around broad promoters histones were highly distributed and aligned in an orderly fashion. This feature was more evident with histones that were methylated or acetylated; moreover, the nucleosome positions around the broad promoters were more stable than those around the peak ones. More strikingly, the overall expression levels of genes associated with broad promoters (but not peak promoters) with modified histones were significantly higher than the levels of genes associated with broad promoters with unmodified histones. Conclusion These results shed light on how epigenetic regulatory networks of histone modifications are associated with promoter architecture

The Case of the Fickle Fingers: How the PRDM9 Zinc Finger Protein Specifies Meiotic Recombination Hotspots in Humans

Recent discoveries have revealed the central role of PRDM9 in mammalian recombination. The precise function of this protein, however, remains poorly understood, as do the causes for its rapid evolution and its role in reproductive isolation

The Anti-Apoptotic Bcl-xL Protein, a New Piece in the Puzzle of Cytochrome C Interactome

A structural model of the adduct between human cytochrome c and the human anti-apoptotic protein Bcl-xL, which defines the protein-protein interaction surface, was obtained from solution NMR chemical shift perturbation data. The atomic level information reveals key intermolecular contacts identifying new potentially druggable areas on cytochrome c and Bcl-xL. Involvement of residues on cytochrome c other than those in its complexes with electron transfer partners is apparent. Key differences in the contact area also exist between the Bcl-xL adduct with the Bak peptide and that with cytochrome c. The present model provides insights to the mechanism by which cytochrome c translocated to cytosol can be intercepted, so that the apoptosome is not assembled

MyoMiner: explore gene co-expression in normal and pathological muscle

International audienceBackground: High-throughput transcriptomics measures mRNA levels for thousands of genes in a biological sample. Most gene expression studies aim to identify genes that are differentially expressed between different biological conditions, such as between healthy and diseased states. However, these data can also be used to identify genes that are co-expressed within a biological condition. Gene co-expression is used in a guilt-by-association approach to prioritize candidate genes that could be involved in disease, and to gain insights into the functions of genes, protein relations, and signaling pathways. Most existing gene co-expression databases are generic, amalgamating data for a given organism regardless of tissue-type.Methods: To study muscle-specific gene co-expression in both normal and pathological states, publicly available gene expression data were acquired for 2376 mouse and 2228 human striated muscle samples, and separated into 142 categories based on species (human or mouse), tissue origin, age, gender, anatomic part, and experimental condition. Co-expression values were calculated for each category to create the MyoMiner database.Results: Within each category, users can select a gene of interest, and the MyoMiner web interface will return all correlated genes. For each co-expressed gene pair, adjusted p-value and confidence intervals are provided as measures of expression correlation strength. A standardized expression-level scatterplot is available for every gene pair r-value. MyoMiner has two extra functions: (a) a network interface for creating a 2-shell correlation network, based either on the most highly correlated genes or from a list of genes provided by the user with the option to include linked genes from the database and (b) a comparison tool from which the users can test whether any two correlation coefficients from different conditions are significantly different.Conclusions: These co-expression analyses will help investigators to delineate the tissue-, cell-, and pathology-specific elements of muscle protein interactions, cell signaling and gene regulation. Changes in co-expression between pathologic and healthy tissue may suggest new disease mechanisms and help define novel therapeutic targets. Thus, MyoMiner is a powerful muscle-specific database for the discovery of genes that are associated with related functions based on their co-expression. MyoMiner is freely available at https://www.sys-myo.com/myominer

Human Macrophages and Dendritic Cells Can Equally Present MART-1 Antigen to CD8+ T Cells after Phagocytosis of Gamma-Irradiated Melanoma Cells

Dendritic cells (DC) can achieve cross-presentation of naturally-occurring tumor-associated antigens after phagocytosis and processing of dying tumor cells. They have been used in different clinical settings to vaccinate cancer patients. We have previously used gamma-irradiated MART-1 expressing melanoma cells as a source of antigens to vaccinate melanoma patients by injecting irradiated cells with BCG and GM-CSF or to load immature DC and use them as a vaccine. Other clinical trials have used IFN-gamma activated macrophage killer cells (MAK) to treat cancer patients. However, the clinical use of MAK has been based on their direct tumoricidal activity rather than on their ability to act as antigen-presenting cells to stimulate an adaptive antitumor response. Thus, in the present work, we compared the fate of MART-1 after phagocytosis of gamma-irradiated cells by clinical grade DC or MAK as well as the ability of these cells to cross present MART-1 to CD8+ T cells. Using a high affinity antibody against MART-1, 2A9, which specifically stains melanoma tumors, melanoma cell lines and normal melanocytes, the expression level of MART-1 in melanoma cell lines could be related to their ability to stimulate IFN-gamma production by a MART-1 specific HLA-A*0201-restricted CD8+ T cell clone. Confocal microscopy with Alexa Fluor®647-labelled 2A9 also showed that MART-1 could be detected in tumor cells attached and/or fused to phagocytes and even inside these cells as early as 1 h and up to 24 h or 48 h after initiation of co-cultures between gamma-irradiated melanoma cells and MAK or DC, respectively. Interestingly, MART-1 was cross-presented to MART-1 specific T cells by both MAK and DC co-cultured with melanoma gamma-irradiated cells for different time-points. Thus, naturally occurring MART-1 melanoma antigen can be taken-up from dying melanoma cells into DC or MAK and both cell types can induce specific CD8+ T cell cross-presentation thereafter