Quantification and discovery of sequence determinants of protein-per-mRNA amount in 29 human tissues.

Despite their importance in determining protein abundance, a comprehensive catalogue of sequence features controlling protein-to-mRNA (PTR) ratios and a quantification of their effects are still lacking. Here, we quantified PTR ratios for 11,575 proteins across 29 human tissues using matched transcriptomes and proteomes. We estimated by regression the contribution of known sequence determinants of protein synthesis and degradation in addition to 45 mRNA and 3 protein sequence motifs that we found by association testing. While PTR ratios span more than 2 orders of magnitude, our integrative model predicts PTR ratios at a median precision of 3.2-fold. A reporter assay provided functional support for two novel UTR motifs, and an immobilized mRNA affinity competition-binding assay identified motif-specific bound proteins for one motif. Moreover, our integrative model led to a new metric of codon optimality that captures the effects of codon frequency on protein synthesis and degradation. Altogether, this study shows that a large fraction of PTR ratio variation in human tissues can be predicted from sequence, and it identifies many new candidate post-transcriptional regulatory elements.

A deep proteome and transcriptome abundance atlas of 29 healthy human tissues.

Genome-, transcriptome- and proteome-wide measurements provide insights into how biological systems are regulated. However, fundamental aspects relating to which human proteins exist, where they are expressed and in which quantities are not fully understood. Therefore, we generated a quantitative proteome and transcriptome abundance atlas of 29 paired healthy human tissues from the Human Protein Atlas project representing human genes by 18,072 transcripts and 13,640 proteins including 37 without prior protein-level evidence. The analysis revealed that hundreds of proteins, particularly in testis, could not be detected even for highly expressed mRNAs, that few proteins show tissue-specific expression, that strong differences between mRNA and protein quantities within and across tissues exist and that protein expression is often more stable across tissues than that of transcripts. Only 238 of 9,848 amino acid variants found by exome sequencing could be confidently detected at the protein level showing that proteogenomics remains challenging, needs better computational methods and requires rigorous validation. Many uses of this resource can be envisaged including the study of gene/protein expression regulation and biomarker specificity evaluation.

Evolution of P2A and P5A ATPases: ancient gene duplications and the red algal connection to green plants revisited.

In a search for slowly evolving nuclear genes that may cast light on the deep evolution of plants, we carried out phylogenetic analyses of two well-characterized subfamilies of P-type pumps (P2A and P5A ATPases) from representative branches of the eukaryotic tree of life. Both P-type ATPase genes were duplicated very early in eukaryotic evolution and before the divergence of the present eukaryotic supergroups. Synapomorphies identified in the sequences provide evidence that green plants and red algae are more distantly related than are green plants and eukaryotic supergroups in which secondary or tertiary plastids are common, such as several groups belonging to the clade that includes Stramenopiles, Alveolata, Rhizaria, Cryptophyta and Haptophyta (SAR). We propose that red algae branched off soon after the first photosynthesizing eukaryote had acquired a primary plastid, while in another lineage that led to SAR, the primary plastid was lost but, in some cases, regained as a secondary or tertiary plastid.

The reversal of human phylogeny: Homo left Africa as erectus, came back as sapiens sapiens.

BACKGROUND: The molecular out of Africa hypothesis, OOAH, has been considered as an established fact amid population geneticists for some 25-30 years despite the early concern with it among phylogeneticists with experience beyond that of Homo. The palaeontological support for the hypothesis is also questionable, a circumstance that in the light of expanding Eurasian palaeontological knowledge has become accentuated through the last decades. RESULTS: The direction of evolution in the phylogenetic tree of modern humans (Homo sapiens sapiens, Hss) was established inter alia by applying progressive phylogenetic analysis to an mtDNA sampling that included a Eurasian, Lund, and the African Mbuti, San and Yoruba. The examination identified the African populations as paraphyletic, thereby compromising the OOAH. The finding, which was consistent with the out of Eurasia hypothesis, OOEH, was corroborated by the mtDNA introgression from Hss into Hsnn (Neanderthals) that demonstrated the temporal and physical Eurasian coexistence of the two lineages. The results are consistent with the palaeontologically established presence of H. erectus in Eurasia, a Eurasian divergence between H. sapiens and H. antecessor ≈ 850,000 YBP, an Hs divergence between Hss and Hsn (Neanderthals + Denisovans) ≈ 800,000 YBP, an mtDNA introgression from Hss into Hsnn* ≈ 500,000 YBP and an Eurasian divergence among the ancestors of extant Hss ≈ 250,000 YBP at the exodus of Mbuti/San into Africa. CONCLUSIONS: The present study showed that Eurasia was not the receiver but the donor in Hss evolution. The findings that Homo left Africa as erectus and returned as sapiens sapiens constitute a change in the understanding of Hs evolution to one that conforms to the extensive Eurasian record of Hs palaeontology and archaeology.

Identification of two abundant Aerococcus urinae cell wall-anchored proteins.

Aerococcus urinae is an emerging pathogen that causes urinary tract infections, bacteremia and infective endocarditis. The mechanisms through which A. urinae cause infection are largely unknown. The aims of this study were to describe the surface proteome of A. urinae and to analyse A. urinae genomes in search for genes encoding surface proteins. Two proteins, denoted Aerococcal surface protein (Asp) 1 and 2, were through the use of mass spectrometry based proteomics found to quantitatively dominate the aerococcal surface. The presence of these proteins on the surface was also shown using ELISA with serum from rabbits immunized with the recombinant Asp. These proteins had a signal sequence in the amino-terminal end and a cell wall-sorting region in the carboxy-terminal end, which contained an LPATG-motif, a hydrophobic domain and a positively charged tail. Twenty-three additional A. urinae genomes were sequenced using Illumina HiSeq technology. Six different variants of asp genes were found (denoted asp1-6). All isolates had either one or two of these asp-genes located in a conserved locus, designated Locus encoding Aerococcal Surface Proteins (LASP). The 25 genomes had in median 13 genes encoding LPXTG-proteins (range 6-24). For other Gram-positive bacteria, cell wall-anchored surface proteins with an LPXTG-motif play a key role for virulence. Thus, it will be of great interest to explore the function of the Asp proteins of A. urinae to establish a better understanding of the molecular mechanisms by which A. urinae cause disease.

Oxidation and cyclization of casbene in the biosynthesis of Euphorbia factors from mature seeds of Euphorbia lathyris L.

The seed oil of Euphorbia lathyris L. contains a series of macrocyclic diterpenoids known as Euphorbia factors. They are the current industrial source of ingenol mebutate, which is approved for the treatment of actinic keratosis, a precancerous skin condition. Here, we report an alcohol dehydrogenase-mediated cyclization step in the biosynthetic pathway of Euphorbia factors, illustrating the origin of the intramolecular carbon-carbon bonds present in lathyrane and ingenane diterpenoids. This unconventional cyclization describes the ring closure of the macrocyclic diterpene casbene. Through transcriptomic analysis of E. lathyris L. mature seeds and in planta functional characterization, we identified three enzymes involved in the cyclization route from casbene to jolkinol C, a lathyrane diterpene. These enzymes include two cytochromes P450 from the CYP71 clan and an alcohol dehydrogenase (ADH). CYP71D445 and CYP726A27 catalyze regio-specific 9-oxidation and 5-oxidation of casbene, respectively. When coupled with these P450-catalyzed monooxygenations, E. lathyris ADH1 catalyzes dehydrogenation of the hydroxyl groups, leading to the subsequent rearrangement and cyclization. The discovery of this nonconventional cyclization may provide the key link to complete elucidation of the biosynthetic pathways of ingenol mebutate and other bioactive macrocyclic diterpenoids.

The terpene synthase gene family in Tripterygium wilfordii harbors a labdane-type diterpene synthase among the monoterpene synthase TPS-b subfamily.

Tripterygium wilfordii (Celastraceae) is a medicinal plant with anti-inflammatory and immunosuppressive properties. Identification of a vast array of unusual sesquiterpenoids, diterpenoids and triterpenoids in T. wilfordii has spurred investigations of their pharmacological properties. The tri-epoxide lactone triptolide was the first of many diterpenoids identified, attracting interest due to the spectrum of bioactivities. To probe the genetic underpinning of diterpenoid diversity, an expansion of the class II diterpene synthase (diTPS) family was recently identified in a leaf transcriptome. Following detection of triptolide and simple diterpene scaffolds in the root, we sequenced and mined the root transcriptome. This allowed identification of the root-specific complement of TPSs and an expansion in the class I diTPS family. Functional characterization of the class II diTPSs established their activities in the formation of four C-20 diphosphate intermediates, precursors of both generalized and specialized metabolism and a novel scaffold for Celastraceae. Functional pairs of the class I and II enzymes resulted in formation of three scaffolds, accounting for some of the terpenoid diversity found in T. wilfordii. The absence of activity-forming abietane-type diterpenes encouraged further testing of TPSs outside the canonical class I diTPS family. TwTPS27, close relative of mono-TPSs, was found to couple with TwTPS9, converting normal-copalyl diphosphate to miltiradiene. The phylogenetic distance to established diTPSs indicates neo-functionalization of TwTPS27 into a diTPS, a function not previously observed in the TPS-b subfamily. This example of evolutionary convergence expands the functionality of TPSs in the TPS-b family and may contribute miltiradiene to the diterpenoids of T. wilfordii.

Personal model-assisted identification of NAD+ and glutathione metabolism as intervention target in NAFLD.

To elucidate the molecular mechanisms underlying non-alcoholic fatty liver disease (NAFLD), we recruited 86 subjects with varying degrees of hepatic steatosis (HS). We obtained experimental data on lipoprotein fluxes and used these individual measurements as personalized constraints of a hepatocyte genome-scale metabolic model to investigate metabolic differences in liver, taking into account its interactions with other tissues. Our systems level analysis predicted an altered demand for NAD+ and glutathione (GSH) in subjects with high HS Our analysis and metabolomic measurements showed that plasma levels of glycine, serine, and associated metabolites are negatively correlated with HS, suggesting that these GSH metabolism precursors might be limiting. Quantification of the hepatic expression levels of the associated enzymes further pointed to altered de novo GSH synthesis. To assess the effect of GSH and NAD+ repletion on the development of NAFLD, we added precursors for GSH and NAD+ biosynthesis to the Western diet and demonstrated that supplementation prevents HS in mice. In a proof-of-concept human study, we found improved liver function and decreased HS after supplementation with serine (a precursor to glycine) and hereby propose a strategy for NAFLD treatment.

Microfluidic screening and whole-genome sequencing identifies mutations associated with improved protein secretion by yeast.

There is an increasing demand for biotech-based production of recombinant proteins for use as pharmaceuticals in the food and feed industry and in industrial applications. Yeast Saccharomyces cerevisiae is among preferred cell factories for recombinant protein production, and there is increasing interest in improving its protein secretion capacity. Due to the complexity of the secretory machinery in eukaryotic cells, it is difficult to apply rational engineering for construction of improved strains. Here we used high-throughput microfluidics for the screening of yeast libraries, generated by UV mutagenesis. Several screening and sorting rounds resulted in the selection of eight yeast clones with significantly improved secretion of recombinant α-amylase. Efficient secretion was genetically stable in the selected clones. We performed whole-genome sequencing of the eight clones and identified 330 mutations in total. Gene ontology analysis of mutated genes revealed many biological processes, including some that have not been identified before in the context of protein secretion. Mutated genes identified in this study can be potentially used for reverse metabolic engineering, with the objective to construct efficient cell factories for protein secretion. The combined use of microfluidics screening and whole-genome sequencing to map the mutations associated with the improved phenotype can easily be adapted for other products and cell types to identify novel engineering targets, and this approach could broadly facilitate design of novel cell factories.

Genomic Variation in IbA10G2 and Other Patient-Derived Cryptosporidium hominis Subtypes.

In order to improve genotyping and epidemiological analysis of Cryptosporidium spp., genomic data need to be generated directly from a broad range of clinical specimens. Utilizing a robust method that we developed for the purification and generation of amplified target DNA, we present its application for the successful isolation and whole-genome sequencing of 14 different Cryptosporidium hominis patient specimens. Six isolates of subtype IbA10G2 were analyzed together with a single representative each of 8 other subtypes: IaA20R3, IaA23R3, IbA9G3, IbA13G3, IdA14, IeA11G3T3, IfA12G1, and IkA18G1. Parasite burden was measured over a range of more than 2 orders of magnitude for all samples, while the genomes were sequenced to mean depths of between 17× and 490× coverage. Sequence homology-based functional annotation identified several genes of interest, including the gene encoding Cryptosporidium oocyst wall protein 9 (COWP9), which presented a predicted loss-of-function mutation in all the sequence subtypes, except for that seen with IbA10G2, which has a sequence identical to the Cryptosporidium parvum reference Iowa II sequence. Furthermore, phylogenetic analysis showed that all the IbA10G2 genomes form a monophyletic clade in the C. hominis tree as expected and yet display some heterogeneity within the IbA10G2 subtype. The current report validates the aforementioned method for isolating and sequencing Cryptosporidium directly from clinical stool samples. In addition, the analysis demonstrates the potential in mining data generated from sequencing multiple whole genomes of Cryptosporidium from human fecal samples, while alluding to the potential for a higher degree of genotyping within Cryptosporidium epidemiology.

Evolutionary engineering reveals divergent paths when yeast is adapted to different acidic environments.

Tolerance of yeast to acid stress is important for many industrial processes including organic acid production. Therefore, elucidating the molecular basis of long term adaptation to acidic environments will be beneficial for engineering production strains to thrive under such harsh conditions. Previous studies using gene expression analysis have suggested that both organic and inorganic acids display similar responses during short term exposure to acidic conditions. However, biological mechanisms that will lead to long term adaptation of yeast to acidic conditions remains unknown and whether these mechanisms will be similar for tolerance to both organic and inorganic acids is yet to be explored. We therefore evolved Saccharomyces cerevisiae to acquire tolerance to HCl (inorganic acid) and to 0.3M L-lactic acid (organic acid) at pH 2.8 and then isolated several low pH tolerant strains. Whole genome sequencing and RNA-seq analysis of the evolved strains revealed different sets of genome alterations suggesting a divergence in adaptation to these two acids. An altered sterol composition and impaired iron uptake contributed to HCl tolerance whereas the formation of a multicellular morphology and rapid lactate degradation was crucial for tolerance to high concentrations of lactic acid. Our findings highlight the contribution of both the selection pressure and nature of the acid as a driver for directing the evolutionary path towards tolerance to low pH. The choice of carbon source was also an important factor in the evolutionary process since cells evolved on two different carbon sources (raffinose and glucose) generated a different set of mutations in response to the presence of lactic acid. Therefore, different strategies are required for a rational design of low pH tolerant strains depending on the acid of interest.

Transcriptomics resources of human tissues and organs.

Quantifying the differential expression of genes in various human organs, tissues, and cell types is vital to understand human physiology and disease. Recently, several large-scale transcriptomics studies have analyzed the expression of protein-coding genes across tissues. These datasets provide a framework for defining the molecular constituents of the human body as well as for generating comprehensive lists of proteins expressed across tissues or in a tissue-restricted manner. Here, we review publicly available human transcriptome resources and discuss body-wide data from independent genome-wide transcriptome analyses of different tissues. Gene expression measurements from these independent datasets, generated using samples from fresh frozen surgical specimens and postmortem tissues, are consistent. Overall, the different genome-wide analyses support a distribution in which many proteins are found in all tissues and relatively few in a tissue-restricted manner. Moreover, we discuss the applications of publicly available omics data for building genome-scale metabolic models, used for analyzing cell and tissue functions both in physiological and in disease contexts.

Gene-specific correlation of RNA and protein levels in human cells and tissues.

An important issue for molecular biology is to establish whether transcript levels of a given gene can be used as proxies for the corresponding protein levels. Here, we have developed a targeted proteomics approach for a set of human non-secreted proteins based on parallel reaction monitoring to measure, at steady-state conditions, absolute protein copy numbers across human tissues and cell lines and compared these levels with the corresponding mRNA levels using transcriptomics. The study shows that the transcript and protein levels do not correlate well unless a gene-specific RNA-to-protein (RTP) conversion factor independent of the tissue type is introduced, thus significantly enhancing the predictability of protein copy numbers from RNA levels. The results show that the RTP ratio varies significantly with a few hundred copies per mRNA molecule for some genes to several hundred thousands of protein copies per mRNA molecule for others. In conclusion, our data suggest that transcriptome analysis can be used as a tool to predict the protein copy numbers per cell, thus forming an attractive link between the field of genomics and proteomics.

Complementing tissue characterization by integrating transcriptome profiling from the Human Protein Atlas and from the FANTOM5 consortium.

Understanding the normal state of human tissue transcriptome profiles is essential for recognizing tissue disease states and identifying disease markers. Recently, the Human Protein Atlas and the FANTOM5 consortium have each published extensive transcriptome data for human samples using Illumina-sequenced RNA-Seq and Heliscope-sequenced CAGE. Here, we report on the first large-scale complex tissue transcriptome comparison between full-length versus 5'-capped mRNA sequencing data. Overall gene expression correlation was high between the 22 corresponding tissues analyzed (R > 0.8). For genes ubiquitously expressed across all tissues, the two data sets showed high genome-wide correlation (91% agreement), with differences observed for a small number of individual genes indicating the need to update their gene models. Among the identified single-tissue enriched genes, up to 75% showed consensus of 7-fold enrichment in the same tissue in both methods, while another 17% exhibited multiple tissue enrichment and/or high expression variety in the other data set, likely dependent on the cell type proportions included in each tissue sample. Our results show that RNA-Seq and CAGE tissue transcriptome data sets are highly complementary for improving gene model annotations and highlight biological complexities within tissue transcriptomes. Furthermore, integration with image-based protein expression data is highly advantageous for understanding expression specificities for many genes.

Solid-phase cloning for high-throughput assembly of single and multiple DNA parts.

We describe solid-phase cloning (SPC) for high-throughput assembly of expression plasmids. Our method allows PCR products to be put directly into a liquid handler for capture and purification using paramagnetic streptavidin beads and conversion into constructs by subsequent cloning reactions. We present a robust automated protocol for restriction enzyme based SPC and its performance for the cloning of >60 000 unique human gene fragments into expression vectors. In addition, we report on SPC-based single-strand assembly for applications where exact control of the sequence between fragments is needed or where multiple inserts are to be assembled. In this approach, the solid support allows for head-to-tail assembly of DNA fragments based on hybridization and polymerase fill-in. The usefulness of head-to-tail SPC was demonstrated by assembly of >150 constructs with up to four DNA parts at an average success rate above 80%. We report on several applications for SPC and we suggest it to be particularly suitable for high-throughput efforts using laboratory workstations.

Cryptosporidium as a testbed for single cell genome characterization of unicellular eukaryotes.

BACKGROUND: Infectious disease involving multiple genetically distinct populations of pathogens is frequently concurrent, but difficult to detect or describe with current routine methodology. Cryptosporidium sp. is a widespread gastrointestinal protozoan of global significance in both animals and humans. It cannot be easily maintained in culture and infections of multiple strains have been reported. To explore the potential use of single cell genomics methodology for revealing genome-level variation in clinical samples from Cryptosporidium-infected hosts, we sorted individual oocysts for subsequent genome amplification and full-genome sequencing. RESULTS: Cells were identified with fluorescent antibodies with an 80 % success rate for the entire single cell genomics workflow, demonstrating that the methodology can be applied directly to purified fecal samples. Ten amplified genomes from sorted single cells were selected for genome sequencing and compared both to the original population and a reference genome in order to evaluate the accuracy and performance of the method. Single cell genome coverage was on average 81 % even with a moderate sequencing effort and by combining the 10 single cell genomes, the full genome was accounted for. By a comparison to the original sample, biological variation could be distinguished and separated from noise introduced in the amplification. CONCLUSIONS: As a proof of principle, we have demonstrated the power of applying single cell genomics to dissect infectious disease caused by closely related parasite species or subtypes. The workflow can easily be expanded and adapted to target other protozoans, and potential applications include mapping genome-encoded traits, virulence, pathogenicity, host specificity and resistance at the level of cells as truly meaningful biological units.

Evolutionary histories of transposable elements in the genome of the largest living marsupial carnivore, the Tasmanian devil.

The largest living carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii), is the sole survivor of a lineage originating about 12 Ma. We set out to investigate the spectrum of transposable elements found in the Tasmanian devil genome, the first high-coverage genome of an Australian marsupial. Marsupial genomes have been shown to have the highest amount of transposable elements among vertebrates. We analyzed the horizontally transmitted DNA transposons OC1 and hAT-1_MEu in the Tasmanian devil genome. OC1 is present in all carnivorous marsupials, while having a very limited distribution among the remaining Australian marsupial orders. In contrast, hAT-1_MEu is present in all Australian marsupial orders, and has so far only been identified in a few placental mammals. We screened 158 introns for phylogenetically informative retrotransposons in the order Dasyuromorphia, and found that the youngest SINE (Short INterspersed Element), WSINE1, is no longer active in the subfamily Dasyuridae. The lack of detectable WSINE1 activity in this group may be due to a retrotransposon inactivation event approximately 30 Ma. We found that the Tasmanian devil genome contains a relatively low number of continuous full-length LINE-1 (Long INterspersed Element 1, L1) retrotransposons compared with the opossum genome. Furthermore, all L1 elements in the Tasmanian devil appeared to be nonfunctional. Hidden Markov Model approaches suggested that other potential sources of functional reverse transcriptase are absent from the genome. We discuss the issues associated with assembling long, highly similar L1 copies from short read Illumina data and describe how assembly artifacts can potentially lead to erroneous conclusions.

Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics.

Global classification of the human proteins with regards to spatial expression patterns across organs and tissues is important for studies of human biology and disease. Here, we used a quantitative transcriptomics analysis (RNA-Seq) to classify the tissue-specific expression of genes across a representative set of all major human organs and tissues and combined this analysis with antibody-based profiling of the same tissues. To present the data, we launch a new version of the Human Protein Atlas that integrates RNA and protein expression data corresponding to ∼80% of the human protein-coding genes with access to the primary data for both the RNA and the protein analysis on an individual gene level. We present a classification of all human protein-coding genes with regards to tissue-specificity and spatial expression pattern. The integrative human expression map can be used as a starting point to explore the molecular constituents of the human body.

The human cardiac and skeletal muscle proteomes defined by transcriptomics and antibody-based profiling.

BACKGROUND: To understand cardiac and skeletal muscle function, it is important to define and explore their molecular constituents and also to identify similarities and differences in the gene expression in these two different striated muscle tissues. Here, we have investigated the genes and proteins with elevated expression in cardiac and skeletal muscle in relation to all other major human tissues and organs using a global transcriptomics analysis complemented with antibody-based profiling to localize the corresponding proteins on a single cell level. RESULTS: Our study identified a comprehensive list of genes expressed in cardiac and skeletal muscle. The genes with elevated expression were further stratified according to their global expression pattern across the human body as well as their precise localization in the muscle tissues. The functions of the proteins encoded by the elevated genes are well in line with the physiological functions of cardiac and skeletal muscle, such as contraction, ion transport, regulation of membrane potential and actomyosin structure organization. A large fraction of the transcripts in both cardiac and skeletal muscle correspond to mitochondrial proteins involved in energy metabolism, which demonstrates the extreme specialization of these muscle tissues to provide energy for contraction. CONCLUSIONS: Our results provide a comprehensive list of genes and proteins elevated in striated muscles. A number of proteins not previously characterized in cardiac and skeletal muscle were identified and localized to specific cellular subcompartments. These proteins represent an interesting starting point for further functional analysis of their role in muscle biology and disease.

Brown and polar bear Y chromosomes reveal extensive male-biased gene flow within brother lineages.

Brown and polar bears have become prominent examples in phylogeography, but previous phylogeographic studies relied largely on maternally inherited mitochondrial DNA (mtDNA) or were geographically restricted. The male-specific Y chromosome, a natural counterpart to mtDNA, has remained underexplored. Although this paternally inherited chromosome is indispensable for comprehensive analyses of phylogeographic patterns, technical difficulties and low variability have hampered its application in most mammals. We developed 13 novel Y-chromosomal sequence and microsatellite markers from the polar bear genome and screened these in a broad geographic sample of 130 brown and polar bears. We also analyzed a 390-kb-long Y-chromosomal scaffold using sequencing data from published male ursine genomes. Y chromosome evidence support the emerging understanding that brown and polar bears started to diverge no later than the Middle Pleistocene. Contrary to mtDNA patterns, we found 1) brown and polar bears to be reciprocally monophyletic sister (or rather brother) lineages, without signals of introgression, 2) male-biased gene flow across continents and on phylogeographic time scales, and 3) male dispersal that links the Alaskan ABC islands population to mainland brown bears. Due to female philopatry, mtDNA provides a highly structured estimate of population differentiation, while male-biased gene flow is a homogenizing force for nuclear genetic variation. Our findings highlight the importance of analyzing both maternally and paternally inherited loci for a comprehensive view of phylogeographic history, and that mtDNA-based phylogeographic studies of many mammals should be reevaluated. Recent advances in sequencing technology render the analysis of Y-chromosomal variation feasible, even in nonmodel organisms.