TILLING-by-Sequencing+ to Decipher Oil Biosynthesis Pathway in Soybeans: A New and Effective Platform for High-Throughput Gene Functional Analysis

Reverse genetic approaches have been widely applied to study gene function in crop species; however, these techniques, including gel-based TILLING, present low efficiency to characterize genes in soybeans due to genome complexity, gene duplication, and the presence of multiple gene family members that share high homology in their DNA sequence. Chemical mutagenesis emerges as a genetically modified-free strategy to produce large-scale soybean mutants for economically important traits improvement. The current study uses an optimized high-throughput TILLING by target capture sequencing technology, or TILLING-by-Sequencing⁺ (TbyS⁺), coupled with universal bioinformatic tools to identify population-wide mutations in soybeans. Four ethyl methanesulfonate mutagenized populations (4032 mutant families) have been screened for the presence of induced mutations in targeted genes. The mutation types and effects have been characterized for a total of 138 soybean genes involved in soybean seed composition, disease resistance, and many other quality traits. To test the efficiency of TbyS⁺ in complex genomes, we used soybeans as a model with a focus on three desaturase gene families, GmSACPD, GmFAD2, and GmFAD3, that are involved in the soybean fatty acid biosynthesis pathway. We successfully isolated mutants from all the six gene family members. Unsurprisingly, most of the characterized mutants showed significant changes either in their stearic, oleic, or linolenic acids. By using TbyS⁺, we discovered novel sources of soybean oil traits, including high saturated and monosaturated fatty acids in addition to low polyunsaturated fatty acid contents. This technology provides an unprecedented platform for highly effective screening of polyploid mutant populations and functional gene analysis. The obtained soybean mutants from this study can be used in subsequent soybean breeding programs for improved oil composition traits.     

Global Identification and Characterization of C2 Domain-Containing Proteins Associated with Abiotic Stress Response in Rice (Oryza sativa L.)

C2 domain-containing proteins (C2DPs) have been identified in different genomes that contain single or multiple C2 domains in their C- or N-terminal. It possesses higher functional activity in the transmembrane regions. The identification of C2 domains were reported in a previous study, such as multiple C2 domains and transmembrane-region proteins (MCTPs) and N-terminal-TM-C2 domain proteins (NTMC2s) of rice, Arabidopsis thaliana, and cotton, whereas the C2DP gene family in rice has not been comprehensively studied, and the role of the C2DP gene in rice in response to abiotic stress is not yet fully understood. In this study, we identified 82 C2DPs in the rice genome and divided them into seven groups through phylogenetic analysis. The synteny analysis revealed that duplication events were either exhibited within the genome of rice or between the genomes of rice and other species. Through the analysis of cis-acting elements in promoters, expression profiles, and qRT-PCR results, the functions of OsC2DPs were found to be widely distributed in diverse tissues and were extensively involved in phytohormones-related and abiotic stresses response in rice. The prediction of the microRNA (miRNA) targets of OsC2DPs revealed the possibility of regulation by consistent miRNAs. Notably, OsC2DP50/51/52 as a co-tandem duplication exhibited similar expression variations and involved the coincident miRNA-regulation pathway. Moreover, the results of the genotypic variation and haplotype analysis revealed that OsC2DP17, OsC2DP29, and OsC2DP49 were associated with cold stress responses. These findings provided comprehensive insights for characterizations of OsC2DPs in rice as well as for their roles for abiotic stress.     

Functional Specialization within the EXO70 Gene Family in Arabidopsis

Localized delivery of plasma-membrane and cell-wall components is a crucial process for plant cell growth. One of the regulators of secretory-vesicle targeting is the exocyst tethering complex. The exocyst mediates first interaction between transport vesicles and the target membrane before their fusion is performed by SNARE proteins. In land plants, genes encoding the EXO70 exocyst subunit underwent an extreme proliferation with 23 paralogs present in the Arabidopsis (Arabidopsis thaliana) genome. These paralogs often acquired specialized functions during evolution. Here, we analyzed functional divergence of selected EXO70 paralogs in Arabidopsis. Performing a systematic cross-complementation analysis of exo70a1 and exo70b1 mutants, we found that EXO70A1 was functionally substituted only by its closest paralog, EXO70A2. In contrast, none of the EXO70 isoforms tested were able to substitute EXO70B1, including its closest relative, EXO70B2, pointing to a unique function of this isoform. The presented results document a high degree of functional specialization within the EXO70 gene family in land plants.     

Genome-Wide Identification of the AGC Protein Kinase Gene Family Related to Photosynthesis in Rice (Oryza sativa)

The cAMP-dependent protein kinase A, cGMP-dependent protein kinase G and phospholipid-dependent protein kinase C (AGC) perform various functions in plants, involving growth, immunity, apoptosis and stress response. AGC gene family is well described in Arabidopsis, however, limited information is provided about AGC genes in rice, an important cereal crop. This research studied the AGC gene family in the AA genome species: Oryza sativa ssp. japonica, Oryza sativa ssp. indica, Oryza nivara, Oryza rufipogon, Oryza glaberrima, Oryza meridionalis, Oryza barthii, Oryza glumaepatula and Oryza longistaminata were searched and classified into six subfamilies, and it was found that these species have similar numbers of members. The analysis of gene duplication and selection pressure indicated that the AGC gene family expanded mainly by segmental or whole genome duplication (WGD), with purifying selection during the long evolutionary period. RNA-seq analysis revealed that OsAGCs of subfamily V were specifically highly expressed in leaves, and the expression patterns of these genes were compared with that of photosynthesis-related genes using qRT-PCR, discovered that OsAGC9, OsAGC20, and OsAGC22 might participate in photosynthesis. These results provide an informative perspective for exploring the evolutionary of AGC gene family and its practical application in rice.     

Extreme Enlargement of the Inverted Repeat Region in the Plastid Genomes of Diatoms from the Genus Climaconeis

We sequenced the plastid genomes of three diatoms from the genus Climaconeis, including two strains formerly designated as Climaconeis scalaris. At 208,097 and 216,580 bp, the plastid genomes of the latter strains are the largest ever sequenced among diatoms and their increased size is explained by the massive expansion of the inverted repeat region. Important rearrangements of gene order were identified among the two populations of Climaconeis cf. scalaris. The other sequenced Climaconeis chloroplast genome is 1.5 times smaller compared with those of the Climaconeis cf. scalaris strains and it features an usual quadripartite structure. The extensive structural changes reported here for the genus Climaconeis are compared with those previously observed for other algae and plants displaying large plastid genomes.     

Comparative Mitogenomics of the Genus Odontobutis (Perciformes: Gobioidei: Odontobutidae) Revealed Conserved Gene Rearrangement and High Sequence Variations

To understand the molecular evolution of mitochondrial genomes (mitogenomes) in the genus Odontobutis, the mitogenome of Odontobutis yaluensis was sequenced and compared with those of another four Odontobutis species. Our results displayed similar mitogenome features among species in genome organization, base composition, codon usage, and gene rearrangement. The identical gene rearrangement of trnS-trnL-trnH tRNA cluster observed in mitogenomes of these five closely related freshwater sleepers suggests that this unique gene order is conserved within Odontobutis. Additionally, the present gene order and the positions of associated intergenic spacers of these Odontobutis mitogenomes indicate that this unusual gene rearrangement results from tandem duplication and random loss of large-scale gene regions. Moreover, these mitogenomes exhibit a high level of sequence variation, mainly due to the differences of corresponding intergenic sequences in gene rearrangement regions and the heterogeneity of tandem repeats in the control regions. Phylogenetic analyses support Odontobutis species with shared gene rearrangement forming a monophyletic group, and the interspecific phylogenetic relationships are associated with structural differences among their mitogenomes. The present study contributes to understanding the evolutionary patterns of Odontobutidae species.     

A Novel Mitochondrial Genome Fragmentation Pattern in the Buffalo Louse Haematopinus tuberculatus (Psocodea: Haematopinidae)

Sucking lice are obligate ectoparasites of mammalian hosts, causing serious public health problems and economic losses worldwide. It is well known that sucking lice have fragmented mitochondrial (mt) genomes, but many remain undetermined. To better understand patterns of mt genome fragmentation in the sucking lice, we sequenced the mt genome of the buffalo louse Haematopinus tuberculatus using next-generation sequencing (NGS). The mt genome of H. tuberculatus has ten circular minichromosomes containing a total of 37 genes. Each minichromosome is 2.9–5.0 kb long and carries one to eight genes plus one large non-coding region. The number of mt minichromosomes of H. tuberculatus (ten) is different from those of congeneric species (horse louse H. asini, domestic pig louse H. suis and wild pig louse H. apri) and other sucking lice. Two events (gene translocation and merger of mt minichromosome) are observed in Haematopinus. Compared to other studies, our phylogeny generated from mt genome datasets showed a different topology, suggesting that inclusion of data other than mt genomes would be required to resolve phylogeny of sucking lice. To our knowledge, this is the first report of a ten mt minichromosomes genome in sucking lice, which opens a new outlook into unexplored mt genome fragmentation patterns in sucking lice.     

Complete Plastid and Mitochondrial Genomes of Aeginetia indica Reveal Intracellular Gene Transfer (IGT), Horizontal Gene Transfer (HGT), and Cytoplasmic Male Sterility (CMS)

Orobanchaceae have become a model group for studies on the evolution of parasitic flowering plants, and Aeginetia indica, a holoparasitic plant, is a member of this family. In this study, we assembled the complete chloroplast and mitochondrial genomes of A. indica. The chloroplast and mitochondrial genomes were 56,381 bp and 401,628 bp long, respectively. The chloroplast genome of A. indica shows massive plastid genes and the loss of one IR (inverted repeat). A comparison of the A. indica chloroplast genome sequence with that of a previous study demonstrated that the two chloroplast genomes encode a similar number of proteins (except atpH) but differ greatly in length. The A. indica mitochondrial genome has 53 genes, including 35 protein-coding genes (34 native mitochondrial genes and one chloroplast gene), 15 tRNA (11 native mitochondrial genes and four chloroplast genes) genes, and three rRNA genes. Evidence for intracellular gene transfer (IGT) and horizontal gene transfer (HGT) was obtained for plastid and mitochondrial genomes. ψndhB and ψcemA in the A. indica mitogenome were transferred from the plastid genome of A. indica. The atpH gene in the plastid of A. indica was transferred from another plastid angiosperm plastid and the atpI gene in mitogenome A. indica was transferred from a host plant like Miscanthus siensis. Cox2 (orf43) encodes proteins containing a membrane domain, making ORF (Open Reading Frame) the most likely candidate gene for CMS development in A. indica.     

Genome-Wide Comprehensive Analysis of the GASA Gene Family in Populus

Gibberellic acid-stimulated Arabidopsis (GASA) proteins, as cysteine-rich peptides (CRPs), play roles in development and reproduction and biotic and abiotic stresses. Although the GASA gene family has been identified in plants, the knowledge about GASAs in Populus euphratica, the woody model plant for studying abiotic stress, remains limited. Here, we referenced the well-sequenced Populus trichocarpa genome, and identified the GASAs in the whole genome of P. euphratica and P. trichocarpa. 21 candidate genes in P. trichocarpa and 19 candidate genes in P. euphratica were identified and categorized into three subfamilies by phylogenetic analysis. Most GASAs with signal peptides were located extracellularly. The GASA genes in Populus have experienced multiple gene duplication events, especially in the subfamily A. The evolution of the subfamily A, with the largest number of members, can be attributed to whole-genome duplication (WGD) and tandem duplication (TD). Collinearity analysis showed that WGD genes played a leading role in the evolution of GASA genes subfamily B. The expression patterns of P. trichocarpa and P. euphratica were investigated using the PlantGenIE database and the real-time quantitative PCR (qRT-PCR), respectively. GASA genes in P. trichocarpa and P. euphratica were mainly expressed in young tissues and organs, and almost rarely expressed in mature leaves. GASA genes in P. euphratica leaves were also widely involved in hormone responses and drought stress responses. GUS activity assay showed that PeuGASA15 was widely present in various organs of the plant, especially in vascular bundles, and was induced by auxin and inhibited by mannitol dramatically. In summary, this present study provides a theoretical foundation for further research on the function of GASA genes in P. euphratica.     

Frontiers in Dissecting and Managing Brassica Diseases: From Reference-Based RGA Candidate Identification to Building Pan-RGAomes

The Brassica genus contains abundant economically important vegetable and oilseed crops, which are under threat of diseases caused by fungal, bacterial and viral pathogens. Resistance gene analogues (RGAs) are associated with quantitative and qualitative disease resistance and the identification of candidate RGAs associated with disease resistance is crucial for understanding the mechanism and management of diseases through breeding. The availability of Brassica genome assemblies has greatly facilitated reference-based quantitative trait loci (QTL) mapping for disease resistance. In addition, pangenomes, which characterise both core and variable genes, have been constructed for B. rapa, B. oleracea and B. napus. Genome-wide characterisation of RGAs using conserved domains and motifs in reference genomes and pangenomes reveals their clustered arrangements and presence of structural variations. Here, we comprehensively review RGA identification in important Brassica genome and pangenome assemblies. Comparison of the RGAs in QTL between resistant and susceptible individuals allows for efficient identification of candidate disease resistance genes. However, the reference-based QTL mapping and RGA candidate identification approach is restricted by the under-represented RGA diversity characterised in the limited number of Brassica assemblies. The species-wide repertoire of RGAs make up the pan-resistance gene analogue genome (pan-RGAome). Building a pan-RGAome, through either whole genome resequencing or resistance gene enrichment sequencing, would effectively capture RGA diversity, greatly expanding breeding resources that can be utilised for crop improvement.     

The Genome of the Mitochondrion-Related Organelle in Cepedea longa,a Large Endosymbiotic Opalinid Inhabiting the Recta of Frogs

Mitochondrion-related organelles (MROs) are loosely defined as degenerated mitochondria in anaerobic and microaerophilic lineages. Opalinids are commonly regarded as commensals in the guts of cold-blooded amphibians. It may represent an intermediate adaptation stage between the conventional aerobic mitochondria and derived anaerobic MROs. In the present study, we sequenced and analyzed the MRO genome of Cepedea longa. It has a linear MRO genome with large inverted repeat gene regions at both ends. Compared to Blastocystis and Proteromonas lacertae, the MRO genome of C. longa has a higher G + C content and repeat sequences near the central region. Although three Opalinata species have different morphological characteristics, phylogenetic analyses based on eight concatenated nad genes indicate that they are close relatives. The phylogenetic analysis showed that C. longa clustered with P. lacertae with strong support. The 18S rRNA gene-based phylogeny resolved the Opalinea clade as a sister clade to Karotomorpha, which then further grouped with Proteromonas. The paraphyly of Proteromonadea needs to be verified due to the lack of MRO genomes for key species, such as Karotomorpha, Opalina and Protoopalina. Besides, our dataset and analyses offered slight support for the paraphyly of Bigyra.     

Genome-Wide Analysis of the Cyclin Gene Family and Their Expression Profile in Medicago truncatula

Cyclins, together with highly conserved cyclin-dependent kinases (CDKs), play an important role in the process of cell cycle in plants, but less is known about the functions of cyclins in legume plants, especially Medicago truncatula. Our genome-wide analysis identified 58, 103, and 51 cyclin members in the M. truncatula, Glycine max, and Phaseolus vulgaris genomes. Phylogenetic analysis suggested that these cyclins could be classified into 10 types, and the CycB-like types (CycBL1-BL8) were the specific subgroups in M. truncatula, which was one reason for the expansion of the B-type in M. truncatula. All putative cyclin genes were mapped onto their own chromosomes of each genome, and 9 segmental duplication gene pairs involving 20 genes were identified in M. truncatula cyclins. Determined by quantitative real-time PCR, the expression profiling suggested that 57 cyclins in M. truncatula were differentially expressed in 9 different tissues, while a few genes were expressed in some specific tissues. Using the publicly available RNAseq data, the expression of Mtcyclins in the wild-type strain A17 and three nodule mutants during rhizobial infection showed that 23 cyclins were highly upregulated in the nodulation (Nod) factor-hypersensitive mutant sickle (skl) mutant after 12 h of rhizobium inoculation. Among these cyclins, six cyclin genes were also specifically expressed in roots and nodules, which might play specific roles in the various phases of Nod factor-mediated cell cycle activation and nodule development. Our results provide information about the cyclin gene family in legume plants, serving as a guide for further functional research on plant cyclins.     

Assembly of the Complete Mitochondrial Genome of Gelsemium elegans Revealed the Existence of Homologous Conformations Generated by a Repeat Mediated Recombination

Gelsemium elegans (G. elegans) is a Chinese medicinal plant with substantial economic and feeding values. There is a lack of detailed studies on the mitochondrial genome of G. elegans. In this study, the mitochondrial genome of G. elegans was sequenced and assembled, and its substructure was investigated. The mitochondrial genome of G. elegans is represented by two circular chromosomes of 406,009 bp in length with 33 annotated protein-coding genes, 15 tRNA genes, and three rRNA genes. We detected 145 pairs of repeats and found that four pairs of repeats could mediate the homologous recombination into one major conformation and five minor conformations, and the presence of conformations was verified by PCR amplification and Sanger sequencing. A total of 124 SSRs were identified in the G. elegans mitochondrial genome. The homologous segments between the chloroplast and mitochondrial genomes accounted for 5.85% of the mitochondrial genome. We also predicted 477 RNA potential editing sites and found that the nad4 gene was edited 38 times, which was the most frequent occurrence. Taken together, the mitochondrial genome of G. elegans was assembled and annotated. We gained a more comprehensive understanding on the genome of this medicinal plant, which is vital for its effective utilization and genetic improvement, especially for cytoplasmic male sterility breeding and evolution analysis in G. elegans.     

Diverse Evolution in 111 Plant Genomes Reveals Purifying and Dosage Balancing Selection Models for F-Box Genes

The F-box proteins function as substrate receptors to determine the specificity of Skp1-Cul1-F-box ubiquitin ligases. Genomic studies revealed large and diverse sizes of the F-box gene superfamily across plant species. Our previous studies suggested that the plant F-box gene superfamily is under genomic drift evolution promoted by epigenomic programming. However, how the size of the superfamily drifts across plant genomes is currently unknown. Through a large-scale genomic and phylogenetic comparison of the F-box gene superfamily covering 110 green plants and one red algal species, I discovered four distinct groups of plant F-box genes with diverse evolutionary processes. While the members in Clusters 1 and 2 are species/lineage-specific, those in Clusters 3 and 4 are present in over 46 plant genomes. Statistical modeling suggests that F-box genes from the former two groups are skewed toward fewer species and more paralogs compared to those of the latter two groups whose presence frequency and sizes in plant genomes follow a random statistical model. The enrichment of known Arabidopsis F-box genes in Clusters 3 and 4, along with comprehensive biochemical evidence showing that Arabidopsis members in Cluster 4 interact with the Arabidopsis Skp1-like 1 (ASK1), demonstrates over-representation of active F-box genes in these two groups. Collectively, I propose purifying and dosage balancing selection models to explain the lineage/species-specific duplications and expansions of F-box genes in plant genomes. The purifying selection model suggests that most, if not all, lineage/species-specific F-box genes are detrimental and are thus kept at low frequencies in plant genomes.     

Evolutionary and Expression Analysis of MOV10 and MOV10L1 Reveals Their Origin,Duplication and Divergence

MOV10 and MOV10L1 both encode ATP-dependent RNA helicases. In mammals, MOV10 and MOV10L1 participate in various kinds of biological contexts, such as defense of RNA virus invasion, neuron system, germ cell and early development. However, mov10 and mov10l1 in zebrafish are obscure and the evolutionary relationships of mov10 among different species remain unclear. In this study, we found MOV10 and MOV10L1 had some variations despite they possessed the conserved feature of RNA helicase, however, they may originate from a single ancestor although they shared limited homology. A single MOV10L1 gene existed among all species, while MOV10 gene experienced lineage-specific intra-chromosomal gene duplication in several species. Interestingly, the mov10 gene expanded to three in zebrafish, which originating from a duplication by whole genome specific duplication of teleost lineage followed by a specific intra-chromosome tandem duplication. The mov10 and mov10l1 showed distinct expression profiles in early stages, however, in adult zebrafish, three mov10 genes exhibited similar diverse expression patterns in almost all tissues. We also demonstrated mov10 genes were upregulated upon virus challenge, highlighting they had redundant conserved roles in virus infection. These results provide valuable data for the evolution of MOV10 and MOV10L1 and they are important to the further functional exploration.