Species-specific partial gene duplication in Arabidopsis thaliana evolved novel phenotypic effects on morphological traits under strong positive selection

Lectin diversity and their positive roles in WSSV replication through regulation of calreticulin expression and inhibiting ALFs expression

In biological evolution, gene duplication (GD) generates new genes to facilitate new functions. C-type lectins (CTLs) in crayfish have been extended by GD to expand their family members. In this study, four CTL genes generated by GD were identified from Procambarus clarkii (PcLec1-4). Among these four genes, PcLec1 can also generate new isoforms with different numbers of tandem repeats through DNA slip mispairing. PcLec1-4 was widely expressed in multiple tissues. The expression levels of PcLec1-4 were upregulated in the intestine of P. clarkii upon white spot syndrome virus (WSSV) challenge at multiple time points. Further analysis indicated that GATA transcription factor regulated PcLec1-4 expression. RNA interference and recombinant PcLec1-4 protein injection experiments suggested that PcLec1-4 promoted the expression of calreticulin (PcCRT) and negatively regulated the expression of antimicrobial peptides, thereby promoting WSSV replication. This study contributes to the understanding of the function of CTLs produced by GD during WSSV invasion in crustaceans.

Comparative Analysis of the GATA Transcription Factors in Five Solanaceae Species and Their Responses to Salt Stress in Wolfberry (Lycium barbarum L.)

GATA proteins are a class of zinc-finger DNA-binding proteins that participate in diverse regulatory processes in plants, including the development processes and responses to environmental stresses. However, a comprehensive analysis of the GATA gene family has not been performed in a wolfberry (Lycium barbarum L.) or other Solanaceae species. There are 156 GATA genes identified in five Solanaceae species (Lycium barbarum L., Solanum lycopersicum L., Capsicum annuum L., Solanum tuberosum L., and Solanum melongena L.) in this study. Based on their phylogeny, they can be categorized into four subfamilies (I-IV). Noticeably, synteny analysis revealed that dispersed- and whole-genome duplication contributed to the expansion of the GATA gene family. Purifying selection was a major force driving the evolution of GATA genes. Moreover, the predicted cis-elements revealed the potential roles of wolfberry GATA genes in phytohormone, development, and stress responses. Furthermore, the RNA-seq analysis identified 31 LbaGATA genes with different transcript profiling under salt stress. Nine candidate genes were then selected for further verification using quantitative real-time PCR. The results revealed that four candidate LbaGATA genes (LbaGATA8, LbaGATA19, LbaGATA20, and LbaGATA24) are potentially involved in salt-stress responses. In conclusion, this study contributes significantly to our understanding of the evolution and function of GATA genes among the Solanaceae species, including wolfberry.

Bioinformatics Analysis of MSH1 Genes of Green Plants: Multiple Parallel Length Expansions, Intron Gains and Losses, Partial Gene Duplications, and Alternative Splicing

MutS homolog 1 (MSH1) is involved in the recombining and repairing of organelle genomes and is essential for maintaining their stability. Previous studies indicated that the length of the gene varied greatly among species and detected species-specific partial gene duplications in Physcomitrella patens. However, there are critical gaps in the understanding of the gene size expansion, and the extent of the partial gene duplication of MSH1 remains unclear. Here, we screened MSH1 genes in 85 selected species with genome sequences representing the main clades of green plants (Viridiplantae). We identified the MSH1 gene in all lineages of green plants, except for nine incomplete species, for bioinformatics analysis. The gene is a singleton gene in most of the selected species with conserved amino acids and protein domains. Gene length varies greatly among the species, ranging from 3234 bp in Ostreococcus tauri to 805,861 bp in Cycas panzhihuaensis. The expansion of MSH1 repeatedly occurred in multiple clades, especially in Gymnosperms, Orchidaceae, and Chloranthus spicatus. MSH1 has exceptionally long introns in certain species due to the gene length expansion, and the longest intron even reaches 101,025 bp. And the gene length is positively correlated with the proportion of the transposable elements (TEs) in the introns. In addition, gene structure analysis indicated that the MSH1 of green plants had undergone parallel intron gains and losses in all major lineages. However, the intron number of seed plants (gymnosperm and angiosperm) is relatively stable. All the selected gymnosperms contain 22 introns except for Gnetum montanum and Welwitschia mirabilis, while all the selected angiosperm species preserve 21 introns except for the ANA grade. Notably, the coding region of MSH1 in algae presents an exceptionally high GC content (47.7% to 75.5%). Moreover, over one-third of the selected species contain species-specific partial gene duplications of MSH1, except for the conserved mosses-specific partial gene duplication. Additionally, we found conserved alternatively spliced MSH1 transcripts in five species. The study of MSH1 sheds light on the evolution of the long genes of green plants.

A robust SNP-haplotype assay for Bct gene region conferring resistance to beet curly top virus in common bean (Phaseolus vulgaris L.)

Beet curly top virus (BCTV), which is synonymous with curly top virus (CTV), causes significant yield loss in common bean (snap and dry beans) cultivars and several other important crops. Common bean cultivars have been found to be resistant to CTV, but screening for resistance is challenging due to the cyclical nature of epidemics and spotty feeding by the leafhopper that vectors the virus. We used an SNP dataset for the Snap Bean Association Panel (SnAP) agro-inoculated with CTV-Logan (CA/Logan) strain to locate the Bct gene region to a 1.7-Mb interval on chromosome Pv07 using genome-wide association study (GWAS) analysis. Recombinant lines from the SnAP were used to further narrow the Bct region to a 58.0-kb interval. A missense SNP (S07_2970381) in candidate gene Phvul.007G036300 Exonuclease V (EXO5) was identified as the most likely causal mutation, and it was the most significant SNP detected by GWAS in a dry bean population (DBP) naturally infected by the CTV-Worland (Wor) strain. Tm-shift assay markers developed for SNP S07_2970381 and two linked SNPs, S07_2970276 and S07_2966197, were useful for tracking different origins of the Bct EXO5 candidate gene resistance to CTV in common bean. The three SNPs identified four haplotypes, with haplotype 3-1 (Haplo3-1) of Middle American origin associated with the highest levels of CTV resistance. This SNP-haplotype assay will enable breeders to track resistance sources and to develop cultivars with better CTV resistance.

An exonuclease V homologue is expressed predominantly during early megasporogenesis in apomictic Brachiaria brizantha

MAIN CONCLUSION

An exonuclease V homologue from apomictic Brachiaria brizantha is expressed and localized in nucellar cells at key moments when these cells differentiate to give rise to unreduced gametophytes. Brachiaria is a genus of forage grasses with economical and agricultural importance to Brazil. Brachiaria reproduces by aposporic apomixis, in which unreduced embryo sacs, derived from nucellar cells, other than the megaspore mother cell (MMC), are formed. The unreduced embryo sacs produce an embryo without fertilization resulting in clones of the mother plant. Comparative gene expression analysis in ovaries of sexual and apomictic Brachiaria spp. revealed a sequence from B. brizantha that showed a distinct pattern of expression in ovaries of sexual and apomictic plants. In this work, we describe a gene named BbrizExoV with strong identity to exonuclease V (Exo V) genes from other grasses. Sequence analysis in signal prediction tools showed that BbrizExoV might have dual localization, depending on the translation point. A longer form to the nucleus and a shorter form which would be directed to the chloroplast. This is also the case for monocot sequences analyzed from other species. The long form of BbrizExoV protein localizes to the nucleus of onion epidermal cells. Analysis of ExoV proteins from dicot species, with exception of Arabidopsis thaliana ExoVL protein, showed only one localization. Using a template-based AlphaFold 2 modelling approach the structure of BbrizExoV in complex with metal and ssDNA was predicted based on the holo structure of the human counterpart. Features predicted to define ssDNA binding but a lack of sequence specificity are shared between the human enzyme and BbrizExoV. Expression analyses indicated the precise site and timing of transcript accumulation during ovule development, which coincides with the differentiation of nucelar cells to form the typical aposporic four-celled unreduced gametophyte. A putative function for this protein is proposed based on its homology and expression pattern.

Interspecies genome divergence is predominantly due to frequent small scale rearrangements in Eucalyptus

Synteny, the ordering of sequences within homologous chromosomes, must be maintained within the genomes of sexually reproducing species for the sharing of alleles and production of viable, reproducing offspring. However, when the genomes of closely related species are compared, a loss of synteny is often observed. Unequal homologous recombination is the primary mechanism behind synteny loss, occurring more often in transposon rich regions, and resulting in the formation of chromosomal rearrangements. To examine patterns of synteny among three closely related, interbreeding, and wild Eucalyptus species, we assembled their genomes using long-read DNA sequencing and de novo assembly. We identify syntenic and rearranged regions between these genomes and estimate that ~48% of our genomes remain syntenic while ~36% is rearranged. We observed that rearrangements highly fragment microsynteny. Our results suggest that synteny between these species is primarily lost through small-scale rearrangements, not through sequence loss, gain, or sequence divergence. Further examination of identified rearrangements suggests that rearrangements may be altering the phenotypes of Eucalyptus species. Our study also underscores that the use of single reference genomes in genomic variation studies could lead to reference bias, especially given the scale at which we show potentially adaptive loci have highly diverged, deleted, duplicated and/or rearranged. This study provides an unbiased framework to look at potential speciation and adaptive loci among a rapidly radiating foundation species of woodland trees that are free from selective breeding seen in most crop species.

Genome-wide analysis of LysM gene family members and their expression in response to Colletotrichum fructicola infection in Octoploid strawberry(Fragaria × ananassa)

The cultivated octoploid strawberry (Fragaria × ananassa) is an economically important fruit that is planted worldwide. The lysin motif (LysM) protein family is composed of the major class of plant pattern recognition receptors, which play important roles in sensing pathogen-associated molecular patterns (PAMPs), and subsequently triggers downstream plant immunity. In the present study, a comprehensive, genome-wide analysis of F. × ananassa LysM (FaLysM) genes was performed to investigate gene structures, phylogenic relationships, chromosome location, collinear relationships, transcription factor binding sites, and protein model analysis. We aimed to identify the LysM genes involved in the defense against plant pathogens. A total of 14 FaLysM genes were identified in the F. × ananassa genome and divided into 2 subgroups (LYP and LYK) on the basis of the phylogenetic analysis. The Ka/Ks ratio for the duplicated pair of most FaLysM genes was less than 1, which indicates that the selection pressure was mostly subject to the purifying selection during evolution. The protein model analysis revealed that FaLysM2-10 contain conserved mode of chitin binding, which suggest the potential role of FaLysM2-10 in pathogen perception and plant immunity. The RNA-Seq results showed the differential regulation of 14 FaLysM genes in response to Colletotrichum fructicola infection, implying the complex interaction between C. fructicola and strawberry. Knockout of candidate effector gene CfLysM2, which was previously proved to be highly expressed during C. fructicola infection, resulted in the up-regulation of six FaLysM genes (FaLysM1, FaLysM2, FaLysM3, FaLysM7, FaLysM8, and FaLysM12), indicating the competitive relations between CfLysM2 and FaLysM genes. Overall, this study provides fundamental information on the roles of LysM proteins in octoploid strawberry and its interaction with C. fructicola, laying useful information for further investigation on the C. fructicola-strawberry interaction and strawberry resistance breeding.

Genome-wide identification, evolution analysis of LysM gene family members and their expression analysis in response to biotic and abiotic stresses in banana (Musa L.)

LysM (Lysin motif), in response to pathogenic molecular stresses, is a crucial signal recognition gene. To understand the molecular characteristics of banana LysM gene family members, we used a series of bioinformatics methods. Based on the genomic databases of Musa acuminata, Musa balbisiana and Musa itinerans, a total of 53 genes and 55 proteins were identified, with 21 genes and 23 proteins in the M.acuminata, 16 genes and 16 proteins in each of M.balbisiana and M.itinerans, respectively. According to the conserved structural domains, LysM can be divided into five classes, namely LysM&MltD, LYK, LYP, LysMn, and LysMe. The LysM gene was relatively highly conserved in the evolution of the three genomes of banana, and some differences occurred. Expression analysis revealed that MaLysM4-5 was relatively highly expressed under high-temperature stress, low-temperature stress and pathogen infection; at the same time, about one-third of the members were down-regulated under low-temperature stress and high-temperature stress, while the expression of MaLysM10-1 and MaLysM4-5 were up-regulated. After the banana wilt fungus FocTR4 infected the banana roots, MaLysM1 was down-regulated and MaLysM11-1 was up-regulated. In conclusion, our study suggests that MaLysMs may be necessary in the response to high- and low-temperature stresses, as well as the banana wilt fungus infestation. Overall, this paper found that LysM genes may be involved in biotic and abiotic stresses in banana, and provided helpful information about LysM's evolution, expression and properties, which will provide theoretical references for further studies on the functions of LysM genes and resistance breeding in the future.

Gene fusion as an important mechanism to generate new genes in the genus Oryza

Background

Events of gene fusion have been reported in several organisms. However, the general role of gene fusion as part of new gene origination remains unknown.

Results

We conduct genome-wide interrogations of four Oryza genomes by designing and implementing novel pipelines to detect fusion genes. Based on the phylogeny of ten plant species, we detect 310 fusion genes across four Oryza species. The estimated rate of origination of fusion genes in the Oryza genus is as high as 63 fusion genes per species per million years, which is fixed at 16 fusion genes per species per million years and much higher than that in flies. By RNA sequencing analysis, we find more than 44% of the fusion genes are expressed and 90% of gene pairs show strong signals of purifying selection. Further analysis of CRISPR/Cas9 knockout lines indicates that newly formed fusion genes regulate phenotype traits including seed germination, shoot length and root length, suggesting the functional significance of these genes.

Conclusions

We detect new fusion genes that may drive phenotype evolution in Oryza. This study provides novel insights into the genome evolution of Oryza.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02696-w.