ABSTRACT
In various organisms, sequencing of selectively bred lines at apparent selection limits has demonstrated that genetic variation can remain at many loci, implying that evolution at the genetic level may continue even if the population mean phenotype remains constant. We compared selection signatures at generations 22 and 61 of the "High Runner" mouse experiment, which includes 4 replicate lines bred for voluntary wheel-running behavior (HR) and 4 non-selected control (C) lines. Previously, we reported multiple regions of differentiation between the HR and C lines, based on whole-genome sequence data for 10 mice from each line at generation 61, which was >31 generations after selection limits had been reached in all HR lines. Here, we analyzed pooled sequencing data from ~20 mice for each of the 8 lines at generation 22, around when HR lines were reaching limits. Differentiation analyses of allele frequencies at ~4.4 million SNP loci used the regularized T-test and detected 258 differentiated regions with FDR = 0.01. Comparable analyses involving pooling generation 61 individual mouse genotypes into allele frequencies by line produced only 11 such regions, with almost no overlap among the largest and most statistically significant peaks between the two generations. These results implicate a sort of "genetic churn" that continues at loci relevant for running. Simulations indicate that loss of statistical power due to random genetic drift and sampling error are insufficient to explain the differences in selection signatures. The 13 differentiated regions at generation 22 with strict culling measures include 79 genes related to a wide variety of functions. Gene ontology identified pathways related to olfaction and vomeronasal pathways as being overrepresented, consistent with generation 61 analyses, despite those specific regions differing between generations. Genes Dspp and Rbm24 are also identified as potentially explaining known bone and skeletal muscle differences, respectively, between the linetypes.
Subject(s)
Polymorphism, Single Nucleotide , Selection, Genetic , Animals , Mice , Gene Frequency , Running , Male , Phenotype , Genotype , Female , Behavior, Animal/physiology , Selective Breeding/geneticsABSTRACT
Apes possess two sex chromosomes-the male-specific Y chromosome and the X chromosome, which is present in both males and females. The Y chromosome is crucial for male reproduction, with deletions being linked to infertility1. The X chromosome is vital for reproduction and cognition2. Variation in mating patterns and brain function among apes suggests corresponding differences in their sex chromosomes. However, owing to their repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the methodology developed for the telomere-to-telomere (T2T) human genome, we produced gapless assemblies of the X and Y chromosomes for five great apes (bonobo (Pan paniscus), chimpanzee (Pan troglodytes), western lowland gorilla (Gorilla gorilla gorilla), Bornean orangutan (Pongo pygmaeus) and Sumatran orangutan (Pongo abelii)) and a lesser ape (the siamang gibbon (Symphalangus syndactylus)), and untangled the intricacies of their evolution. Compared with the X chromosomes, the ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements-owing to the accumulation of lineage-specific ampliconic regions, palindromes, transposable elements and satellites. Many Y chromosome genes expand in multi-copy families and some evolve under purifying selection. Thus, the Y chromosome exhibits dynamic evolution, whereas the X chromosome is more stable. Mapping short-read sequencing data to these assemblies revealed diversity and selection patterns on sex chromosomes of more than 100 individual great apes. These reference assemblies are expected to inform human evolution and conservation genetics of non-human apes, all of which are endangered species.
Subject(s)
Hominidae , X Chromosome , Y Chromosome , Animals , Female , Male , Gorilla gorilla/genetics , Hominidae/genetics , Hominidae/classification , Hylobatidae/genetics , Pan paniscus/genetics , Pan troglodytes/genetics , Phylogeny , Pongo abelii/genetics , Pongo pygmaeus/genetics , Telomere/genetics , X Chromosome/genetics , Y Chromosome/genetics , Evolution, Molecular , DNA Copy Number Variations/genetics , Humans , Endangered Species , Reference StandardsABSTRACT
Apes possess two sex chromosomes-the male-specific Y and the X shared by males and females. The Y chromosome is crucial for male reproduction, with deletions linked to infertility. The X chromosome carries genes vital for reproduction and cognition. Variation in mating patterns and brain function among great apes suggests corresponding differences in their sex chromosome structure and evolution. However, due to their highly repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the state-of-the-art experimental and computational methods developed for the telomere-to-telomere (T2T) human genome, we produced gapless, complete assemblies of the X and Y chromosomes for five great apes (chimpanzee, bonobo, gorilla, Bornean and Sumatran orangutans) and a lesser ape, the siamang gibbon. These assemblies completely resolved ampliconic, palindromic, and satellite sequences, including the entire centromeres, allowing us to untangle the intricacies of ape sex chromosome evolution. We found that, compared to the X, ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements. This divergence on the Y arises from the accumulation of lineage-specific ampliconic regions and palindromes (which are shared more broadly among species on the X) and from the abundance of transposable elements and satellites (which have a lower representation on the X). Our analysis of Y chromosome genes revealed lineage-specific expansions of multi-copy gene families and signatures of purifying selection. In summary, the Y exhibits dynamic evolution, while the X is more stable. Finally, mapping short-read sequencing data from >100 great ape individuals revealed the patterns of diversity and selection on their sex chromosomes, demonstrating the utility of these reference assemblies for studies of great ape evolution. These complete sex chromosome assemblies are expected to further inform conservation genetics of nonhuman apes, all of which are endangered species.
ABSTRACT
In most animals, mitochondrial DNA is strictly maternally inherited and non-recombining. One exception to this pattern is called doubly uniparental inheritance (DUI), a phenomenon involving the independent transmission of female and male mitochondrial genomes. DUI is known only from the molluskan class Bivalvia. The phylogenetic distribution of male-transmitted mitochondrial DNA (M mtDNA) in bivalves is consistent with several evolutionary scenarios, including multiple independent gains, losses, and varying degrees of recombination with female-transmitted mitochondrial DNA (F mtDNA). In this study, we use phylogenetic methods to test M mtDNA origination hypotheses and infer the prevalence of mitochondrial recombination in bivalves with DUI. Phylogenetic modeling using site concordance factors supported a single origin of M mtDNA in bivalves coupled with recombination acting over long evolutionary timescales. Ongoing mitochondrial recombination is present in Mytilida and Venerida, which results in a pattern of concerted evolution of F mtDNA and M mtDNA. Mitochondrial recombination could be favored to offset the deleterious effects of asexual inheritance and maintain mitonuclear compatibility across tissues. Cardiida and Unionida have gone without recent recombination, possibly due to an extension of the COX2 gene in male mitochondrial DNA. The loss of recombination could be connected to the role of M mtDNA in sex determination or sexual development. Our results support that recombination events may occur throughout the mitochondrial genomes of DUI species. Future investigations may reveal more complex patterns of inheritance of recombinants, which could explain the retention of signal for a single origination of M mtDNA in protein-coding genes.
Subject(s)
Bivalvia , Genome, Mitochondrial , Animals , Female , Male , Phylogeny , Mitochondria/genetics , Bivalvia/genetics , DNA, Mitochondrial/genetics , Inheritance Patterns , Recombination, GeneticABSTRACT
Restriction-site-associated DNA sequencing (RADseq) has become an accessible way to obtain genome-wide data in the form of single-nucleotide polymorphisms (SNPs) for phylogenetic inference. Nonetheless, how differences in RADseq methods influence phylogenetic estimation is poorly understood because most comparisons have largely relied on conceptual predictions rather than empirical tests. We examine how differences in ddRAD and 2bRAD data influence phylogenetic estimation in two non-model frog groups. We compare the impact of method choice on phylogenetic information, missing data, and allelic dropout, considering different sequencing depths. Given that researchers must balance input (funding, time) with output (amount and quality of data), we also provide comparisons of laboratory effort, computational time, monetary costs, and the repeatability of library preparation and sequencing. Both 2bRAD and ddRAD methods estimated well-supported trees, even at low sequencing depths, and had comparable amounts of missing data, patterns of allelic dropout, and phylogenetic signal. Compared to ddRAD, 2bRAD produced more repeatable datasets, had simpler laboratory protocols, and had an overall faster bioinformatics assembly. However, many fewer parsimony-informative sites per SNP were obtained from 2bRAD data when using native pipelines, highlighting a need for further investigation into the effects of each pipeline on resulting datasets. Our study underscores the importance of comparing RADseq methods, such as expected results and theoretical performance using empirical datasets, before undertaking costly experiments.
ABSTRACT
Fructose (C6H12O6) is acutely obesogenic and is a risk factor for hypertension, cardiovascular disease, and nonalcoholic fatty liver disease. However, the possible long-lasting effects of early-life fructose consumption have not been studied. We tested for effects of early-life fructose and/or wheel access (voluntary exercise) in a line of selectively bred High Runner (HR) mice and a non-selected Control (C) line. Exposures began at weaning and continued for 3 weeks to sexual maturity, followed by a 23-week "washout" period (equivalent to â¼17 human years). Fructose increased total caloric intake, body mass, and body fat during juvenile exposure, but had no effect on juvenile wheel running and no important lasting effects on adult physical activity or body weight/composition. Interestingly, adult maximal aerobic capacity (VO2max) was reduced in mice that had early-life fructose and wheel access. Consistent with previous studies, early-life exercise promoted adult wheel running. In a 3-way interaction, C mice that had early-life fructose and no wheel access gained body mass in response to 2 weeks of adult wheel access, while all other groups lost mass. Overall, we found some long-lasting positive effects of early-life exercise, but minimal effects of early-life fructose, regardless of the mouse line.
Subject(s)
Motor Activity , Physical Conditioning, Animal , Humans , Mice , Animals , Motor Activity/physiology , Body Composition/physiology , Adipose Tissue , Energy Intake , Physical Conditioning, Animal/physiologyABSTRACT
Replicate lines under uniform selection often evolve in different ways. Previously, analyses using whole-genome sequence data for individual mice (Mus musculus) from 4 replicate High Runner lines and 4 nonselected control lines demonstrated genomic regions that have responded consistently to selection for voluntary wheel-running behavior. Here, we ask whether the High Runner lines have evolved differently from each other, even though they reached selection limits at similar levels. We focus on 1 High Runner line (HR3) that became fixed for a mutation at a gene of major effect (Myh4Minimsc) that, in the homozygous condition, causes a 50% reduction in hindlimb muscle mass and many pleiotropic effects. We excluded HR3 from SNP analyses and identified 19 regions not consistently identified in analyses with all 4 lines. Repeating analyses while dropping each of the other High Runner lines identified 12, 8, and 6 such regions. (Of these 45 regions, 37 were unique.) These results suggest that each High Runner line indeed responded to selection somewhat uniquely, but also that HR3 is the most distinct. We then applied 2 additional analytical approaches when dropping HR3 only (based on haplotypes and nonstatistical tests involving fixation patterns). All 3 approaches identified 7 new regions (as compared with analyses using all 4 High Runner lines) that include genes associated with activity levels, dopamine signaling, hippocampus morphology, heart size, and body size, all of which differ between High Runner and control lines. Our results illustrate how multiple solutions and "private" alleles can obscure general signatures of selection involving "public" alleles.
Subject(s)
Genomics , Selective Breeding , Mice , Animals , Body Size , LocomotionABSTRACT
With limited sampling, geographic variation within a single species can be difficult to distinguish from interspecific variation, confounding our ability to draw accurate species boundaries. We argue that thorough sampling and analysis of contact zones between putative taxa can determine if assortative mating or selection against hybrids exists (supporting the presence of two distinct species), or alternatively if mating is random among genotypes and admixture among adjacent populations is gradual and continuous (supporting geographic variation within a single species). Here, we test two alternative hypotheses for two pairs of named taxa at contact zones within the American milksnake (Lampropeltis triangulum) complex. A prior morphological analysis found areas of gradual intergradation among named taxa, and concluded that the taxa represented geographical races of a single polytypic species. In contrast, a subsequent analysis of gene sequence data, but with limited sampling near the contact zones, hypothesized distinct boundaries between species at the contact zones. At the contact zone between proposed species L. triangulum and Lampropeltis gentilis, we examined a $\sim$700 km-wide transect across the states of Kansas and Missouri, with thorough sampling and reduced-representation genomic-level sequencing, to test the two opposing taxonomic hypotheses. Our transect analyses included examinations of population structure, fixed differences, cline-fitting, and an admixture index analysis. These analyses all supported a gradual and continuous geographic cline across a broad intergrade zone between two geographic forms of L. triangulum, thus providing strong support for a single species in this region (and no support for the recognition of L. gentilis as a distinct species). At a second contact zone between proposed species L. triangulum and Lampropeltis elapsoides (but variously treated as species or subspecies by different researchers) in Kentucky and Tennessee, we re-evaluated morphological data. In this case, the contact zone analysis indicated sympatry and reproductive isolation of the two taxa, and thus strongly supported L. triangulum and L. elapsoides as distinct species. We conclude that detailed studies of contact zones, based on either genetic or morphological data, are essential for distinguishing intraspecific from interspecific variation in the case of widely and continuously distributed taxa. [Contact zones; speciation; species concepts; species delimitation; taxonomy.].
Subject(s)
Genetics, Population , Reproductive Isolation , Phylogeny , Genotype , Geography , Hybridization, GeneticABSTRACT
Defense against ultraviolet (UV) radiation exposure is essential for survival, especially in high-elevation species. Although some specific genes involved in UV response have been reported, the full view of UV defense mechanisms remains largely unexplored. Herein, we used integrated approaches to analyze UV responses in the highest-elevation frog, Nanorana parkeri. We show less damage and more efficient antioxidant activity in skin of this frog than those of its lower-elevation relatives after UV exposure. We also reveal genes related to UV defense and a corresponding temporal expression pattern in N. parkeri. Genomic and metabolomic analysis along with large-scale transcriptomic profiling revealed a time-dependent coordinated defense mechanism in N. parkeri. We also identified several microRNAs that play important regulatory roles, especially in decreasing the expression levels of cell cycle genes. Moreover, multiple defense genes (i.e., TYR for melanogenesis) exhibit positive selection with function-enhancing substitutions. Thus, both expression shifts and gene mutations contribute to UV adaptation in N. parkeri. Our work demonstrates a genetic framework for evolution of UV defense in a natural environment.
Subject(s)
Anura , Ultraviolet Rays , Animals , Anura/genetics , Skin , Gene Expression Profiling , AntioxidantsABSTRACT
Although physical activity and sedentary behavior are moderately heritable, little is known about the mechanisms that influence these traits. Combining data for up to 703,901 individuals from 51 studies in a multi-ancestry meta-analysis of genome-wide association studies yields 99 loci that associate with self-reported moderate-to-vigorous intensity physical activity during leisure time (MVPA), leisure screen time (LST) and/or sedentary behavior at work. Loci associated with LST are enriched for genes whose expression in skeletal muscle is altered by resistance training. A missense variant in ACTN3 makes the alpha-actinin-3 filaments more flexible, resulting in lower maximal force in isolated type IIA muscle fibers, and possibly protection from exercise-induced muscle damage. Finally, Mendelian randomization analyses show that beneficial effects of lower LST and higher MVPA on several risk factors and diseases are mediated or confounded by body mass index (BMI). Our results provide insights into physical activity mechanisms and its role in disease prevention.
Subject(s)
Genome-Wide Association Study , Sedentary Behavior , Actinin/genetics , Cross-Sectional Studies , Exercise/physiology , Humans , Leisure ActivitiesABSTRACT
The gut microbiome can affect various aspects of both behavior and physiology, including exercise ability, but effects on voluntary exercise have rarely been studied. We studied females from a selection experiment in which 4 replicate High Runner (HR) lines of mice are bred for voluntary exercise and compared with 4 non-selected control (C) lines. HR and C mice differ in several traits that likely interact with the gut microbiome, including higher daily running distance, body temperatures when running, spontaneous physical activity when housed without wheels, and food consumption. After two weeks of wheel access to reach a stable plateau in daily running, mice were administered broad-spectrum antibiotics for 10 days. Antibiotic treatment caused a significant reduction in daily wheel-running distance in the HR mice (-21%) but not in the C mice. Antibiotics did not affect body mass or food consumption in either HR or C mice, and we did not observe sickness behavior. Wheel running by HR mice did not recover during the 12 days following cessation of antibiotics. The decreased wheel-running in HR but not C mice, with no apparent negative side effects of antibiotics, suggests that the HR microbiome is an important component of their high-running phenotype.
Subject(s)
Physical Conditioning, Animal , Sports , Animals , Anti-Bacterial Agents/pharmacology , Female , Mice , Motor Activity/physiology , Phenotype , Physical Conditioning, Animal/physiologyABSTRACT
SignificanceTo adapt to arboreal lifestyles, treefrogs have evolved a suite of complex traits that support vertical movement and gliding, thus presenting a unique case for studying the genetic basis for traits causally linked to vertical niche expansion. Here, based on two de novo-assembled Asian treefrog genomes, we determined that genes involved in limb development and keratin cytoskeleton likely played a role in the evolution of their climbing systems. Behavioral and morphological evaluation and time-ordered gene coexpression network analysis revealed the developmental patterns and regulatory pathways of the webbed feet used for gliding in Rhacophorus kio.
Subject(s)
Locomotion , Trees , Adaptation, Physiological/genetics , Animals , Anura , Biological Evolution , Biomechanical Phenomena , Genomics , Humans , Locomotion/geneticsABSTRACT
The Himalaya are among the youngest and highest mountains in the world, but the exact timing of their uplift and origins of their biodiversity are still in debate. The Himalayan region is a relatively small area but with exceptional diversity and endemism. One common hypothesis to explain the rich montane diversity is uplift-driven diversification-that orogeny creates conditions favoring rapid in situ speciation of resident lineages. We test this hypothesis in the Himalayan region using amphibians and reptiles, two environmentally sensitive vertebrate groups. In addition, analysis of diversification of the herpetofauna provides an independent source of information to test competing geological hypotheses of Himalayan orogenesis. We conclude that the origins of the Himalayan herpetofauna date to the early Paleocene, but that diversification of most groups was concentrated in the Miocene. There was an increase in both rates and modes of diversification during the early to middle Miocene, together with regional interchange (dispersal) between the Himalaya and adjacent regions. Our analyses support a recently proposed stepwise geological model of Himalayan uplift beginning in the Paleocene, with a subsequent rapid increase of uplifting during the Miocene, finally giving rise to the intensification of the modern South Asian Monsoon.
ABSTRACT
While nursing, mammals progress through critical developmental periods for the cardiovascular, musculoskeletal, and central nervous systems. The suckling period in mammals is therefore especially vulnerable to environmental factors that may affect the "developmental programming" of many complex traits. As a result, various aspects of maternal behavior and physiology can influence offspring in ways that have lasting effects into adulthood. Several recent studies of animal models have shown that maternal effects can partially program adult activity behaviors, which has important implications for health and locomotor performance. Here, we used cross-fostering to test for possible maternal effects on adult wheel-running behavior (voluntary exercise), maximal aerobic capacity during forced exercise (VO2max), body mass and composition, and organ masses. Subjects were from a line of mice that has been selectively bred for â¼90 generations for high voluntary wheel-running behavior (High Runner; HR) and a non-selected Control (C) line. Adult HR mice run â¼3-fold the daily distances of C mice and have evolved other differences associated with exercise capacity, including elevated VO2max, reduced body mass and fat mass, and larger hearts. At birth, we fostered offspring to create 4 experimental groups: C pups to other C dams (in-foster), HR pups to other HR dams (in-foster), C pups to HR dams (cross-foster), HR pups to C dams (cross-foster). Thus, all pups were fostered to a different mother. Mice were weaned 3 weeks later, and adult testing began at â¼6 weeks of age. At weaning, pups raised by HR dams were smaller than those raised by C dams for both sexes and as expected, HR pups raised by HR dams weighed less than C pups raised by C dams. As adults, mice raised by HR dams continued to have reduced body masses. As expected, adult HR mice ran approximately 3-fold more than their C counterparts and females ran more than males. However, cross-fostering did not statistically affect any aspect of wheel-running behavior (distance, duration, speed). Similarly, with body mass as a covariate, HR mice had higher VO2max than C mice, and males had higher VO2max than females, but cross-fostering had no effect. With body mass as a covariate, cross-fostering had variable effects on adult organ masses in a sex-specific manner. Overall, our results indicate that development of the adult High Runner phenotype does not require rearing by an HR dam, suggesting that high adult activity in humans may be independent of high maternal activity.
Subject(s)
Motor Activity , Physical Conditioning, Animal , Adult , Animals , Female , Humans , Male , Maternal Behavior , Mice , PhenotypeABSTRACT
As DNA sequencing technologies and methods for delimiting species with genomic data become more accessible and numerous, researchers have more tools than ever to investigate questions in systematics and phylogeography. However, easy access to sophisticated computational tools is not without its drawbacks. Choosing the right approach for one's question can be challenging when presented with multitudinous options, some of which fail to distinguish between species and intraspecific population structure. Here, we employ a methodology that emphasizes intensive geographic sampling, particularly at contact zones between populations, with a focus on differentiating intraspecific genetic clusters from species in the Pantherophis guttatus complex, a group of North American ratsnakes. Using a mitochondrial marker as well as ddRADseq data, we find evidence of mitonuclear discordance which has contributed to historical confusion about the relationships within this group. Additionally, we identify geographically and genetically structured populations within the species Pantherophis emoryi that are congruent with previously described morphological variation. Importantly, we find that these structured populations within P. emoryi are highly admixed throughout the range of the species and show no evidence of any reproductive isolation. Our data support a revision of the taxonomy of this group, and we recognize two species within the complex and three subspecies within P. emoryi. This study illustrates the importance of thorough sampling of contact zones and consideration of gene flow when delimiting species in widespread complexes containing parapatric lineages.
Subject(s)
Cell Nucleus/genetics , DNA, Mitochondrial/genetics , Genetic Variation , Genomics , Phylogeny , Phylogeography , Snakes/genetics , Animals , Sequence Analysis, DNAABSTRACT
Ethologically relevant chemical senses and behavioral habits are likely to coadapt in response to selection. As olfaction is involved in intrinsically motivated behaviors in mice, we hypothesized that selective breeding for a voluntary behavior would enable us to identify novel roles of the chemosensory system. Voluntary wheel running (VWR) is an intrinsically motivated and naturally rewarding behavior, and even wild mice run on a wheel placed in nature. We have established 4 independent, artificially evolved mouse lines by selectively breeding individuals showing high VWR activity (High Runners; HRs), together with 4 non-selected Control lines, over 88 generations. We found that several sensory receptors in specific receptor clusters were differentially expressed between the vomeronasal organ (VNO) of HRs and Controls. Moreover, one of those clusters contains multiple single-nucleotide polymorphism loci for which the allele frequencies were significantly divergent between the HR and Control lines, i.e., loci that were affected by the selective breeding protocol. These results indicate that the VNO has become genetically differentiated between HR and Control lines during the selective breeding process. Although the role of the vomeronasal chemosensory receptors in VWR activity remains to be determined, the current results suggest that these vomeronasal chemosensory receptors are important quantitative trait loci for voluntary exercise in mice. We propose that olfaction may play an important role in motivation for voluntary exercise in mammals.
Subject(s)
Behavior, Animal , Physical Conditioning, Animal , Vomeronasal Organ/metabolism , Animals , Female , Gene Frequency , Genetic Loci , Male , Mice , Mice, Inbred ICR , Polymorphism, Single Nucleotide , Selection, GeneticABSTRACT
The biological basis of exercise behavior is increasingly relevant for maintaining healthy lifestyles. Various quantitative genetic studies and selection experiments have conclusively demonstrated substantial heritability for exercise behavior in both humans and laboratory rodents. In the "High Runner" selection experiment, four replicate lines of Mus domesticus were bred for high voluntary wheel running (HR), along with four nonselected control (C) lines. After 61 generations, the genomes of 79 mice (9-10 from each line) were fully sequenced and single nucleotide polymorphisms (SNPs) were identified. We used nested ANOVA with MIVQUE estimation and other approaches to compare allele frequencies between the HR and C lines for both SNPs and haplotypes. Approximately 61 genomic regions, across all somatic chromosomes, showed evidence of differentiation; 12 of these regions were differentiated by all methods of analysis. Gene function was inferred largely using Panther gene ontology terms and KO phenotypes associated with genes of interest. Some of the differentiated genes are known to be associated with behavior/motivational systems and/or athletic ability, including Sorl1, Dach1, and Cdh10 Sorl1 is a sorting protein associated with cholinergic neuron morphology, vascular wound healing, and metabolism. Dach1 is associated with limb bud development and neural differentiation. Cdh10 is a calcium ion binding protein associated with phrenic neurons. Overall, these results indicate that selective breeding for high voluntary exercise has resulted in changes in allele frequencies for multiple genes associated with both motivation and ability for endurance exercise, providing candidate genes that may explain phenotypic changes observed in previous studies.
Subject(s)
Directed Molecular Evolution , Polymorphism, Single Nucleotide , Running , Selection, Genetic , Animals , Cadherins/genetics , Chromosomes/genetics , Eye Proteins/genetics , Female , Hybridization, Genetic , Male , Membrane Transport Proteins/genetics , Mice , Mice, Inbred ICR , Multifactorial Inheritance , Receptors, LDL/geneticsABSTRACT
Sexually selected traits can be expected to increase in importance when the period of sexual behavior is constrained, such as in seasonally restricted breeders. Anolis lizard male dewlaps are classic examples of multifaceted signaling traits, with demonstrated intraspecific reproductive function reflected in courtship behavior. Fitch and Hillis found a correlation between dewlap size and seasonality in mainland Anolis using traditional statistical methods and suggested that seasonally restricted breeding seasons enhanced the differentiation of this signaling trait. Here, we present two tests of the Fitch-Hillis Hypothesis using new phylogenetic and morphological data sets for 44 species of Mexican Anolis. A significant relationship between dewlap size and seasonality is evident in phylogenetically uncorrected analyses but erodes once phylogeny is accounted for. This loss of strong statistical support for a relationship between a key aspect of dewlap morphology and seasonality also occurs within a species complex (A. sericeus group) that inhabits seasonal and aseasonal environments. Our results fail to support seasonality as a strong driver of evolution of Anolis dewlap size. We discuss the implications of our results and the difficulty of disentangling the strength of single mechanisms on trait evolution when multiple selection pressures are likely at play.
ABSTRACT
Many recent species delimitation studies rely exclusively on limited analyses of genetic data analyzed under the multispecies coalescent (MSC) model, and results from these studies often are regarded as conclusive support for taxonomic changes. However, most MSC-based species delimitation methods have well-known and often unmet assumptions. Uncritical application of these genetic-based approaches (without due consideration of sampling design, the effects of a priori group designations, isolation by distance, cytoplasmic-nuclear mismatch, and population structure) can lead to over-splitting of species. Here, we argue that in many common biological scenarios, researchers must be particularly cautious regarding these limitations, especially in cases of well-studied, geographically variable, and parapatrically distributed species complexes. We consider these points with respect to a historically controversial species group, the American milksnakes (Lampropeltis triangulum complex), using genetic data from a recent analysis (Ruane et al. 2014). We show that over-reliance on the program Bayesian Phylogenetics and Phylogeography, without adequate consideration of its assumptions and of sampling limitations, resulted in over-splitting of species in this study. Several of the hypothesized species of milksnakes instead appear to represent arbitrary slices of continuous geographic clines. We conclude that the best available evidence supports three, rather than seven, species within this complex. More generally, we recommend that coalescent-based species delimitation studies incorporate thorough analyses of geographic variation and carefully examine putative contact zones among delimited species before making taxonomic changes.
Subject(s)
Classification/methods , Phylogeography , Animals , Phylogeny , Snakes/classification , Snakes/geneticsABSTRACT
Viviparous (live-bearing) vertebrates have evolved repeatedly within otherwise oviparous (egg-laying) clades. Over two-thirds of these changes in vertebrate reproductive parity mode happened in squamate reptiles, where the transition has happened between 98 and 129 times. The transition from oviparity to viviparity requires numerous physiological, morphological, and immunological changes to the female reproductive tract, including eggshell reduction, delayed oviposition, placental development for supply of water and nutrition to the embryo by the mother, enhanced gas exchange, and suppression of maternal immune rejection of the embryo. We performed genomic and transcriptomic analyses of a closely related oviparous-viviparous pair of lizards (Phrynocephalus przewalskii and Phrynocephalus vlangalii) to examine these transitions. Expression patterns of maternal oviduct through reproductive development of the egg and embryo differ markedly between the two species. We found changes in expression patterns of appropriate genes that account for each of the major aspects of the oviparity to viviparity transition. In addition, we compared the gene sequences in transcriptomes of four oviparous-viviparous pairs of lizards in different genera (Phrynocephalus, Eremias, Scincella, and Sphenomorphus) to look for possible gene convergence at the sequence level. We discovered low levels of convergence in both amino acid replacement and evolutionary rate shift. This suggests that most of the changes that produce the oviparity-viviparity transition are changes in gene expression, so occasional reversals to oviparity from viviparity may not be as difficult to achieve as has been previously suggested.