Tetraploid interspecific hybrids between Asian and African rice species restore fertility depending on killer-protector loci for hybrid sterility.

Interspecific F1 hybrids between Asian (Oryza sativa) and African rice (Oryza glaberrima) exhibit severe sterility caused by the accumulation of hybrid sterility genes/loci at 15 or more loci. The mechanisms underlying the hybrid sterility genes are largely unknown; however, a few genes associated with the killer-protector system, which is the system most frequently associated with hybrid sterility genes, have been identified. We previously produced fertile plants as tetraploids derived from diploid interspecific F1 hybrids through anther culture; therefore, it was suggested that hybrid sterility could be overcome following tetraploidization. We investigated whether tetraploid interspecific plants produced by crossing are fertile and tested the involvement of hybrid sterility genes in the process. Fertile tetraploid interspecific F1 hybrid plants were obtained by crossing two tetraploids of Oryza sativa and Oryza glaberrima. To elucidate the relationships between pollen fertility and the hybrid sterility loci in the tetraploid F1 microspores, we performed genetic analyses of the tetraploid F2 hybrids and diploid plants obtained from the microspores of tetraploid interspecific hybrids by anther culture. The result suggested that the tetraploid interspecific hybrids overcame pollen and seed infertility based on the proportion of loci with the killer-protector system present in the tetraploids. The heterozygous hybrid sterility loci with the killer-protector system in the tetraploid segregate the homozygous killed allele (16.7-21.4%), with more than three-quarters of the gametes surviving. We theoretically and experimentally demonstrated that fertile rice progenies can be grown from tetraploid interspecific hybrids.

OlCHR, encoding a chromatin remodeling factor, is a killer causing hybrid sterility between rice species Oryza sativa and O. longistaminata.

The genetic mechanisms of reproductive isolation have been widely investigated within Asian cultivated rice (Oryza sativa); however, relevant genes between diverged species have been in sighted rather less. Herein, a gene showing selfish behavior was discovered in hybrids between the distantly related rice species Oryza longistaminata and O. sativa. The selfish allele S13l in the S13 locus impaired male fertility, discriminately eliminating pollens containing the allele S13s from O. sativa in heterozygotes (S13s/S13l). Genetic analysis revealed that a gene encoding a chromatin-remodeling factor (CHR) is involved in this phenomenon and a variety of O. sativa owns the truncated gene OsCHR745, whereas its homologue OlCHR has a complete structure in O. longistaminata. CRISPR-Cas9-mediated loss of function mutants restored fertility in hybrids. African cultivated rice, which naturally lacks the OlCHR homologue, is compatible with both S13s and S13l carriers. These results suggest that OlCHR is a Killer gene, which leads to reproductive isolation.

Different contributions of PROG1 and TAC1 to the angular kinematics of the main culm and tillers of wild rice (Oryza rufipogon).

MAIN CONCLUSION: PROG1 is necessary but insufficient for the main culm inclination while TAC1 partially takes part in it, and both genes promote tiller inclination in Asian wild rice. Asian wild rice (Oryza rufipogon), the ancestor of cultivated rice (O. sativa), has a prostrate architecture, with tillers branching from near the ground. The main culm of each plant grows upward and then tilts during the vegetative stage. Genes controlling tiller angle have been reported; however, their genetic contributions to the culm movement have not been quantified. Here, we quantified their genetic contributions to angular kinematics in the main culm and tillers. For the main culm inclination, one major QTL surrounding the PROG1 region was found. In cultivated rice, tillers firstly inclined and lately rose, while it kept inclining in wild rice. It was suggested that PROG1 affected the tiller elevation angle in the later kinematics, whereas TAC1 was weakly associated with the tiller angle in the whole vegetative stage. Micro-computed tomography (micro-CT) suggested that these angular changes are produced by the bending of culm bases. Because near-isogenic lines (NILs) of wild rice-type Prog1 and Tac1 alleles in the genetic background of cultivated rice did not show the prostrate architecture, the involvement of another gene(s) for inclination of the main culm was suggested. Our findings will not only contribute to the understanding of the morphological transition during domestication but also be used in plant breeding to precisely reproduce the ideal plant architecture by combining the effects of multiple genes.

Phylogenetic analysis of endogenous viral elements in the rice genome reveals local chromosomal evolution in Oryza AA-genome species.

Introduction: Rice genomes contain endogenous viral elements homologous to rice tungro bacilliform virus (RTBV) from the pararetrovirus family Caulimoviridae. These viral elements, known as endogenous RTBV-like sequences (eRTBVLs), comprise five subfamilies, eRTBVL-A, -B, -C, -D, and -X. Four subfamilies (A, B, C, and X) are present to a limited degree in the genomes of the Asian cultivated rice Oryza sativa (spp. japonica and indica) and the closely related wild species Oryza rufipogon. Methods: The eRTBVL-D sequences are widely distributed within these and other Oryza AA-genome species. Fifteen eRTBVL-D segments identified in the japonica (Nipponbare) genome occur mostly at orthologous chromosomal positions in other AA-genome species. The eRTBVL-D sequences were inserted into the genomes just before speciation of the AA-genome species. Results and discussion: Ten eRTBVL-D segments are located at six loci, which were used for our evolutionary analyses during the speciation of the AA-genome species. The degree of genetic differentiation varied among the eRTBVL-D segments. Of the six loci, three showed phylogenetic trees consistent with the standard speciation pattern (SSP) of the AA-genome species (Type A), and the other three represented phylogenies different from the SSP (Type B). The atypical phylogenetic trees for the Type B loci revealed chromosome region-specific evolution among the AA-genome species that is associated with phylogenetic incongruences: complex genome rearrangements between eRTBVL-D segments, an introgression between the distant species, and low genetic diversity of a shared eRTBVL-D segment. Using eRTBVL-D as an indicator, this study revealed the phylogenetic incongruence of local chromosomal regions with different topologies that developed during speciation.

Genetics of chilling response at early growth stage in rice: a recessive gene for tolerance and importance of acclimation.

Low-temperature adaptation in rice is mediated by the ability of a genotype to tolerate chilling temperatures. A genetic locus on chromosome 11 was analysed for chilling tolerance at the plumule stage in rice. The tolerant allele of A58, a japonica landrace in Japan, was inherited as a recessive gene (ctp-1A58), whereas the susceptible alleles from wild rice (Ctp-1W107) and modern variety (Ctp-1HY) were the dominant genes. Another recessive tolerant allele (ctp-1Silewah) was found in a tropical japonica variety (Silewah). Fine-mapping revealed that a candidate gene for the ctp-1 locus encoded a protein similar to the nucleotide-binding domain and leucine-rich repeat (NLR) protein, in which frameshift mutation by a 73 bp-deletion might confer chilling tolerance in ctp-1A58. Analysis of near-isogenic lines demonstrated that ctp-1A58 imparted tolerance effects only at severe chilling temperatures of 0.5 °C and 2 °C, both at plumule and seedling stages. Chilling acclimation treatments at a wide range of temperatures (8 °C-16 °C) for 72 h concealed the susceptible phenotype of Ctp-1W107 and Ctp-1HY. Furthermore, short-term acclimation treatment of 12 h at 8 °C was enough to be fully acclimated. These results suggest that the NLR gene induces a susceptible response upon exposure to severe chilling stress, however, another interacting gene(s) for acclimation response could suppress the maladaptive phenotype caused by the Ctp-1 allele. This study provides new insights for the adaptation and breeding of rice in a low-temperature environment.

Identification of a Saltol-Independent Salinity Tolerance Polymorphism in Rice Mekong Delta Landraces and Characterization of a Promising Line, Doc Phung.

The Mekong Delta River in Vietnam is facing salinity intrusion caused by climate change and sea-level rise that is severely affecting rice cultivation. Here, we evaluated salinity responses of 97 rice accessions (79 landraces and 18 improved accessions) from the Mekong Delta population by adding 100 mM NaCl to the nutrient solution for up to 20 days. We observed a wide distribution in salinity tolerance/sensitivity, with two major peaks across the 97 accessions when using the standard evaluation system (SES) developed by the International Rice Research Institute. SES scores revealed strong negative correlations (ranging from - 0.68 to - 0.83) with other phenotypic indices, such as shoot elongation length, root elongation length, shoot dry weight, and root dry weight. Mineral concentrations of Na+ in roots, stems, and leaves and Ca2+ in roots and stems were positively correlated with SES scores, suggesting that tolerant accessions lower their cation exchange capacity in the root cell wall. The salinity tolerance of Mekong Delta accessions was independent from the previously described salinity tolerance-related locus Saltol, which encodes an HKT1-type transporter in the salinity-tolerant cultivars Nona Bokra and Pokkali. Indeed, genome-wide association studies using SES scores and shoot dry weight ratios of the 79 accessions as traits identified a single common peak located on chromosome 1. This SNP did not form a linkage group with other nearby SNPs and mapped to the 3' untranslated region of gene LOC_Os01g32830, over 6.5 Mb away from the Saltol locus. LOC_Os01g32830 encodes chloroplast glycolate/glycerate translocator 1 (OsPLGG1), which is responsible for photorespiration and growth. SES and shoot dry weight ratios differed significantly between the two possible haplotypes at the causal SNP. Through these analyses, we characterize Doc Phung, one of the most salinity-tolerant varieties in the Mekong Delta population and a promising new genetic resource.

How to establish a mutually beneficial relationship between a transposon and its host: lessons from Tam3 in Antirrhinum.

The transposon Tam3 of Antirrhinum (snapdragon) has acquired properties that distinguish it from other transposons. Mobile DNA, commonly referred to as a transposable element or transposon, is considered to be synonymous with a selfish factor. That is, a transposable element increases in copy number and moves copies of itself independently of the survival of the host organism. Therefore, the host collectively regulates the transposition activities of most transposable elements in its genome by epigenetic means. However, our analyses of the structure and behavior of Tam3, as shown by the following five results, provide evidence that it does not behave in a selfish manner in relation to the host. 1) Active transposable elements normally increase the abundance of their non-autonomous elements, whereas Tam3 is known to have no non-autonomous elements, and a limited number of around 10 copies of autonomous elements present in the genome have been isolated as active copies. 2) Tam3 does not transpose at 25 ℃, which is the optimal growth temperature for Antirrhinum. Transposition of Tam3 occurs only at low temperatures of about 15 ℃, which is stressful for Antirrhinum. 3) Few strains of Antirrhinum have been found to contain genes that specifically suppress Tam3 transposition. 4) Most of the Tam3 insertions found in Antirrhinum genes do not affect the host genome, and the expression of these host genes is not completely suppressed. 5) Transcription and translation of the Tam3 transposase gene are not epigenetically regulated by the host. These five experimental results constitute evidence that Tam3 retains features that are dissimilar to those of many other transposons and that it does not behave in a selfish manner that is detrimental to the survival of the host. In this review, we consider what kinds of behavior are required if transposons are to establish a mutually beneficial relationship with their hosts, with reference to Tam3.

Temporal changes in transcripts of miniature inverted-repeat transposable elements during rice endosperm development.

The repression of transcription from transposable elements (TEs) by DNA methylation is necessary to maintain genome integrity and prevent harmful mutations. However, under certain circumstances, TEs may escape from the host defense system and reactivate their transcription. In Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), DNA demethylases target the sequences derived from TEs in the central cell, the progenitor cell for the endosperm in the female gametophyte. Genome-wide DNA demethylation is also observed in the endosperm after fertilization. In the present study, we used a custom microarray to survey the transcripts generated from TEs during rice endosperm development and at selected time points in the embryo as a control. The expression patterns of TE transcripts are dynamically up- and downregulated during endosperm development, especially those of miniature inverted-repeat TEs (MITEs). Some TE transcripts were directionally controlled, whereas the other DNA transposons and retrotransposons were not. We also discovered the NUCLEAR FACTOR Y binding motif, CCAAT, in the region near the 5' terminal inverted repeat of Youren, one of the transcribed MITEs in the endosperm. Our results uncover dynamic changes in TE activity during endosperm development in rice.

Seed management using NGS technology to rapidly eliminate a deleterious allele from rice breeder seeds.

Spontaneous mutations are stochastic phenomena that occur in every population. However, deleterious mutated allele present in seeds distributed to farmers must be detected and removed. Here, we eliminated undesirable mutations from the parent population in one generation through a strategy based on next-generation sequencing (NGS). This study dealt with a spontaneous albino mutant in the 'Hinohikari' rice variety grown at the Miyazaki Comprehensive Agricultural Experiment Station, Japan. The incidence of albinism in the population was 1.36%. NGS analysis revealed the genomic basis for differences between green and albino phenotypes. Every albino plant had a C insertion in the Snow-White Leaf1 (SWL1) gene on chromosome 4 causing a frameshift mutation. Selfing plants heterozygous for the mutant allele, swl1-R332P, resulted in a 3:1 green/albino ratio, confirming that a single recessive gene controls albinism. Ultrastructural leaf features in the swl1-R332P mutants displayed deformed chlorophyll-associated organelles in albino plants that were similar to those of previously described swl1 mutants. Detection of the causative gene and its confirmation using heterozygous progenies were completed within a year. The NGS technique outlined here facilitates rapid identification of spontaneous mutations that can occur in breeder seeds.

Effectiveness of Create ML in microscopy image classifications: a simple and inexpensive deep learning pipeline for non-data scientists.

Observing chromosomes is a time-consuming and labor-intensive process, and chromosomes have been analyzed manually for many years. In the last decade, automated acquisition systems for microscopic images have advanced dramatically due to advances in their controlling computer systems, and nowadays, it is possible to automatically acquire sets of tiling-images consisting of large number, more than 1000, of images from large areas of specimens. However, there has been no simple and inexpensive system to efficiently select images containing mitotic cells among these images. In this paper, a classification system of chromosomal images by deep learning artificial intelligence (AI) that can be easily handled by non-data scientists was applied. With this system, models suitable for our own samples could be easily built on a Macintosh computer with Create ML. As examples, models constructed by learning using chromosome images derived from various plant species were able to classify images containing mitotic cells among samples from plant species not used for learning in addition to samples from the species used. The system also worked for cells in tissue sections and tetrads. Since this system is inexpensive and can be easily trained via deep learning using scientists' own samples, it can be used not only for chromosomal image analysis but also for analysis of other biology-related images.

Revision of the relationship between anther morphology and pollen sterility by cold stress at the booting stage in rice.

BACKGROUND AND AIMS: Cold stress in rice (Oryza sativa) plants at the reproductive stage prevents normal anther development and causes pollen sterility. Tapetum hypertrophy in anthers has been associated with pollen sterility in response to cold at the booting stage. Here, we re-examined whether the relationships between anther abnormality and pollen sterility caused by cold stress at the booting stage in rice can be explained by a monovalent factor such as tapetum hypertrophy. METHODS: After exposing plants to a 4-d cold treatment at the booting stage, we collected and processed anthers for transverse sectioning immediately and at the flowering stage. We anatomically evaluated the effect of cold treatment on anther internal morphologies, pollen fertilities and pollen numbers in the 13 cultivars with various cold sensitivities. KEY RESULTS: We observed four types of morphological anther abnormalities at each stage. Pollen sterility was positively correlated with the frequency of undeveloped locules, but not with tapetum hypertrophy as commonly believed. In cold-sensitive cultivars grown at low temperatures, pollen sterility was more frequent than anther morphological abnormalities, and some lines showed remarkably high pollen sterility without any anther morphological alterations. Most morphological anomalies occurred only in specific areas within large and small locules. Anther length tended to shorten in response to cold treatment and was positively correlated with pollen number. One cultivar showed a considerably reduced pollen number, but fertile pollen grains under cold stress. We propose three possible relationships to explain anther structure and pollen sterility and reduction due to cold stress. CONCLUSIONS: The pollen sterility caused by cold stress at the booting stage was correlated with the frequency of entire locule-related abnormalities, which might represent a phenotypic consequence, but not a direct cause of pollen abortion. Multivalent factors might underlie the complicated relationships between anther abnormality and pollen sterility in rice.

Diploid Male Gametes Circumvent Hybrid Sterility Between Asian and African Rice Species.

In F1 hybrids of Oryza sativa (Asian rice) and Oryza glaberrima (African rice), heterozygosity leads to a complete gamete abortion because of allelic conflict at each of the 13 hybrid sterility (HS) loci. We systematically produced 19 plants from the F1 hybrids of both the rice species by the anther culture (AC) method. Five of the 19 interspecific hybrid plants were partially fertile and able to produce seeds. Unlike ordinal doubled haploid plants resulting from AC, these regenerated plants showed various ploidy levels (diploid to pentaploid) and different zygosities (completely homozygous, completely heterozygous, and a combination). These properties were attributable to meiotic anomalies in the interspecific hybrid F1 plants. Examination of the genetic structures of the regenerated plants suggested meiotic non-reduction took place in the interspecific hybrid F1 plants. The centromeric regions in the regenerated plants revealed that the abnormal first and/or second divisions of meiosis, namely the first division restitution (FDR) and/or second division restitution (SDR), had occurred in the interspecific hybrid. Immunohistochemical observations also verified these phenomena. FDR and SDR occurrences at meiosis might strongly lead to the formation of diploid microspores. The results demonstrated that meiotic anomalies functioned as a reproductive barrier occurred before the HS genes acted in gamete of the interspecific hybrid. Although such meiotic anomalies are detrimental to pollen development, the early rescue of microspores carrying the diploid gamete resulted in the fertile regenerated plants. The five partially fertile plants carrying tetraploid genomes with heterozygous alleles of the HS loci produced fertile diploid pollens, implying that the diploid gametes circumvented the allelic conflicts at the HS loci. We also proposed how diploid male gametes avoid HS with the killer-protector model.

Fertile Tetraploids: New Resources for Future Rice Breeding?

Ploidy manipulation is an efficient technique for the development of novel phenotypes in plant breeding. However, in rice (Oryza sativa L.), severe seed sterility has been considered a barrier preventing cultivation of autotetraploids since the 1930s. Recently, a series of studies identified two fertile autotetraploids, identified herein as the PMeS (Polyploid Meiosis Stability) and Neo-Tetraploid lines. Here, we summarize their characteristics, focusing on the recovery of seed fertility, and discuss potential future directions of study in this area, providing a comprehensive understanding of current progress in the study of fertile tetraploid rice, a classical, but promising, concept for rice breeding.

Infection with an asymptomatic virus in rice results in a delayed drought response.

Infection of viruses in plants often modifies plant responses to biotic and abiotic stresses. In the present study we examined the effects of Rice tungro spherical virus (RTSV) infection on drought response in rice. RTSV infection delayed the onset of leaf rolling by 1-2 days. During the delay in drought response, plants infected with RTSV showed higher stomatal conductance and less negative leaf water potential under drought than those of uninfected plants, indicating that RTSV-infected leaves were more hydrated. Other growth and physiological traits of plants under drought were not altered by infection with RTSV. An expression analysis of genes for drought response-related transcription factors showed that the expression of OsNAC6 and OsDREB2a was less activated by drought in RTSV-infected plants than in uninfected plants, further suggesting improved water status of the plants due to RTSV infection. RTSV accumulated more in plants under drought than in well-watered plants, indicating the increased susceptibility of rice plants to RTSV infection by drought. Collectively, these results indicated that infection with RTSV can transiently mitigate the influence of drought stress on rice plants by increasing leaf hydration, while drought increased the susceptibility of rice plants to RTSV.

Genetic Properties Responsible for the Transgressive Segregation of Days to Heading in Rice.

Transgressive segregation produces hybrid progeny phenotypes that exceed the parental phenotypes. Unlike heterosis, extreme phenotypes caused by transgressive segregation are heritably stable. We examined transgressive phenotypes of flowering time in rice, and revealed transgressive segregation in F2 populations derived from a cross between parents with similar (proximal) days to heading (DTH). The DTH phenotypes of the A58 × Kitaake F2 progenies were frequently more extreme than those of either parent. These transgressive phenotypes were maintained in the F3 and F4 populations. Both A58 and Kitaake are japonica rice cultivars adapted to Hokkaido, Japan, which is a high-latitude region, and have a short DTH. Among the four known loci required for a short DTH, three loci had common alleles in A58 and Kitaake, implying there is a similar genetic basis for DTH between the two varieties. A genome-wide single nucleotide polymorphism (SNP) analysis based on the F4 population identified five new quantitative trait loci (QTL) associated with transgressive DTH phenotypes. Each of these QTL had different degrees of additive effects on DTH, and two QTL had an epistatic effect on each other. Thus, a genome-wide SNP analysis facilitated the detection of genetic loci associated with extreme DTH phenotypes, and revealed that the transgressive phenotypes were produced by exchanging the complementary alleles of a few minor QTL in the similar parental phenotypes.

Profiling SNP and Nucleotide Diversity to Characterize Mekong Delta Rice Landraces in Southeast Asian Populations.

CORE IDEAS: Single nucleotide polymorphism (SNP) analyses are a powerful tool to examine structure of local rice population. 3000 dataset of IRRI facilitates SNP profiling of Southeast Asian rice populations. Mekong Delta population is featured by comparisons with the other populations. The low π-value SNPs well-profile unique genetic regions in their genomes. Recent analyses using single nucleotide polymorphism (SNP) are a feasible mean for local collections which potentially possess useful, but not large, genetic variations. Genomic sequences of more than 3000 accessions released by the International Rice Research Institute (IRRI) can be used to characterize various local rice (Oryza sativa) populations. The aim of this study was to develop a method to facilitate genomic characterization of local rice populations. We mainly used 99 indica rice accessions (81 landraces and 18 improved varieties) from the Mekong Delta Development Research Institute (MDI). We obtained 2301 SNPs after a genomic sequencing analysis of the 99 rice accessions and subsequent filtering. Within the IRRI's dataset, the landraces fell into a cluster consisting of accessions from Southeast Asian countries (Ind3 cluster), and the MDI improved varieties were grouped in a cluster containing IRRI improved varieties (Ind1B cluster). A principal component analysis suggested that geographical location strongly affects phylogenetic relationships, and the MDI landraces were placed into a Vietnam+Cambodia group. To detect the nucleotide diversity within a population, π-value is commonly used. We think that whole genome distribution of π-values representing the nucleotide diversity of each population can be used to characterize local populations. Our simple profiling using low π-value genomic regions was able to reveal regional characteristics of rice genomes and should be useful for identifying local rice populations.

Anther culture in rice proportionally rescues microspores according to gametophytic gene effect and enhances genetic study of hybrid sterility.

BACKGROUND: To investigate plant hybrid sterility, we studied interspecific hybrids of two cultivated rice species, Asian rice (Oryza sativa) and African rice (O. glaberrima). Male gametes of these hybrids display complete sterility owing to a dozen of hybrid sterility loci, termed HS loci, but this complicated genetic system remains poorly understood. RESULTS: Microspores from these interspecific hybrids form sterile pollen but are viable at the immature stage. Application of the anther culture (AC) method caused these immature microspores to induce callus. The segregation distortion of 11 among 13 known HS loci was assessed in the callus population. Using many individual calli, fine mapping of the HS loci was attempted based on heterozygotes produced from chromosome segment substitution lines (CSSLs). Transmission ratio distortion (TRD) from microspores was detected at 6 of 11 HS loci in the callus population. The fine mapping of S 1 and S 19 loci using CSSLs revealed precise distances of markers from the positions of HS loci exhibiting excessive TRD. CONCLUSIONS: We demonstrated that AC to generate callus populations derived from immature microspores is a useful methodology for genetic study. The callus population facilitated detection of TRD at multiple HS loci and dramatically shortened the process for mapping hybrid sterility genes.

Ancient Endogenous Pararetroviruses in Oryza Genomes Provide Insights into the Heterogeneity of Viral Gene Macroevolution.

Endogenous viral sequences in eukaryotic genomes, such as those derived from plant pararetroviruses (PRVs), can serve as genomic fossils to study viral macroevolution. Many aspects of viral evolutionary rates are heterogeneous, including substitution rate differences between genes. However, the evolutionary dynamics of this viral gene rate heterogeneity (GRH) have been rarely examined. Characterizing such GRH may help to elucidate viral adaptive evolution. In this study, based on robust phylogenetic analysis, we determined an ancient endogenous PRV group in Oryza genomes in the range of being 2.41-15.00 Myr old. We subsequently used this ancient endogenous PRV group and three younger groups to estimate the GRH of PRVs. Long-term substitution rates for the most conserved gene and a divergent gene were 2.69 × 10-8 to 8.07 × 10-8 and 4.72 × 10-8 to 1.42 × 10-7 substitutions/site/year, respectively. On the basis of a direct comparison, a long-term GRH of 1.83-fold was identified between these two genes, which is unexpectedly low and lower than the short-term GRH (>3.40-fold) of PRVs calculated using published data. The lower long-term GRH of PRVs was due to the slightly faster rate decay of divergent genes than of conserved genes during evolution. To the best of our knowledge, we quantified for the first time the long-term GRH of viral genes using paleovirological analyses, and proposed that the GRH of PRVs might be heterogeneous on time scales (time-dependent GRH). Our findings provide special insights into viral gene macroevolution and should encourage a more detailed examination of the viral GRH.

Genomic fossils reveal adaptation of non-autonomous pararetroviruses driven by concerted evolution of noncoding regulatory sequences.

The interplay of different virus species in a host cell after infection can affect the adaptation of each virus. Endogenous viral elements, such as endogenous pararetroviruses (PRVs), have arisen from vertical inheritance of viral sequences integrated into host germline genomes. As viral genomic fossils, these sequences can thus serve as valuable paleogenomic data to study the long-term evolutionary dynamics of virus-virus interactions, but they have rarely been applied for this purpose. All extant PRVs have been considered autonomous species in their parasitic life cycle in host cells. Here, we provide evidence for multiple non-autonomous PRV species with structural defects in viral activity that have frequently infected ancient grass hosts and adapted through interplay between viruses. Our paleogenomic analyses using endogenous PRVs in grass genomes revealed that these non-autonomous PRV species have participated in interplay with autonomous PRVs in a possible commensal partnership, or, alternatively, with one another in a possible mutualistic partnership. These partnerships, which have been established by the sharing of noncoding regulatory sequences (NRSs) in intergenic regions between two partner viruses, have been further maintained and altered by the sequence homogenization of NRSs between partners. Strikingly, we found that frequent region-specific recombination, rather than mutation selection, is the main causative mechanism of NRS homogenization. Our results, obtained from ancient DNA records of viruses, suggest that adaptation of PRVs has occurred by concerted evolution of NRSs between different virus species in the same host. Our findings further imply that evaluation of within-host NRS interactions within and between populations of viral pathogens may be important.