Human genetic variation in the Baltic Sea region: Features of population history and natural selection

In this thesis, the genetic variation of human populations from the Baltic Sea region was studied in order to elucidate population history as well as evolutionary adaptation in this region. The study provided novel understanding of how the complex population level processes of migration, genetic drift, and natural selection have shaped genetic variation in North European populations. Results from genome-wide, mitochondrial DNA and Y-chromosomal analyses suggested that the genetic background of the populations of the Baltic Sea region lies predominantly in Continental Europe, which is consistent with earlier studies and archaeological evidence. The late settlement of Fennoscandia after the Ice Age and the subsequent small population size have led to pronounced genetic drift, especially in Finland and Karelia but also in Sweden, evident especially in genome-wide and Y-chromosomal analyses. Consequently, these populations show striking genetic differentiation, as opposed to much more homogeneous pattern of variation in Central European populations. Additionally, the eastern side of the Baltic Sea was observed to have experienced eastern influence in the genome-wide data as well as in mitochondrial DNA and Y-chromosomal variation – consistent with linguistic connections. However, Slavic influence in the Baltic Sea populations appears minor on genetic level. While the genetic diversity of the Finnish population overall was low, genome-wide and Y-chromosomal results showed pronounced regional differences. The genetic distance between Western and Eastern Finland was larger than for many geographically distant population pairs, and provinces also showed genetic differences. This is probably mainly due to the late settlement of Eastern Finland and local isolation, although differences in ancestral migration waves may contribute to this, too. In contrast, mitochondrial DNA and Y-chromosomal analyses of the contemporary Swedish population revealed a much less pronounced population structure and a fusion of the traces of ancient admixture, genetic drift, and recent immigration. Genome-wide datasets also provide a resource for studying the adaptive evolution of human populations. This study revealed tens of loci with strong signs of recent positive selection in Northern Europe. These results provide interesting targets for future research on evolutionary adaptation, and may be important for understanding the background of disease-causing variants in human populations.

Molecular Alterations in Asbestos-Related Lung Cancer

Background: Asbestos is a well known cancer-causing mineral fibre, which has a synergistic effect on lung cancer risk in combination with tobacco smoking. Several in vitro and in vivo experiments have demonstrated that asbestos can evoke chromosomal damage and cause alterations as well as gene expression changes. Lung tumours, in general, have very complex karyotypes with several recurrently gained and lost chromosomal regions and this has made it difficult to identify specific molecular changes related primarily to asbestos exposure. The main aim of these studies has been to characterize asbestos-related lung cancer at a molecular level. Methods: Samples from asbestos-exposed and non-exposed lung cancer patients were studied using array comparative genomic hybridization (aCGH) and fluorescent in situ hybridization (FISH) to detect copy number alterations (CNA) as well as microsatellite analysis to detect allelic imbalance (AI). In addition, asbestos-exposed cell lines were studied using gene expression microarrays. Results: Eighteen chromosomal regions showing differential copy number in the lung tumours of asbestos-exposed patients compared to those of non-exposed patients were identified. The most significant differences were detected at 2p21-p16.3, 5q35.3, 9q33.3-q34.11, 9q34.13-q34.3, 11p15.5, 14q11.2 and 19p13.1-p13.3 (p<0.005). The alterations at 2p and 9q were validated and characterized in detail using AI and FISH analysis in a larger study population. Furthermore, in vitro studies were performed to examine the early gene expression changes induced by asbestos in three different lung cell lines. The results revealed specific asbestos-associated gene expression profiles and biological processes as well as chromosomal regions enriched with genes believed to contribute to the common asbestos-related responses in the cell lines. Interestingly, the most significant region enriched with asbestos-response genes was identified at 2p22, close to the previously identified region showing asbestos-related CNA in lung tumours. Additionally, in this thesis, the dysregulated biological processes (Gene Ontology terms) detected in the cell line experiment were compared to dysregulated processes identified in patient samples in a later study (Ruosaari et al., 2008a). Commonly affected processes such as those related to protein ubiquitination, ion transport and surprisingly sensory perception of smell were identified. Conclusions: The identification of specific CNA and dysregulated biological processes shed some light on the underlying genes acting as mediators in asbestos-related lung carcinogenesis. It is postulated that the combination of several asbestos-specific molecular alterations could be used to develop a diagnostic method for the identification of asbestos-related lung cancer.

Phylogeography and hybrid swarms: history of brackish water bivalve diversity in North European marginal seas

This study addressed the large-scale molecular zoogeography in two brackish water bivalve molluscs, Macoma balthica and Cerastoderma glaucum, and genetic signatures of the postglacial colonization of Northern Europe by them. The traditional view poses that M. balthica in the Baltic, White and Barents seas (i.e. marginal seas) represent direct postglacial descendants of the adjacent Northeast Atlantic populations, but this has recently been challenged by observations of close genetic affinities between these marginal populations and those of the Northeast Pacific. The primary aim of the thesis was to verify, quantify and characterize the Pacific genetic contribution across North European populations of M. balthica and to resolve the phylogeographic histories of the two bivalve taxa in range-wide studies using information from mitochondrial DNA (mtDNA) and nuclear allozyme polymorphisms. The presence of recent Pacific genetic influence in M. balthica of the Baltic, White and Barents seas, along with an Atlantic element, was confirmed by mtDNA sequence data. On a broader temporal and geographical scale, altogether four independent trans-Arctic invasions of Macoma from the Pacific since the Miocene seem to have been involved in generating the current North Atlantic lineage diversity. The latest trans-Arctic invasion that affected the current Baltic, White and Barents Sea populations probably took place in the early post-glacial. The nuclear genetic compositions of these marginal sea populations are intermediate between those of pure Pacific and Atlantic subspecies. In the marginal sea populations of mixed ancestry (Barents, White and Northern Baltic seas), the Pacific and Atlantic components are now randomly associated in the genomes of individual clams, which indicates both pervasive historical interbreeding between the previously long-isolated lineages (subspecies), and current isolation of these populations from the adjacent pure Atlantic populations. These mixed populations can be characterized as self-supporting hybrid swarms, and they arguably represent the most extensive marine animal hybrid swarms so far documented. Each of the three swarms still has a distinct genetic composition, and the relative Pacific contributions vary from 30 to 90 % in local populations. This diversity highlights the potential of introgressive hybridization to rapidly give rise to new evolutionarily and ecologically significant units in the marine realm. In the south of the Danish straits and in the Southern Baltic Sea, a broad genetic transition zone links the pure North Sea subspecies M. balthica rubra to the inner Baltic hybrid swarm, which has about 60 % of Pacific contribution in its genome. This transition zone has no regular smooth clinal structure, but its populations show strong genotypic disequilibria typical of a hybrid zone maintained by the interplay of selection and gene flow by dispersing pelagic larvae. The structure of the genetic transition is partly in line with features of Baltic water circulation and salinity stratification, with greater penetration of Atlantic genes on the Baltic south coast and in deeper water populations. In all, the scenarios of historical isolation and secondary contact that arise from the phylogeographic studies of both Macoma and Cerastoderma shed light to the more general but enigmatic patterns seen in marine phylogeography, where deep genetic breaks are often seen in species with high dispersal potential.

Avian diversification patterns across the K-Pg boundary: Influence of calibrations, datasets, and model misspecification

Birds represent the most diverse extant tetrapod clade, with ca. 10,000 extant species, and the timing of the crown avian radiation remains hotly debated. The fossil record supports a primarily Cenozoic radiation of crown birds, whereas molecular divergence dating analyses generally imply that this radiation was well underway during the Cretaceous. Furthermore, substantial differences have been noted between published divergence estimates. These have been variously attributed to clock model, calibration regime, and gene type. One underappreciated phenomenon is that disparity between fossil ages and molecular dates tends to be proportionally greater for shallower nodes in the avian Tree of Life. Here, we explore potential drivers of disparity in avian divergence dates through a set of analyses applying various calibration strategies and coding methods to a mitochondrial genome dataset and an 18-gene nuclear dataset, both sampled across 72 taxa. Our analyses support the occurrence of two deep divergences (i.e., the Palaeognathae/Neognathae split and the Galloanserae/Neoaves split) well within the Cretaceous, followed by a rapid radiation of Neoaves near the K-Pg boundary. However, 95% highest posterior density intervals for most basal divergences in Neoaves cross the boundary, and we emphasize that, barring unreasonably strict prior distributions, distinguishing between a rapid Early Paleocene radiation and a Late Cretaceous radiation may be beyond the resolving power of currently favored divergence dating methods. In contrast to recent observations for placental mammals, constraining all divergences within Neoaves to occur in the Cenozoic does not result in unreasonably high inferred substitution rates. Comparisons of nuclear DNA (nDNA) versus mitochondrial DNA (mtDNA) datasets and NT- versus RY-coded mitochondrial data reveal patterns of disparity that are consistent with substitution model misspecifications that result in tree compression/tree extension artifacts, which may explain some discordance between previous divergence estimates based on different sequence types. Comparisons of fully calibrated and nominally calibrated trees support a correlation between body mass and apparent dating error. Overall, our results are consistent with (but do not require) a Paleogene radiation for most major clades of crown birds.

HER2 positive gastric and gastroesophageal adenocarcinoma; An Irish tertiary center experience

Background: Trastuzumab has been approved for patients with human epidermal growth factor receptor 2 (HER2) over expression and gene amplification metastatic gastric cancer. Here we present the prevalence of HER2 positive gastric cancer in an Irish population, the use of Trastuzumab in first line and beyond progression. Methods: The study was conducted in St James's Hospital, Dublin. A retrospective analysis of the date of patients with HER2 positive gastric cancer over a period of 3 years was carried out. Her2 positive was defined as immunohistochemistry (IHC) score of +3, of IHC score of +2 and increased gene copy number by fluorescence in situ hybridization (FISH). Overall survival was calculated from the day of initiation of treatment with Trastuzumab until death. Results: During the study period 140 patients with gastric and gastro-esophageal junction adenocarcinoma were treated. Out of those, 30 (21.4%) had HER2 positive disease. Among HER2 positive disease patients 18 (12.8%) were treated with first line Trastuzumab containing regimen with a median overall survival of 13 months. Nine (50%) developed progressive disease while on Trastuzumab and of those, 4 (22.2%) patients continued on Trastuzumab beyond progression, two (11.1%) of whom achieved stable disease and a prolonged survival. Conclusion: HER2 positivity rate in an Irish population with advanced gastric and gastro-esophageal junction adenocarcinoma is 21.4%. Treatment with Trastuzumab in the first line in combination with chemotherapy is a reasonable approach. Continuation of Trastuzumab beyond progression is a feasible strategy that requires further exploration.

Data integration from two microarray platforms identifies bi-allelic genetic inactivation of RIC8A in a breast cancer cell line

Background: Using array comparative genomic hybridization (aCGH), a large number of deleted genomic regions have been identified in human cancers. However, subsequent efforts to identify target genes selected for inactivation in these regions have often been challenging. Methods: We integrated here genome-wide copy number data with gene expression data and non-sense mediated mRNA decay rates in breast cancer cell lines to prioritize gene candidates that are likely to be tumour suppressor genes inactivated by bi-allelic genetic events. The candidates were sequenced to identify potential mutations. Results: This integrated genomic approach led to the identification of RIC8A at 11p15 as a putative candidate target gene for the genomic deletion in the ZR-75-1 breast cancer cell line. We identified a truncating mutation in this cell line, leading to loss of expression and rapid decay of the transcript. We screened 127 breast cancers for RIC8A mutations, but did not find any pathogenic mutations. No promoter hypermethylation in these tumours was detected either. However, analysis of gene expression data from breast tumours identified a small group of aggressive tumours that displayed low levels of RIC8A transcripts. qRT-PCR analysis of 38 breast tumours showed a strong association between low RIC8A expression and the presence of TP53 mutations (P = 0.006). Conclusion: We demonstrate a data integration strategy leading to the identification of RIC8A as a gene undergoing a classical double-hit genetic inactivation in a breast cancer cell line, as well as in vivo evidence of loss of RIC8A expression in a subgroup of aggressive TP53 mutant breast cancers.

Genomic alterations in Myeloproliferative Neoplasms and Myelodysplastic Syndromes

Myeloproliferative neoplasms (MPN) and myelodysplastic syndromes (MDS) are a heterogeneous group of clonal hematopoietic disorders whose etiology and molecular pathogenesis are poorly understood. During the past decade, enormous developments in microarray technology and bioinformatics methods have made it possible to mine novel molecular alterations in a large number of malignancies, including MPN and MDS, which has facilitated the detection of new prognostic, predictive and therapeutic biomarkers for disease stratification. By applying novel microarray techniques, we profiled copy number alterations and microRNA (miRNA) expression changes in bone marrow aspirate and blood samples. In addition, we set up and validated an miRNA expression test for bone marrow core biopsies in order to utilize the large archive material available in many laboratories. We also tested JAK2 mutation status and compare it with the in vitro growth pattern of hematologic progenitors cells. In the study focusing on 100 MPN cases, we detected a Janus kinase 2 (JAK2) mutation in 71 cases. We observed spontaneous erythroid colony growth in all mutation-positive cases in addition to nine mutation negative cases. Interestingly, seven JAK2V167F negative ET cases showed spontaneous megakaryocyte colony formation, one case of which also harbored a myeloproliferative leukemia virus oncogene (MPL) mutation. We studied copy number alterations in 35 MPN and 37 MDS cases by using oligonucleotide-based array comparative hybridization (array CGH). Only one essential thrombocythemia (ET) case presented copy number alterations in chromosomes 1q and 13q. In contrast, MDS cases were characterized by numerous novel cryptic chromosomal aberrations with the most common copy number losses at 5q21.3q33.1 and 7q22.1q33, while the most common copy number gain was trisomy 8. As for the study of the bone marrow core biopsy samples, we showed that even though these samples were embedded in paraffin and underwent decalcification, they were reliable sources of miRNA and suitable for array expression analysis. Further, when studying the miRNA expression profiles of the 19 MDS cases, we found that, compared to controls, two miRNAs (one human Epstein-Barr virus (miR-BART13) miRNA and one human (has-miR-671-5p) miRNA) were downregulated, whereas two other miRNAs (hsa-miR-720 and hsa-miR-21) were upregulated. However, we could find no correlation between copy number alterations and microRNA expression when integrating these two data. This thesis brings to light new information about genomic changes implicated in the development of MPN and MDS, and also underlines the power of applying genome-wide array screening techniques in neoplasias. Rapid advances in molecular techniques and the integration of different genomic data will enable the discovery of the biological contexts of many complex disorders, including myeloid neoplasias.

Combination antiretroviral therapy -associated lipodystrophy : insights into pathogenesis and treatment

Introduction: Combination antiretroviral therapy (cART) has decreased morbidity and mortality of individuals infected with human immunodeficiency virus type 1 (HIV-1). Its use, however, is associated with adverse effects which increase the patients risk of conditions such as diabetes and coronary heart disease. Perhaps the most stigmatizing side effect is lipodystrophy, i.e., the loss of subcutaneous adipose tissue (SAT) in the face, limbs and trunk while fat accumulates intra-abdominally and dorsocervically. The pathogenesis of cART-associated lipodystrophy is obscure. Nucleoside reverse transcriptase inhibitors (NRTI) have been implicated to cause lipoatrophy via mitochondrial toxicity. There is no known effective treatment for cART-associated lipodystrophy during unchanged antiretroviral regimen in humans, but in vitro data have shown uridine to abrogate NRTI-induced toxicity in adipocytes. Aims: To investigate whether i) cART or lipodystrophy associated with its use affect arterial stiffness; ii) lipoatrophic SAT is inflamed compared to non-lipoatrophic SAT; iii) abdominal SAT from patients with compared to those without cART-associated lipoatrophy differs with respect to mitochondrial DNA (mtDNA) content, adipose tissue inflammation and gene expression, and if NRTIs stavudine and zidovudine are associated with different degree of changes; iv) lipoatrophic abdominal SAT differs from preserved dorsocervical SAT with respect to mtDNA content, adipose tissue inflammation and gene expression in patients with cART-associated lipodystrophy and v) whether uridine can revert lipoatrophy and the associated metabolic disturbances in patients on stavudine or zidovudine based cART. Subjects and methods: 64 cART-treated patients with (n=45) and without lipodystrophy/-atrophy (n=19) were compared cross-sectionally. A marker of arterial stiffness, heart rate corrected augmentation index (AgIHR), was measured by pulse wave analysis. Body composition was measured by magnetic resonance imaging and dual-energy X-ray absorptiometry, and liver fat content by proton magnetic resonance spectroscopy. Gene expression and mtDNA content in SAT were assessed by real-time polymerase chain reaction and microarray. Adipose tissue composition and inflammation were assessed by histology and immunohistochemistry. Dorsocervical and abdominal SAT were studied. The efficacy and safety of uridine for the treatment of cART-associated lipoatrophy were evaluated in a randomized, double-blind, placebo-controlled 3-month trial in 20 lipoatrophic cART-treated patients. Results: Duration of antiretroviral treatment and cumulative exposure to NRTIs and protease inhibitors, but not the presence of cART-associated lipodystrophy, predicted AgIHR independent of age and blood pressure. Gene expression of inflammatory markers was increased in SAT of lipodystrophic as compared to non-lipodystrophic patients. Expression of genes involved in adipogenesis, triglyceride synthesis and glucose disposal was lower and of those involved in mitochondrial biogenesis, apoptosis and oxidative stress higher in SAT of patients with than without cART-associated lipoatrophy. Most changes were more pronounced in stavudine-treated than in zidovudine-treated individuals. Lipoatrophic SAT had lower mtDNA than SAT of non-lipoatrophic patients. Expression of inflammatory genes was lower in dorsocervical than in abdominal SAT. Neither depot had characteristics of brown adipose tissue. Despite being spared from lipoatrophy, dorsocervical SAT of lipodystrophic patients had lower mtDNA than the phenotypically similar corresponding depot of non-lipodystrophic patients. The greatest difference in gene expression between dorsocervical and abdominal SAT, irrespective of lipodystrophy status, was in expression of homeobox genes that regulate transcription and regionalization of organs during embryonal development. Uridine increased limb fat and its proportion of total fat, but had no effect on liver fat content and markers of insulin resistance. Conclusions: Long-term cART is associated with increased arterial stiffness and, thus, with higher cardiovascular risk. Lipoatrophic abdominal SAT is characterized by inflammation, apoptosis and mtDNA depletion. As mtDNA is depleted even in non-lipoatrophic dorsocervical SAT, lipoatrophy is unlikely to be caused directly by mtDNA depletion. Preserved dorsocervical SAT of patients with cART-associated lipodystrophy is less inflamed than their lipoatrophic abdominal SAT, and does not resemble brown adipose tissue. The greatest difference in gene expression between dorsocervical and abdominal SAT is in expression of transcriptional regulators, homeobox genes, which might explain the differential susceptibility of these adipose tissue depots to cART-induced toxicity. Uridine is able to increase peripheral SAT in lipoatrophic patients during unchanged cART.

Genetic transformation of Brassica nigra by agrobacterium based vector and direct plasmid uptake

Genetic transformation systems have been established for Brassica nigra (cv. IC 257) by using an Agrobacterium binary vector as well as by direct DNA uptake of a plasmid vector. Both the type of vectors carried nptII gene and gus gene. For Agrobacterium mediated transformation, hypocotyl tissue explants were used, and up to 33% of the explants produced calli on selection medium. All of these expressed B-glucuronidase gene on histochemical staining. Protoplasts isolated from hypocotyl tissues of seedlings could be transformed with a plasmid vector by FEG mediated uptake of vector DNA. A number of fertile kanamycin resistant plants were obtained using both the methods, and their transformed nature was confirmed by Southern blot analysis and histochemical staining for GUS. Backcrossed and selfed progenies of these transformed plants showed the presence of npt and gus genes.

Proposed role of w chromosome inactivation and the absence of dosage compensation in avian sex determination

Three features of avian sex chromosomes - female heterogamety (ZZ male, ZW female), the apparently inactive state of the W chromosome, and dose-dependent expression of Z-linked genes - are examined in regard to their possible relation to sex determination. It is proposed that the W chromosome is facultatively heterochromatic and that the Z and W chromosomes carry one or more homologous sex-determination genes. The absence of dosage compensation in ZZ embryos, and W inactivation in ZW embryos, would then bring about a 2n(ZZ)-n(ZW) inequality in the effective copy number of such genes. The absence of dosage compensation of Z-linked genes in ZZ embryos is viewed as a means by which two copies of Z-W homologous sex determination genes are kept active to meet the requirements of testis determination. W inactivation may promote ovarian development by reducing the effective copy number of these genes from 2n to n. If there is a W-specific gene for femaleness, spread of heterochromatization to this gene in cells forming the right gonadal primordium may explain the latter's normally undifferentiated state; reversal of heterochromatization may similarly explain the development of the right gonad into a testis following left ovariectomy.

Population genetic structure and conservation of Asian elephants (Elephas maximus) across India

This study examines the population genetic structure of Asian elephants (Elephas maximus) across India, which harbours over half the world's population of this endangered species. Mitochondrial DNA control region sequences and allele frequencies at six nuclear DNA microsatellite markers obtained from the dung of free-ranging elephants reveal low mtDNA and typical microsatellite diversity. Both known divergent clades of mtDNA haplotypes in the Asian elephant are present in India, with southern and central India exhibiting exclusively the β clade of Fernando et al. (2000), northern India exhibiting exclusively the α clade and northeastern India exhibiting both, but predominantly the α clade. A nested clade analysis revealed isolation by distance as the principal mechanism responsible for the observed haplotype distributions within the α and β clades. Analyses of molecular variance and pairwise population FST tests based on both mitochondrial and microsatellite DNA suggest that northern-northeastern India, central India, Nilgiris (in southern India) and Anamalai-Periyar (in southern India) are four demographically autonomous population units and should be managed separately. In addition, evidence for female philopatry, male-mediated gene flow and two possible historical biogeographical barriers is described.

Multiple molecular alterations in phosphodegrons 1-3 within AAV2 capsid demonstrates higher hepatic gene transfer efficiency

The success of AAV2 mediated hepatic gene transfer in human trials for diseases such as hemophilia has been hampered by a combination of low transduction efficiency and a robust immune response directed against these vectors. We have previously shown that AAV2 is targeted for destruction in the cytoplasm by the host-cellular kinase/ubiquitination/proteasomal degradation machinery and modification of the serine(S)/threonine(T) kinase and lysine(K) targets on AAV capsid is beneficial. Thus targeted single mutations of S/T>A(S489A, S498A, T251A) and K>R (K532R) improved the efficiency of gene transfer in vivo as compared to wild type (WT)-AAV2 vectors (∼6-14 fold). In the present study, we evaluated if combined alteration of the phosphodegrons (PD), which are the phosphorylation sites recognized as degradation signals by ubiquitin ligases, improves further the gene transfer efficiency. Thus, we generated four multiple mutant vectors (PD: 1+3, S489A+K532R, PD: 1+3, S489A+K532R together with T251 residue which did not lie in any of the phosphodegrons but had shown increased transduction efficiency compared to the WT-AAV2 vector (∼6 fold) and was also conserved in 9 out of 10 AAV serotypes (AAV 1 to 10), PD: 1+3, S489A+K532R+S498A and a fourth combination of PD: 3, K532R+T251. We then evaluated them in vitro and in vivo and compared their gene transfer efficiency with either the WT-AAV2 or the best single mutant S489A-AAV2 vector. The novel multiple mutations on the AAV2 capsid did not affect the overall vector packaging efficiency. All the multiple AAV2 mutants showed superior transduction efficiency in HeLa cells in vitro when compared to either the WT (62-72% Vs 21%) or the single mutant S489A (62-72% Vs 50%) AAV2 vectors as demonstrated by FACS analysis (Fig. 1A). On hepatic gene transfer with 5x10^10 vgs per animal in C57BL/6 mice, all the multiple mutants showed increased transgene expression compared to either the WT-AAV2 (∼15-23 fold) or the S489A single mutant vector (∼2-3 fold) (Fig.1B and C). These novel multiple mutant AAV2 vectors also showed higher vector copy number in murine hepatocytes 4 weeks post transduction, as compared to the WT-AAV2 (∼5-6 Vs 1.4 vector copies/diploid genome) and further higher when compared to the single mutant S489A(∼5-6 fold Vs 3.8 fold) (Fig.1D). Further ongoing studies will demonstrate the therapeutic benefit of one or more of the multiple mutants vectors in preclinical models of hemophilia.

How Many Initiator tRNA Genes Does Escherichia coli Need?

Multiple copies of a gene require enhanced investment on the part of the cell and, as such, call for an explanation. The observation that Escherichia coli has four copies of initiator tRNA (tRNA(i)) genes, encoding a special tRNA (tRNA(fMet)) required to start protein synthesis, is puzzling particularly because the cell appears to be unaffected by the removal of one copy. However, the fitness of an organism has both absolute and relative connotations. Thus, we carried out growth competition experiments between E. coli strains that differ in the number of tRNA(i) genes they contain. This has enabled us to uncover an unexpected link between the number of tRNA(i) genes and protein synthesis, nutritional status, and fitness. Wild-type strains with the canonical four tRNA(i) genes are favored in nutrient-rich environments, and those carrying fewer are favored in nutrient-poor environments. Auxotrophs behave as if they have a nutritionally poor internal environment. A heuristic model that links tRNA(i) gene copy number, genetic stress, and growth rate accounts for the findings. Our observations provide strong evidence that natural selection can work through seemingly minor quantitative variations in gene copy number and thereby impact organismal fitness.