Nuclear Genotypes Research Articles

Mitochondrial-nuclear coadaptation revealed through mtDNA replacements in Saccharomyces cerevisiae

BackgroundMitochondrial function requires numerous genetic interactions between mitochondrial- and nuclear- encoded genes. While selection for optimal mitonuclear interactions should result in coevolution between both genomes, evidence for mitonuclear coadaptation is challenging to document. Genetic models where mitonuclear interactions can be explored are needed.ResultsWe systematically exchanged mtDNAs between 15 Saccharomyces cerevisiae isolates from a variety of ecological niches to create 225 unique mitochondrial-nuclear genotypes. Analysis of phenotypic profiles confirmed that environmentally-sensitive interactions between mitochondrial and nuclear genotype contributed to growth differences. Exchanges of mtDNAs between strains of the same or different clades were just as likely to demonstrate mitonuclear epistasis although epistatic effect sizes increased with genetic distances. Strains with their original mtDNAs were more fit than strains with synthetic mitonuclear combinations when grown in media that resembled isolation habitats.ConclusionsThis study shows that natural variation in mitonuclear interactions contributes to fitness landscapes. Multiple examples of coadapted mitochondrial-nuclear genotypes suggest that selection for mitonuclear interactions may play a role in helping yeasts adapt to novel environments and promote coevolution.

Nuclear Dynamics in the Arbuscular Mycorrhizal Fungi.

Arbuscular mycorrhizal fungi (AMF) are plant root symbionts that continuously carry thousands of nuclei in their spores and hyphae. This unique cellular biology raises fundamental questions regarding their nuclear dynamics. This review aims to address these by synthesizing current knowledge of nuclear content and behavior in these ubiquitous soil fungi. Overall, we find that that nuclear counts, as well as the nuclei shape and organization, vary drastically both within and among species in this group. By comparing these features with those of other fungi, we highlight unique aspects of the AMF nuclear biology that require further attention. The potential implications of the observed nuclear variability for the biology and evolution of these widespread plant symbionts are discussed.

Patterns of mtDNA introgression suggest population replacement in Palaearctic whiskered bat species.

Secondary contacts can play a major role in the evolutionary histories of species. Various taxa diverge in allopatry and later on come into secondary contact during range expansions. When they meet, their interactions and the extent of gene flow depend on the level of their ecological differentiation and the strength of their reproductive isolation. In this study, we present the multilocus phylogeography of two cryptic whiskered bat species, Myotis mystacinus and M. davidii, with a particular focus on their putative sympatric zone. Our findings suggest that M. mystacinus and M. davidii evolved in allopatry and came into secondary contact during range expansions. Individuals in the area of secondary contact, in Anatolia and the Balkans, have discordant population assignments based on the mitochondrial and the nuclear datasets. These observed patterns suggest that the local M. mystacinus populations hybridized with expanding M. davidii populations, which resulted in mitochondrial introgression from the former. In the introgression area, M. mystacinus individuals with concordant nuclear and mitochondrial genotypes were identified in relatively few locations, suggesting that the indigenous populations might have been largely replaced by invading M. davidii. Changing environmental conditions coupled with ecological competition is the likely reason for this replacement. Our study presents one possible example of a historical population replacement that was captured in phylogeographic patterns.

Phenological shifts and genetic differentiation between sympatric populations of Sargassum horneri (Fucales, Phaeophyceae) in Japan

Marked seasonality, especially in sexual reproduction, is common among seaweed species along temperate coasts and increases the possibility of successful fertilization in outcrossing species. A phenological shift in reproductive seasons, therefore, could be an effective isolation barrier between conspecific seasonal populations, although its power has not been verified in algae.Sargassum horneri, a major component of seaweed beds along the temperate coast of Japan, is known for variability in its reproductive phenology. To understand the significance of phenological shift as an isolation barrier in seaweed species, phenological investigations ofS. horneriseasonal populations on the Sea of Japan coast of central Honshu, Japan, were combined with Bayesian cluster analysis based on a nuclear simple sequence repeat genotype. Results from these analyses concordantly suggest a genetic differentiation between the seasonal populations, although almost 20% of field-collected plants were estimated to be hybrids or have a hybrid origin based on results of Bayesian cluster analyses using experimental hybrids. A collapse of seasonal isolation was also detected at the site of the field investigation, and a high percentage of putative hybrids in the following generation at the site (41%) suggested significant seasonal isolation in the differentiation observed in this study.

Genetic Diversity and Connectivity of Southern Right Whales (Eubalaena australis) Found in the Brazil and Chile–Peru Wintering Grounds and the South Georgia (Islas Georgias del Sur) Feeding Ground

As species recover from exploitation, continued assessments of connectivity and population structure are warranted to provide information for conservation and management. This is particularly true in species with high dispersal capacity, such as migratory whales, where patterns of connectivity could change rapidly. Here we build on a previous long-term, large-scale collaboration on southern right whales (Eubalaena australis) to combine new (nnew) and published (npub) mitochondrial (mtDNA) and microsatellite genetic data from all major wintering grounds and, uniquely, the South Georgia (Islas Georgias del Sur: SG) feeding grounds. Specifically, we include data from Argentina (npub mtDNA/microsatellite = 208/46), Brazil (nnew mtDNA/microsatellite = 50/50), South Africa (nnew mtDNA/microsatellite = 66/77, npub mtDNA/microsatellite = 350/47), Chile–Peru (nnew mtDNA/microsatellite = 1/1), the Indo-Pacific (npub mtDNA/microsatellite = 769/126), and SG (npub mtDNA/microsatellite = 8/0, nnew mtDNA/microsatellite = 3/11) to investigate the position of previously unstudied habitats in the migratory network: Brazil, SG, and Chile–Peru. These new genetic data show connectivity between Brazil and Argentina, exemplified by weak genetic differentiation and the movement of 1 genetically identified individual between the South American grounds. The single sample from Chile–Peru had an mtDNA haplotype previously only observed in the Indo-Pacific and had a nuclear genotype that appeared admixed between the Indo-Pacific and South Atlantic, based on genetic clustering and assignment algorithms. The SG samples were clearly South Atlantic and were more similar to the South American than the South African wintering grounds. This study highlights how international collaborations are critical to provide context for emerging or recovering regions, like the SG feeding ground, as well as those that remain critically endangered, such as Chile–Peru.

Mitonuclear mismatch alters performance and reproductive success in naturally introgressed populations of a montane leaf beetle.

Coordination between nuclear and mitochondrial genomes is critical to metabolic processes underlying animals' ability to adapt to local environments, yet consequences of mitonuclear interactions have rarely been investigated in populations where individuals with divergent mitochondrial and nuclear genomes naturally interbreed. Genetic variation in the leaf beetle Chrysomela aeneicollis was assessed along a latitudinal thermal gradient in California's Sierra Nevada. Variation at mitochondrial cytochrome oxidase II (COII) and the nuclear gene phosphoglucose isomerase (PGI) shows concordance and was significantly greater along a 65 km transect than 10 other loci. STRUCTURE analyses using neutral loci identified a southern and northern subpopulation, which interbreed in the central drainage Bishop Creek. COII and PGI were used as indicators of mitochondrial and nuclear genetic variation in field and laboratory experiments conducted on beetles from this admixed population. Fecundity, larval development rate, running speed and male mating frequency were higher for beetles with geographically "matched" than "mismatched" mitonuclear genotypes. Effects of mitonuclear mismatch were largest for individuals with northern nuclear genotypes possessing southern mitochondria and were most pronounced after heat treatment or at high elevation. These findings suggest that mitonuclear incompatibility diminishes performance and reproductive success in nature, effects that could intensify at environmental extremes.

Species-specific genetic variation in response to deep-sea environmental variation amongst Vulnerable Marine Ecosystem indicator taxa

Understanding the ecological processes that shape spatial genetic patterns of population structure is critical for understanding evolutionary dynamics and defining significant evolutionary and management units in the deep sea. Here, the role of environmental factors (topographic, physico-chemical and biological) in shaping the population genetic structure of four deep-sea habitat-forming species (one sponge - Poecillastra laminaris, three corals - Goniocorella dumosa, Madrepora oculata, Solenosmilia variabilis) was investigated using seascape genetics. Genetic data (nuclear and mitochondrial sequences and microsatellite multilocus genotypes) and environmental variables were employed to build individual-based and population-level models. The results indicated that environmental factors affected genetic variation differently amongst the species, as well as at different geographic scales. For individual-based analyses, different environmental variables explained genetic variation in P. laminaris (dissolved oxygen), G. dumosa (dynamic topography), M. oculata (sea surface temperature and surface water primary productivity), and S. variabilis (tidal current speed). At the population level, factors related to current and food source explained the regional genetic structure in all four species, whilst at the geomorphic features level, factors related to food source and topography were most important. Environmental variation in these parameters may be acting as barriers to gene flow at different scales. This study highlights the utility of seascape genetic studies to better understand the processes shaping the genetic structure of organisms, and to identify environmental factors that can be used to locate sites for the protection of deep-sea Vulnerable Marine Ecosystems.

Mitochondrial genome-wide association study of migraine - the HUNT Study.

BackgroundVariation in mitochondrial DNA (mtDNA) has been indicated in migraine pathogenesis, but genetic studies to date have focused on candidate variants, with sparse findings. We aimed to perform the first mitochondrial genome-wide association study of migraine, examining both single variants and mitochondrial haplogroups.MethodsIn total, 71,860 participants from the population-based Nord-Trøndelag Health Study were genotyped. We excluded samples not passing quality control for nuclear genotypes, in addition to samples with low call rate and closely maternally related. We analysed 775 mitochondrial DNA variants in 4021 migraine cases and 14,288 headache-free controls, using logistic regression. In addition, we analysed 3831 cases and 13,584 controls who could be reliably assigned to a mitochondrial haplogroup. Lastly, we attempted to replicate previously reported mitochondrial DNA candidate variants.ResultsNeither of the mitochondrial variants or haplogroups were associated with migraine. In addition, none of the previously reported mtDNA candidate variants replicated in our data.ConclusionsOur findings do not support a major role of mitochondrial genetic variation in migraine pathophysiology, but a larger sample is needed to detect rare variants and future studies should also examine heteroplasmic variation, epigenetic changes and copy-number variation.

The cuticle inward barrier in Drosophila melanogaster is shaped by mitochondrial and nuclear genotypes and a sex-specific effect of diet.

An important role of the insect cuticle is to prevent wetting (i.e., permeation of water) and also to prevent penetration of potentially harmful substances. This barrier function mainly depends on the hydrophobic cuticle surface composed of lipids including cuticular hydrocarbons (CHCs). We investigated to what extent the cuticle inward barrier function depends on the genotype, comprising mitochondrial and nuclear genes in the fruit fly Drosophila melanogaster, and investigated the contribution of interactions between mitochondrial and nuclear genotypes (mito-nuclear interactions) on this function. In addition, we assessed the effects of nutrition and sex on the cuticle barrier function. Based on a dye penetration assay, we find that cuticle barrier function varies across three fly lines that were captured from geographically separated regions in three continents. Testing different combinations of mito-nuclear genotypes, we show that the inward barrier efficiency is modulated by the nuclear and mitochondrial genomes independently. We also find an interaction between diet and sex. Our findings provide new insights into the regulation of cuticle inward barrier function in nature.

Multiscale spatial genetic structure within and between populations of wild cherry trees in nuclear genotypes and chloroplast haplotypes.

Spatial genetic structure (SGS) of plants mainly depends on the effective population size and gene dispersal. Maternally inherited loci are expected to have higher genetic differentiation between populations and more intensive SGS within populations than biparentally inherited loci because of smaller effective population sizes and fewer opportunities of gene dispersal in the maternally inherited loci. We investigated biparentally inherited nuclear genotypes and maternally inherited chloroplast haplotypes of microsatellites in 17 tree populations of three wild cherry species under different conditions of tree distribution and seed dispersal. As expected, interpopulation genetic differentiation was 6–9 times higher in chloroplast haplotypes than in nuclear genotypes. This difference indicated that pollen flow 4–7 times exceeded seed flow between populations. However, no difference between nuclear and chloroplast loci was detected in within‐population SGS intensity due to their substantial variation among the populations. The SGS intensity tended to increase as trees became more aggregated, suggesting that tree aggregation biased pollen and seed dispersal distances toward shorter. The loss of effective seed dispersers, Asian black bears, did not affect the SGS intensity probably because of mitigation of the bear loss by other vertebrate dispersers and too few tree generations after the bear loss to alter SGS. The findings suggest that SGS is more variable in smaller spatial scales due to various ecological factors in local populations.

Long-term isolation at a low effective population size greatly reduced genetic diversity in Gulf of California fin whales

The Gulf of California, Mexico is home to many cetacean species, including a presumed resident population of fin whales, Balaenoptera physalus. Past studies reported very low levels of genetic diversity among Gulf of California fin whales and a significant level of genetic differentiation from con-specifics in the eastern North Pacific. The aim of the present study was to assess the degree and timing of the isolation of Gulf of California fin whales in a population genetic analysis of 18 nuclear microsatellite genotypes from 402 samples and 565 mitochondrial control region DNA sequences (including mitochondrial sequences retrieved from NCBI). The analyses revealed that the Gulf of California fin whale population was founded ~2.3 thousand years ago and has since remained at a low effective population size (~360) and isolated from the eastern North Pacific (Nem between 0.89–1.4). The low effective population size and high degree of isolation implied that Gulf of California fin whales are vulnerable to the negative effects of genetic drift, human-caused mortality and habitat change.

Nuclear and mitochondrial marker sequences reveal close relationship between Coronocyclus coronatus and a potential Cylicostephanus calicatus cryptic species complex

The Cyathostominae (Nematoda, Strongyloidea) parasitising equines represent a diverse group currently including 50 species. However, their taxonomy has been repeatedly revised and occasionally the presence of cryptic genospecies was suggested. Moreover, molecular- and morphology-based phylogenetic analyses give divergent results. For instance, molecular data have suggested close relationship between Coronocyclus coronatus and Cylicostephanus calicatus, although morphology-based taxonomy places them in different genera. Here, nuclear (internal transcribed spacer 2, ITS-2) and mitochondrial (cytochrome oxidase I, COI) sequences were obtained from the same individual, morphologically identified worms. In both morphospecies, two ITS-2 sequences types were observed: In Cor. coronatus, a small PCR product of 278 bp (nuclear haplotype group nHGBco) was always present but often in combination with a larger 369–370 bp fragment (nHGAco). In Cys. calicatus, either a large 370 bp product (nHGAca) or a short 281 bp amplicon (nHGBca) were found, but never both. Sequence identity between morphospecies was up to 100%. The smaller differed from the larger fragments by deletion of the region 110–198 in Cor. coronatus and 112–203 in Cys. calicatus. In COI, three and five mitochondrial haplotype groups (HGs), mtHG1co-mtHG3co and mtHG1ca-mtHG5ca were identified for Cor. coronatus and Cys. calicatus, respectively. In Cor. coronatus, there was no particular association of mtHG with nuclear genotypes (only nHGBco vs. both nHGBco plus nHGAco). In Cys. calicatus the nHGAca was always associated with the mtHG1ca, mtHG2ca or mtHG5ca whereas nHGBca was exclusively associated with mtHG3ca or mtHG4ca. Despite up to 100% identity in the nHGs, no mixing of mtHGs was observed between both species. Clear separation of certain nHGs with particular mtHGs in Cys. calicatus, despite the fact that the same host individuals were infected with both groups simultaneously, suggests presence of two non-interbreeding genospecies within Cys. calicatus, which needs further confirmation using additional samples from diverse geographical origins.

Lethal Interaction of Nuclear and Mitochondrial Genotypes in Drosophila melanogaster.

Drosophila melanogaster, like most animal species, displays considerable genetic variation in both nuclear and mitochondrial DNA (mtDNA). Here we tested whether any of four natural mtDNA variants was able to modify the effect of the phenotypically mild, nuclear tko25t mutation, affecting mitochondrial protein synthesis. When combined with tko25t, the mtDNA from wild strain KSA2 produced pupal lethality, accompanied by the presence of melanotic nodules in L3 larvae. KSA2 mtDNA, which carries a substitution at a conserved residue of cytochrome b that is predicted to be involved in subunit interactions within respiratory complex III, conferred drastically decreased respiratory capacity and complex III activity in the tko25t but not a wild-type nuclear background. The complex III inhibitor antimycin A was able to phenocopy effects of the tko25t mutation in the KSA2 mtDNA background. This is the first report of a lethal, nuclear-mitochondrial interaction within a metazoan species, representing a paradigm for understanding genetic interactions between nuclear and mitochondrial genotype relevant to human health and disease.

Our Unique Cell Lines Derived Serially After Injection With Human T24 HRAS Transformed NIH/3T3 Mouse Cells (GhrasT‐NIH/3T3) Showed The Highly Tumorigenic T1‐A Cells Cultured From The First Primary Tumor Had A Six‐Fold Higher LD50 For Resveratrol Killing Than The LD50 For Highly Metastatic Lung T4‐PA Cells Derived After Three Subsequent In Vivo / In Vitro Serial Manipulations

Our new cell lines demonstrating progressively increasingly metastatic behavior as they were derived in a series of in vivo/in vitro manipulations (Figure 1) provide the opportunity to evaluate in vivo and in vitro compounds and methods to kill primary tumor cells and various forms of increasingly metastatic cells derived from the same primary tumor (Figure 2) (Ray, D.B., et. al., Experimental Cell Research, 2016 Jan; 340 (1): pp. 1–11). We determined the dose of the anti‐cancer natural phytoalexin resveratrol to kill 50% (LD50s) of our various cells that represent progressive cancer stages. Cell counts were recorded daily, over the course of 3 to 6 days, using gridded plates and photomicrographs of six pre‐determined 2 mm square sections of each plate. An average percent change in cell number from each individual grid was then calculated. P‐values were determined between each treatment time and dose and significance was determined. The NIH Swiss and NIH/3T3 cells were not significantly affected by 3–300 μM after 3 days. After 3 or 4 days the LD50s for the tumorigenic cells were GhrasT‐NIH/3T3 cells: LD50 = 3.5+/− 1.0 μM; T1‐A cells: LD50 = 45 +/− 5.0 μM; T4‐PA cells: LD50 = 7.0 +/− 1.5 μM. Thus, unexpectantly, the highly tumorigenic T1‐A cells were significantly more difficult to kill than the highly metastatic T4‐PA cells. This pattern may mimic the epithelial to mesenchymal transition (EMT) seen in breast and other cancer progression. EMT enhances mobility, invasion, and resistance to apoptotic stimuli as tumor cells acquire stem cell properties and exhibit marked therapeutic resistance. RNA‐Seq, genome sequencing (WGS) and spectral karyotyping (SKY) analysis of these cells demonstrating increasingly aggressive metastatic cancer tumorigenesis that progressed from a single nuclear genotype and mitochondrial haplotype should help contribute to the understanding of cancer progression. These cells are available by contacting our laboratory.Support or Funding InformationSupported in part by the National Science Foundation grant # 8602787, the Grove City College Research Fund, the Grove City College Swezey Research Fund and the Grove City College Janicki Foundation Research Fund.This abstract is from the Experimental Biology 2019 Meeting. There is no full text article associated with this abstract published in The FASEB Journal.

An admixture of Quercus dentata in the coastal ecotype of Q. mongolica var. crispula in northern Hokkaido and genetic and environmental effects on their traits.

In northern Japan, coastal oak forests consist of Quercus dentata (Qd) on the coastal side and Q. mongolica var. crispula (Qc) on the inland side. In the forests of northern Hokkaido, Qd is rare, and a coastal ecotype of Qc with some Qd-like traits grows on the coastal side. To reveal the genetic background of this ecotype, nuclear microsatellite genotypes in closely related oak taxa were obtained from the Eurasian continent, Sakhalin, and Hokkaido. The clustering of these genotypes suggests an admixture of Qd in the coastal ecotype of Qc. Next, we evaluated the effects of admixture and coastal stress on the leaf and shoot traits of Qc and Qd along coastal-inland gradients in northern Hokkaido. The admixture of Qd in Qc was quantified by the Qd ancestry proportions. Coastal stress causes bud mortality in the upper parts of shoots and was quantified by the survival patterns of buds in shoots. The genetic and environmental effects on the traits at Qd-abundant and Qd-rare sites were estimated using linear mixed models. The genetic effect was detected in all traits. Both genetic and environmental effects were detected in most traits. Some traits differed between Qd-abundant and Qd-rare sites in addition to these effects, indicating more Qd-like traits at Qd-rare sites. The findings suggest that an admixture of Qd characterizes the genetic background of the coastal ecotype of Qc and that not only the coastal stress but also the genetic background is responsible for the leaf and shoot traits of Qc and Qd in northern Hokkaido.

Genetic variability of the mtDNA control region and the intron of the nuclear gene tumor necrosis factor alpha in graylings, Thymallus

Genetic variability of the mtDNA control region and the intron of the nuclear gene tumor necrosis factor alpha in graylings, Thymallus

Genetic identification of native populations of Mediterranean brown trout Salmo trutta L. complex (Osteichthyes: Salmonidae) in central Italy

Italian native populations of Mediterranean brown trout belong to the Salmo trutta complex. This species complex includes many mitochondrial lineages and phenotypes that have caused taxonomic controversies over time. The spatial distribution and the genetic diversity of these fishes are threatened by habitat destruction, global warming and, mainly, by the introduction of domestic trout of Atlantic origin. Indeed allochthonous trouts were massively restocked in Italian rivers for a century and they admixed with native populations. In order to identify residual native populations of Mediterranean brown trout, a genetic analysis of specimens collected within Latium region, on the Tyrrhenian slope of central Italy, was undertaken. To this purpose, 210 trout specimens were collected from six different rivers and analyzed for the identification of their nuclear (LDH-C1* RFLP) and mitochondrial (Control Region sequences) genotypes. Genetic characterization with these molecular markers allowed a quantitative estimate of allochthonous genotypes, which are present in all brown trout populations of the six sites, even if not equally distributed across the sampling area. At least three populations, inhabiting diverse lotic ecosystems (mountain, hilly and coastal streams respectively), are characterized by a high percentage of native nuclear allele *100 at locus LDH-C1* and typical Mediterranean haplotypes (of AD and ME lineages), which can be considered as different management units (MUs). This finding highlighted the aquatic ecosystems of the Latium Region as an important hotspot of salmonid biodiversity within the Italian peninsula, with important implications from a conservation perspective.

Assessing the fitness consequences of mitonuclear interactions in natural populations.

Metazoans exist only with a continuous and rich supply of chemical energy from oxidative phosphorylation in mitochondria. The oxidative phosphorylation machinery that mediates energy conservation is encoded by both mitochondrial and nuclear genes, and hence the products of these two genomes must interact closely to achieve coordinated function of core respiratory processes. It follows that selection for efficient respiration will lead to selection for compatible combinations of mitochondrial and nuclear genotypes, and this should facilitate coadaptation between mitochondrial and nuclear genomes (mitonuclear coadaptation). Herein, we outline the modes by which mitochondrial and nuclear genomes may coevolve within natural populations, and we discuss the implications of mitonuclear coadaptation for diverse fields of study in the biological sciences. We identify five themes in the study of mitonuclear interactions that provide a roadmap for both ecological and biomedical studies seeking to measure the contribution of intergenomic coadaptation to the evolution of natural populations. We also explore the wider implications of the fitness consequences of mitonuclear interactions, focusing on central debates within the fields of ecology and biomedicine.

Deep Intraspecific Divergence in the Endemic Herb Lancea tibetica (Mazaceae) Distributed Over the Qinghai-Tibetan Plateau.

Qinghai-Tibetan Plateau (QTP) is an important biodiversity hub, which is very sensitive to climate change. Here in this study, we investigated genetic diversity and past population dynamics of Lancea tibetica (Mazaceae), an endemic herb to QTP and adjacent highlands. We sequenced chloroplast and nuclear ribosomal DNA fragments for 429 individuals, collected from 29 localities, covering their major distribution range at the QTP. A total of 19 chloroplast haplotypes and 13 nuclear genotypes in two well-differentiated lineages, corresponding to populations into two groups isolated by Tanggula and Bayangela Mountains. Meanwhile, significant phylogeographical structure was detected among sampling range of L. tibetica, and 61.50% of genetic variations was partitioned between groups. Gene flow across the whole region appears to be restricted by high mountains, suggesting a significant role of geography in the genetic differences between the two groups. Divergence time between the two lineages dated to 8.63 million years ago, which corresponded to the uplifting of QTP during the late Miocene and Pliocene. Ecological differences were found between both the lineages represent species-specific characteristics, sufficient to keep the lineages separated to a high degree. The simulated distribution from the last interglacial period to the current period showed that the distribution of L. tibetica experienced shrinkage and expansion. Climate changes during the Pleistocene glacial-interglacial cycles had a dramatic effect on L. tibetica distribution ranges. Multiple refugia of L. tibetica might have remained during the species history, to south of the Tanggula and north of Bayangela Mountains, both appeared as topological barrier and contributed to restricting gene flow between the two lineages. Together, geographic isolation and climatic factors have played a fundamental role in promoting diversification and evolution of L. tibetica.

Red Fox Ancestry and Connectivity Assessments Reveal Minimal Fur Farm Introgression in Greater Yellowstone Ecosystem

AbstractRocky Mountain red foxes Vulpes vulpes macroura potentially encounter other red fox Vulpes vulpes lineages at lower elevations, which may include nonindigenous red foxes derived from fur farms. Introgression from nonindigenous red foxes could have negative evolutionary consequences for the rare Rocky Mountain red fox subspecies. Red foxes at high elevations in the Greater Yellowstone Ecosystem exhibit lighter coat colors than those at lower elevations, potentially indicating that they represent the indigenous subspecies and that gene flow across the elevational gradient is restricted. We collected tissue samples across a 1,750-m elevation range and examined mitochondrial DNA sequences and nuclear DNA microsatellite genotypes to assess the ancestry and genetic population structure of red foxes in the northern Greater Yellowstone Ecosystem. We also used reference samples from fur farm red foxes and indigenous red foxes of the western United States to assess the extent of nonindigenous introgression across the ecosystem. We found little overlap in the elevational distribution of maternally inherited mitochondrial DNA haplotypes: above 1,600 m, we only found indigenous Rocky Mountain haplotypes (n = 4), whereas below 1,600 m, we found haplotypes not indigenous to the Rocky Mountains (n = 5) that were associated with fur farms or indigenous to the Great Plains. In contrast, biparentally inherited microsatellite variation showed little population structure across the elevational gradient. Despite this evidence of nuclear gene flow across the elevational gradient, we found little fur farm introgression in the microsatellite genotypes. It is possible that long-standing nuclear (but apparently not mitochondrial) gene flow between Rocky Mountain red foxes and indigenous red foxes on the Great Plains explained the low nuclear differentiation of these populations. Importantly, our results suggested that high elevations of the northern Greater Yellowstone Ecosystem remained free of significant fur farm introgression. Mitonuclear discordance could reflect sex-biased dispersal, which we hypothesize could be the effect of elevational differences in reproductive phenology.