Sugarcane (Saccharum spp.) is one of the most important economic crops because of its high sugar production and biofuel potential. Due to the high polyploid level and complex genome of sugarcane, it has been a huge challenge to investigate genomic sequence variations, which are critical for identifying alleles contributing to important agronomic traits. In order to mine the genetic variations in sugarcane, genotyping by sequencing (GBS), was used to genotype 14 representative Saccharum complex accessions. GBS is a method to generate a large number of markers, enabled by next generation sequencing (NGS) and the genome complexity reduction using restriction enzymes. To use GBS for high throughput genotyping highly polyploid sugarcane, the GBS analysis pipelines in 14 Saccharum complex accessions were established by evaluating different alignment methods, sequence variants callers, and sequence depth for single nucleotide polymorphism (SNP) filtering. By using the established pipeline, a total of 76,251 non-redundant SNPs, 5642 InDels, 6380 presence/absence variants (PAVs), and 826 copy number variations (CNVs) were detected among the 14 accessions. In addition, non-reference based universal network enabled analysis kit and Stacks de novo called 34,353 and 109,043 SNPs, respectively. In the 14 accessions, the percentages of single dose SNPs ranged from 38.3% to 62.3% with an average of 49.6%, much more than the portions of multiple dosage SNPs. Concordantly called SNPs were used to evaluate the phylogenetic relationship among the 14 accessions. The results showed that the divergence time between the Erianthus genus and the Saccharum genus was more than 10 million years ago (MYA). The Saccharum species separated from their common ancestors ranging from 0.19 to 1.65 MYA. The GBS pipelines including the reference sequences, alignment methods, sequence variant callers, and sequence depth were recommended and discussed for the Saccharum complex and other related species. A large number of sequence variations were discovered in the Saccharum complex, including SNPs, InDels, PAVs, and CNVs. Genome-wide SNPs were further used to illustrate sequence features of polyploid species and demonstrated the divergence of different species in the Saccharum complex. The results of this study showed that GBS was an effective NGS-based method to discover genomic sequence variations in highly polyploid and heterozygous species.

Thioredoxin z (TRX z) plays a significant role in chloroplast development by regulating the transcription of chloroplast genes. In this study, we identified a pentatricopeptide repeat (PPR) protein, rice albino seedling-lethal (RAS), that interacts with OsTRX z. This interaction was initially discovered by using a yeast two-hybrid (Y2H) screening technique and was further validated through Y2H and bimolecular fluorescence complementation (BiFC) experiments. RAS contains 16 PPR motifs and features a small MutS-related (SMR) domain at its C-terminus. CRISPR/Cas9-generated ras mutants exhibited an albino seedling-lethal phenotype characterized by abnormal chloroplast structures and a significantly reduced chlorophyll content. RAS localizes to the chloroplast and is predominantly expressed in young leaves. Mutations in RAS affect RNA editing at the rpl2, rps14, and ndhA sites, as well as RNA splicing at the rpl2, atpF, and ndhA transcripts within the chloroplast. Furthermore, the expression levels of genes associated with chloroplast formation are altered in the ras mutant. Both OsTRX z and RAS were found to interact with chloroplast signal recognition particle (cpSRP) proteins, indicating that their proper localization within the chloroplast may be dependent on the SRP pathway. Collectively, our findings highlight the critical role of RAS in chloroplast development, as it is involved in RNA processing and the regulation of chloroplast gene expression.

The pathogen-related protein 1 (PR-1) family plays an important role in plant response to biotic and abiotic stresses. PR-1 proteins have been studied in many plant species; however, they were not systematically studied in brown algae, which are important components of coastal ecosystems and have great economic value in the aquaculture industry. In the present study, we characterized the structure, evolution and expression of PR-1 proteins in brown algal genomes. A total of 141 PR-1s were identified in the 19 brown algal genomes, with an average of 7 genes in each species. Most PR-1s are acidic, while only 18 PR-1s are basic. Phylogenetic analysis showed that PR-1s in brown algae clustered into five clades, and showed no strong relationship with other lineages, suggesting an ancient origin. All the PR-1s contain a conserved CAP superfamily domain. Some PR-1s contain distinct functional domains, such as the WSC, Blect, and Bulb-type lectin domains, which are involved in carbohydrate binding. Their promoter regions were enriched in stress-response elements, hormone-response elements, growth and development elements. GO and KEGG annotation showed that brown algal PR-1 proteins may be involved in diverse roles and pathways. Moreover, expression analysis shows that some PR-1s, especially basic proteins are responsive to abiotic stress conditions and life stage development, further suggesting they participate in multiple functional pathways. Our results provide important data for future research on the function of brown algal PR-1 family genes.