A non-synonymous SNP within the isopentenyl transferase 2 locus is associated with kernel weight in Chinese maize inbreds (Zea mays L.)

Functional characterization of transcriptional activator gene SIARRI in tomato reveals its role in fruit growth and ripening

Auxins regulate several characteristics of plant development and growth. Here, we characterized a new transcriptional activator SIARRI which binds specific DNA sequences and was revealed in Arabidopsis (ARR1). SIARRI acts as a two-component response regulator and its Arabidopsis homologous gene is AT3G16857. It belongs to the subfamily of type-B response regulators in the cytokinin signaling pathway. The study aimed to characterize the transgenic Micro-Tom plants by the overexpression of Solanum lycopersicum two-component response regulator ARR1. Overexpression of SIARRI results in a pleiotropic phenotype during fruit development and ripening. This study indicates that SIARRI is a primary regulator of leaf morphology and fruit development. Moreover, overexpressed plants showed variations in growth related to auxin as well as shorter hypocotyl elongation, enlarged leaf vascularization, and decreased apical dominance. The qRT-PCR investigation revealed that expression was downregulated at the breaker stage and high at Br+6 at various stages of fruit growth and ripening. In contrast to the fruit color, lycopene and β-carotene concentrations in red-yellow overexpression line fruits were reduced significantly, and also slightly reduced in some red fruits. The quantity of β-carotene in the transgenic fruits was lower than that of lycopene. This study showed that this gene might be a new transcriptional activator in fruit development and ripening. Furthermore, this study will provide new insights into tomato fruit ripening.

Overexpression of ZmIPT2 gene delays leaf senescence and improves grain yield in maize

Cytokinins (CTKs) are a major phytohormone group that are significant in the promotion of cellular division, growth, and divergence. Isopentenyl transferase (IPT) regulates a rate-limiting step in plant CTK synthesis, promotes the synthesis of isopentenyl adenonucleotides from 5-AMP and isopentenyl pyrophosphate, and then converts both these chemicals into various CTKs. Here, the full-length cDNA of ZmIPT2, which encodes 322 amino acids, was isolated and was introduced into a maize inbred line by Agrobacterium-mediated transformation. In both controlled environments and field experiments, the overexpression of ZmIPT2 gene in the transformed plants delayed leaf senescence. Compared to the receptor line, the transgenic maize lines retained higher chlorophyll levels, photosynthetic rates, and cytokinin content for an extended period of time, and produced significantly higher grain yield by a margin of 17.71-20.29% under normal field planting conditions. Subsequently, ten possible genes that interacted with ZmIPT2 were analyzed by qRT-PCR, showing that the expression pattern of GRMZM2G022904 was consistent with ZmIPT2 expression. Through comprehensive analysis, we screened for transgenic lines with stable inheritance of ZmIPT2 gene, clear functional efficiency, and significant yield improvement, in order to provide theoretical basis and material support for the breeding of new high-yield transgenic maize varieties.

Using Transcriptome Analysis to Explore Gray Mold Resistance-Related Genes in Onion (Allium cepa L.)

Gray mold disease caused by Botrytis in onions (Allium cepa L.) during growth and storage negatively affects their yield and quality. Exploring the genes related to gray mold resistance in onion and their application to the breeding of resistant onion lines will support effective and ecological control methods of the disease. Here, the genetic relationship of 54 onion lines based on random amplified polymorphic DNA (RAPD) and in vitro-cultured onion lines infected with gray mold were used for screening resistance and susceptibility traits. Two genetically related onion lines were selected, one with a resistant and one with a susceptible phenotype. In vitro gray mold infection was repeated with these two lines, and leaf samples were collected for gene expression studies in time series. Transcript sequences obtained by RNA sequencing were subjected to DEG analysis, variant analysis, and KEGG mapping. Among the KEGG pathways, ‘α-linoleic acid metabolism’ was selected because the comparison of the time series expression pattern of Jasmonate resistant 1 (JAR1), Coronatine-insensitive protein 1 (COI 1), and transcription factor MYC2 (MYC2) genes between the resistant and susceptible lines revealed its significant relationship with gray-mold-resistant phenotypes. Expression pattern and SNP of the selected genes were verified by quantitative real-time PCR and high-resolution melting (HRM) analysis, respectively. The results of this study will be useful for the development of molecular marker and finally breeding of gray-mold-resistant onions.

Genetic and molecular control of grain yield in maize

Understanding the genetic and molecular basis of grain yield is important for maize improvement. Here, we identified 49 consensus quantitative trait loci (cQTL) controlling maize yield-related traits using QTL meta-analysis. Then, we collected yield-related traits associated SNPs detected by association mapping and identified 17 consensus significant loci. Comparing the physical positions of cQTL with those of significant SNPs revealed that 47 significant SNPs were located within 20 cQTL regions. Furthermore, intensive reviews of 31 genes regulating maize yield-related traits found that the functions of many genes were conservative in maize and other plant species. The functional conservation indicated that some of the 575 maize genes (orthologous to 247 genes controlling yield or seed traits in other plant species) might be functionally related to maize yield-related traits, especially the 49 maize orthologous genes in cQTL regions, and 41 orthologous genes close to the physical positions of significant SNPs. In the end, we prospected on the integration of the public sources for exploring the genetic and molecular mechanisms of maize yield-related traits, and on the utilization of genetic and molecular mechanisms for maize improvement.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01214-3.

TaCKX gene family, at large, is associated with thousand-grain weight and plant height in common wheat

KEY MESSAGE

We used SMRT sequencing and explored the haplotypes of TaCKX genes, linked with thousand-grain weight and plant height, and developed the functionally validated markers, which can be used in the marker-assisted breeding program. Cytokinin oxidase/dehydrogenase (CKX) enzymes catalyze the permanent degradation of cytokinins. Identification of the TaCKX alleles associated with yield traits and the development of functional markers is the first step in using these alleles in marker-assisted breeding program. To identify the alleles, we sequenced the genome fragments, containing TaCKX genes from 48 wheat genotypes, by PacBio® sequencing. Six out of 22 TaCKX genes were found polymorphic, forming 14 distinct haplotypes. Functional markers were developed and validated for all the polymorphic TaCKX genes. Four specific haplotypes, i.e., TaCKX2A_2, TaCKX4A_2, TaCKX5A_3, and TaCKX9A_2, were found significantly associated with high thousand-grain weight (TGW) and short plant height (PH) in Chinese wheat micro-core collection (MCC) and GWAS open population (GWAS-OP), whereas TaCKX1B_2 in GWAS-OP and TaCKX11A_3 in MCC were significantly associated with high TGW and short PH. The mean values of TGW and PH for cumulative favorable haplotypes from chromosome 3A, i.e., TaCKX2A_2, TaCKX4A_2, and TaCKX5A_3, were significantly higher as compared to the cumulative unfavored haplotypes, and the change was additive in manner. Frequency distribution analysis revealed that since the 1960s, the frequency of the favorable haplotypes and TGW has gradually increased in Chinese wheat cultivars. Expression profiling in the seed tissue excised at 2, 4, 6, and 8 days after anthesis depicted that the favorable haplotypes are significantly less expressive as compared to the unfavored haplotypes. We conclude that the functional markers developed in this study can be used to select the favorable haplotypes of TaCKX genes in wheat marker-assisted breeding programs.

Molecular characteristics of segment 5, a unique fragment encoding two partially overlapping ORFs in the genome of rice black-streaked dwarf virus

Rice black-streaked dwarf virus (RBSDV), a ds-RNA virus in Fijivirus genus with family Reoviridae, which is transmitted by the small brown planthopper, is responsible for incidence of maize rough dwarf disease (MRDD) and rice black-streaked dwarf disease (RBSDD). To understand the variation and evolution of S5, a unique fragment in the genome of RBSDV which encodes two partially overlapping ORFs (ORF5-1 and ORF5-2), we analyzed 127 sequences from maize and rice exhibiting symptoms of dwarfism. The nucleotide diversity of both ORF5-1 (π = 0.039) and ORF5-2 (π = 0.027) was higher than that of the overlapping region (π = 0.011) (P < 0.05). ORF5-2 was under the greatest selection pressure based on codon bias analysis, and its activation was possibly influenced by the overlapping region. The recombinant fragments of three recombinant events (14NM23, 14BM20, and 14NM17) cross the overlapping region. Based on neighbor-joining tree analysis, the overlapping region could represent the evolutionary basis of the full-length S5, which was classified into three main groups. RBSDV populations were expanding and haplotype diversity resulted mainly from the overlapping region. The genetic differentiation of combinations (T127-B35, T127-J34, A58-B35, A58-J34, and B35-J34) reached significant or extremely significant levels. Gene flow was most frequent between subpopulations A58 and B35, with the smallest |Fst| (0.02930). We investigated interactions between 13 RBSDV proteins by two-hybrid screening assays and identified interactions between P5-1/P6, P6/P9-1, and P3/P6. We also observed self-interactive effects of P3, P6, P7-1, and P10. In short, we have proven that RBSDV populations were expanding and the overlapping region plays an important role in the genetic variation and evolution of RBSDV S5. Our results enable ongoing research into the evolutionary history of RBSDV-S5 with two partly overlapping ORFs.

Natural variation in the THICK TASSEL DWARF1 (TD1) gene in the regulation of maize (Zea mays L.) ear-related traits

THICK TASSEL DWARF1 (TD1) is a critical gene that negatively modulates meristem size during maize inflorescence development and may also regulate ear-related traits. In the present study, the whole genomic DNA sequences and the promoter regions of TD1 were analyzed in 165 diverse maize inbred lines. Polymorphism analysis identified 39 SNPs and five InDels in the genic region of TD1 and allowed 23 haplotypes to be classified. Among these sites, eight SNPs and one InDel were significantly associated with kernel number (KN) (p < 0.05), seven SNPs and one InDel were significantly associated with kernel number per row (KNPR) (p < 0.05), and three SNPs were significantly associated with kernel row number (KRN) (p < 0.05). In addition, 21 SNPs and 14 InDels were identified in the promoter regions of TD1, and two SNPs and seven InDels of these sites were found to be significantly associated with KRN (p < 0.05). The results denote that Hap_7 was the favorable haplotype for both KN and KNPR, and Hap_2 was the elite haplotype for KRN. These favorable haplotypes could be utilized in molecular marker-assisted selection (MAS) to improve KN, KNPR, or KRN, and thereby further increase grain yield in maize breeding programs.

Hormonal regulation of cereal endosperm development with a focus on rice (Oryza sativa)

The endosperm of cereal grain forms the staple diet for most of the world's population, and feeds much of their stock. Grain size and quality are determined largely by events taking place during coenocytic nuclear division, endosperm cellularisation and cell differentiation, and the production of storage molecules. Thus, understanding the complex signalling processes occurring at each of these steps is essential for maintaining and improving our food supply. Here, we critically review evidence for the effects of phytohormones on grain size, as well as hormone homeostasis, signalling and crosstalk. We focus on rice endosperm due to the importance of rice as a food crop and a model grass, as well as its relative neglect in recent reviews; however, data from other cereals are also discussed due to strong evidence for conserved signalling networks operating during grain development. Discussion is restricted to auxin, cytokinin, ethylene, abscisic acid and gibberellin. Our review highlights the need for accurate hormone determinations combined with information on gene expression. We present evidence for separate, localised signalling roles for auxin at different stages of grain development and highlight key research questions for other hormones where much less data are available.

Differences in Molecular Characteristics of Segment 8 in Rice black-streaked dwarf virus and Southern rice black-streaked dwarf virus

Rice black-streaked dwarf virus (RBSDV) and Southern rice black-streaked dwarf virus (SRBSDV) cause maize rough dwarf disease (MRDD) and rice black-streaked dwarf disease (RBSDD) in China. RBSDV segment 8 (S8) contains the only deletion mutation in the genomes of these viruses, which are both members of the genus Fijivirus. To illuminate the molecular differences between the RBSDV and SRBSDV genomes and better understand the evolution of these viruses, and to determine which virus is specifically associated with MRDD and RBSDD in each region, S8 was analyzed in 66 virus isolates collected from 10 geographic locations in China and 14 S8 sequences obtained from the National Center for Biotechnology Information GenBank. Phylogenetic analysis showed that the pathogen associated with MRDD and RBSDD in the Yellow and Huai River valleys was RBSDV, whereas the pathogen associated with these diseases in Sanya was SRBSDV. Codon usage bias in S8 differed significantly between RBSDV and SRBSDV, as indicated by effective number of codons used by a gene (Nc) and GC values, Nc plots, and variation explained by the first axis in correspondence analysis. The nucleotide identities among these 66 RBSDV and SRBSDV isolates ranged from 66.2 to 68.2%, and were considerably lower than the nucleotide identities within RBSDV (from 94.1 to 99.9%) or SRBSDV (from 93.9 to 100%) isolates. Most S8 polymorphisms were identified in the region from 1,000 to 1,200 bp in RBSDV and in the region from 500 to 700 bp in SRBSDV. The difference in the lengths of RBSDV (1,936 bp) and SRBSDV (1,928 bp) was due to an 8-bp deletion in the 3'-untranslated region of SRBSDV. Six recombination events were detected in S8 in RBSDV and two recombination events were detected in S8 in SRBSDV. Recombination breakpoints were found within the region containing the deletion mutation in nine isolates. However, no recombination events were detected between RBSDV and SRBSDV. Both of these viruses were under negative and purifying selection, although the ratio of nonsynonymous mutations to synonymous mutations (Ka/Ks) for RBSDV S8 (0.0530) was not significantly lower than that of SRBSDV S8 (0.0823, P = 0.1550). We found that SRBSDV was more highly genetically differentiated (product of effective population size and the migration rate among populations < 1; values for the among-populations component of genetic variation or normalized variation > 0.33; and P values of the sequence statistic, the rank statistic, and the nearest-neighbor statistic < 0.01) than RBSDV. However, gene flow between RBSDV and SRBSDV was not frequent.

Diverse Roles of Ethylene in Regulating Agronomic Traits in Rice

Gaseous hormone ethylene has diverse effects in various plant processes. These processes include seed germination, plant growth, senescence, fruit ripening, biotic and abiotic stresses responses, and many other aspects. The biosynthesis and signaling of ethylene have been extensively studied in model Arabidopsis in the past two decades. However, knowledge about the ethylene signaling mechanism in crops and roles of ethylene in regulation of crop agronomic traits are still limited. Our recent findings demonstrate that rice possesses both conserved and diverged mechanism for ethylene signaling compared with Arabidopsis. Here, we mainly focused on the recent advances in ethylene regulation of important agronomic traits. Of special emphasis is its impact on rice growth, flowering, grain filling, and grain size control. Similarly, the influence of ethylene on other relevant crops will be compared. Additionally, interactions of ethylene with other hormones will also be discussed in terms of crop growth and development. Increasing insights into the roles and mechanisms of ethylene in regulating agronomic traits will contribute to improvement of crop production through precise manipulation of ethylene actions in crops.

Rice black-streaked dwarf virus Genome in China: Diversification, Phylogeny, and Selection

Rice black-streaked dwarf virus (RBSDV), a Fijivirus, causes maize rough dwarf disease and rice black-streaked dwarf disease in the summer maize-growing regions of the Yellow and Huai rivers, respectively, in China. Nevertheless, the diversification and selection of the entire genome from S1 to S10 have not been illuminated. Molecular variation, evolution, conserved regions, and other genomic properties were analyzed in 21 RBSDV isolates from maize (Zea mays L.) and rice (Oryza sativa) hosts sampled from nine geographic locations in China. Low codon adaptation index values ranging from 0.1878 to 0.2918 indicated a low degree of codon-usage bias and low potential expression for all 13 RBSDV open reading frames (ORFs). ORF9-2 showed a stronger effect of codon usage bias than did other ORFs, as the majority of points for this ORF lay close to the standard curve in the Nc plot (the effective number of codons [Nc] versus the frequency of G+C at synonymous third-base positions [GC3]). A 9-bp deletion mutation was detected in the RBSDV genome in the 3' UTR of S8. Nucleotide diversity analysis indicated that the structural proteins of RBSDV, such as S2 and S4, were all more conserved than nonstructural proteins such as S9. Nucleotide diversity (π) was highest among S9 sequences (0.0656), and was significantly higher than among S4 sequences (0.0225, P < 0.01). The number of conserved regions among the 10 segments varied substantially. The highest number of conserved regions (5) was found in S5, whereas no conserved regions were identified in S9. Nucleotide diversity and the number of conserved regions were independent of the lengths of segments. Nucleotide diversity was also not correlated with the number of conserved regions in segments. Ten recombination events in 21 isolates were found in seven segments with breakpoint positions in UTRs, intergenic spacer regions, and gene coding regions. The number of recombination events was also independent of the lengths of segments. RBSDV isolates from China could be phylogenetically classified into two groups using either 10 segment sequences or the concatenated sequence of S1 through S10, regardless of host or geographical location. The phylogenetic tree generated from pairwise nucleotide identities of individual RBSDV segments such as S9 and S3, with nucleotide identity values of 93.74% and 95.86%, respectively, is similar to the tree constructed from the concatenated sequences of the entire RBSDV genome. The 13 RBSDV ORFs were under negative and purifying selection (Ka/Ks < 1). ORF5-2 was under the greatest selection pressure; however, ORF2, which encodes the core protein of RBSDV, was under the lowest selection pressure.

Analysis of the genetic architecture of maize ear and grain morphological traits by combined linkage and association mapping

Using combined linkage and association mapping, 26 stable QTL and six stable SNPs were detected across multiple environments for eight ear and grain morphological traits in maize. One QTL, PKS2, might play an important role in maize yield improvement. In the present study, one bi-parental population and an association panel were used to identify quantitative trait loci (QTL) for eight ear and grain morphological traits. A total of 108 QTL related to these traits were detected across four environments using an ultra-high density bin map constructed using recombinant inbred lines (RILs) derived from a cross between Ye478 and Qi319, and 26 QTL were identified in more than two environments. Furthermore, 64 single nucleotide polymorphisms (SNPs) were found to be significantly associated with the eight ear and grain morphological traits (-log₁₀(P) > 4) in an association panel of 240 maize inbred lines. Combining the two mapping populations, a total of 17 pleiotropic QTL/SNPs (pQTL/SNPs) were associated with various traits across multiple environments. PKS2, a stable locus influencing kernel shape identified on chromosome 2 in a genome-wide association study (GWAS), was within the QTL confidence interval defined by the RILs. The candidate region harbored a short 13-Kb LD block encompassing four SNPs (SYN11386, PHM14783.16, SYN11392, and SYN11378). In the association panel, 13 lines derived from the hybrid PI78599 possessed the same allele as Qi319 at the PHM14783.16 (GG) locus, with an average value of 0.21 for KS, significantly lower than that of the 34 lines derived from Ye478 that carried a different allele (0.25, P < 0.05). Therefore, further fine mapping of PKS2 will provide valuable information for understanding the genetic components of grain yield and improving molecular marker-assisted selection (MAS) in maize.

Whole-Genome Analysis of Candidate genes Associated with Seed Size and Weight in Sorghum bicolor Reveals Signatures of Artificial Selection and Insights into Parallel Domestication in Cereal Crops

Seed size and seed weight are major quality attributes and important determinants of yield that have been strongly selected for during crop domestication. Limited information is available about the genetic control and genes associated with seed size and weight in sorghum. This study identified sorghum orthologs of genes with proven effects on seed size and weight in other plant species and searched for evidence of selection during domestication by utilizing resequencing data from a diversity panel. In total, 114 seed size candidate genes were identified in sorghum, 63 of which exhibited signals of purifying selection during domestication. A significant number of these genes also had domestication signatures in maize and rice, consistent with the parallel domestication of seed size in cereals. Seed size candidate genes that exhibited differentially high expression levels in seed were also found more likely to be under selection during domestication, supporting the hypothesis that modification to seed size during domestication preferentially targeted genes for intrinsic seed size rather than genes associated with physiological factors involved in the carbohydrate supply and transport. Our results provide improved understanding of the complex genetic control of seed size and weight and the impact of domestication on these genes.

Whole-Genome Analysis of Candidate genes Associated with Seed Size and Weight in Sorghum bicolor Reveals Signatures of Artificial Selection and Insights into Parallel Domestication in Cereal Crops

Seed size and seed weight are major quality attributes and important determinants of yield that have been strongly selected for during crop domestication. Limited information is available about the genetic control and genes associated with seed size and weight in sorghum. This study identified sorghum orthologs of genes with proven effects on seed size and weight in other plant species and searched for evidence of selection during domestication by utilizing resequencing data from a diversity panel. In total, 114 seed size candidate genes were identified in sorghum, 63 of which exhibited signals of purifying selection during domestication. A significant number of these genes also had domestication signatures in maize and rice, consistent with the parallel domestication of seed size in cereals. Seed size candidate genes that exhibited differentially high expression levels in seed were also found more likely to be under selection during domestication, supporting the hypothesis that modification to seed size during domestication preferentially targeted genes for intrinsic seed size rather than genes associated with physiological factors involved in the carbohydrate supply and transport. Our results provide improved understanding of the complex genetic control of seed size and weight and the impact of domestication on these genes.

Molecular evolution and functional characterisation of haplotypes of an important rubber biosynthesis gene in Hevea brasiliensis

Hydroxy-methylglutaryl coenzyme-A synthase (HMGS) is a rate-limiting enzyme in the cytoplasmic isoprenoid biosynthesis pathway leading to natural rubber production in Hevea brasiliensis (rubber). Analysis of the structural variants of this gene is imperative to understand their functional significance in rubber biosynthesis so that they can be properly utilised for ongoing crop improvement programmes in Hevea. We report here allele richness and diversity of the HMGS gene in selected popular rubber clones. Haplotypes consisting of single nucleotide polymorphisms (SNPs) from the coding and non-coding regions with a high degree of heterozygosity were identified. Segregation and linkage disequilibrium analysis confirmed that recombination is the major contributor to the generation of allelic diversity, rather than point mutations. The evolutionarily conserved nature of some SNPs was identified by comparative DNA sequence analysis of HMGS orthologues from diverse taxa, demonstrating the molecular evolution of rubber biosynthesis genes in general. In silico three-dimensional structural studies highlighting the structural positioning of non-synonymous SNPs from different HMGS haplotypes revealed that the ligand-binding site on the enzyme remains impervious to the reported sequence variations. In contrast, gene expression results indicated the possibility of association between specific haplotypes and HMGS expression in Hevea clones, which may have a downstream impact up to the level of rubber production. Moreover, haplotype diversity of the HMGS gene and its putative association with gene expression can be the basis for further genetic association studies in rubber. Furthermore, the data also show the role of SNPs in the evolution of candidate genes coding for functional traits in plants.

A comprehensive meta-analysis of plant morphology, yield, stay-green, and virus disease resistance QTL in maize (Zea mays L.)

The meta-QTL and candidate genes will facilitate the elucidation of molecular bases underlying agriculturally important traits and open new avenues for functional markers development and elite alleles introgression in maize breeding program. A large number of QTLs attributed to grain productivity and other agriculturally important traits have been identified and deposited in public repositories. The integration of fruitful QTL becomes a major issue in current plant genomics. To this end, we first collected QTL for six agriculturally important traits in maize, including yield, plant height, ear height, leaf angle, stay-green, and maize rough dwarf disease resistance. The meta-analysis method was then employed to retrieve 113 meta-QTL. Additionally, we also isolated candidate genes for target traits by the bioinformatic technique. Several candidates, including some well-characterized genes, GA3ox2 for plant height, lg1 and lg4 for leaf angle, zfl1 and zfl2 for flowering time, were co-localized with established meta-QTL intervals. Intriguingly, in a relatively narrow meta-QTL region, the maize ortholog of rice yield-related gene GW8/OsSPL16 was believed to be a candidate for yield. Leveraging results presented in this study will provide further insights into the genetic architecture of maize agronomic traits. Moreover, the meta-QTL and candidate genes reported here could be harnessed for the enhancement of stress tolerance and yield performance in maize and translation to other crops.

Molecular Genetic Analysis and Evolution of Segment 7 in Rice Black-Streaked Dwarf Virus in China

Rice black-streaked dwarf virus (RBSDV) causes maize rough dwarf disease or rice black-streaked dwarf disease and can lead to severe yield losses in maize and rice. To analyse RBSDV evolution, codon usage bias and genetic structure were investigated in 111 maize and rice RBSDV isolates from eight geographic locations in 2013 and 2014. The linear dsRNA S7 is A+U rich, with overall codon usage biased toward codons ending with A (A3s, S7-1: 32.64%, S7-2: 29.95%) or U (U3s, S7-1: 44.18%, S7-2: 46.06%). Effective number of codons (Nc) values of 45.63 in S7-1 (the first open reading frame of S7) and 39.96 in S7-2 (the second open reading frame of S7) indicate low degrees of RBSDV-S7 codon usage bias, likely driven by mutational bias regardless of year, host, or geographical origin. Twelve optimal codons were detected in S7. The nucleotide diversity (π) of S7 sequences in 2013 isolates (0.0307) was significantly higher than in 2014 isolates (0.0244, P = 0.0226). The nucleotide diversity (π) of S7 sequences in isolates from Jinan (0.0391) was higher than that from the other seven locations (P < 0.01). Only one S7 recombinant was detected in Baoding. RBSDV isolates could be phylogenetically classified into two groups according to S7 sequences, and further classified into two subgroups. S7-1 and S7-2 were under negative and purifying selection, with respective Ka/Ks ratios of 0.0179 and 0.0537. These RBSDV populations were expanding (P < 0.01) as indicated by negative values for Tajima's D, Fu and Li's D, and Fu and Li's F. Genetic differentiation was detected in six RBSDV subpopulations (P < 0.05). Absolute Fst (0.0790) and Nm (65.12) between 2013 and 2014, absolute Fst (0.1720) and Nm (38.49) between maize and rice, and absolute Fst values of 0.0085-0.3069 and Nm values of 0.56-29.61 among these eight geographic locations revealed frequent gene flow between subpopulations. Gene flow between 2013 and 2014 was the most frequent.

Phylogenetic and recombination analysis of rice black-streaked dwarf virus segment 9 in China

Rice black-streaked dwarf virus (RBSDV) is an economically important virus that causes maize rough dwarf disease and rice black-streaked dwarf disease in East Asia. To study RBSDV variation and recombination, we examined the segment 9 (S9) sequences of 49 RBSDV isolates from maize and rice in China. Three S9 recombinants were detected in Baoding, Jinan, and Jining, China. Phylogenetic analysis showed that Chinese RBSDV isolates could be classified into two groups based on their S9 sequences, regardless of host or geographical origin. Further analysis suggested that S9 has undergone negative and purifying selection.