A positive correlation between GC content and growth temperature in prokaryotes.

BACKGROUND: GC pairs are generally more stable than AT pairs; GC-rich genomes were proposed to be more adapted to high temperatures than AT-rich genomes. Previous studies consistently showed positive correlations between growth temperature and the GC contents of structural RNA genes. However, for the whole genome sequences and the silent sites of the codons in protein-coding genes, the relationship between GC content and growth temperature is in a long-lasting debate. RESULTS: With a dataset much larger than previous studies (681 bacteria and 155 archaea with completely assembled genomes), our phylogenetic comparative analyses showed positive correlations between optimal growth temperature (Topt) and GC content both in bacterial and archaeal structural RNA genes and in bacterial whole genome sequences, chromosomal sequences, plasmid sequences, core genes, and accessory genes. However, in the 155 archaea, we did not observe a significant positive correlation of Topt with whole-genome GC content (GCw) or GC content at four-fold degenerate sites. We randomly drew 155 samples from the 681 bacteria for 1000 rounds. In most cases (> 95%), the positive correlations between Topt and genomic GC contents became statistically nonsignificant (P > 0.05). This result suggested that the small sample sizes might account for the lack of positive correlations between growth temperature and genomic GC content in the 155 archaea and the bacterial samples of previous studies. Comparing the GC content among four categories (psychrophiles/psychrotrophiles, mesophiles, thermophiles, and hyperthermophiles) also revealed a positive correlation between GCw and growth temperature in bacteria. By including the GCw of incompletely assembled genomes, we expanded the sample size of archaea to 303. Positive correlations between GCw and Topt appear especially after excluding the halophilic archaea whose GC contents might be strongly shaped by intense UV radiation. CONCLUSIONS: This study explains the previous contradictory observations and ends a long debate. Prokaryotes growing in high temperatures have higher GC contents. Thermal adaptation is one possible explanation for the positive association. Meanwhile, we propose that the elevated efficiency of DNA repair in response to heat mutagenesis might have the by-product of increasing GC content like that happens in intracellular symbionts and marine bacterioplankton.

Desiccation does not drastically increase the accessibility of exogenous DNA to nuclear genomes: evidence from the frequency of endosymbiotic DNA transfer.

BACKGROUND: Although horizontal gene transfer (HGT) is a widely accepted force in the evolution of prokaryotic genomes, its role in the evolution of eukaryotic genomes remains hotly debated. Some bdelloid rotifers that are resistant to extreme desiccation and radiation undergo a very high level of HGT, whereas in another desiccation-resistant invertebrate, the tardigrade, the pattern does not exist. Overall, the DNA double-strand breaks (DSBs) induced by prolonged desiccation have been postulated to open a gateway to the nuclear genome for exogenous DNA integration and thus to facilitate the HGT process, thereby enhancing the rate of endosymbiotic DNA transfer (EDT). RESULTS: We first surveyed the abundance of nuclear mitochondrial DNAs (NUMTs) and nuclear plastid DNAs (NUPTs) in five eukaryotes that are highly resistant to desiccation: the bdelloid rotifers Adineta vaga and Adineta ricciae, the tardigrade Ramazzottius varieornatus, and the resurrection plants Dorcoceras hygrometricum and Selaginella tamariscina. Excessive NUMTs or NUPTs were not detected. Furthermore, we compared 24 groups of desiccation-tolerant organisms with their relatively less desiccation-tolerant relatives but did not find a significant difference in NUMT/NUPT contents. CONCLUSIONS: Desiccation may induce DSBs, but it is unlikely to dramatically increase the frequency of exogenous sequence integration in most eukaryotes. The capture of exogenous DNA sequences is possible only when DSBs are repaired through a subtype of non-homologous end joining, named alternative end joining (alt-EJ). Due to the deleterious effects of the resulting insertion mutations, alt-EJ is less frequently initiated than other mechanisms.

Aerobic prokaryotes do not have higher GC contents than anaerobic prokaryotes, but obligate aerobic prokaryotes have.

BACKGROUND: Among the four bases, guanine is the most susceptible to damage from oxidative stress. Replication of DNA containing damaged guanines results in G to T mutations. Therefore, the mutations resulting from oxidative DNA damage are generally expected to predominantly consist of G to T (and C to A when the damaged guanine is not in the reference strand) and result in decreased GC content. However, the opposite pattern was reported 16 years ago in a study of prokaryotic genomes. Although that result has been widely cited and confirmed by nine later studies with similar methods, the omission of the effect of shared ancestry requires a re-examination of the reliability of the results. RESULTS: When aerobic and obligate aerobic prokaryotes were mixed together and anaerobic and obligate anaerobic prokaryotes were mixed together, phylogenetic controlled analyses did not detect significant difference in GC content between aerobic and anaerobic prokaryotes. This result is consistent with two generally neglected studied that had accounted for the phylogenetic relationship. However, when obligate aerobic prokaryotes were compared with aerobic prokaryotes, anaerobic prokaryotes, and obligate anaerobic prokaryotes separately using phylogenetic regression analysis, a significant positive association was observed between aerobiosis and GC content, no matter it was calculated from whole genome sequences or the 4-fold degenerate sites of protein-coding genes. Obligate aerobes have significantly higher GC content than aerobes, anaerobes, and obligate anaerobes. CONCLUSIONS: The positive association between aerobiosis and GC content could be attributed to a mutational force resulting from incorporation of damaged deoxyguanosine during DNA replication rather than oxidation of the guanine nucleotides within DNA sequences. Our results indicate a grade in the aerobiosis-associated mutational force, strong in obligate aerobes, moderate in aerobes, weak in anaerobes and obligate anaerobes.

Evaluation of the mechanisms of intron loss and gain in the social amoebae Dictyostelium.

BACKGROUND: Spliceosomal introns are a common feature of eukaryotic genomes. To approach a comprehensive understanding of intron evolution on Earth, studies should look beyond repeatedly studied groups such as animals, plants, and fungi. The slime mold Dictyostelium belongs to a supergroup of eukaryotes not covered in previous studies. RESULTS: We found 441 precise intron losses in Dictyostelium discoideum and 202 precise intron losses in Dictyostelium purpureum. Consistent with these observations, Dictyostelium discoideum was found to have significantly more copies of reverse transcriptase genes than Dictyostelium purpureum. We also found that the lost introns are significantly further from the 5' end of genes than the conserved introns. Adjacent introns were prone to be lost simultaneously in Dictyostelium discoideum. In both Dictyostelium species, the exonic sequences flanking lost introns were found to have a significantly higher GC content than those flanking conserved introns. Together, these observations support a reverse-transcription model of intron loss in which intron losses were caused by gene conversion between genomic DNA and cDNA reverse transcribed from mature mRNA. We also identified two imprecise intron losses in Dictyostelium discoideum that may have resulted from genomic deletions. Ninety-eight putative intron gains were also observed. Consistent with previous studies of other lineages, the source sequences were found in only a small number of cases, with only two instances of intron gain identified in Dictyostelium discoideum. CONCLUSIONS: Although they diverged very early from animals and fungi, Dictyostelium species have similar mechanisms of intron loss.

Isolation and complete nucleotide sequence of a Batai virus strain in Inner Mongolia, China.

BACKGROUND: Batai virus (BATV) is a member of the Orthobunyavirus genus of the family Bunyaviridae, and a vector-borne pathogen. Genomic variations of BATV exist in different regions of the world, due to genetic reassortment. Whole-genome sequencing of any isolate is necessary for a phylogenetic analysis. In 1998, a BATV strain was isolated from an Anopheles philippines mosquito in Yunnan Province, China. This strain has not been found to infect any other host. We investigated BATV infection in cattle in Inner Mongolia, China and performed deep sequencing of the genome of the BATV isolate. FINDINGS: Ninety-five blood samples were collected from cattle in Inner Mongolia, China in 2012. The BATV infection rate was 2.1%. Previously, BATV strain NM/12 was isolated from two cattle in Inner Mongolia, China, and the whole genomic sequence of the strain has been available. We determined the complete genomic nucleotide sequences of the small (S), medium (M), and large (L) genome segments using bovine blood obtained in 2012, and the nucleotide homologies of these segments with those from GenBank were 88.7%-97%, 84%-95.4%, and 72.6%-95.8%, respectively. The deduced amino acid identities were 87.2-99.7%, 64.2-96.8%, and 81.1-98.6%. Phylogenetic analyses based on full-length genomic sequences indicated that the M and L segments, and a portion of the S segment, of NM/12 are most closely related to the BATV strains isolated in Asia. The S and M segments of NM/12 were independent of phylogenetic lineages. The L segment was the most closely related to Chittoor/IG-20217 (isolated in India), and distantly related to isolated strains in Italy. Nucleotide substitution rates in the nucleotide sequences that code for the nucleocapsid, envelope glycoprotein, and polymerase protein of NM/12 strain were 2.56%, 4.69%, and 4.21%, respectively, relative to the original strain of MM2222. CONCLUSION: A novel BATV NM/12 strain from bovine serum collected in Inner Mongolia was isolated from cattle in China for the first time. Our findings elucidate the evolutionary status of the BATV NM/12 strain among different orthobunyavirus strains and may provide some clues to prevent the emergence of BATV infection in cattle and humans.

Frequency of intron loss correlates with processed pseudogene abundance: a novel strategy to test the reverse transcriptase model of intron loss.

BACKGROUND: Although intron loss in evolution has been described, the mechanism involved is still unclear. Three models have been proposed, the reverse transcriptase (RT) model, genomic deletion model and double-strand-break repair model. The RT model, also termed mRNA-mediated intron loss, suggests that cDNA molecules reverse transcribed from spliced mRNA recombine with genomic DNA causing intron loss. Many studies have attempted to test this model based on its predictions, such as simultaneous loss of adjacent introns, 3'-side bias of intron loss, and germline expression of intron-lost genes. Evidence either supporting or opposing the model has been reported. The mechanism of intron loss proposed in the RT model shares the process of reverse transcription with the formation of processed pseudogenes. If the RT model is correct, genes that have produced more processed pseudogenes are more likely to undergo intron loss. RESULTS: In the present study, we observed that the frequency of intron loss is correlated with processed pseudogene abundance by analyzing a new dataset of intron loss obtained in mice and rats. Furthermore, we found that mRNA molecules of intron-lost genes are mostly translated on free cytoplasmic ribosomes, a feature shared by mRNA molecules of the parental genes of processed pseudogenes and long interspersed elements. This feature is likely convenient for intron-lost gene mRNA molecules to be reverse transcribed. Analyses of adjacent intron loss, 3'-side bias of intron loss, and germline expression of intron-lost genes also support the RT model. CONCLUSIONS: Compared with previous evidence, the correlation between the abundance of processed pseudogenes and intron loss frequency more directly supports the RT model of intron loss. Exploring such a correlation is a new strategy to test the RT model in organisms with abundant processed pseudogenes.

Characteristics and factors of repeated influenza vaccination among elderly individuals in Shanghai, China from 2020 to 2022.

Elderly individuals face a high risk of hospitalization and death related to influenza, thus prioritizing them for influenza vaccination. Due to variations in the influenza virus and waning protective antibodies, annual influenza vaccination is recommended. However, research on repeated influenza vaccination among elderly individuals in China is limited. From 2020 to 2022, the average influenza vaccination coverage among registered elderly individuals in Shanghai was 4.1%, showing a declining trend over time. In 2020, the rate of repeated influenza vaccination among elderly individuals was 28.35%, which rose to almost two-thirds both in 2021 and 2022. No increased risk of adverse events following immunization was observed after repeated influenza vaccination during this period. Our study also found that elderly individuals with Shanghai household registration, managed by community clinics, and older age tended to receive more doses of repeated influenza vaccination throughout the period from 2020 to 2022. Increasing influenza vaccine coverage among elderly individuals in Shanghai is both urgent and challenging. Health authorities should intensify educational and promotional campaigns to encourage uptake of annual repeated influenza vaccination among elderly individuals.

Can ENCODE tell us how much junk DNA we carry in our genome?

One of the large, unsolved problems in human genetics is the proportion of functional sequences in genomes. Recently, the encyclopedia of DNA elements consortium revealed that the majority of the genome is biochemically active, which were described as biochemical functions. This has been used as evidence to pronounce the death of the junk DNA concept. In evolutionary biology, junk DNAs are sequences whose gain or loss does not seriously affect fitness of the host organism. In the human genome, a large amount of biochemical activity should be to repress the sequences so as to avoid their harmful expression. The biochemical activity is very different from functionality in the light of evolution. The single nucleotide polymorphism sites associated with disease and other phenotypes may be functional, but their abundance in the active genome regions is not reliable evidence of functionality. Because of sequence-independent functions, the proportion of functional regions would be underestimated when sequence constraints are used alone. Knockout may be the most effective means of distinguishing functional sequences from junk DNA.

Investigating the Relationship between CRISPR-Cas Content and Growth Rate in Bacteria.

CRISPR-Cas systems provide adaptive immunity for prokaryotic cells by recognizing and eliminating the recurrent genetic invaders whose sequences had been captured in a prior infection and stored in the CRISPR arrays as spacers. However, the biological/environmental factors determining the efficiency of this immune system have yet to be fully characterized. Recent studies in cultured bacteria showed that slowing the growth rate of bacterial cells could promote their acquisition of novel spacers. This study examined the relationship between the CRISPR-Cas content and the minimal doubling time across the bacteria and the archaea domains. Every completely sequenced genome could be used to predict a minimal doubling time. With a large data set of 4,142 bacterial samples, we found that the predicted minimal doubling times are positively correlated with spacer number and other parameters of the CRISPR-Cas systems, like array number, Cas gene cluster number, and Cas gene number. Different data sets gave different results. Weak results were obtained in analyzing bacterial empirical minimal doubling times and the archaea domain. Still, the conclusion of more spacers in slowly grown prokaryotes was supported. In addition, we found that the minimal doubling times are negatively correlated with the occurrence of prophages, and the spacer numbers per array are negatively associated with the number of prophages. These observations support the existence of an evolutionary trade-off between bacterial growth and adaptive defense against virulent phages. IMPORTANCE Accumulating evidence indicates that slowing the growth of cultured bacteria could stimulate their CRISPR spacer acquisition. We observed a positive correlation between CRISPR-Cas content and cell cycle duration across the bacteria domain. This observation extends the physiological conclusion to an evolutionary one. In addition, the correlation provides evidence supporting the existence of a trade-off between bacterial growth/reproduction and antiviral resistance.

Precipitous Increase of Bacterial CRISPR-Cas Abundance at Around 45°C.

Although performing adaptive immunity, CRISPR-Cas systems are present in only 40% of bacterial genomes. We observed an abrupt increase of bacterial CRISPR-Cas abundance at around 45°C. Phylogenetic comparative analyses confirmed that the abundance correlates with growth temperature only at the temperature range around 45°C. From the literature, we noticed that the diversities of cellular predators (like protozoa, nematodes, and myxobacteria) have a steep decline at this temperature range. The grazing risk faced by bacteria reduces substantially at around 45°C and almost disappears above 60°C. We propose that viral lysis would become the dominating factor of bacterial mortality, and antivirus immunity has a higher priority at higher temperatures. In temperature ranges where the abundance of cellular predators does not change with temperature, the growth temperatures of bacteria would not significantly affect their CRISPR-Cas contents. The hypothesis predicts that bacteria should also be rich in CRISPR-Cas systems if they live in other extreme conditions inaccessible to grazing predators.

Intron losses and gains in the nematodes.

BACKGROUND: The evolution of spliceosomal introns has been widely studied among various eukaryotic groups. Researchers nearly reached the consensuses on the pattern and the mechanisms of intron losses and gains across eukaryotes. However, according to previous studies that analyzed a few genes or genomes, Nematoda seems to be an eccentric group. RESULTS: Taking advantage of the recent accumulation of sequenced genomes, we extensively analyzed the intron losses and gains using 104 nematode genomes across all the five Clades of the phylum. Nematodes have a wide range of intron density, from less than one to more than nine per kbp coding sequence. The rates of intron losses and gains exhibit significant heterogeneity both across different nematode lineages and across different evolutionary stages of the same lineage. The frequency of intron losses far exceeds that of intron gains. Five pieces of evidence supporting the model of cDNA-mediated intron loss have been observed in ten Caenorhabditis species, the dominance of the precise intron losses, frequent loss of adjacent introns, high-level expression of the intron-lost genes, preferential losses of short introns, and the preferential losses of introns close to 3'-ends of genes. Like studies in most eukaryotic groups, we cannot find the source sequences for the limited number of intron gains detected in the Caenorhabditis genomes. CONCLUSIONS: These results indicate that nematodes are a typical eukaryotic group rather than an outlier in intron evolution.

Evaluation of the effect of a three-color ladder management model for knee osteoarthritis in the community.

OBJECTIVE: To explore the application effect of a three-color ladder management system for knee osteoarthritis in the community. METHODS: Eighty-six patients with knee osteoarthritis in our community were obtained for study and randomly grouped. The control group received routine management, while the research group received three-color ladder management for 12 months. The knee joint function (WOMAC score), pain degree (VAS score), joint flexibility, health-related behavior score, self-care ability scale (exercise of self-care agency scale (ESCA) score), quality of life (knee osteoarthritis quality of life scale (AIMS2-SF) score) and knee replacement rate were compared between the two groups before and after management, and the changes of patients' visits and treatment costs before and after management were observed. RESULTS: After 12 months, the scores of WOMAC and VAS in the research group were significantly lower than those of the control group (P<0.05), while the scores of joint flexibility and extension, cognition, behavior and condition of Omaha System health-related behaviors, ESCA and AIMS2-SF were significantly higher than those of the control group (P<0.05). After 12 months, the monthly visits and expenses of green cards, yellow cards and red cards in the research group were significantly lower than those before entering the group (P<0.05). After 12 months, the knee replacement rate was 20.93% (9/43) in the research group, while it was 27.91% (12/43) in the control group, with no significant difference between the two groups (P>0.05). CONCLUSION: The three-color ladder management system for knee osteoarthritis patients in the community can reduce the number of doctor visits and overall expenses, improve knee joint function, reduce pain, improve self-management ability and quality of life, and it has high community popularization.

Clinical effect of warm needle penetration in treating knee osteoarthritis.

OBJECTIVE: Our aim was to explore the clinical effect of warm needle penetration in treating knee osteoarthritis. METHODS: We randomly divided 118 patients with knee osteoarthritis into the observation group (n=59) and the control group (n=59). The observation group was treated with warm needle penetration combined with western medicine therapy, and the control group was only given western medicine therapy. The clinical effective rate, Visual analogue scale pain score, knee joint score, Western Ontario and McMaster Universities osteoarthritis index and subjective satisfaction were observed and compared between the two groups. RESULTS: The clinical effective rate, knee joint score and subjective satisfaction rate were higher (all P<0.05), while the Visual analogue scale pain score and Western Ontario and McMaster Universities osteoarthritis indices score were lower in the observation group than in the control group (both P<0.05). CONCLUSION: Warm needle penetration can ameliorate the efficiency of treatment, reduce pain, and improve the range of knee joint motion and subjective satisfaction of patients with knee osteoarthritis.

Evaluation of models of the mechanisms underlying intron loss and gain in Aspergillus fungi.

Although intron loss and gain have been widely observed, their mechanisms are still to be determined. In four Aspergillus genomes, we found 204 cases of intron loss and 84 cases of intron gain. Using this data, we tested common hypotheses of intron loss or gain. Statistical analysis showed that adjacent introns tend to be lost simultaneously and small introns were preferentially lost, supporting the model of mRNA-mediated intron loss. The lost introns reside in internal regions of genes, which is inconsistent with the traditional version of the model (partial length cDNAs are reverse transcribed from 3' ends of mRNAs), but consistent with an alternate version (partial length cDNAs are produced by self-primed reverse transcription). The latter version was not supported by examination of the abundance of T-rich segments in mRNAs. Preferential loss of internal introns might be explained by highly efficient recombination at internal regions of genes. Among the 84 cases of intron gain, we found a significantly higher frequency of short direct repeats near exon-intron boundary than in conserved introns, supporting the double-strand break repair model. We also found possible source sequences for two cases of intron gain, one by gene conversion and one by insertion of a mitochondrial sequence during double-strand break repair. Source sequences for most gained introns could not be identified and the possible reasons were discussed. In the four Aspergillus genomes studied, we did not find evidence of frequent parallel intron gains.

Low contents of carbon and nitrogen in highly abundant proteins: evidence of selection for the economy of atomic composition.

Proteins that assimilate particular elements were found to avoid using amino acids containing the element, which indicates that the metabolic constraints of amino acids may influence the evolution of proteins. We suspected that low contents of carbon, nitrogen, and sulfur may also be selected for economy in highly abundant proteins that consume large amounts of the resources of cells. By analyzing recently available proteomic data in Escherichia coli, Saccharomyces cerevisiae, and Schizosaccharomyces pombe, we found that at least the carbon and nitrogen contents in amino acid side chains are negatively correlated with protein abundance. An amino acid with a high number of carbon atoms in its side chain generally requires relatively more energy for its synthesis. Thus, it may be selected against in highly abundant proteins either because of economy in building blocks or because of economy in energy. Previous studies showed that highly abundant proteins preferentially use cheap (in terms of energy) amino acids. We found that the carbon content is still negatively correlated with protein abundance after controlling for the energetic cost of the amino acids. However, the negative correlation between protein abundance and energetic cost disappeared after controlling for carbon content. Building blocks seem to be more restricted than energy. It seems that the amino acid sequences of highly abundant proteins have to compromise between optimization for their biological functions and reducing the consumption of limiting resources. By contrast, the amino acid sequences of weakly expressed proteins are more likely to be optimized for their biological functions.

Characterization of Human Dosage-Sensitive Transcription Factor Genes.

Copy number changes in protein-coding genes are detrimental if the consequent changes in protein concentrations disrupt essential cellular functions. The dosage sensitivity of transcription factor (TF) genes is particularly interesting because their products are essential in regulating the expression of genetic information. From four recently curated data sets of dosage-sensitive genes (genes with conserved copy numbers across mammals, ohnologs, and two data sets of haploinsufficient genes), we compiled a data set of the most reliable dosage-sensitive (MRDS) genes and a data set of the most reliable dosage-insensitive (MRDIS) genes. The MRDS genes were those present in all four data sets, while the MRDIS genes were those absent from any one of the four data sets and with the probability of being loss of function-intolerant (pLI) values < 0.5 in both of the haploinsufficient gene data sets. Enrichment analysis of TF genes among the MRDS and MRDIS gene data sets showed that TF genes are more likely to be dosage-sensitive than other genes in the human genome. The nuclear receptor family was the most enriched TF family among the dosage-sensitive genes. TF families with very few members were also deemed more likely to be dosage-sensitive than TF families with more members. In addition, we found a certain number of dosage-insensitive TFs. The most typical were the Krüppel-associated box domain-containing zinc-finger proteins (KZFPs). Gene ontology (GO) enrichment analysis showed that the MRDS TFs were enriched for many more terms than the MRDIS TFs; however, the proteins interacting with these two groups of TFs did not show such sharp differences. Furthermore, we found that the MRDIS KZFPs were not significantly enriched for any GO terms, whereas their interacting proteins were significantly enriched for thousands of GO terms. Further characterizations revealed significant differences between MRDS TFs and MRDIS TFs in the lengths and nucleotide compositions of DNA-binding sites as well as in expression level, protein size, and selective force.

Evidence against the energetic cost hypothesis for the short introns in highly expressed genes.

BACKGROUND: In animals, the moss Physcomitrella patens and the pollen of Arabidopsis thaliana, highly expressed genes have shorter introns than weakly expressed genes. A popular explanation for this is selection for transcription efficiency, which includes two sub-hypotheses: to minimize the energetic cost or to minimize the time cost. RESULTS: In an individual human, different organs may differ up to hundreds of times in cell number (for example, a liver versus a hypothalamus). Considered at the individual level, a gene specifically expressed in a large organ is actually transcribed tens or hundreds of times more than a gene with a similar expression level (a measure of mRNA abundance per cell) specifically expressed in a small organ. According to the energetic cost hypothesis, the former should have shorter introns than the latter. However, in humans and mice we have not found significant differences in intron length between large-tissue/organ-specific genes and small-tissue/organ-specific genes with similar expression levels. Qualitative estimation shows that the deleterious effect (that is, the energetic burden) of long introns in highly expressed genes is too negligible to be efficiently selected against in mammals. CONCLUSION: The short introns in highly expressed genes should not be attributed to energy constraint. We evaluated evidence for the time cost hypothesis and other alternatives.

Association between the availability of environmental resources and the atomic composition of organismal proteomes: evidence from Prochlorococcus strains living at different depths.

The cyanobacteria Prochlorococcus is a cyanbacterial genus, with some strains adapted to sea surface environments, which are poor in nutrients and have high-light intensity, and some strains adapted to deep sea conditions, which have relatively higher concentrations of nitrogen and phosphorus and lower light intensity. Here, we report pairwise comparisons between strains isolated from different depths of the same sea, which reveal a close association between atomic composition of the proteome and the availability nitrogen and phosphorus in the environment. The atomic composition of proteomes differs significantly among Prochlorococcus strains with different supplies of nitrogen in vivo; these different supplies result from different capacities for nitrogen assimilation. We repeated our whole-proteome analysis with the core proteomes of Prochlorococcus and obtained similar results. Our findings indicate that the elemental composition of proteomes is shaped by the availability of resources in the environment.

Genotypic variations between wild-type and small colony variant of Staphylococcus aureus in prosthetic valve infectious endocarditis: a comparative genomic and transcriptomic analysis.

Staphylococcus aureus small colony variants (SCVs) can cause persistent infections. However, the genomes and transcriptomes of S. aureus SCVs remain poorly understood. A pair of isogenic wild-type and SCV methicillin-resistant S. aureus (MRSA) strains (IE1 and IE2, respectively) were isolated from a patient with prosthetic valve infectious endocarditis. The SCV strain IE2 grew more slowly than the wild-type strain, and serum killing and mouse lethality assays revealed that the virulence of SCV strain IE2 was decreased. Whole-genome sequencing of the SCV and wild-type strains revealed 15 mutations in nine genes associated with metabolism, virulence and DNA repair, including serine/threonine-protein kinase PrkC (prkC), glycerol-3-phosphate acyltransferase (plsY), 2-deoxyribose-5-phosphate aldolase (deoC), extracellular adherence protein (eap), iron compound ABC uptake transporter substrate-binding protein (sstD), RecU Holliday junction resolvase (recU), excinuclease ABC subunit B (uvrB), type I restriction-modification system, M subunit (hsdM) and smooth muscle caldesmon. Sequencing of RNA transcripts revealed that expression levels of 321 genes were upregulated and 582 genes were downregulated in SCV strain IE2 compared with IE1. Most of the differentially expressed genes were involved in metabolism. Expression levels of several genes involved in the pathways to which plsY, deoC, eap and sstD belonged were changed, associated with the metabolism and virulence of S. aureus. In conclusion, the reduced growth rate and decreased virulence of MRSA SCV strains may be related to mutations in and downregulation of genes associated with metabolism and virulence, especially plsY, deoC, eap and sstD.