An integrated spatio-temporal view of riverine biodiversity using environmental DNA metabarcoding.

Anthropogenically forced changes in global freshwater biodiversity demand more efficient monitoring approaches. Consequently, environmental DNA (eDNA) analysis is enabling ecosystem-scale biodiversity assessment, yet the appropriate spatio-temporal resolution of robust biodiversity assessment remains ambiguous. Here, using intensive, spatio-temporal eDNA sampling across space (five rivers in Europe and North America, with an upper range of 20-35 km between samples), time (19 timepoints between 2017 and 2018) and environmental conditions (river flow, pH, conductivity, temperature and rainfall), we characterise the resolution at which information on diversity across the animal kingdom can be gathered from rivers using eDNA. In space, beta diversity was mainly dictated by turnover, on a scale of tens of kilometres, highlighting that diversity measures are not confounded by eDNA from upstream. Fish communities showed nested assemblages along some rivers, coinciding with habitat use. Across time, seasonal life history events, including salmon and eel migration, were detected. Finally, effects of environmental conditions were taxon-specific, reflecting habitat filtering of communities rather than effects on DNA molecules. We conclude that riverine eDNA metabarcoding can measure biodiversity at spatio-temporal scales relevant to species and community ecology, demonstrating its utility in delivering insights into river community ecology during a time of environmental change.

Long-read genome sequencing and variant reanalysis increase diagnostic yield in neurodevelopmental disorders.

Variant detection from long-read genome sequencing (lrGS) has proven to be considerably more accurate and comprehensive than variant detection from short-read genome sequencing (srGS). However, the rate at which lrGS can increase molecular diagnostic yield for rare disease is not yet precisely characterized. We performed lrGS using Pacific Biosciences "HiFi" technology on 96 short-read-negative probands with rare disease that were suspected to be genetic. We generated hg38-aligned variants and de novo phased genome assemblies, and subsequently annotated, filtered, and curated variants using clinical standards. New disease-relevant or potentially relevant genetic findings were identified in 16/96 (16.7%) probands, eight of which (8/96, 8.33%) harbored pathogenic or likely pathogenic variants. Newly identified variants were visible in both srGS and lrGS in nine probands (~9.4%) and resulted from changes to interpretation mostly from recent gene-disease association discoveries. Seven cases included variants that were only interpretable in lrGS, including copy-number variants, an inversion, a mobile element insertion, two low-complexity repeat expansions, and a 1 bp deletion. While evidence for each of these variants is, in retrospect, visible in srGS, they were either: not called within srGS data, were represented by calls with incorrect sizes or structures, or failed quality-control and filtration. Thus, while reanalysis of older data clearly increases diagnostic yield, we find that lrGS allows for substantial additional yield (7/96, 7.3%) beyond srGS. We anticipate that as lrGS analysis improves, and as lrGS datasets grow allowing for better variant frequency annotation, the additional lrGS-only rare disease yield will grow over time.

Biosynthesis of Haloterpenoids in Red Algae via Microbial-like Type I Terpene Synthases.

Red algae or seaweeds produce highly distinctive halogenated terpenoid compounds, including the pentabromochlorinated monoterpene halomon that was once heralded as a promising anticancer agent. The first dedicated step in the biosynthesis of these natural product molecules is expected to be catalyzed by terpene synthase (TS) enzymes. Recent work has demonstrated an emerging class of type I TSs in red algal terpene biosynthesis. However, only one such enzyme from a notoriously haloterpenoid-producing red alga (Laurencia pacifica) has been functionally characterized and the product structure is not related to halogenated terpenoids. Herein, we report 10 new type I TSs from the red algae Portieria hornemannii, Plocamium pacificum, L. pacifica, and Laurencia subopposita that produce a diversity of halogenated mono- and sesquiterpenes. We used a combination of genome sequencing, terpenoid metabolomics, in vitro biochemistry, and bioinformatics to establish red algal TSs in all four species, including those associated with the selective production of key halogenated terpene precursors myrcene, trans-ß-ocimene, and germacrene D-4-ol. These results expand on a small but growing number of characterized red algal TSs and offer insight into the biosynthesis of iconic halogenated algal compounds that are not without precedence elsewhere in biology.

Convergent evolution and horizontal gene transfer in Arctic Ocean microalgae.

Microbial communities in the world ocean are affected strongly by oceanic circulation, creating characteristic marine biomes. The high connectivity of most of the ocean makes it difficult to disentangle selective retention of colonizing genotypes (with traits suited to biome specific conditions) from evolutionary selection, which would act on founder genotypes over time. The Arctic Ocean is exceptional with limited exchange with other oceans and ice covered since the last ice age. To test whether Arctic microalgal lineages evolved apart from algae in the global ocean, we sequenced four lineages of microalgae isolated from Arctic waters and sea ice. Here we show convergent evolution and highlight geographically limited HGT as an ecological adaptive force in the form of PFAM complements and horizontal acquisition of key adaptive genes. Notably, ice-binding proteins were acquired and horizontally transferred among Arctic strains. A comparison with Tara Oceans metagenomes and metatranscriptomes confirmed mostly Arctic distributions of these IBPs. The phylogeny of Arctic-specific genes indicated that these events were independent of bacterial-sourced HGTs in Antarctic Southern Ocean microalgae.

A generalist-specialist trade-off between switchgrass cytotypes impacts climate adaptation and geographic range.

Polyploidy results from whole-genome duplication and is a unique form of heritable variation with pronounced evolutionary implications. Different ploidy levels, or cytotypes, can exist within a single species, and such systems provide an opportunity to assess how ploidy variation alters phenotypic novelty, adaptability, and fitness, which can, in turn, drive the development of unique ecological niches that promote the coexistence of multiple cytotypes. Switchgrass, Panicum virgatum, is a widespread, perennial C4 grass in North America with multiple naturally occurring cytotypes, primarily tetraploids (4×) and octoploids (8×). Using a combination of genomic, quantitative genetic, landscape, and niche modeling approaches, we detect divergent levels of genetic admixture, evidence of niche differentiation, and differential environmental sensitivity between switchgrass cytotypes. Taken together, these findings support a generalist (8×)­specialist (4×) trade-off. Our results indicate that the 8× represent a unique combination of genetic variation that has allowed the expansion of switchgrass' ecological niche and thus putatively represents a valuable breeding resource.

A draft genome provides hypotheses on drought tolerance in a keystone plant species in Western North America threatened by climate change.

Climate change presents distinct ecological and physiological challenges to plants as extreme climate events become more common. Understanding how species have adapted to drought, especially ecologically important nonmodel organisms, will be crucial to elucidate potential biological pathways for drought adaptation and inform conservation strategies. To aid in genome-to-phenome research, a draft genome was assembled for a diploid individual of Artemisia tridentata subsp. tridentata, a threatened keystone shrub in western North America. While this taxon has few genetic resources available and genetic/genomics work has proven difficult due to genetic heterozygosity in the past, a draft genome was successfully assembled. Aquaporin (AQP) genes and their promoter sequences were mined from the draft genome to predict mechanisms regulating gene expression and generate hypotheses on key genes underpinning drought response. Fifty-one AQP genes were fully assembled within the draft genome. Promoter and phylogenetic analyses revealed putative duplicates of A. tridentata subsp. tridentata AQPs which have experienced differentiation in promoter elements, potentially supporting novel biological pathways. Comparison with nondrought-tolerant congener supports enrichments of AQP genes in this taxon during adaptation to drought stress. Differentiation of promoter elements revealed that paralogues of some genes have evolved to function in different pathways, highlighting these genes as potential candidates for future research and providing critical hypotheses for future genome-to-phenome work.

Long-read genome sequencing for the molecular diagnosis of neurodevelopmental disorders.

Exome and genome sequencing have proven to be effective tools for the diagnosis of neurodevelopmental disorders (NDDs), but large fractions of NDDs cannot be attributed to currently detectable genetic variation. This is likely, at least in part, a result of the fact that many genetic variants are difficult or impossible to detect through typical short-read sequencing approaches. Here, we describe a genomic analysis using Pacific Biosciences circular consensus sequencing (CCS) reads, which are both long (>10 kb) and accurate (>99% bp accuracy). We used CCS on six proband-parent trios with NDDs that were unexplained despite extensive testing, including genome sequencing with short reads. We identified variants and created de novo assemblies in each trio, with global metrics indicating these datasets are more accurate and comprehensive than those provided by short-read data. In one proband, we identified a likely pathogenic (LP), de novo L1-mediated insertion in CDKL5 that results in duplication of exon 3, leading to a frameshift. In a second proband, we identified multiple large de novo structural variants, including insertion-translocations affecting DGKB and MLLT3, which we show disrupt MLLT3 transcript levels. We consider this extensive structural variation likely pathogenic. The breadth and quality of variant detection, coupled to finding variants of clinical and research interest in two of six probands with unexplained NDDs, support the hypothesis that long-read genome sequencing can substantially improve rare disease genetic discovery rates.

Assessing the impact of the threatened crucian carp (Carassius carassius) on pond invertebrate diversity: A comparison of conventional and molecular tools.

Fishes stocked for recreation and angling can damage freshwater habitats and negatively impact biodiversity. The pond-associated crucian carp (Carassius carassius) is rare across Europe and is stocked for conservation management in England, but its impacts on pond biota are understudied. Freshwater invertebrates contribute substantially to aquatic biodiversity, encompassing many rare and endemic species, but their small size and high abundance complicate their assessment. Practitioners have employed sweep-netting and kick-sampling with microscopy (morphotaxonomy), but specimen size/quality and experience can bias identification. DNA and environmental DNA (eDNA) metabarcoding offer alternative means of invertebrate assessment. We compared invertebrate diversity in ponds (N = 18) with and without crucian carp using morphotaxonomic identification, DNA metabarcoding and eDNA metabarcoding. Five 2 L water samples and 3 min sweep-net samples were collected at each pond. Inventories produced by morphotaxonomic identification of netted samples, DNA metabarcoding of bulk tissue samples and eDNA metabarcoding of water samples were compared. Alpha diversity was greatest with DNA or eDNA metabarcoding, depending on whether standard or unbiased methods were considered. DNA metabarcoding reflected morphotaxonomic identification, whereas eDNA metabarcoding produced markedly different communities. These complementary tools should be combined for comprehensive invertebrate assessment. Crucian carp presence minimally reduced alpha diversity in ponds, but positively influenced beta diversity through taxon turnover (i.e., ponds with crucian carp contained different invertebrates to fishless ponds). Crucian carp presence contributes to landscape-scale invertebrate diversity, supporting continued conservation management in England. Our results show that molecular tools can enhance freshwater invertebrate assessment and facilitate development of more accurate and ecologically effective pond management strategies.

Targeted and passive environmental DNA approaches outperform established methods for detection of quagga mussels, Dreissena rostriformis bugensis in flowing water.

The early detection of invasive non-native species (INNS) is important for informing management actions. Established monitoring methods require the collection or observation of specimens, which is unlikely at the beginning of an invasion when densities are likely to be low. Environmental DNA (eDNA) analysis is a highly promising technique for the detection of INNS-particularly during the early stages of an invasion.Here, we compared the use of traditional kick-net sampling with two eDNA approaches (targeted detection using both conventional and quantitative PCR and passive detection via metabarcoding with conserved primers) for detection of quagga mussel, Dreissena rostriformis bugensis, a high priority INNS, along a density gradient on the River Wraysbury, UK.All three molecular tools outperformed traditional sampling in terms of detection. Conventional PCR and qPCR both had 100% detection rate in all samples and outperformed metabarcoding when the target species was at low densities. Additionally, quagga mussel DNA copy number (qPCR) and relative read count (metabarcoding) were significantly influenced by both mussel density and distance from source population, with distance being the most significant predictor. Synthesis and application. All three molecular approaches were more sensitive than traditional kick-net sampling for the detection of the quagga mussel in flowing water, and both qPCR and metabarcoding enabled estimates of relative abundance. Targeted approaches were more sensitive than metabarcoding, but metabarcoding has the advantage of providing information on the wider community and consequently the impacts of INNS.

The Genomic Basis of Color Pattern Polymorphism in the Harlequin Ladybird.

Many animal species comprise discrete phenotypic forms. A common example in natural populations of insects is the occurrence of different color patterns, which has motivated a rich body of ecological and genetic research [1-6]. The occurrence of dark, i.e., melanic, forms displaying discrete color patterns is found across multiple taxa, but the underlying genomic basis remains poorly characterized. In numerous ladybird species (Coccinellidae), the spatial arrangement of black and red patches on adult elytra varies wildly within species, forming strikingly different complex color patterns [7, 8]. In the harlequin ladybird, Harmonia axyridis, more than 200 distinct color forms have been described, which classic genetic studies suggest result from allelic variation at a single, unknown, locus [9, 10]. Here, we combined whole-genome sequencing, population-based genome-wide association studies, gene expression, and functional analyses to establish that the transcription factor Pannier controls melanic pattern polymorphism in H. axyridis. We show that pannier is necessary for the formation of melanic elements on the elytra. Allelic variation in pannier leads to protein expression in distinct domains on the elytra and thus determines the distinct color patterns in H. axyridis. Recombination between pannier alleles may be reduced by a highly divergent sequence of ∼170 kb in the cis-regulatory regions of pannier, with a 50 kb inversion between color forms. This most likely helps maintain the distinct alleles found in natural populations. Thus, we propose that highly variable discrete color forms can arise in natural populations through cis-regulatory allelic variation of a single gene.

Needle in a haystack? A comparison of eDNA metabarcoding and targeted qPCR for detection of the great crested newt (Triturus cristatus).

Environmental DNA (eDNA) analysis is a rapid, cost-effective, non-invasive biodiversity monitoring tool which utilises DNA left behind in the environment by organisms for species detection. The method is used as a species-specific survey tool for rare or invasive species across a broad range of ecosystems. Recently, eDNA and "metabarcoding" have been combined to describe whole communities rather than focusing on single target species. However, whether metabarcoding is as sensitive as targeted approaches for rare species detection remains to be evaluated. The great crested newt Triturus cristatus is a flagship pond species of international conservation concern and the first UK species to be routinely monitored using eDNA. We evaluate whether eDNA metabarcoding has comparable sensitivity to targeted real-time quantitative PCR (qPCR) for T. cristatus detection. Extracted eDNA samples (N = 532) were screened for T. cristatus by qPCR and analysed for all vertebrate species using high-throughput sequencing technology. With qPCR and a detection threshold of 1 of 12 positive qPCR replicates, newts were detected in 50% of ponds. Detection decreased to 32% when the threshold was increased to 4 of 12 positive qPCR replicates. With metabarcoding, newts were detected in 34% of ponds without a detection threshold, and in 28% of ponds when a threshold (0.028%) was applied. Therefore, qPCR provided greater detection than metabarcoding but metabarcoding detection with no threshold was equivalent to qPCR with a stringent detection threshold. The proportion of T. cristatus sequences in each sample was positively associated with the number of positive qPCR replicates (qPCR score) suggesting eDNA metabarcoding may be indicative of eDNA concentration. eDNA metabarcoding holds enormous potential for holistic biodiversity assessment and routine freshwater monitoring. We advocate this community approach to freshwater monitoring to guide management and conservation, whereby entire communities can be initially surveyed to best inform use of funding and time for species-specific surveys.

The effect of filtration method on the efficiency of environmental DNA capture and quantification via metabarcoding.

Environmental DNA (eDNA) is a promising tool for rapid and noninvasive biodiversity monitoring. eDNA density is low in environmental samples, and a capture method, such as filtration, is often required to concentrate eDNA for downstream analyses. In this study, six treatments, with differing filter types and pore sizes for eDNA capture, were compared for their efficiency and accuracy to assess fish community structure with known fish abundance and biomass via eDNA metabarcoding. Our results showed that different filters (with the exception of 20-µm large-pore filters) were broadly consistent in their DNA capture ability. The 0.45-µm filters performed the best in terms of total DNA yield, probability of species detection, repeatability within pond and consistency between ponds. However performance of 0.45-µm filters was only marginally better than for 0.8-µm filters, while filtration time was significantly longer. Given this trade-off, the 0.8-µm filter is the optimal pore size of membrane filter for turbid, eutrophic and high fish density ponds analysed here. The 0.45-µm Sterivex enclosed filters performed reasonably well and are suitable in situations where on-site filtration is required. Finally, prefilters are applied only if absolutely essential for reducing the filtration time or increasing the throughput volume of the capture filters. In summary, we found encouraging similarity in the results obtained from different filtration methods, but the optimal pore size of filter or filter type might strongly depend on the water type under study.

Genetic evidence challenges the native status of a threatened freshwater fish (Carassius carassius) in England.

A fundamental consideration for the conservation of a species is the extent of its native range, that is, regions naturally colonized. However, both natural processes and human-mediated introductions can drive species distribution shifts. Ruling out the human-mediated introduction of a species into a given region is vital for its conservation, but remains a significant challenge in most cases. The crucian carp Carassius carassius (L.) is a threatened freshwater fish thought to be native to much of Europe. However, its native status in England is based only on anecdotal evidence. Here, we devise an approach that can be used to empirically test the native status of English fauna. We use this approach, along with 13 microsatellite loci, population structure analyses, and Approximate Bayesian Computation (ABC), to test hypotheses for the origins of C. carassius in England. Contrary to the current consensus, we find strong support for the human-mediated introduction of C. carassius into England during the 15th century. This result stimulates an interesting and timely debate surrounding motivations for the conservation of species. We discuss this topic, and the potential for continued conservation of C. carassius in England, despite its non-native origins.

Environmental DNA metabarcoding of lake fish communities reflects long-term data from established survey methods.

Organisms continuously release DNA into their environments via shed cells, excreta, gametes and decaying material. Analysis of this 'environmental DNA' (eDNA) is revolutionizing biodiversity monitoring. eDNA outperforms many established survey methods for targeted detection of single species, but few studies have investigated how well eDNA reflects whole communities of organisms in natural environments. We investigated whether eDNA can recover accurate qualitative and quantitative information about fish communities in large lakes, by comparison to the most comprehensive long-term gill-net data set available in the UK. Seventy-eight 2L water samples were collected along depth profile transects, gill-net sites and from the shoreline in three large, deep lakes (Windermere, Bassenthwaite Lake and Derwent Water) in the English Lake District. Water samples were assayed by eDNA metabarcoding of the mitochondrial 12S and cytochrome b regions. Fourteen of the 16 species historically recorded in Windermere were detected using eDNA, compared to four species in the most recent gill-net survey, demonstrating eDNA is extremely sensitive for detecting species. A key question for biodiversity monitoring is whether eDNA can accurately estimate abundance. To test this, we used the number of sequence reads per species and the proportion of sampling sites in which a species was detected with eDNA (i.e. site occupancy) as proxies for abundance. eDNA abundance data consistently correlated with rank abundance estimates from established surveys. These results demonstrate that eDNA metabarcoding can describe fish communities in large lakes, both qualitatively and quantitatively, and has great potential as a complementary tool to established monitoring methods.

Comparing RADseq and microsatellites to infer complex phylogeographic patterns, an empirical perspective in the Crucian carp, Carassius carassius, L.

The conservation of threatened species must be underpinned by phylogeographic knowledge. This need is epitomized by the freshwater fish Carassius carassius, which is in decline across much of its European range. Restriction site-associated DNA sequencing (RADseq) is increasingly used for such applications; however, RADseq is expensive, and limitations on sample number must be weighed against the benefit of large numbers of markers. This trade-off has previously been examined using simulation studies; however, empirical comparisons between these markers, especially in a phylogeographic context, are lacking. Here, we compare the results from microsatellites and RADseq for the phylogeography of C. carassius to test whether it is more advantageous to genotype fewer markers (microsatellites) in many samples, or many markers (SNPs) in fewer samples. These data sets, along with data from the mitochondrial cytochrome b gene, agree on broad phylogeographic patterns, showing the existence of two previously unidentified C. carassius lineages in Europe: one found throughout northern and central-eastern European drainages and a second almost exclusively confined to the Danubian catchment. These lineages have been isolated for approximately 2.15 m years and should be considered separate conservation units. RADseq recovered finer population structure and stronger patterns of IBD than microsatellites, despite including only 17.6% of samples (38% of populations and 52% of samples per population). RADseq was also used along with approximate Bayesian computation to show that the postglacial colonization routes of C. carassius differ from the general patterns of freshwater fish in Europe, likely as a result of their distinctive ecology.

The globalization of naval provisioning: ancient DNA and stable isotope analyses of stored cod from the wreck of the Mary Rose, AD 1545.

A comparison of ancient DNA (single-nucleotide polymorphisms) and carbon and nitrogen stable isotope evidence suggests that stored cod provisions recovered from the wreck of the Tudor warship Mary Rose, which sank in the Solent, southern England, in 1545, had been caught in northern and transatlantic waters such as the northern North Sea and the fishing grounds of Iceland and Newfoundland. This discovery, underpinned by control data from archaeological samples of cod bones from potential source regions, illuminates the role of naval provisioning in the early development of extensive sea fisheries, with their long-term economic and ecological impacts.

Cannibalism in invasive, native and biocontrol populations of the harlequin ladybird.

BACKGROUND: Cannibalism is widespread in both vertebrates and invertebrates but its extent is variable between and within species. Cannibalism depends on population density and nutritional conditions, and could be beneficial during colonisation of new environments. Empirical studies are needed to determine whether this trait might facilitate invasion of a new area in natural systems. We investigated whether the propensity for cannibalism in H. axyridis differs both between native and invasive populations and between invasive populations from the core and from the front of the invasive area in Western Europe. We also compared the propensity for cannibalism of these natural populations with that of laboratory-reared biocontrol populations. We measured the cannibalism rates of eggs by first instar larvae and adult females at two different individual densities of ladybirds from three types of population (invasive, native and biocontrol), in laboratory-controlled conditions. RESULTS: Cannibalism was significantly greater in larvae from invasive populations compared to native or biocontrol populations, but there was no difference in cannibalism rates between populations from the core or front of the invaded range. Cannibalism was significantly lower in larvae from biocontrol populations compared to wild (invasive and native) populations. No differences in cannibalism rates of adult females were found between any populations. While high population density significantly increased cannibalism in both larvae and adults, the norm of reaction of cannibalism to individual density did not change significantly during the invasion and/or laboratory rearing processes. CONCLUSION: This study is the first to provide evidence for a higher propensity for cannibalism in invasive populations compared to native ones. Our experiments also shed light on the difference in cannibalism evolution with respect to life stages. However, we are still at an early stage in understanding the underlying mechanisms and several different research perspectives are needed to determine whether the higher propensity for cannibalism is a general feature of the invasion process.

Characteristics and drivers of high-altitude ladybird flight: insights from vertical-looking entomological radar.

Understanding the characteristics and drivers of dispersal is crucial for predicting population dynamics, particularly in range-shifting species. Studying long-distance dispersal in insects is challenging, but recent advances in entomological radar offer unique insights. We analysed 10 years of radar data collected at Rothamsted Research, U.K., to investigate characteristics (altitude, speed, seasonal and annual trends) and drivers (aphid abundance, air temperature, wind speed and rainfall) of high-altitude flight of the two most abundant U.K. ladybird species (native Coccinella septempunctata and invasive Harmonia axyridis). These species cannot be distinguished in the radar data since their reflectivity signals overlap, and they were therefore analysed together. However, their signals do not overlap with other, abundant insects so we are confident they constitute the overwhelming majority of the analysed data. The target species were detected up to ∼1100 m above ground level, where displacement speeds of up to ∼60 km/h were recorded, however most ladybirds were found between ∼150 and 500 m, and had a mean displacement of 30 km/h. Average flight time was estimated, using tethered flight experiments, to be 36.5 minutes, but flights of up to two hours were observed. Ladybirds are therefore potentially able to travel 18 km in a "typical" high-altitude flight, but up to 120 km if flying at higher altitudes, indicating a high capacity for long-distance dispersal. There were strong seasonal trends in ladybird abundance, with peaks corresponding to the highest temperatures of mid-summer, and warm air temperature was the key driver of ladybird flight. Climatic warming may therefore increase the potential for long-distance dispersal in these species. Low aphid abundance was a second significant factor, highlighting the important role of aphid population dynamics in ladybird dispersal. This research illustrates the utility of radar for studying high-altitude insect flight and has important implications for predicting long-distance dispersal.

Climate shaped the worldwide distribution of human mitochondrial DNA sequence variation.

There is an ongoing discussion in the literature on whether human mitochondrial DNA (mtDNA) evolves neutrally. There have been previous claims for natural selection on human mtDNA based on an excess of non-synonymous mutations and higher evolutionary persistence of specific mitochondrial mutations in Arctic populations. However, these findings were not supported by the reanalysis of larger datasets. Using a geographical framework, we perform the first direct test of the relative extent to which climate and past demography have shaped the current spatial distribution of mtDNA sequences worldwide. We show that populations living in colder environments have lower mitochondrial diversity and that the genetic differentiation between pairs of populations correlates with difference in temperature. These associations were unique to mtDNA; we could not find a similar pattern in any other genetic marker. We were able to identify two correlated non-synonymous point mutations in the ND3 and ATP6 genes characterized by a clear association with temperature, which appear to be plausible targets of natural selection producing the association with climate. The same mutations have been previously shown to be associated with variation in mitochondrial pH and calcium dynamics. Our results indicate that natural selection mediated by climate has contributed to shape the current distribution of mtDNA sequences in humans.