Genomic Incompatibilities Research Articles

Asymmetric mate preference and reproductive interference mediate climate-induced changes in mate availability in a small mammal hybrid zone.

Range expansion and contraction are among the most common biotic responses to changing environmental conditions, yet much is to be learned about the mechanisms that underlie range-edge population dynamics, especially when those areas are points of secondary contact between closely related species. Here, we present field-measured parentage data that document the reproductive outcomes of changes in mate availability at a secondary contact zone between two species of woodrat in the genus Neotoma. Changes in mate availability resulted from drought-driven differential survival between the species and their hybrids. As the availability of conspecific mates declined, rates of hybridization increased, leading to the accumulation of admixed individuals in the zone of contact. Patterns of reproductive success in the wild appear to be the result of a combination of both pre-mating isolation and post-zygotic selection resulting from genomic incompatibilities between the parental lineages. Evidence of asymmetric mate preference between the parental lineages came from both skewed reproductive output in the field and laboratory preference trials. Moreover, partial genomic incompatibility was evident from the near-zero reproductive success of F1 males and because nearly all surviving hybrids had one pure parent. Nonetheless, the high reproductive success of F1 females and backcrossing in both parental directions allow for introgression between the parental species. These findings reveal how climate change may alter evolutionary outcomes for species at the edge of their ranges through an interplay of behavioral, demographic, and genetic mechanisms.

A Review on Role of Wide Hybridization in Crop Improvement

Wide hybridization in crop plants involves crossing two distantly related species to initiate desirable traits, such as disease resistance, into a cultivated crop. This technique expands the genetic diversity accessible for breeding programs, enhancing the adaptability and resilience of crops. Wide hybrids often surpass limitations of conventional breeding by incorporating unique genetic material. However, difficulties arise due to reproductive barriers and genomic incompatibilities between divergent species. Researchers employ techniques like embryo rescue and tissue culture to overcome these hurdles. Despite its potentiality, wide hybridization requires careful selection and extensive backcrossing to stabilize desired traits. The resulting crops may show improved resistance to pests, diseases, or environmental stressors, contributing to sustainable agriculture and food security. Continuous advancements in molecular biology and genomics facilitate the identification and transfer of specific genes, accelerating the development of wide hybrid varieties with enhanced agronomic traits.

Functional characterization of archaic-specific variants in mitonuclear genes: insights from comparative analysis in S. cerevisiae.

Neanderthal and Denisovan hybridisation with modern humans has generated a non-random genomic distribution of introgressed regions, the result of drift and selection dynamics. Cross-species genomic incompatibility and more efficient removal of slightly deleterious archaic variants have been proposed as selection-based processes involved in the post-hybridisation purge of archaic introgressed regions. Both scenarios require the presence of functionally different alleles across Homo species onto which selection operated differently according to which populations hosted them, but only a few of these variants have been pinpointed so far. In order to identify functionally divergent archaic variants removed in humans, we focused on mitonuclear genes, which are underrepresented in the genomic landscape of archaic humans. We searched for non-synonymous, fixed, archaic-derived variants present in mitonuclear genes, rare or absent in human populations. We then compared the functional impact of archaic and human variants in the model organism Saccharomyces cerevisiae. Notably, a variant within the mitochondrial tyrosyl-tRNA synthetase 2 (YARS2) gene exhibited a significant decrease in respiratory activity and a substantial reduction of Cox2 levels, a proxy for mitochondrial protein biosynthesis, coupled with the accumulation of the YARS2 protein precursor and a lower amount of mature enzyme. Our work suggests that this variant is associated with mitochondrial functionality impairment, thus contributing to the purging of archaic introgression in YARS2. While different molecular mechanisms may have impacted other mitonuclear genes, our approach can be extended to the functional screening of mitonuclear genetic variants present across species and populations.

Mapping restricted introgression across the genomes of admixed indigenous African cattle breeds

BackgroundThe genomes of indigenous African cattle are composed of components with Middle Eastern (taurine) and South Asian (indicine) origins, providing a valuable model to study hybridization and to identify genetic barriers to gene flow. In this study, we analysed indigenous African cattle breeds as models of hybrid zones, considering taurine and indicine samples as ancestors. In a genomic cline analysis of whole-genome sequence data, we considered over 8 million variants from 144 animals, which allows for fine-mapping of potential genomic incompatibilities at high resolution across the genome.ResultsWe identified several thousand variants that had significantly steep clines (‘SCV’) across the whole genome, indicating restricted introgression. Some of the SCV were clustered into extended regions, with the longest on chromosome 7, spanning 725 kb and including 27 genes. We found that variants with a high phenotypic impact (e.g. indels, intra-genic and missense variants) likely represent greater genetic barriers to gene flow. Furthermore, our findings provide evidence that a large proportion of breed differentiation in African cattle could be linked to genomic incompatibilities and reproductive isolation. Functional evaluation of genes with SCV suggest that mitonuclear incompatibilities and genes associated with fitness (e.g. resistance to paratuberculosis) could account for restricted gene flow in indigenous African cattle.ConclusionsTo our knowledge, this is the first time genomic cline analysis has been applied to identify restricted introgression in the genomes of indigenous African cattle and the results provide extended insights into mechanisms (e.g. genomic incompatibilities) contributing to hybrid differentiation. These results have important implications for our understanding of genetic incompatibilities and reproductive isolation and provide important insights into the impact of cross-breeding cattle with the aim of producing offspring that are both hardy and productive.

Genomic mechanisms and consequences of diverse postzygotic barriers between monkeyflower species.

The evolution of genomic incompatibilities causing postzygotic barriers to hybridization is a key step in species divergence. Incompatibilities take two general forms - structural divergence between chromosomes leading to severe hybrid sterility in F1 hybrids and epistatic interactions between genes causing reduced fitness of hybrid gametes or zygotes (Dobzhansky-Muller incompatibilities). Despite substantial recent progress in understanding the molecular mechanisms and evolutionary origins of both types of incompatibility, how each behaves across multiple generations of hybridization remains relatively unexplored. Here, we use genetic mapping in F2 and RIL hybrid populations between the phenotypically divergent but naturally hybridizing monkeyflowers Mimulus cardinalis and M. parishii to characterize the genetic basis of hybrid incompatibility and examine its changing effects over multiple generations of experimental hybridization. In F2s, we found severe hybrid pollen inviability (< 50% reduction vs. parental genotypes) and pseudolinkage caused by a reciprocal translocation between Chromosomes 6 and 7 in the parental species. RILs retained excess heterozygosity around the translocation breakpoints, which caused substantial pollen inviability when interstitial crossovers had not created compatible heterokaryotypic configurations. Strong transmission ratio distortion and inter-chromosomal linkage disequilibrium in both F2s and RILs identified a novel two-locus genic incompatibility causing sex-independent gametophytic (haploid) lethality. The latter interaction eliminated three of the expected nine F2 genotypic classes via F1 gamete loss without detectable effects on the pollen number or viability of F2 double heterozygotes. Along with the mapping of numerous milder incompatibilities, these key findings illuminate the complex genetics of plant hybrid breakdown and are an important step toward understanding the genomic consequences of natural hybridization in this model system.

Evidence of hybridization between genetically distinct Baltic cod stocks during peak population abundance(s).

Range expansions can lead to increased contact of divergent populations, thus increasing the potential of hybridization events. Whether viable hybrids are produced will most likely depend on the level of genomic divergence and associated genomic incompatibilities between the different entities as well as environmental conditions. By taking advantage of historical Baltic cod (Gadus morhua) otolith samples combined with genotyping and whole genome sequencing, we here investigate the genetic impact of the increased spawning stock biomass of the eastern Baltic cod stock in the mid 1980s. The eastern Baltic cod is genetically highly differentiated from the adjacent western Baltic cod and locally adapted to the brackish environmental conditions in the deeper Eastern basins of the Baltic Sea unsuitable for its marine counterparts. Our genotyping results show an increased proportion of eastern Baltic cod in western Baltic areas (Mecklenburg Bay and Arkona Basin)-indicative of a range expansion westwards-during the peak population abundance in the 1980s. Additionally, we detect high frequencies of potential hybrids (including F1, F2 and backcrosses), verified by whole genome sequencing data for a subset of individuals. Analysis of mitochondrial genomes further indicates directional gene flow from eastern Baltic cod males to western Baltic cod females. Our findings unravel that increased overlap in distribution can promote hybridization between highly divergent populations and that the hybrids can be viable and survive under specific and favourable environmental conditions. However, the observed hybridization had seemingly no long-lasting impact on the continuous separation and genetic differentiation between the unique Baltic cod stocks.

Asymmetric gene expression in grain development of reciprocal crosses between tetraploid and hexaploid wheats

Production of viable progeny from interploid crosses requires precise regulation of gene expression from maternal and paternal chromosomes, yet the transcripts contributed to hybrid seeds from polyploid parent species have rarely been explored. To investigate the genome-wide maternal and paternal contributions to polyploid grain development, we analyzed the transcriptomes of developing embryos, from zygote to maturity, alongside endosperm in two stages of development, using reciprocal crosses between tetraploid and hexaploid wheats. Reciprocal crosses between species with varied levels of ploidy displayed broad impacts on gene expression, including shifts in alternative splicing events in select crosses, as illustrated by active splicing events, enhanced protein synthesis and chromatin remodeling. Homoeologous gene expression was repressed on the univalent D genome in pentaploids, but this suppression was attenuated in crosses with a higher ploidy maternal parent. Imprinted genes were identified in endosperm and early embryo tissues, supporting predominant maternal effects on early embryogenesis. By systematically investigating the complex transcriptional networks in reciprocal-cross hybrids, this study presents a framework for understanding the genomic incompatibility and transcriptome shock that results from interspecific hybridization and uncovers the transcriptional impacts on hybrid seeds created from agriculturally-relevant polyploid species.

Analysis of the Whole-Genome Sequences from an Equus Parent-Offspring Trio Provides Insight into the Genomic Incompatibilities in the Hybrid Mule.

Interspecific hybridization often shows negative effects on hybrids. However, only a few multicellular species, limited to a handful of plants and animals, have shown partial genetic mechanisms by which hybridization leads to low fitness in hybrids. Here, to explore the outcome of combining the two genomes of a horse and donkey, we analyzed the whole-genome sequences from an Equus parent-offspring trio using Illumina platforms. We generated 41.39× and 46.21× coverage sequences for the horse and mule, respectively. For the donkey, a 40.38× coverage sequence was generated and stored in our laboratory. Approximately 24.86 million alleles were discovered that varied from the reference genome. Single nucleotide polymorphisms were used as polymorphic markers for assigning alleles to their parental genomic inheritance. We identified 25,703 Mendelian inheritance error single nucleotide polymorphisms in the mule genome that were not inherited from the parents through Mendelian inheritance. A total of 555 de novo single nucleotide polymorphisms were also identified. The rate of de novo single nucleotide polymorphisms was 2.21 × 10-7 in the mule from the Equus parent-offspring trio. This rate is obviously higher than the natural mutation rate for Equus, which is also consistent with the previous hypothesis that interracial crosses may have a high mutation rate. The genes associated with these single nucleotide polymorphisms are mainly involved in immune processes, DNA repair, and cancer processes. The results of the analysis of three genomes from an Equus parent-offspring trio improved our knowledge of the consequences of the integration of parental genomes in mules.

Contrasting continental patterns of adaptive population divergence in the holarctic ectomycorrhizal fungus Boletus edulis.

In the hyperdiverse fungi, the process of speciation is virtually unknown, including for the > 20 000 species of ectomycorrhizal mutualists. To understand this process, we investigated patterns of genome-wide differentiation in the ectomycorrhizal porcini mushroom, Boletus edulis, a globally distributed species complex with broad ecological amplitude. By whole-genome sequencing 160 individuals from across the Northern Hemisphere, we genotyped 792 923 single nucleotide polymorphismsto characterize patterns of genome-wide differentiation and to identify the adaptive processes shaping global population structure. We show that B. edulis exhibits contrasting patterns of genomic divergence between continents, with multiple lineages present across North America, while a single lineage dominates Europe. These geographical lineages are inferred to have diverged 1.62-2.66 million years ago, during a period of climatic upheaval and the onset of glaciation in the Pliocene-Pleistocene boundary. High levels of genomic differentiation were observed among lineages despite evidence of substantial and ongoing introgression. Genome scans, demographic inference, and ecological niche models suggest that genomic differentiation is maintained by environmental adaptation, not physical isolation. Our study uncovers striking patterns of genome-wide differentiation on a global scale and emphasizes the importance of local adaptation and ecologically mediated divergence, rather than prezygotic barriers such as allopatry or genomic incompatibility, in fungal population differentiation.

Evidence that genomic incompatibilities and other multilocus processes impact hybrid fitness in a rattlesnake hybrid zone.

Hybrid zones provide valuable opportunities to understand the genomic mechanisms that promote speciation by providing insight into factors involved in intermediate stages of speciation. Here, we investigate introgression in a hybrid zone between two rattlesnake species (Crotalus viridis and Crotalus oreganus concolor) that have undergone historical allopatric divergence and recent range expansion and secondary contact. We use Bayesian genomic cline models to characterize genomic patterns of introgression between these lineages and identify loci potentially subject to selection in hybrids. We find evidence for a large number of genomic regions with biased ancestry that deviate from the genomic background in hybrids (i.e., excess ancestry loci), which tend to be associated with genomic regions with higher recombination rates. We also identify suites of excess ancestry loci that show highly correlated allele frequencies (including conspecific and heterospecific combinations) across physically unlinked genomic regions in hybrids. Our findings provide evidence for multiple multilocus evolutionary processes impacting hybrid fitness in this system.

Symmetric subgenomes and balanced homoeolog expression stabilize the establishment of allopolyploidy in cyprinid fish

BackgroundInterspecific postzygotic reproduction isolation results from large genetic divergence between the subgenomes of established hybrids. Polyploidization immediately after hybridization may reset patterns of homologous chromosome pairing and ameliorate deleterious genomic incompatibility between the subgenomes of distinct parental species in plants and animals. However, the observation that polyploidy is less common in vertebrates raises the question of which factors restrict its emergence. Here, we perform analyses of the genome, epigenome, and gene expression in the nascent allotetraploid lineage (2.95 Gb) derived from the intergeneric hybridization of female goldfish (Carassius auratus, 1.49 Gb) and male common carp (Cyprinus carpio, 1.42 Gb), to shed light on the changes leading to the stabilization of hybrids.ResultsWe firstly identify the two subgenomes derived from the parental lineages of goldfish and common carp. We find variable unequal homoeologous recombination in somatic and germ cells of the intergeneric F1 and allotetraploid (F22 and F24) populations, reflecting high plasticity between the subgenomes, and rapidly varying copy numbers between the homoeolog genes. We also find dynamic changes in transposable elements accompanied by genome merger and duplication in the allotetraploid lineage. Finally, we observe the gradual decreases in cis-regulatory effects and increases in trans-regulatory effects along with the allotetraploidization, which contribute to increases in the symmetrical homoeologous expression in different tissues and developmental stages, especially in early embryogenesis.ConclusionsOur results reveal a series of changes in transposable elements, unequal homoeologous recombination, cis- and trans-regulations (e.g. DNA methylation), and homoeologous expression, suggesting their potential roles in mediating adaptive stabilization of regulatory systems of the nascent allotetraploid lineage. The symmetrical subgenomes and homoeologous expression provide a novel way of balancing genetic incompatibilities, providing a new insight into the early stages of allopolyploidization in vertebrate evolution.

The molecular biology of grafting: Recent research may provide new applications for a millennia-old agricultural technology.

Plant biologists have begun to unravel the molecular mechanisms behind grafting, yielding insights into plant evolution and expanding the scope of this technique.

Genetic Architecture Underlying Nascent Speciation-The Evolution of Eurasian Pigs under Domestication.

Speciation is a process whereby the evolution of reproductive barriers leads to isolated species. Although many studies have addressed large-effect genetic footprints in the advanced stages of speciation, the genetics of reproductive isolation in nascent stage of speciation remains unclear. Here, we show that pig domestication offers an interesting model for studying the early stages of speciation in great details. Pig breeds have not evolved the large X-effect of hybrid incompatibility commonly observed between “good species.” Instead, deleterious epistatic interactions among multiple autosomal loci are common. These weak Dobzhansky–Muller incompatibilities confer partial hybrid inviability with sex biases in crosses between European and East Asian domestic pigs. The genomic incompatibility is enriched in pathways for angiogenesis, androgen receptor signaling and immunity, with an observation of many highly differentiated cis-regulatory variants. Our study suggests that partial hybrid inviability caused by pervasive but weak interactions among autosomal loci may be a hallmark of nascent speciation in mammals.

Asymmetric Fitness of Second-Generation Interspecific Hybrids Between Ciona robusta and Ciona intestinalis.

Reproductive isolation is central to speciation, but interspecific crosses between two closely related species can produce viable and fertile hybrids. Two different species of tunicates in the same ascidian genus, Ciona robusta and Ciona intestinalis, can produce hybrids. However, wild sympatric populations display limited gene flow, suggesting the existence of obstacles to interspecific reproduction that remain unknown. Here, we took advantage of a closed culture system to cross C. robusta with C. intestinalis and established F1 and F2 hybrids. We monitored post-embryonic development, survival, and sexual maturation to characterize the genetic basis of simple traits, and further probe the physiological mechanisms underlying reproductive isolation. Partial viability of first and second generation hybrids suggested that both pre- and postzygotic mechanisms contributed to genomic incompatibilities in hybrids. We observed asymmetric fitness, whereby the C. intestinalis maternal lines fared more poorly in our system, pointing to maternal origins of species-specific sensitivity. We discuss the possibility that asymmetrical second generation inviability and infertility emerge from interspecific incompatibilities between the nuclear and mitochondrial genomes, or other maternal effect genes. This work paves the way to quantitative genetic approaches to study the mechanisms underlying genomic incompatibilities and other complex traits in the genome-enabled Ciona model.

Safeguarding mitochondrial genomes in higher eukaryotes.

Mitochondria respond to DNA damage and preserve their own genetic material in a manner distinct from that of the nucleus but that requires organized mito-nuclear communication. Failure to resolve mtDNA breaks leads to mitochondrial dysfunction and affects host cells and tissues. Here, we review the pathways that safeguard mitochondrial genomes and examine the insights gained from studies of cellular and tissue-wide responses to mtDNA damage and mito-nuclear genome incompatibility.

Effect of Whole-Genome Duplication on the Evolutionary Rescue of Sterile Hybrid Monkeyflowers

Hybridization is a creative evolutionary force, increasing genomic diversity and facilitating adaptation and even speciation. Hybrids often face significant challenges to establishment, including reduced fertility that arises from genomic incompatibilities between their parents. Whole-genome duplication in hybrids (allopolyploidy) can restore fertility, cause immediate phenotypic changes, and generate reproductive isolation. Yet the survival of polyploid lineages is uncertain, and few studies have compared the performance of recently formed allopolyploids and their parents under field conditions. Here, we use natural and synthetically produced hybrid and polyploid monkeyflowers (Mimulus spp.) to study how polyploidy contributes to the fertility, reproductive isolation, phenotype, and performance of hybrids in the field. We find that polyploidization restores fertility and that allopolyploids are reproductively isolated from their parents. The phenotype of allopolyploids displays the classic gigas effect of whole-genome duplication, in which plants have larger organs and are slower to flower. Field experiments indicate that survival of synthetic hybrids before and after polyploidization is intermediate between that of the parents, whereas natural hybrids have higher survival than all other taxa. We conclude that hybridization and polyploidy can act as sources of genomic novelty, but adaptive evolution is key in mediating the establishment of young allopolyploid lineages.

Fitness costs associated with ancestry to isolated populations of an endangered species

Habitat fragmentation from urban development leaves species vulnerable to inbreeding depression and genomic erosion. Restoring gene flow can reduce inbreeding and preserve genetic diversity, but a common concern is that genomic incompatibilities may lead to outbreeding depression. The introduction of deleterious genetic load is less commonly considered. The endangered Pacific pocket mouse (Perognathus longimembris pacificus) persists in three isolated populations in southern California. Mitochondrial and microsatellite data indicated that effective population sizes were extremely small (Ne< 50), and continued declines prompted a conservation breeding program founded by individuals from each population. We tracked genetic diversity and individual fitness (measured by reproductive success) in a captive setting over six generations of admixture. Although we observed an increase in fitness in the F1 and F2 generations relative to the founding populations, inbreeding depression alone did not explain the fitness patterns, and outbreeding depression was difficult to rule out as reproductive success waned after F2. However, reproductive success was consistently higher in admixed individuals than founders from Dana Point, the smallest population with the lowest heterozygosity. Across generations, we saw a strong negative correlation between individual reproductive success and ancestry to Dana Point, leading to a rapid decrease of Dana Point alleles. Although the genomic underpinnings remain to be determined, reduced fitness associated with Dana Point ancestry is consistent with high deleterious genetic load in this population, and thus any facilitated migration should be unidirectional. Our findings highlight that, even in the absence of outbreeding depression, there may be a cost to restoration of gene flow if populations harbor high genetic load.

High-resolution linkage and quantitative trait locus mapping using an interspecific cross between Argopecten irradians irradians (♀) and A. purpuratus (♂)

The bay scallop and Peruvian scallop are economically important species. Interspecific hybrids of these two scallops outperformed both of their parent species in multiple growth traits but exhibited decreased fertility, which provides good models for the study of heterosis and species divergence. Genetic mapping serves as a chromosomal-level framework to investigate the molecular mechanisms of hybridization and introgression. In this study, high-resolution linkage maps were constructed for the bay and Peruvian scallops with an interspecific hybrid family. The linkage map of the bay scallop covered over 98.9% of the whole genome with 2994 mapped markers and the average marker interval of 0.32 cM. For the Peruvian scallop, 1585 markers were mapped with the average maker interval of 0.51 cM, covering 97.7% of the genome. Both the two linkage maps have 16 linkage groups, corresponding to the haploid chromosome number of the two species. Approximately, 54.5% of markers exhibited significant deviation from the expected Mendelian ratio of segregation, lending insights into the intrinsic incompatibility between the two species. QTLs related to growth and shell coloration were detected, which could explain 13.1% and 74.9% of the phenotypic variance, respectively. This represents important information for further evaluation. These findings are an important addition to the genomic resources for scallop genetic studies, and are especially useful for investigations on genomic incompatibility for hybridization, genome evolution of closely related species, and genetic enhancement programs in aquaculture.

The Effects of Genetically Modified Organisms (GMO) on Environment and Molecular Techniques to Minimize Its Risk

Biotechnology is the application of scientific techniques to modify and improve plants, animals, and microorganisms to enhance their value. Genetically modified organism (GMO) is any organism whose genetic material has been altered using genetic engineering techniques. Transgenic technology deals with the integration of exogenous DNA into the plant genome using gene transfer technologies. Genetically modified (GM) crop plants contain artificially inserted gene (s) or “transgenes” from another unrelated plant or from a completely different species via advanced genetic engineering techniques. Transgenic crops do not present new categories of environmental risk compared to conventional methods of crop improvement. Transgenic plants can be generated using (1) biolistic technique where in the desired gene is coated on to either gold or tungsten particles are shot into plant cells using a gene-gun. The necessary criterion for this is that cells or plant tissues should be suitable for transformation permit regeneration of a whole plant thereafter. (2) The other major gene transfer technique makes use of the soil bacterium Agrobacterium tumefaciens, containing a tumor-inducing (Ti) plasmid including virulence (vir) gene (s) and a transferred-DNA (T-DNA) region, in which genes of interest can be inserted. Generally the gene that transfer from other organism to the desired crop (GM) can cause different environmental risks like, generation of super weed, development of tolerance to target herbicide, loss of biodiversity and sustainable resistance in insect pests, through gene flow. These risks can be minimized using different molecular techniques such as, maternal inheritance, male sterility, cleistogamy and apomixes, genetic use restriction technologies, genome incompatibility, controlling gene expression and transgenic mitigation.

Mitochondrial haplotypes affect metabolic phenotypes in the Drosophila Genetic Reference Panel.

The nature and extent of mitochondrial DNA variation in a population and how it affects traits is poorly understood. Here we resequence the mitochondrial genomes of 169 Drosophila Genetic Reference Panel lines, identifying 231 variants that stratify along 12 mitochondrial haplotypes. We identify 1,845 cases of mitonuclear allelic imbalances, thus implying that mitochondrial haplotypes are reflected in the nuclear genome. However, no major fitness effects are associated with mitonuclear imbalance, suggesting that such imbalances reflect population structure at the mitochondrial level rather than genomic incompatibilities. Although mitochondrial haplotypes have no direct impact on mitochondrial respiration, some haplotypes are associated with stress- and metabolism-related phenotypes, including food intake in males. Finally, through reciprocal swapping of mitochondrial genomes, we demonstrate that a mitochondrial haplotype associated with high food intake can rescue a low food intake phenotype. Together, our findings provide new insight into population structure at the mitochondrial level and point to the importance of incorporating mitochondrial haplotypes in genotype-phenotype relationship studies.