Rates Of Sequence Evolution Research Articles

CURRENT PERSPECTIVES OF DNA BARCODING IN BOSNIA AND HERZEGOVINA THROUGH BOLD DATABASE

UDK: 577.212:574]:007.5(497.6)
 DNA barcoding is a method designed to provide rapid and precise species identifications by using one or more of short gene sequences called barcodes. In most plant and fungi studies, the standard barcodes of choice are three plastid (rbcL, matK and trnH-psbA) and one nuclear (ITS) gene regions. The relatively high, but comparatively conserved rate of sequence evolution of mtDNA has made COI the marker of choice in animals. BOLD is a freely available cloud-based data storage and analysis platform developed with the aim to advance biodiversity science through DNA barcoding species identification. To date, over 6 million barcodes have been deposited in BOLD with 196,000 animal species, 68,000 plant species and 22,000 species of fungi and other organism entries. In this database, there are currently 447 entries for organisms from Bosnia and Herzegovina, which makes 0.0067% of the total number of BOLD entries. According to BOLD statistics, only 1.11% of all organism entries from B&H were submitted by B&H institutions.
 Despite the fact that Bosnia and Herzegovina has valuable natural resources with a high percentage of endemic and autochthonous species, BOLD statistics elucidated the lack of coordinated and systematic DNA barcoding research so far. It is necessary to establish continuous progress of molecular-genetic characterization of these resources in the future. It is up to B&H institutions to decide if they want to continue the practice of submitting the data sporadically or if they will animate the research community to actively participate in this global project.

Sex-biased gene expression in an hemimetabolous insect: pattern during development, extent, functions involved, rate of sequence evolution, and comparison with an holometabolous insect

A recommendation of: Jelisaveta Djordjevic, Zoé Dumas, Marc Robinson-Rechavi, Tanja Schwander, Darren James Parker Dynamics of sex-biased gene expression over development in the stick insect <i>Timema californicum</i> https://doi.org/10.1101/2021.01.23.427895

Faster Rates of Molecular Sequence Evolution in Reproduction-Related Genes and in Species with Hypodermic Sperm Morphologies.

Sexual selection drives the evolution of many striking behaviors and morphologies and should leave signatures of selection at loci underlying these phenotypes. However, although loci thought to be under sexual selection often evolve rapidly, few studies have contrasted rates of molecular sequence evolution at such loci across lineages with different sexual selection contexts. Furthermore, work has focused on separate sexed animals, neglecting alternative sexual systems. We investigate rates of molecular sequence evolution in hermaphroditic flatworms of the genus Macrostomum. Specifically, we compare species that exhibit contrasting sperm morphologies, strongly associated with multiple convergent shifts in the mating strategy, reflecting different sexual selection contexts. Species donating and receiving sperm in every mating have sperm with bristles, likely to prevent sperm removal. Meanwhile, species that hypodermically inject sperm lack bristles, potentially as an adaptation to the environment experienced by hypodermic sperm. Combining functional annotations from the model, Macrostomum lignano, with transcriptomes from 93 congeners, we find genus-wide faster sequence evolution in reproduction-related versus ubiquitously expressed genes, consistent with stronger sexual selection on the former. Additionally, species with hypodermic sperm morphologies had elevated molecular sequence evolution, regardless of a gene’s functional annotation. These genome-wide patterns suggest reduced selection efficiency following shifts to hypodermic mating, possibly due to higher selfing rates in these species. Moreover, we find little evidence for convergent amino acid changes across species. Our work not only shows that reproduction-related genes evolve rapidly also in hermaphroditic animals, but also that well-replicated contrasts of different sexual selection contexts can reveal underappreciated genome-wide effects.

Variables Influencing Differences in Sequence Conservation in the Fission Yeast Schizosaccharomyces pombe

Which variables determine the constraints on gene sequence evolution is one of the most central questions in molecular evolution. In the fission yeast Schizosaccharomyces pombe, an important model organism, the variables influencing the rate of sequence evolution have yet to be determined. Previous studies in other single celled organisms have generally found gene expression levels to be most significant, with numerous other variables such as gene length and functional importance identified as having a smaller impact. Using publicly available data, we used partial least squares regression, principal components regression, and partial correlations to determine the variables most strongly associated with sequence evolution constraints. We identify centrality in the protein–protein interactions network, amino acid composition, and cellular location as the most important determinants of sequence conservation. However, each factor only explains a small amount of variance, and there are numerous variables having a significant or heterogeneous influence. Our models explain more than half of the variance in dN, raising the possibility that future refined models could quantify the role of stochastics in evolutionary rate variation.

Genome assembly of primitive cultivated potato Solanum stenotomum provides insights into potato evolution.

Genetic diversity is the raw material for germplasm enhancement. Landraces and wild species relatives of potato, which contain a rich gene pool of valuable agronomic traits, can provide insights into the genetic diversity behind the adaptability of the common potato. The diploid plant, Solanum stenotomum (Sst), is believed to have an ancestral relationship with modern potato cultivars and be a potential source of resistance against disease. Sequencing of the Sst genome generated an assembly of 852.85 Mb (N50 scaffold size, 3.7 Mb). Pseudomolecule construction anchored 788.75 Mb of the assembly onto 12 pseudochromosomes, with an anchor rate of 92.4%. Genome annotation yielded 41,914 high-confidence protein-coding gene models and comparative analyses with closely related Solanaceae species identified 358 Sst-specific gene families, 885 gene families with expansion along the Sst lineage, and 149 genes experiencing accelerated rates of protein sequence evolution in Sst, the functions of which were mainly associated with defense responses, particularly against bacterial and fungal infection. Insights into the Sst genome and the genomic variation of cultivated potato taxa are valuable in elaborating the impact of potato evolution in early landrace diploid and facilitate modern potato breeding.

Evolutionary dynamics of sex-biased genes expressed in cricket brains and gonads.

Sex-biased gene expression, particularly sex-biased expression in the gonad, has been linked to rates of protein sequence evolution (nonsynonymous to synonymous substitutions, dN/dS) in animals. However, in insects, sex-biased expression studies remain centred on a few holometabolous species. Moreover, other major tissue types such as the brain remain underexplored. Here, we studied sex-biased gene expression and protein evolution in a hemimetabolous insect, the cricket Gryllus bimaculatus. We generated novel male and female RNA-seq data for two sexual tissue types, the gonad and somatic reproductive system, and for two core components of the nervous system, the brain and ventral nerve cord. From a genome-wide analysis, we report several core findings. Firstly, testis-biased genes had accelerated evolution, as compared to ovary-biased and unbiased genes, which was associated with positive selection events. Secondly, although sex-biased brain genes were much less common than for the gonad, they exhibited a striking tendency for rapid protein sequence evolution, an effect that was stronger for the female than male brain. Further, some sex-biased brain genes were linked to sexual functions and mating behaviours, which we suggest may have accelerated their evolution via sexual selection. Thirdly, a tendency for narrow cross-tissue expression breadth, suggesting low pleiotropy, was observed for sex-biased brain genes, suggesting relaxed purifying selection, which we speculate may allow enhanced freedom to evolve adaptive protein functional changes. The findings of rapid evolution of testis-biased genes and male and female-biased brain genes are discussed with respect to pleiotropy, positive selection and the mating biology of this cricket.

Phylogenetic Inference and the Misplaced Premise of Substitution Rates.

Three competing 'methods' have been endorsed for inferring phylogenetic hypotheses: parsimony, likelihood, and Bayesianism. The latter two have been claimed superior because they take into account rates of sequence substitution. Can rates of substitution be justified on its own accord in inferences of explanatory hypotheses? Answering this question requires addressing four issues: (1) the aim of scientific inquiry, (2) the nature of why-questions, (3) explanatory hypotheses as answers to why-questions, and (4) acknowledging that neither parsimony, likelihood, nor Bayesianism are inferential actions leading to explanatory hypotheses. The aim of scientific inquiry is to acquire causal understanding of effects. Observation statements of organismal characters lead to implicit or explicit why-questions. Those questions, conveyed in data matrices, assume the truth of observation statements, which is contrary to subsequently invoking substitution rates within inferences to phylogenetic hypotheses. Inferences of explanatory hypotheses are abductive in form, such that some version of an evolutionary theory(ies) is/are included or implied. If rates of sequence evolution are to be considered, it must be done prior to, rather than within abduction, which requires renaming those putatively-shared nucleotides subject to substitution rates. There are, however, no epistemic grounds for renaming characters to accommodate rates, calling into question the legitimacy of causally accounting for sequence data.

The Missing Expression Level-Evolutionary Rate Anticorrelation in Viruses Does Not Support Protein Function as a Main Constraint on Sequence Evolution.

One of the central goals in molecular evolutionary biology is to determine the sources of variation in the rate of sequence evolution among proteins. Gene expression level is widely accepted as the primary determinant of protein evolutionary rate, because it scales with the extent of selective constraints imposed on a protein, leading to the well-known negative correlation between expression level and protein evolutionary rate (the E–R anticorrelation). Selective constraints have been hypothesized to entail the maintenance of protein function, the avoidance of cytotoxicity caused by protein misfolding or nonspecific protein–protein interactions, or both. However, empirical tests evaluating the relative importance of these hypotheses remain scarce, likely due to the nontrivial difficulties in distinguishing the effect of a deleterious mutation on a protein’s function versus its cytotoxicity. We realized that examining the sequence evolution of viral proteins could overcome this hurdle. It is because purifying selection against mutations in a viral protein that result in cytotoxicity per se is likely relaxed, whereas purifying selection against mutations that impair viral protein function persists. Multiple analyses of SARS-CoV-2 and nine other virus species revealed a complete absence of any E–R anticorrelation. As a control, the E–R anticorrelation does exist in human endogenous retroviruses where purifying selection against cytotoxicity is present. Taken together, these observations do not support the maintenance of protein function as the main constraint on protein sequence evolution in cellular organisms.

Detecting de novo mitochondrial mutations in angiosperms with highly divergent evolutionary rates.

Although plant mitochondrial genomes typically show low rates of sequence evolution, levels of divergence in certain angiosperm lineages suggest anomalously high mitochondrial mutation rates. However, de novo mutations have never been directly analyzed in such lineages. Recent advances in high-fidelity DNA sequencing technologies have enabled detection of mitochondrial mutations when still present at low heteroplasmic frequencies. To date, these approaches have only been performed on a single plant species (Arabidopsis thaliana). Here, we apply a high-fidelity technique (Duplex Sequencing) to multiple angiosperms from the genus Silene, which exhibits extreme heterogeneity in rates of mitochondrial sequence evolution among close relatives. Consistent with phylogenetic evidence, we found that Silene latifolia maintains low mitochondrial variant frequencies that are comparable with previous measurements in Arabidopsis. Silene noctiflora also exhibited low variant frequencies despite high levels of historical sequence divergence, which supports other lines of evidence that this species has reverted to lower mitochondrial mutation rates after a past episode of acceleration. In contrast, S. conica showed much higher variant frequencies in mitochondrial (but not in plastid) DNA, consistent with an ongoing bout of elevated mitochondrial mutation rates. Moreover, we found an altered mutational spectrum in S. conica heavily biased towards AT→GC transitions. We also observed an unusually low number of mitochondrial genome copies per cell in S. conica, potentially pointing to reduced opportunities for homologous recombination to accurately repair mismatches in this species. Overall, these results suggest that historical fluctuations in mutation rates are driving extreme variation in rates of plant mitochondrial sequence evolution.

Another blow to the conserved gene order in Annelida: Evidence from mitochondrial genomes of the calcareous tubeworm genus Hydroides

Mitochondrial genomes are frequently applied in phylogenetic and evolutionary studies across metazoans, yet they are still poorly represented in many groups of invertebrates, including annelids. Here, we report ten mitochondrial genomes from the annelid genus Hydroides (Serpulidae) and compare them with all available annelid mitogenomes. We detected all 13 protein coding genes in Hydroides spp., including the atp8 which was reported as a missing gene in the Christmas Tree worm Spirobranchus giganteus, another annelid of the family Serpulidae. All available mitochondrial genomes of Hydroides show a highly positive GC skew combined with a highly negative AT skew – a feature consistent with that found only in the mitogenome of S. giganteus. In addition, amino acid sequences of the 13 protein-coding genes showed a high genetic distance between the Hydroides clade and S. giganteus, suggesting a fast rate of mitochondrial sequence evolution in Serpulidae. The gene order of protein-coding genes within Hydroides exhibited extensive rearrangements at species level, and were different from the arrangement patterns of other annelids, including S. giganteus. Phylogenetic analyses based on protein-coding genes recovered Hydroides as a monophyletic group sister to Spirobranchus with a long branch, and sister to the fan worm Sabellidae. Yet the Serpulidae + Sabellidae clade was unexpectedly grouped with Sipuncula, suggesting that mitochondrial genomes alone are insufficient to resolve the phylogenetic position of Serpulidae within Annelida due to its high base substitution rates. Overall, our study revealed a high variability in the gene order arrangement of mitochondrial genomes within Serpulidae, provided evidence to question the conserved pattern of the mitochondrial gene order in Annelida and called for caution when applying mitochondrial genes to infer their phylogenetic relationships.

Genome-wide signatures of plastid-nuclear coevolution point to repeated perturbations of plastid proteostasis systems across angiosperms.

Nuclear and plastid (chloroplast) genomes experience different mutation rates, levels of selection, and transmission modes, yet key cellular functions depend on their coordinated interactions. Functionally related proteins often show correlated changes in rates of sequence evolution across a phylogeny [evolutionary rate covariation (ERC)], offering a means to detect previously unidentified suites of coevolving and cofunctional genes. We performed phylogenomic analyses across angiosperm diversity, scanning the nuclear genome for genes that exhibit ERC with plastid genes. As expected, the strongest hits were highly enriched for genes encoding plastid-targeted proteins, providing evidence that cytonuclear interactions affect rates of molecular evolution at genome-wide scales. Many identified nuclear genes functioned in post-transcriptional regulation and the maintenance of protein homeostasis (proteostasis), including protein translation (in both the plastid and cytosol), import, quality control, and turnover. We also identified nuclear genes that exhibit strong signatures of coevolution with the plastid genome, but their encoded proteins lack organellar-targeting annotations, making them candidates for having previously undescribed roles in plastids. In sum, our genome-wide analyses reveal that plastid-nuclear coevolution extends beyond the intimate molecular interactions within chloroplast enzyme complexes and may be driven by frequent rewiring of the machinery responsible for maintenance of plastid proteostasis in angiosperms.

Hybrid Incompatibilities and Transgressive Gene Expression Between Two Closely Related Subspecies of Drosophila.

Genome-wide assays of expression between species and their hybrids have identified genes that become either over- or underexpressed relative to the parental species (i.e., transgressive). Transgressive expression in hybrids is of interest because it highlights possible changes in gene regulation linked to hybrid dysfunction. Previous studies in Drosophila that used long-diverged species pairs with complete or nearly complete isolation (i.e., full sterility and partial inviability of hybrids) and high-levels of genome misregulation have found correlations between expression and coding sequence divergence. The work highlighted the possible effects of directional selection driving sequence divergence and transgressive expression. Whether the same is true for taxa at early stages of divergence that have only achieved partial isolation remains untested. Here, we reanalyze previously published genome expression data and available genome sequence reads from a pair of partially isolated subspecies of Drosophila to compare expression and sequence divergence. We find a significant correlation in rates of expression and sequence evolution, but no support for directional selection driving transgressive expression in hybrids. We find that most transgressive genes in hybrids show no differential expression between parental subspecies and used SNP data to explore the role of stabilizing selection through compensatory mutations. We also examine possible misregulation through cascade effects that could be driven by interacting gene networks or co-option of off-target cis-regulatory elements.

Mutational Pressure Drives Differential Genome Conservation in Two Bacterial Endosymbionts of Sap-Feeding Insects

Compared with free-living bacteria, endosymbionts of sap-feeding insects have tiny and rapidly evolving genomes. Increased genetic drift, high mutation rates, and relaxed selection associated with host control of key cellular functions all likely contribute to genome decay. Phylogenetic comparisons have revealed massive variation in endosymbiont evolutionary rate, but such methods make it difficult to partition the effects of mutation versus selection. For example, the ancestor of Auchenorrhynchan insects contained two obligate endosymbionts, Sulcia and a betaproteobacterium (BetaSymb; called Nasuia in leafhoppers) that exhibit divergent rates of sequence evolution and different propensities for loss and replacement in the ensuing ∼300 Ma. Here, we use the auchenorrhynchan leafhopper Macrosteles sp. nr. severini, which retains both of the ancestral endosymbionts, to test the hypothesis that differences in evolutionary rate are driven by differential mutagenesis. We used a high-fidelity technique known as duplex sequencing to measure and compare low-frequency variants in each endosymbiont. Our direct detection of de novode novo mutations reveals that the rapidly evolving endosymbiont (Nasuia) has a much higher frequency of single-nucleotide variants than the more stable endosymbiont (Sulcia) and a mutation spectrum that is potentially even more AT-biased than implied by the 83.1% AT content of its genome. We show that indels are common in both endosymbionts but differ substantially in length and distribution around repetitive regions. Our results suggest that differences in long-term rates of sequence evolution in Sulcia versus BetaSymb, and perhaps the contrasting degrees of stability of their relationships with the host, are driven by differences in mutagenesis.

MSH1 is required for maintenance of the low mutation rates in plant mitochondrial and plastid genomes

Mitochondrial and plastid genomes in land plants exhibit some of the slowest rates of sequence evolution observed in any eukaryotic genome, suggesting an exceptional ability to prevent or correct mutations. However, the mechanisms responsible for this extreme fidelity remain unclear. We tested seven candidate genes involved in cytoplasmic DNA replication, recombination, and repair (POLIA, POLIB, MSH1, RECA3, UNG, FPG, and OGG1) for effects on mutation rates in the model angiosperm Arabidopsis thaliana by applying a highly accurate DNA sequencing technique (duplex sequencing) that can detect newly arisen mitochondrial and plastid mutations even at low heteroplasmic frequencies. We find that disrupting MSH1 (but not the other candidate genes) leads to massive increases in the frequency of point mutations and small indels and changes to the mutation spectrum in mitochondrial and plastid DNA. We also used droplet digital PCR to show transmission of de novo heteroplasmies across generations in msh1 mutants, confirming a contribution to heritable mutation rates. This dual-targeted gene is part of an enigmatic lineage within the mutS mismatch repair family that we find is also present outside of green plants in multiple eukaryotic groups (stramenopiles, alveolates, haptophytes, and cryptomonads), as well as certain bacteria and viruses. MSH1 has previously been shown to limit ectopic recombination in plant cytoplasmic genomes. Our results point to a broader role in recognition and correction of errors in plant mitochondrial and plastid DNA sequence, leading to greatly suppressed mutation rates perhaps via initiation of double-stranded breaks and repair pathways based on faithful homologous recombination.

Heterogeneity in the rate of molecular sequence evolution substantially impacts the accuracy of detecting shifts in diversification rates.

As species richness varies along the tree of life, there is a great interest in identifying factors that affect the rates by which lineages speciate or go extinct. To this end, theoretical biologists have developed a suite of phylogenetic comparative methods that aim to identify where shifts in diversification rates had occurred along a phylogeny and whether they are associated with some traits. Using these methods, numerous studies have predicted that speciation and extinction rates vary across the tree of life. In this study, we show that asymmetric rates of sequence evolution lead to systematic biases in the inferred phylogeny, which in turn lead to erroneous inferences regarding lineage diversification patterns. The results demonstrate that as the asymmetry in sequence evolution rates increases, so does the tendency to select more complicated models that include the possibility of diversification rate shifts. These results thus suggest that any inference regarding shifts in diversification pattern should be treated with great caution, at least until any biases regarding the molecular substitution rate have been ruled out.

Variation in the rates of evolution of the insulin and glucagon hormone and receptor genes in rodents

Insulin and glucagon are important hormones for regulating blood glucose levels. Rodents are useful models for understanding human physiology, however, differences exist between rodents and humans. Here I examined the evolution of the genes encoding insulin (Ins) and glucagon (Gcg, which also encodes GLP-1 and GLP-2) and the receptors for these hormones (Insr, Gcgr, Glp1r, and Glp2r). Our results show that the insulin 1 gene (Ins1) that originated by retroposition in some rodents such as mice, experienced selective constraints that are as strong as those acting upon the Ins2 gene found in the locus-of-origin. Previous studies had shown that the insulin hormones and genes in hystricomorph rodents, such as the guinea pig, have altered function and selective constraints, respectively. Here I show that the insulin receptor genes in hystricomorph rodents also experienced changes in evolutionary rates, but that these changes did not alter sites involved in hormone binding. While glucagon, but not GLP-1 and GLP-2, in hystricomorph rodents also show increased rates of sequence evolution, no changes in the evolution of the glucagon receptor gene (Gcgr) was seen. Intriguingly, the GLP2 receptor gene (Glp2r) in mice-like rodents evolved more rapidly than those in hystricomorph rodents. When the rates of evolution of the genes encoding the receptors for proglucagon-derived peptides, which are all G-protein coupled receptors, were compared, the GLP-1 receptor gene (Glp1r) was found to display increased levels of sequence constraint compared to the Gcgr and Glp2r genes.

Local adaptation drives the diversification of effectors in the fungal wheat pathogen Parastagonospora nodorum in the United States.

Filamentous fungi rapidly evolve in response to environmental selection pressures in part due to their genomic plasticity. Parastagonospora nodorum, a fungal pathogen of wheat and causal agent of septoria nodorum blotch, responds to selection pressure exerted by its host, influencing the gain, loss, or functional diversification of virulence determinants, known as effector genes. Whole genome resequencing of 197 P. nodorum isolates collected from spring, durum, and winter wheat production regions of the United States enabled the examination of effector diversity and genomic regions under selection specific to geographically discrete populations. 1,026,859 SNPs/InDels were used to identify novel loci, as well as SnToxA and SnTox3 as factors in disease. Genes displaying presence/absence variation, predicted effector genes, and genes localized on an accessory chromosome had significantly higher pN/pS ratios, indicating a higher rate of sequence evolution. Population structure analyses indicated two P. nodorum populations corresponding to the Upper Midwest (Population 1) and Southern/Eastern United States (Population 2). Prevalence of SnToxA varied greatly between the two populations which correlated with presence of the host sensitivity gene Tsn1 in the most prevalent cultivars in the corresponding regions. Additionally, 12 and 5 candidate effector genes were observed to be under diversifying selection among isolates from Population 1 and 2, respectively, but under purifying selection or neutrally evolving in the opposite population. Selective sweep analysis revealed 10 and 19 regions that had recently undergone positive selection in Population 1 and 2, respectively, involving 92 genes in total. When comparing genes with and without presence/absence variation, those genes exhibiting this variation were significantly closer to transposable elements. Taken together, these results indicate that P. nodorum is rapidly adapting to distinct selection pressures unique to spring and winter wheat production regions by rapid adaptive evolution and various routes of genomic diversification, potentially facilitated through transposable element activity.

Recent and old duplications in crustaceans "Internal Transcribed Spacer 1″: structural and phylogenetic implications.

Internal Transcribed Spacer structures are important in preserving accessibility to specific enzymes for the maturation of rRNAs. ITS1 sequences reported in the literature in Crustaceans range between 182and 820bp and are characterized by the absence of repeats or the presence of only a limited number of microsatellites. Here, we sequenced ITS1 for a range of shrimp families (infraorder Caridea) and show that most taxa have much larger ITS1 sequences. We find a high number of microsatellites in Alpheushebes and Crangoncrangon and we report repeat units in Pandalidae, Palaemonidae and mainly in Alpheidae species. Up to four repeats were found in A.vanderbilti (1915bp), A.rostratus (1635bp) and A.lottini (1625bp). In general, four helices were found in ITS1. Repeat units led to extra hairpins and loops. No conserved positions occurred except in helix 4. Three clades were defined in A. lottini for the first time. We estimated the ITS1 divergence rate for the three clades of A.lottini collected in French Polynesia using existing calibrations of substitution rates. Rates of sequence evolution are largely influenced by repeat units, which likely evolve separately. By comparison with COI marker, we estimated the divergence rate of the whole ITS1 sequence to range from 0.5 to 1.4% Pmy and between 0.12 and 0.5% for the 3' end of ITS1 located outside the repeat units. Given the degree of identity between repeats, we suggest that a duplication event recently occurred in A. floridanus (98% identity) whereas an ancient duplication happened in A.sulcatus (50% identity) early at the origination of the group Alpheidae, approximately 50mya ago. In conclusion, our results highlight an over representation of shorter ITS1 sequences in public repositories, and underlines the importance to further understand patterns of molecular evolution of this functionally important gene.

Mitochondrial genomes of the bird genus Piranga: rates of sequence evolution, and discordance between mitochondrial and nuclear markers

We report the characteristics of the mitochondrial genomes of 22 individuals in the bird genus Piranga, including all currently recognized species in the genus (n = 11). Elements follow the standard avian mitogenome series, including two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, 13 protein coding genes, and the mitochondrial control region. Excluding tRNA sequences, sequence divergence rate was lowest in rRNA genes and highest in genes encoding NADH (specifically ND1, ND2, ND3) and the control region. Gene trees assembled from 16 elements (non-tRNAs) varied greatly in topological concordance compared to the recognized species tree (based on thousands of nuclear loci), with no one gene tree precisely recovering all relationships in the genus. We also investigated patterns of concordance between the mitogenome tree and the nuclear species tree and found some discrepancies. Across non-tRNA gene trees (n = 16), the species tree topology was recovered by as few as three elements at a particular node and complete concordance (i.e. 16/16 gene trees matched the species tree topology) was recovered at only one node. We found mitochondrial gene regions that are often used in vertebrate systematics (e.g. CytB, ND2) recovered nearly the exact same topology as the nuclear species tree topology.

Parallels between natural selection in the cold-adapted crop-wild relative Tripsacum dactyloides and artificial selection in temperate adapted maize.

Artificial selection has produced varieties of domesticated maize that thrive in temperate climates around the world. However, the direct progenitor of maize, teosinte, is indigenous only to a relatively small range of tropical and subtropical latitudes and grows poorly or not at all outside of this region. Tripsacum, a sister genus to maize and teosinte, is naturally endemic to the majority of areas in the western hemisphere where maize is cultivated. A full-length reference transcriptome for Tripsacum dactyloides generated using long-read Iso-Seq data was used to characterize independent adaptation to temperate climates in this clade. Genes related to phospholipid biosynthesis, a critical component of cold acclimation in other cold-adapted plant lineages, were enriched among those genes experiencing more rapid rates of protein sequence evolution in T.dactyloides. In contrast with previous studies of parallel selection, we find that there is a significant overlap between the genes that were targets of artificial selection during the adaptation of maize to temperate climates and those that were targets of natural selection in temperate-adapted T.dactyloides. Genes related to growth, development, response to stimulus, signaling, and organelles were enriched in the set of genes identified as both targets of natural and artificial selection.