Natural selection in a postglacial range expansion: the case of the colour cline in the European barn owl.

Gradients of variation-or clines-have always intrigued biologists. Classically, they have been interpreted as the outcomes of antagonistic interactions between selection and gene flow. Alternatively, clines may also establish neutrally with isolation by distance (IBD) or secondary contact between previously isolated populations. The relative importance of natural selection and these two neutral processes in the establishment of clinal variation can be tested by comparing genetic differentiation at neutral genetic markers and at the studied trait. A third neutral process, surfing of a newly arisen mutation during the colonization of a new habitat, is more difficult to test. Here, we designed a spatially explicit approximate Bayesian computation (ABC) simulation framework to evaluate whether the strong cline in the genetically based reddish coloration observed in the European barn owl (Tyto alba) arose as a by-product of a range expansion or whether selection has to be invoked to explain this colour cline, for which we have previously ruled out the actions of IBD or secondary contact. Using ABC simulations and genetic data on 390 individuals from 20 locations genotyped at 22 microsatellites loci, we first determined how barn owls colonized Europe after the last glaciation. Using these results in new simulations on the evolution of the colour phenotype, and assuming various genetic architectures for the colour trait, we demonstrate that the observed colour cline cannot be due to the surfing of a neutral mutation. Taking advantage of spatially explicit ABC, which proved to be a powerful method to disentangle the respective roles of selection and drift in range expansions, we conclude that the formation of the colour cline observed in the barn owl must be due to natural selection.

Population differentiation and its consequences in the white campion "Silene latifolia" (Caryophyllaceae): a functional analysis

Summary In his theory On the Origin of Species by Means of Natural Selection (1859), Darwin describes evolution as a gradual change in population over time and that natural selection is a process that caused evolution. Because quantitative variation in species is partly influenced by several genes and thus heritable, association between levels of genetic variation at neutral markers and at quantitative traits and their partitioning within and among populations are important to study mechanisms that drive evolution in populations. Most studies addressing quantitative variation in plants focused on morphological and life history traits but not in traits affecting reproductive success. The aim of this thesis is to better understand how patterns of variation for neutral molecular markers and phenotypic traits drive the evolution of reproduction and defensive mechanisms in six European populations of Silene latifolia, a dioecious plant species. We found evidence for extremely high within and between population variation at six microsatellite loci and at most quantitative traits studied in plants grown under standardized conditions (morphology, life history and reproductive traits). Interestingly, there was clinal variation between age at first flowering and latitude. This pattern is likely due to natural selection since differentiation of this trait was high, heritable and probably higher than differentiation at neutral markers. Our study focused on sex specific selective pressures: mechanisms of intersexual coadaptation and defence mechanism against the seed predator Hadena bicruris. To address divergence at reproductive traits, we studied male and female population of origin effects and in particular pollen competitive ability on male post-pollination success in the study populations with within and between populations crosses. We crossed the same female plant with pollen from a male within the same population of origin and pollen from two males from two distinct populations, using a fixed tester male as a competitor. Additionally, we conducted control crosses with pollen from each male as a single donor. We analysed paternity success of each competitor with two microsatellite loci, seed set and offspring fitness. Male population of origin showed significant among-population variation for siring success at pollen competition. In vitro pollen germination rate showed heritable variation among populations and was positively correlated to siring success. Local or foreign pollen did not have a consistent advantage. Furthermore, female population of origin affected the outcome of pollen competition in some populations. There was no difference of seed set or offspring fitness in within/ between population crosses. This suggests that reproductive divergence may occur via pollen competition in Silene latifolia. The specialist seed predator Hadena bicruris may also induce divergence between populations. We tested potential constitutive and induced defence mechanisms against the specialist predator Hadena bicruris. Because fruit wall thickness is smaller in the invasive range (Northern America) were the moth is absent, this suggests that a thicker fruit wall is a potentially defensive trait against larval attack, and that relaxed selection in the absence of the seed predator has resulted in an evolutionary loss of this defence in the invasive range. Fruit wall thickness was different among three populations. Experimental exposure to moth eggs increased fruit abortion. Fruits built after attack on exposed plants did not have thicker fruit walls compared to fruits on non-exposed plants. Furthermore, fruits with thicker fruit walls were not less profitable, nor did they require longer handling time when exposed to larvae, suggesting no defensive role of fruit wall thickness. Our results show that there is high molecular and phenotypic variation in Silene latifolia and that traits potentially involved in reproductive success both for intra-specific (between sexes) and inter-specific interactions are heritable. Different selective forces may thus interact and cause differential evolution of geographically separated Silene latifolia populations in Europe, leading to the observed differentiation. Résumé Dans sa théorie de l'évolution, L'origine des espèces, ch. 4 (1859), Darwin décrit l'évolution comme un processus continu au cours du temps à l'intérieur de populations et que la sélection naturelle en est le moteur. La variation quantitative est en partie déterminée par plusieurs gènes, donc transmissible à la descendance. Associer le niveau de variation génétique à des marqueurs neutres au niveau de la variation à des traits quantitatifs, ainsi que la répartition à l'intérieur et entre les populations d'une espèce donnée de cette variation, sont importants dans la compréhension des forces évolutives. La plupart des études scientifiques sur la variation quantitative chez les plantes se sont intéressées à la morphologie et à la phénologie mais pas aux caractères impliqués dans le succès reproducteur. L'objectif de cette thèse est de mieux comprendre comment la répartition de la variation à des marqueurs neutres et des caractères quantitatifs influence l'évolution de la reproduction et des mécanismes de défense dans six populations Européennes de l'espèce dioïque Silene latifolia. Nous avons mis en évidence une grande diversité intra et inter-population à six loci microsatellites ainsi qu'à la plupart des caractères quantitatifs mesurés (morphologie, phénologie et traits reproducteurs) sur des plantes cultivées dans des conditions standardisées. Un résultat intéressant est la présence d'un cline latitudinal pour l'âge à la floraison. Ceci est probablement une conséquence de la sélection naturelle, puisque ce caractère est différencié entre les populations étudiées, héritable et que la différenciation de ce trait est supérieure à la différenciation des marqueurs neutres étudiés. Notre étude a ensuite porté plus précisément sur les pressions de sélection spécifiques aux sexes : la coadaptation entre les sexes et les mécanismes de défense contre l'insecte granivore Hadena bicruris. Afin d'évaluer la divergence sur les traits reproducteurs, nous avons étudié les effets des populations d'origine des mâles et des femelles et en particulier le succès reproducteur des mâles après pollinisation à l'aide de croisements inter et intra-population. Nous avons pollinisé la même femelle avec du pollen provenant d'un mâle de la même population ainsi qu'avec le pollen de deux mâles provenant de deux autres populations en situation de compétition avec un pollen provenant d'une population test. Des croisements contrôle ont été réalisés avec les mêmes mâles en pollinisation pure. Nous avons évalué le succès reproducteur de chaque mâle à l'aide d'analyses de paternité ainsi que la production de graines et la fitness de la descendance. L'origine du mâle avait un effet sur la paternité. Le taux de croissance in vitro du pollen est un caractère héritable et a eu un effet positif sur le succès reproducteur. De plus, l'origine de la femelle avait un effet sur le succès des mâles en compétition dans certaines populations. Nos résultats suggèrent qu'une divergence reproductive chez Silene latifolia pourrait apparaître suite à la compétition pollinique. Nous avons ensuite testé des mécanismes potentiels de défense constitutive et induite contre l'herbivore spécialiste Hadena bicruris, un papillon nocturne qui pourrait aussi jouer un rôle dans la différenciation des populations. L'épaisseur des fruits étant plus faible dans les régions où la plante est invasive (Amérique du Nord) et où l'insecte est absent, ce trait pourrait jouer un rôle défensif. Une pression de sélection plus faible causée par l'absence de l'herbivore aurait abouti à une perte de cette défense dans ces régions. Nous avons montré que l'épaisseur du fruit est variable selon les populations. L'infestation artificielle de fruit par l'insecte induit l'abscission sélective des fruits. Les fruits produits après une infestation n'étaient pas plus épais que les fruits issus de plantes non infestées. De plus, les fruits épais n'étaient pas moins nutritifs et ne causaient pas de perte de temps pour la prédation pour les larves, ce qui suggère que l'épaisseur des fruits ne joue pas un rôle défensif. Nos résultats montrent que plusieurs pressions de sélection interviennent et interagissent dans l'évolution de populations distantes, provoquant la divergence des populations Européennes de l'espèce Silene latifolia.

Geographical and altitudinal population genetic structure of two dung fly species with contrasting mobility and temperature preference.

Local adaptation of populations requires some degree of spatio-temporal isolation. Previous studies of the two dung fly species Scathophaga stercoraria and Sepsis cynipsea have revealed low levels of geographic and altitudinal genetic differentiation in quantitative life history and morphological traits, but instead high degrees of phenotypic plasticity. These patterns suggest that gene flow is extensive despite considerable geographic barriers and large spatio-temporal variation in selection on body size and related traits. In this study we addressed this hypothesis by investigating genetic differentiation of dung fly populations throughout Switzerland based on the same 10 electrophoretic loci in each species. Overall, we found no significant geographic differentiation of populations for either species. This is inconsistent with the higher rates of gene flow expected due to better flying capacity of the larger S. stercoraria. However, heterozygote deficiencies within populations indicated structuring on a finer scale, seen for several loci in S. cynipsea, and for the locus PGM (Phosphoglucomutase) in S. stercoraria. Additionally, S. cynipsea showed a tendency towards a greater gene diversity at higher altitudes, mediated primarily by the locus MDH (malate dehydrogenase), at which a second allele was only present in populations above 1000 m. This may be caused by increased environmental stress at higher altitudes in this warm-adapted species. MDH might thus be a candidate locus subject to thermal selection in this species, but this remains to be corroborated by direct evidence. In S. stercoraria, no altitudinal variation was found.

Bioinformatics-driven identification and examination of candidate genes for non-alcoholic fatty liver disease.

ObjectiveCandidate genes for non-alcoholic fatty liver disease (NAFLD) identified by a bioinformatics approach were examined for variant associations to quantitative traits of NAFLD-related phenotypes.Research Design and MethodsBy integrating public database text mining, trans-organism protein-protein interaction transferal, and information on liver protein expression a protein-protein interaction network was constructed and from this a smaller isolated interactome was identified. Five genes from this interactome were selected for genetic analysis. Twenty-one tag single-nucleotide polymorphisms (SNPs) which captured all common variation in these genes were genotyped in 10,196 Danes, and analyzed for association with NAFLD-related quantitative traits, type 2 diabetes (T2D), central obesity, and WHO-defined metabolic syndrome (MetS).Results273 genes were included in the protein-protein interaction analysis and EHHADH, ECHS1, HADHA, HADHB, and ACADL were selected for further examination. A total of 10 nominal statistical significant associations (P<0.05) to quantitative metabolic traits were identified. Also, the case-control study showed associations between variation in the five genes and T2D, central obesity, and MetS, respectively. Bonferroni adjustments for multiple testing negated all associations.ConclusionsUsing a bioinformatics approach we identified five candidate genes for NAFLD. However, we failed to provide evidence of associations with major effects between SNPs in these five genes and NAFLD-related quantitative traits, T2D, central obesity, and MetS.

Genome-wide association analysis identifies variants associated with nonalcoholic fatty liver disease that have distinct effects on metabolic traits.

Nonalcoholic fatty liver disease (NAFLD) clusters in families, but the only known common genetic variants influencing risk are near PNPLA3. We sought to identify additional genetic variants influencing NAFLD using genome-wide association (GWA) analysis of computed tomography (CT) measured hepatic steatosis, a non-invasive measure of NAFLD, in large population based samples. Using variance components methods, we show that CT hepatic steatosis is heritable (∼26%-27%) in family-based Amish, Family Heart, and Framingham Heart Studies (n = 880 to 3,070). By carrying out a fixed-effects meta-analysis of genome-wide association (GWA) results between CT hepatic steatosis and ∼2.4 million imputed or genotyped SNPs in 7,176 individuals from the Old Order Amish, Age, Gene/Environment Susceptibility-Reykjavik study (AGES), Family Heart, and Framingham Heart Studies, we identify variants associated at genome-wide significant levels (p<5×10(-8)) in or near PNPLA3, NCAN, and PPP1R3B. We genotype these and 42 other top CT hepatic steatosis-associated SNPs in 592 subjects with biopsy-proven NAFLD from the NASH Clinical Research Network (NASH CRN). In comparisons with 1,405 healthy controls from the Myocardial Genetics Consortium (MIGen), we observe significant associations with histologic NAFLD at variants in or near NCAN, GCKR, LYPLAL1, and PNPLA3, but not PPP1R3B. Variants at these five loci exhibit distinct patterns of association with serum lipids, as well as glycemic and anthropometric traits. We identify common genetic variants influencing CT-assessed steatosis and risk of NAFLD. Hepatic steatosis associated variants are not uniformly associated with NASH/fibrosis or result in abnormalities in serum lipids or glycemic and anthropometric traits, suggesting genetic heterogeneity in the pathways influencing these traits.

Models that include supercoiling of topological domains reproduce several known features of interphase chromosomes.

Understanding the structure of interphase chromosomes is essential to elucidate regulatory mechanisms of gene expression. During recent years, high-throughput DNA sequencing expanded the power of chromosome conformation capture (3C) methods that provide information about reciprocal spatial proximity of chromosomal loci. Since 2012, it is known that entire chromatin in interphase chromosomes is organized into regions with strongly increased frequency of internal contacts. These regions, with the average size of ∼1 Mb, were named topological domains. More recent studies demonstrated presence of unconstrained supercoiling in interphase chromosomes. Using Brownian dynamics simulations, we show here that by including supercoiling into models of topological domains one can reproduce and thus provide possible explanations of several experimentally observed characteristics of interphase chromosomes, such as their complex contact maps.

A genomewide linkage study on suicidality in major depressive disorder confirms evidence for linkage to 2p12.

Suicidal behavior is commonly associated with depression. Twin studies indicate that both suicidality and major depressive disorder (MDD) are heritable. However, epidemiological evidence suggests that the inheritance of suicidality is likely to be independent of the underlying psychiatric disorder, implying a distinct genetic contribution to suicidality. We conducted a genomewide linkage search aiming to detect genomic loci that may harbor susceptibility genes contributing to risk for suicidality in recurrent MDD. Affected sibling pair (ASP) variance components analysis was performed using the Depression Network cohort of 971 ASPs. The quantitative trait measuring suicidality as a broad phenotype, encompassing ideation and suicide attempts, was established from Schedules for Clinical Assessment in Neuropsychiatry interview items. We examined 1,060 genotyped microsatellite markers with an average spacing of 3.3 cM. Empirical thresholds for linkage evidence were set by whole-genome simulations (LOD = 2.71 for genomewide significance, 1.71 for suggestive linkage). No genomewide significant findings were found. Marker D3S1234 on 3p14 achieved suggestive linkage and yielded a maximum LOD of 1.853 (P = 0.0017), loci 9p24.3 and 18q22-q23 achieved LOD scores >1.5. We found some support for linkage to 2p12 (LOD = 1.2, P = 0.0087) which was previously implicated in linkage studies of suicidality. Our follow-up meta-analysis of five studies showed strong linkage to this region (P = 2 × 10(-6) ). In conclusion, this study analyzed suicidality as a continuous trait in MDD. We found modest evidence for linkage on 3p14. Our meta-analysis supports previous evidence of linkage to suicidality on 2p12. Some candidate genes in these regions may plausibly be implicated in suicidality.

Dyslipidemia in HIV-infected individuals: from pharmacogenetics to pharmacogenomics.

HIV-infected individuals may have accelerated atherogenesis and an increased risk for premature coronary artery disease. Dyslipidemia represents a key pro-atherogenic mechanism. In HIV-infected patients, dyslipidemia is typically attributed to the adverse effects of antiretroviral therapy. Nine recent genome-wide association studies have afforded a comprehensive, unbiased inventory of common SNPs at 36 genetic loci that are reproducibly associated with dyslipidemia in the general population. Genome-wide association study-validated SNPs have now been demonstrated to contribute to dyslipidemia in the setting of HIV infection and antiretroviral therapy. In a Swiss HIV-infected study population, a similar proportion of serum lipid variability was explained by antiretroviral therapy and by genetic background. In the individual patient, both antiretroviral therapy and the cumulative effect of SNPs contribute to the risk of high low-density lipoprotein cholesterol, low high-density lipoprotein cholesterol and hypertriglyceridemia. Genetic variants presumably contribute to additional major metabolic complications in HIV-infected individuals, including diabetes mellitus and coronary artery disease. In an effort to explain an increasing proportion of the heritability of complex metabolic traits, ongoing large-scale gene resequencing studies are focusing on the effects of rare SNPs and structural genetic variants.

Do riverine barriers, history or introgression shape the genetic structuring of a common shrew (Sorex araneus) population?

The common shrew (Sorex araneus) is subdivided into numerous chromosome races. The Valais and Cordon chromosome races meet and hybridize at a mountain river in Les Houches (French Alps). Significant genetic structuring was recently reported among populations found on the Valais side of this hybrid zone. In this paper, a phylogenetic analysis and partial Mantel tests are used to investigate the patterns and causes of this structuring. A total of 185 shrews were trapped at 12 localities. All individuals were typed for nine microsatellite loci. Although several mountain rivers are found in the study area, riverine barriers do not have a significant influence on gene flow. Partial Mantel tests show that our result is caused by the influence of the hybrid zone with the Cordon race. The geographical patterns of this structuring are discussed in the context of the contact zone, which appears to extend up to a group of two rivers. The glacier they originate from is known to have cut the Arve valley as recently as 1818. The recent history of this glacier, its moraine and possibly rivers, may therefore be linked to the history of this hybrid zone.

Is the Gibraltar strait a barrier to gene flow for the bat Myotis myotis (Chiroptera: Vespertilionidae)?

Because of their role in limiting gene flow, geographical barriers like mountains or seas often coincide with intraspecific genetic discontinuities. Although the Strait of Gibraltar represents such a potential barrier for both plants and animals, few studies have been conducted on its impact on gene flow. Here we test this effect on a bat species (Myotis myotis) which is apparently distributed on both sides of the strait. Six colonies of 20 Myotis myotis each were sampled in southern Spain and northern Morocco along a linear transect of 1350 km. Results based on six nuclear microsatellite loci reveal no significant population structure within regions, but a complete isolation between bats sampled on each side of the strait. Variability at 600 bp of a mitochondrial gene (cytochrome b) confirms the existence of two genetically distinct and perfectly segregating clades, which diverged several million years ago. Despite the narrowness of the Gibraltar Strait (14 km), these molecular data suggest that neither males, nor females from either region have ever reproduced on the opposite side of the strait. Comparisons of molecular divergence with bats from a closely related species (M. blythii) suggest that the North African clade is possibly a distinct taxon warranting full species rank. We provisionally refer to it as Myotis cf punicus Felten 1977, but a definitive systematic understanding of the whole Mouse-eared bat species complex awaits further genetic sampling, especially in the Eastern Mediterranean areas.

Human metabolic individuality in biomedical and pharmaceutical research.

Genome-wide association studies (GWAS) have identified many risk loci for complex diseases, but effect sizes are typically small and information on the underlying biological processes is often lacking. Associations with metabolic traits as functional intermediates can overcome these problems and potentially inform individualized therapy. Here we report a comprehensive analysis of genotype-dependent metabolic phenotypes using a GWAS with non-targeted metabolomics. We identified 37 genetic loci associated with blood metabolite concentrations, of which 25 show effect sizes that are unusually high for GWAS and account for 10-60% differences in metabolite levels per allele copy. Our associations provide new functional insights for many disease-related associations that have been reported in previous studies, including those for cardiovascular and kidney disorders, type 2 diabetes, cancer, gout, venous thromboembolism and Crohn's disease. The study advances our knowledge of the genetic basis of metabolic individuality in humans and generates many new hypotheses for biomedical and pharmaceutical research.

Comparison of Strategies to Detect Epistasis from eQTL Data.

Genome-wide association studies have been instrumental in identifying genetic variants associated with complex traits such as human disease or gene expression phenotypes. It has been proposed that extending existing analysis methods by considering interactions between pairs of loci may uncover additional genetic effects. However, the large number of possible two-marker tests presents significant computational and statistical challenges. Although several strategies to detect epistasis effects have been proposed and tested for specific phenotypes, so far there has been no systematic attempt to compare their performance using real data. We made use of thousands of gene expression traits from linkage and eQTL studies, to compare the performance of different strategies. We found that using information from marginal associations between markers and phenotypes to detect epistatic effects yielded a lower false discovery rate (FDR) than a strategy solely using biological annotation in yeast, whereas results from human data were inconclusive. For future studies whose aim is to discover epistatic effects, we recommend incorporating information about marginal associations between SNPs and phenotypes instead of relying solely on biological annotation. Improved methods to discover epistatic effects will result in a more complete understanding of complex genetic effects.

Association of genetic variation with systolic and diastolic blood pressure among African Americans: the Candidate Gene Association Resource study.

The prevalence of hypertension in African Americans (AAs) is higher than in other US groups; yet, few have performed genome-wide association studies (GWASs) in AA. Among people of European descent, GWASs have identified genetic variants at 13 loci that are associated with blood pressure. It is unknown if these variants confer susceptibility in people of African ancestry. Here, we examined genome-wide and candidate gene associations with systolic blood pressure (SBP) and diastolic blood pressure (DBP) using the Candidate Gene Association Resource (CARe) consortium consisting of 8591 AAs. Genotypes included genome-wide single-nucleotide polymorphism (SNP) data utilizing the Affymetrix 6.0 array with imputation to 2.5 million HapMap SNPs and candidate gene SNP data utilizing a 50K cardiovascular gene-centric array (ITMAT-Broad-CARe [IBC] array). For Affymetrix data, the strongest signal for DBP was rs10474346 (P= 3.6 × 10(-8)) located near GPR98 and ARRDC3. For SBP, the strongest signal was rs2258119 in C21orf91 (P= 4.7 × 10(-8)). The top IBC association for SBP was rs2012318 (P= 6.4 × 10(-6)) near SLC25A42 and for DBP was rs2523586 (P= 1.3 × 10(-6)) near HLA-B. None of the top variants replicated in additional AA (n = 11 882) or European-American (n = 69 899) cohorts. We replicated previously reported European-American blood pressure SNPs in our AA samples (SH2B3, P= 0.009; TBX3-TBX5, P= 0.03; and CSK-ULK3, P= 0.0004). These genetic loci represent the best evidence of genetic influences on SBP and DBP in AAs to date. More broadly, this work supports that notion that blood pressure among AAs is a trait with genetic underpinnings but also with significant complexity.

Two high throughput technologies to detect segmental aneuploidies identify new Williams-Beuren syndrome patients with atypical deletions.

OBJECTIVE: To develop and compare two new technologies for diagnosing a contiguous gene syndrome, the Williams-Beuren syndrome (WBS). METHODS: The first proposed method, named paralogous sequence quantification (PSQ), is based on the use of paralogous sequences located on different chromosomes and quantification of specific mismatches present at these loci using pyrosequencing technology. The second exploits quantitative real time polymerase chain reaction (QPCR) to assess the relative quantity of an analysed locus. RESULTS: A correct and unambiguous diagnosis was obtained for 100% of the analysed samples with either technique (n = 165 and n = 155, respectively). These methods allowed the identification of two patients with atypical deletions in a cohort of 182 WBS patients. Both patients presented with mild facial anomalies, mild mental retardation with impaired visuospatial cognition, supravalvar aortic stenosis, and normal growth indices. These observations are consistent with the involvement of GTF2IRD1 or GTF2I in some of the WBS facial features. CONCLUSIONS: Both PSQ and QPCR are robust, easy to interpret, and simple to set up. They represent a competitive alternative for the diagnosis of segmental aneuploidies in clinical laboratories. They have advantages over fluorescence in situ hybridisation or microsatellites/SNP genotyping for detecting short segmental aneuploidies as the former is costly and labour intensive while the latter depends on the informativeness of the polymorphisms.

A multiplicity of factors contributes to selective RNA polymerase III occupancy of a subset of RNA polymerase III genes in mouse liver.

The genomic loci occupied by RNA polymerase (RNAP) III have been characterized in human culture cells by genome-wide chromatin immunoprecipitations, followed by deep sequencing (ChIP-seq). These studies have shown that only ∼40% of the annotated 622 human tRNA genes and pseudogenes are occupied by RNAP-III, and that these genes are often in open chromatin regions rich in active RNAP-II transcription units. We have used ChIP-seq to characterize RNAP-III-occupied loci in a differentiated tissue, the mouse liver. Our studies define the mouse liver RNAP-III-occupied loci including a conserved mammalian interspersed repeat (MIR) as a potential regulator of an RNAP-III subunit-encoding gene. They reveal that synteny relationships can be established between a number of human and mouse RNAP-III genes, and that the expression levels of these genes are significantly linked. They establish that variations within the A and B promoter boxes, as well as the strength of the terminator sequence, can strongly affect RNAP-III occupancy of tRNA genes. They reveal correlations with various genomic features that explain the observed variation of 81% of tRNA scores. In mouse liver, loci represented in the NCBI37/mm9 genome assembly that are clearly occupied by RNAP-III comprise 50 Rn5s (5S RNA) genes, 14 known non-tRNA RNAP-III genes, nine Rn4.5s (4.5S RNA) genes, and 29 SINEs. Moreover, out of the 433 annotated tRNA genes, half are occupied by RNAP-III. Transfer RNA gene expression levels reflect both an underlying genomic organization conserved in dividing human culture cells and resting mouse liver cells, and the particular promoter and terminator strengths of individual genes.