Sequence-Based Analysis of Thermal Adaptation and Protein Energy Landscapes in an Invasive Blue Mussel (Mytilus galloprovincialis).

Adaptive responses to thermal stress in poikilotherms plays an important role in determining competitive ability and species distributions. Amino acid substitutions that affect protein stability and modify the thermal optima of orthologous proteins may be particularly important in this context. Here, we examine a set of 2,770 protein-coding genes to determine if proteins in a highly invasive heat tolerant blue mussel (Mytilus galloprovincialis) contain signals of adaptive increases in protein stability relative to orthologs in a more cold tolerant M. trossulus. Such thermal adaptations might help to explain, mechanistically, the success with which the invasive marine mussel M. galloprovincialis has displaced native species in contact zones in the eastern (California) and western (Japan) Pacific. We tested for stabilizing amino acid substitutions in warm tolerant M. galloprovincialis relative to cold tolerant M. trossulus with a generalized linear model that compares in silico estimates of recent changes in protein stability among closely related congeners. Fixed substitutions in M. galloprovincialis were 3,180.0 calories per mol per substitution more stabilizing at genes with both elevated dN/dS ratios and transcriptional responses to heat stress, and 705.8 calories per mol per substitution more stabilizing across all 2,770 loci investigated. Amino acid substitutions concentrated in a small number of genes were more stabilizing in M. galloprovincialis compared with cold tolerant M. trossulus. We also tested for, but did not find, enrichment of a priori GO terms in genes with elevated dN/dS ratios in M. galloprovincialis. This might indicate that selection for thermodynamic stability is generic across all lineages, and suggests that the high change in estimated protein stability that we observed in M. galloprovincialis is driven by selection for extra stabilizing substitutions, rather than by higher incidence of selection in a greater number of genes in this lineage. Nonetheless, our finding of more stabilizing amino acid changes in the warm adapted lineage is important because it suggests that adaption for thermal stability has contributed to M. galloprovincialis' superior tolerance to heat stress, and that pairing tests for positive selection and tests for transcriptional response to heat stress can identify candidates of protein stability adaptation

Petőfi Sándor levele Arany Jánosnak

BackgroundComparative genomics studies investigating the signals of positive selection among groups of closely related species are still rare and limited in taxonomic breadth. Such studies show great promise in advancing our knowledge about the proportion and the identity of genes experiencing diversifying selection. However, methodological challenges have led to high levels of false positives in past studies. Here, we use the well-annotated genome of the purple sea urchin, Strongylocentrotus purpuratus, as a reference to investigate the signals of positive selection at 6520 single-copy orthologs from nine sea urchin species belonging to the family Strongylocentrotidae paying careful attention to minimizing false positives.ResultsWe identified 1008 (15.5%) candidate positive selection genes (PSGs). Tests for positive selection along the nine terminal branches of the phylogeny identified 824 genes that showed lineage-specific adaptive diversification (1.67% of branch-sites tests performed). Positively selected codons were not enriched at exon borders or near regions containing missing data, suggesting a limited contribution of false positives caused by alignment or annotation errors. Alignments were validated at 10 loci with re-sequencing using Sanger methods. No differences were observed in the rates of synonymous substitution (d S), GC content, and codon bias between the candidate PSGs and those not showing positive selection. However, the candidate PSGs had 68% higher rates of nonsynonymous substitution (d N) and 33% lower levels of heterozygosity, consistent with selective sweeps and opposite to that expected by a relaxation of selective constraint. Although positive selection was identified at reproductive proteins and innate immunity genes, the strongest signals of adaptive diversification were observed at extracellular matrix proteins, cell adhesion molecules, membrane receptors, and ion channels. Many candidate PSGs have been widely implicated as targets of pathogen binding, inactivation, mimicry, or exploitation in other groups (notably mammals).ConclusionsOur study confirmed the widespread action of positive selection across sea urchin genomes and allowed us to reject the possibility that annotation and alignment errors (including paralogs) were responsible for creating false signals of adaptive molecular divergence. The candidate PSGs identified in our study represent promising targets for future research into the selective agents responsible for their adaptive diversification and their contribution to speciation

Whole-genome sequencing reveals host factors underlying critical COVID-19

Critical COVID-19 is caused by immune-mediated inflammatory lung injury. Host genetic variation influences the development of illness requiring critical care1 or hospitalization2–4 after infection with SARS-CoV-2. The GenOMICC (Genetics of Mortality in Critical Care) study enables the comparison of genomes from individuals who are critically ill with those of population controls to find underlying disease mechanisms. Here we use whole-genome sequencing in 7,491 critically ill individuals compared with 48,400 controls to discover and replicate 23 independent variants that significantly predispose to critical COVID-19. We identify 16 new independent associations, including variants within genes that are involved in interferon signalling (IL10RB and PLSCR1), leucocyte differentiation (BCL11A) and blood-type antigen secretor status (FUT2). Using transcriptome-wide association and colocalization to infer the effect of gene expression on disease severity, we find evidence that implicates multiple genes—including reduced expression of a membrane flippase (ATP11A), and increased expression of a mucin (MUC1)—in critical disease. Mendelian randomization provides evidence in support of causal roles for myeloid cell adhesion molecules (SELE, ICAM5 and CD209) and the coagulation factor F8, all of which are potentially druggable targets. Our results are broadly consistent with a multi-component model of COVID-19 pathophysiology, in which at least two distinct mechanisms can predispose to life-threatening disease: failure to control viral replication; or an enhanced tendency towards pulmonary inflammation and intravascular coagulation. We show that comparison between cases of critical illness and population controls is highly efficient for the detection of therapeutically relevant mechanisms of disease

Genetic mechanisms of critical illness in COVID-19.

Host-mediated lung inflammation is present1, and drives mortality2, in the critical illness caused by coronavirus disease 2019 (COVID-19). Host genetic variants associated with critical illness may identify mechanistic targets for therapeutic development3. Here we report the results of the GenOMICC (Genetics Of Mortality In Critical Care) genome-wide association study in 2,244 critically ill patients with COVID-19 from 208 UK intensive care units. We have identified and replicated the following new genome-wide significant associations: on chromosome 12q24.13 (rs10735079, P = 1.65 × 10-8) in a gene cluster that encodes antiviral restriction enzyme activators (OAS1, OAS2 and OAS3); on chromosome 19p13.2 (rs74956615, P = 2.3 × 10-8) near the gene that encodes tyrosine kinase 2 (TYK2); on chromosome 19p13.3 (rs2109069, P = 3.98 ×  10-12) within the gene that encodes dipeptidyl peptidase 9 (DPP9); and on chromosome 21q22.1 (rs2236757, P = 4.99 × 10-8) in the interferon receptor gene IFNAR2. We identified potential targets for repurposing of licensed medications: using Mendelian randomization, we found evidence that low expression of IFNAR2, or high expression of TYK2, are associated with life-threatening disease; and transcriptome-wide association in lung tissue revealed that high expression of the monocyte-macrophage chemotactic receptor CCR2 is associated with severe COVID-19. Our results identify robust genetic signals relating to key host antiviral defence mechanisms and mediators of inflammatory organ damage in COVID-19. Both mechanisms may be amenable to targeted treatment with existing drugs. However, large-scale randomized clinical trials will be essential before any change to clinical practice

Whole-genome sequencing reveals host factors underlying critical COVID-19

Critical COVID-19 is caused by immune-mediated inflammatory lung injury. Host genetic variation influences the development of illness requiring critical care1 or hospitalization2,3,4 after infection with SARS-CoV-2. The GenOMICC (Genetics of Mortality in Critical Care) study enables the comparison of genomes from individuals who are critically ill with those of population controls to find underlying disease mechanisms. Here we use whole-genome sequencing in 7,491 critically ill individuals compared with 48,400 controls to discover and replicate 23 independent variants that significantly predispose to critical COVID-19. We identify 16 new independent associations, including variants within genes that are involved in interferon signalling (IL10RB and PLSCR1), leucocyte differentiation (BCL11A) and blood-type antigen secretor status (FUT2). Using transcriptome-wide association and colocalization to infer the effect of gene expression on disease severity, we find evidence that implicates multiple genes—including reduced expression of a membrane flippase (ATP11A), and increased expression of a mucin (MUC1)—in critical disease. Mendelian randomization provides evidence in support of causal roles for myeloid cell adhesion molecules (SELE, ICAM5 and CD209) and the coagulation factor F8, all of which are potentially druggable targets. Our results are broadly consistent with a multi-component model of COVID-19 pathophysiology, in which at least two distinct mechanisms can predispose to life-threatening disease: failure to control viral replication; or an enhanced tendency towards pulmonary inflammation and intravascular coagulation. We show that comparison between cases of critical illness and population controls is highly efficient for the detection of therapeutically relevant mechanisms of disease

Biochemical studies on the expression of overdominance at the phosphoglucomutase-2 locus in the Pacific oyster, Crassostrea gigas (Thunberg)

Numerous studies have documented significant associations between multiple-locus heterozygosity and fitness-related traits in natural populations, but the explanations for these patterns remain unknown. The objective of the present study was to examine the merits of the overdominance hypothesis as the mechanism responsible for a positive correlation between adult body weight and heterozygosity involving the phosphoglucomutase-2 (Pgm-2) locus in the Pacific oyster, Crassostrea gigas. The kinetic and structural properties of seven Pgm-2 genotypes were examined over physiological ranges of temperature and pH. Significant differences were detected between Pgm-2 genotypes in a variety of enzymic parameters, but these were largely confined to genotypes possessing the Pgm-2-92 allele, and heterozygotes displayed strict intermediacy for all functional and structural properties examined. The expression of marginal overdominance at the Pgm-2 locus was considered unlikely because of the limited scope of the observed variation between allozymes, and its incompatibility with allelic frequencies in natural populations. The three most common heterozygotes at the Pgm-2 locus displayed the extremely unusual property of overdominant enzyme activities. The magnitude of this overdominance was similar in the mantle and adductor muscle tissues, and was consistently observed in population samples from two intertidal positions in three different seasons. A physiological impact of the Pgm-2 polymorphism was demonstrated on the metabolism of glycogen, the biochemical pathway in which PGM functions. Pgm-2 genotypes exhibited different concentrations of glycogen in their mantle, but not their adductor muscle tissues, -which were directly associated with variation in their PGM activity levels. It was suggested that Pgm-2 genotype-dependent enzyme activity variation may affect rates of glycogen synthesis by a partitioning effect at the glucose-6-phosphate branch point. Non-random associations were detected between the PGM activities of Pgm-2 genotypic groups and the activities of adjacent glycogen synthesis pathway enzymes, but none that could clearly account for the differing glycogen concentrations observed between genotypes. The expression of overdominance for PGM activity, and its impact on mantle glycogen levels, provided direct evidence favoring the overdominance explanation as the cause of the larger body weights of heterozygotes at the Pgm-2 locus in Crassostrea gigas.Science, Faculty ofZoology, Department ofGraduat

INTROGRESS_v.1.22_input_file_(parent_2_data)

This is file 3 (parental 2 data) of 4 input files needed for running the R package INTROGRESS v1.22 (Gompert & Berkle 2010). These files together allow maximum-likelihood estimate of the ancestry (multilocus hybrid score) of each individual, and genomic cline analysis, which compares admixture at a single locus to average admixture across the rest of the genome. We used the parametric approach, first estimating a multilocus hybrid index and then fitting clines in genotype frequencies at individual SNPs as a function of the neutral expectation derived from the genome-wide hybrid index

Data from: Introgression between invasive and native blue mussels (genus Mytilus) in the central California hybrid zone.

The ecological and genetic factors determining the extent of introgression between species in secondary contact zones remain poorly understood. Here, we investigate the relative importance of isolating barriers and the demographic expansion of invasive Mytilus galloprovincialis on the magnitude and the direction of introgression with the native Mytilus trossulus in a hybrid zone in central California. We use double-digest restriction-site-associated DNA sequencing (ddRADseq) to genotype 1337 randomly selected single nucleotide polymorphisms and accurately distinguish early and advanced generation hybrids for the first time in the central California Mytilus spp. hybrid zone. Weak levels of introgression were observed in both directions but were slightly more prevalent from the native M. trossulus into the invasive M. galloprovincialis. Few early and advanced backcrossed individuals were observed across the hybrid zone confirming the presence of strong barriers to interbreeding. Heterogeneous patterns of admixture across the zone of contact were consistent with the colonization history of M. galloprovincialis with more extensive introgression in northern localities furthest away from the putative site of introduction in southern California. These observations reinforce the importance of dynamic spatial and demographic expansions in determining patterns of introgression between close congeners, even in those with high dispersal potential and well-developed reproductive barriers. Our results suggest that the threat posed by invasive M. galloprovincialis is more ecological than genetic as it has displaced, and continues to displace the native M. trossulus from much of central and southern California

Genome-wide patterns of codon bias are shaped by natural selection in the purple sea urchin, Strongylocentrotus purpuratus.

Codon usage bias has been documented in a wide diversity of species, but the relative contributions of mutational bias and various forms of natural selection remain unclear. Here, we describe for the first time genome-wide patterns of codon bias at 4623 genes in the purple sea urchin, Strongylocentrotus purpuratus. Preferred codons were identified at 18 amino acids that exclusively used G or C at third positions, which contrasted with the strong AT bias of the genome (overall GC content is 36.9%). The GC content of third positions and coding regions exhibited significant correlations with the magnitude of codon bias. In contrast, the GC content of introns and flanking regions was indistinguishable from the genome-wide background, which suggested a limited contribution of mutational bias to synonymous codon usage. Five distinct clusters of genes were identified that had significantly different synonymous codon usage patterns. A significant correlation was observed between codon bias and mRNA expression supporting translational selection, but this relationship was driven by only one highly biased cluster that represented only 8.6% of all genes. In all five clusters preferred codons were evolutionarily conserved to a similar degree despite differences in their synonymous codon usage distributions and magnitude of codon bias. The third positions of preferred codons in two codon usage groups also paired significantly more often in stems than in loops of mRNA secondary structure predictions, which suggested that codon bias might also affect mRNA stability. Our results suggest that mutational bias has played a minor role in determining codon bias in S. purpuratus and that preferred codon usage may be heterogeneous across different genes and subject to different forms of natural selection

Molecular Population Genetics of the Male and Female Mitochondrial DNA Molecules of the California Sea Mussel, Mytilus californianus

The presence of two gender-associated mitochondrial genomes in marine mussels provides a unique opportunity to investigate the dynamics of mtDNA evolution without complications inherent in interspecific comparisons. Here, we assess the relative importance of selection, mutation, and differential constraint in shaping the patterns of polymorphism within and divergence between the male (M) and female (F) mitochondrial genomes of the California sea mussel, Mytilus californianus. Partial sequences were obtained from homologous regions of four genes (nad2, cox1, atp6, and nad5) totaling 2307 bp in length. The M and F mtDNA molecules of M. californianus exhibited extensive levels of nucleotide polymorphism and were more highly diverged than observed in other mytilids (overall Tamura–Nei distances >40%). Consistent with previous studies, the M molecule had significantly higher levels of silent and replacement polymorphism relative to F. Both genomes possessed large numbers of singleton and low-frequency mutations that gave rise to significantly negative Tajima's D values. Mutation-rate scalars estimated for silent and replacement mutations were elevated in the M genome but were not sufficient to account for its higher level of polymorphism. McDonald–Kreitman tests were highly significant at all loci due to excess numbers of fixed replacement mutations between molecules. Strong purifying selection was evident in both genomes in keeping the majority of replacement mutations at low population frequencies but appeared to be slightly relaxed in M. Our results suggest that a reduction in selective constraint acting on the M genome remains the best explanation for its greater levels of polymorphism and faster rate of evolution