Search for: All records

Abstract Mitochondrial RNA editing has evolved independently in numerous eukaryotic lineages, where it generally restores conserved sequences and functional reading frames in mRNA transcripts derived from altered or disrupted mitochondrial protein-coding genes. In contrast to this “restorative” RNA editing in mitochondria, most editing of nuclear mRNAs introduces novel sequence variants and diversifies the proteome. This Perspective addresses the hypothesis that these completely opposite effects of mitochondrial vs. nuclear RNA editing arise from the enormous difference in gene number between the respective genomes. Because mitochondria produce a much smaller transcriptome, they likely create less opportunity for off-target editing, which has been supported by recent experimental work expressing mitochondrial RNA editing machinery in foreign contexts. In addition, there is recent evidence that the size and complexity of RNA targets may slow the kinetics and reduce efficiency of on-target RNA editing. These findings suggest that efficient targeting and a low risk of off-target editing have facilitated the repeated emergence of disrupted mitochondrial genes and associated restorative RNA editing systems via (potentially non-adaptive) evolutionary pathways that are not feasible in larger nuclear transcriptomes due to lack of precision. 

Plant mitochondrial genomes (mitogenomes) experience remarkable levels of horizontal gene transfer, including the recent discovery that orchids anciently acquired DNA from fungal mitogenomes. Thus far, however, there is no evidence that any of the genes from this interkingdom horizontal gene transfer are functional in orchid mitogenomes. Here, we applied a specialized sequencing approach to the orchid Corallorhiza maculata and found that some fungal-derived tRNA genes in the transferred region are transcribed, post-transcriptionally modified, and aminoacylated. In contrast, all the transferred protein-coding sequences appear to be pseudogenes. These findings show that fungal horizontal gene transfer has altered the composition of the orchid mitochondrial tRNA pool and suggest that these foreign tRNAs function in translation. The exceptional capacity of tRNAs for horizontal gene transfer and functional replacement is further illustrated by the diversity of tRNA genes in the C. maculata mitogenome, which also include genes of plastid and bacterial origin in addition to their native mitochondrial counterparts. 

SUMMARY The function of transfer RNAs (tRNAs) depends on enzymes that cleave primary transcript ends, add a 3′ CCA tail, introduce post‐transcriptional base modifications, and charge (aminoacylate) mature tRNAs with the correct amino acid. Maintaining an available pool of the resulting aminoacylated tRNAs is essential for protein synthesis. High‐throughput sequencing techniques have recently been developed to provide a comprehensive view of aminoacylation state in a tRNA‐specific fashion. However, these methods have never been applied to plants. Here, we treatedArabidopsis thalianaRNA samples with periodate and then performed tRNA‐seq to distinguish between aminoacylated and uncharged tRNAs. This approach successfully captured every tRNA isodecoder family and detected expression of additional tRNA‐like transcripts. We found that estimated aminoacylation rates and CCA tail integrity were significantly higher on average for organellar (mitochondrial and plastid) tRNAs than for nuclear/cytosolic tRNAs. Reanalysis of previously published human cell line data showed a similar pattern. Base modifications result in nucleotide misincorporations and truncations during reverse transcription, which we quantified and used to test for relationships with aminoacylation levels. We also determined that the Arabidopsis tRNA‐like sequences (t‐elements) that are cleaved from the ends of some mitochondrial messenger RNAs have post‐transcriptionally modified bases and CCA‐tail addition. However, these t‐elements are not aminoacylated, indicating that they are only recognized by a subset of tRNA‐interacting enzymes and do not play a role in translation. Overall, this work provides a characterization of the baseline landscape of plant tRNA aminoacylation rates and demonstrates an approach for investigating environmental and genetic perturbations to plant translation machinery. 

Assigning gene function from genome sequences is a rate-limiting step in molecular biology research. A protein's position within an interaction network can potentially provide insights into its molecular mechanisms. Phylogenetic analysis of evolutionary rate covariation (ERC) in protein sequence has been shown to be effective for large-scale prediction of functional relationships and interactions. However, gene duplication, gene loss, and other sources of phylogenetic incongruence are barriers for analyzing ERC on a genome-wide basis. Here, we developed ERCnet, a bioinformatic program designed to overcome these challenges, facilitating efficient all-versus-all ERC analyses for large protein sequence datasets. We simulated proteome datasets and found that ERCnet achieves combined false positive and negative error rates well below 10% and that our novel “branch-by-branch” length measurements outperforms “root-to-tip” approaches in most cases, offering a valuable new strategy for performing ERC. We also compiled a sample set of 35 angiosperm genomes to test the performance of ERCnet on empirical data, including its sensitivity to user-defined analysis parameters such as input dataset size and branch-length measurement strategy. We investigated the overlap between ERCnet runs with different species samples to understand how species number and composition affect predicted interactions and to identify the protein sets that consistently exhibit ERC across angiosperms. Our systematic exploration of the performance of ERCnet provides a roadmap for design of future ERC analyses to predict functional interactions in a wide array of genomic datasets. ERCnet code is freely available at https://github.com/EvanForsythe/ERCnet. 

Abstract PremiseA complicating factor in analyzing allopolyploid genomes is the possibility of physical interactions between homoeologous chromosomes during meiosis, resulting in either crossover (homoeologous exchanges) or non‐crossover products (homoeologous gene conversion). Homoeologous gene conversion was first described in cotton by comparing SNP patterns in sequences from two diploid progenitors with those from the allopolyploid subgenomes. These analyses, however, did not explicitly consider other evolutionary scenarios that may give rise to similar SNP patterns as homoeologous gene conversion, creating uncertainties about the reality of the inferred gene conversion events. MethodsHere, we use an expanded phylogenetic sampling of high‐quality genome assemblies from seven allopolyploidGossypiumspecies (all derived from the same polyploidy event), four diploid species (two closely related to each subgenome), and a diploid outgroup to derive a robust method for identifying potential genomic regions of gene conversion and homoeologous exchange. ResultsWe found little evidence for homoeologous gene conversion in allopolyploid cottons, and that only two of the 40 best‐supported events were shared by more than one species. We did, however, reveal a single, shared homoeologous exchange event at one end of chromosome 1, which occurred shortly after allopolyploidization but prior to divergence of the descendant species. ConclusionsOverall, our analyses demonstrated that homoeologous gene conversion and homoeologous exchanges are uncommon inGossypium, affecting between zero and 24 genes per subgenome (0.0–0.065%) across the seven species. More generally, we highlighted the potential problems of using simple four‐taxon tests to investigate patterns of homoeologous gene conversion in established allopolyploids. 

Eukaryotic nuclear genomes often encode distinct sets of translation machinery for function in the cytosol vs. organelles (mitochondria and plastids). This raises questions about why multiple translation systems are maintained even though they are capable of comparable functions and whether they evolve differently depending on the compartment where they operate. These questions are particularly interesting in plants because translation machinery, including aminoacyl-transfer RNA (tRNA) synthetases (aaRS), is often dual-targeted to the plastids and mitochondria. These organelles have different functions, with much higher rates of translation in plastids to supply the abundant, rapid-turnover proteins required for photosynthesis. Previous studies have indicated that plant organellar aaRS evolve more slowly compared to mitochondrial aaRS in eukaryotes that lack plastids. Thus, we investigated the evolution of nuclear-encoded organellar and cytosolic aaRS and tRNA maturation enzymes across a broad sampling of angiosperms, including nonphotosynthetic (heterotrophic) plant species with reduced plastid gene expression, to test the hypothesis that translational demands associated with photosynthesis constrain the evolution of enzymes involved in organellar tRNA metabolism. Remarkably, heterotrophic plants exhibited wholesale loss of many organelle-targeted aaRS and other enzymes, even though translation still occurs in their mitochondria and plastids. These losses were often accompanied by apparent retargeting of cytosolic enzymes and tRNAs to the organelles, sometimes preserving aaRS–tRNA charging relationships but other times creating surprising mismatches between cytosolic aaRS and mitochondrial tRNA substrates. Our findings indicate that the presence of a photosynthetic plastid drives the retention of specialized systems for organellar tRNA metabolism. 

Abstract Hybridization in plants is often accompanied by nuclear genome doubling (allopolyploidy), which has been hypothesized to perturb interactions between nuclear and organellar (mitochondrial and plastid) genomes by creating imbalances in the relative copy number of these genomes and producing genetic incompatibilities between maternally derived organellar genomes and the half of the allopolyploid nuclear genome from the paternal progenitor. Several evolutionary responses have been predicted to ameliorate these effects, including selection for changes in protein sequences that restore cytonuclear interactions; biased gene retention/expression/conversion favoring maternal nuclear gene copies; and fine-tuning of relative cytonuclear genome copy numbers and expression levels. Numerous recent studies, however, have found that evolutionary responses are inconsistent and rarely scale to genome-wide generalities. The apparent robustness of plant cytonuclear interactions to allopolyploidy may reflect features that are general to allopolyploids such as the lack of F2 hybrid breakdown under disomic inheritance, and others that are more plant-specific, including slow sequence divergence in organellar genomes and pre-existing regulatory responses to changes in cell size and endopolyploidy during development. Thus, cytonuclear interactions may only rarely act as the main barrier to establishment of allopolyploid lineages, perhaps helping to explain why allopolyploidy is so pervasive in plant evolution. 

Abstract The angiosperm genus Silene has been the subject of extensive study in the field of ecology and evolution, but the availability of high-quality reference genome sequences has been limited for this group. Here, we report a chromosome-level assembly for the genome of Silene conica based on Pacific Bioscience HiFi, Hi-C, and Bionano technologies. The assembly produced 10 scaffolds (1 per chromosome) with a total length of 862 Mb and only ∼1% gap content. These results confirm previous observations that S. conica and its relatives have a reduced base chromosome number relative to the genus's ancestral state of 12. Silene conica has an exceptionally large mitochondrial genome (>11 Mb), predominantly consisting of sequence of unknown origins. Analysis of shared sequence content suggests that it is unlikely that transfer of nuclear DNA is the primary driver of this mitochondrial genome expansion. More generally, this assembly should provide a valuable resource for future genomic studies in Silene, including comparative analyses with related species that recently evolved sex chromosomes. 

Abstract Eukaryotes maintain separate protein translation systems for nuclear and organellar genes, including distinct sets of tRNAs and aminoacyl-tRNA synthetases (aaRSs). In animals, mitochondrial-targeted aaRSs are expressed at lower levels and are less conserved in sequence than cytosolic aaRSs involved in translation of nuclear mRNAs, likely reflecting lower translational demands in mitochondria. In plants, translation is further complicated by the presence of plastids, which share most aaRSs with mitochondria. In addition, plant mitochondrial tRNA pools have a dynamic history of gene loss and functional replacement by tRNAs from other compartments. To investigate the consequences of these distinctive features of translation in plants, we analyzed sequence evolution in angiosperm aaRSs. In contrast to previously studied eukaryotic systems, we found that plant organellar and cytosolic aaRSs exhibit only a small difference in expression levels, and organellar aaRSs are slightly more conserved than cytosolic aaRSs. We hypothesize that these patterns result from high translational demands associated with photosynthesis in mature chloroplasts. We also investigated aaRS evolution in Sileneae, an angiosperm lineage with extensive mitochondrial tRNA replacement and aaRS retargeting. We predicted positive selection for changes in aaRS sequence resulting from these recent changes in subcellular localization and tRNA substrates but found little evidence for accelerated sequence divergence. Overall, the complex tripartite translation system in plant cells appears to have imposed more constraints on the long-term evolutionary rates of organellar aaRSs compared with other eukaryotic lineages, and plant aaRS protein sequences appear largely robust to more recent perturbations in subcellular localization and tRNA interactions.