Tandem Gene Duplication Research Articles

Genomic and structural basis for evolution of tropane alkaloid biosynthesis

The tropane alkaloids (TAs) cocaine and hyoscyamine have been used medicinally for thousands of years. To understand the evolutionary origins and trajectories of serial biosynthetic enzymes of TAs and especially the characteristic tropane skeletons, we generated the chromosome-level genome assemblies of cocaine-producing Erythroxylum novogranatense (Erythroxylaceae, rosids clade) and hyoscyamine-producing Anisodus acutangulus (Solanaceae, asterids clade). Comparative genomic and phylogenetic analysis suggested that the lack of spermidine synthase/N-methyltransferase (EnSPMT1) in ancestral asterids species contributed to the divergence of polyamine (spermidine or putrescine) methylation in cocaine and hyoscyamine biosynthesis. Molecular docking analysis and key site mutation experiments suggested that ecgonone synthases CYP81AN15 and CYP82M3 adopt different active-site architectures to biosynthesize the same product ecgonone from the same substrate in Erythroxylaceae and Solanaceae. Further synteny analysis showed different evolutionary origins and trajectories of CYP81AN15 and CYP82M3, particularly the emergence of CYP81AN15 through the neofunctionalization of ancient tandem duplication genes. The combination of structural biology and comparative genomic analysis revealed that ecgonone methyltransferase, which is responsible for the biosynthesis of characteristic 2-substituted carboxymethyl group in cocaine, evolved from the tandem copies of salicylic acid methyltransferase by the mutations of critical E216 and S153 residues. Overall, we provided strong evidence for the independent origins of serial TA biosynthetic enzymes on the genomic and structural level, underlying the chemotypic convergence of TAs in phylogenetically distant species.

Selective modes affect gene feature and function differentiation of tetraploid Brassica species in their evolution and domestication.

The genus Brassica contains a diverse group of important vegetables and oilseed crops. Genome sequencing has been completed for the six species (B. rapa, B. oleracea, B. nigra, B. carinata, B. napus, and B. juncea) in U's triangle model. The purpose of the study is to investigate whether positively and negatively selected genes (PSGs and NSGs) affect gene feature and function differentiation of Brassica tetraploids in their evolution and domestication. A total of 9,701 PSGs were found in the A, B and C subgenomes of the three tetraploids, of which, a higher number of PSGs were identified in the C subgenome as comparing to the A and B subgenomes. The PSGs of the three tetraploids had more tandem duplicated genes, higher single copy, lower multi-copy, shorter exon length and fewer exon number than the NSGs, suggesting that the selective modes affected the gene feature of Brassica tetraploids. The PSGs of all the three tetraploids enriched in a few common KEGG pathways relating to environmental adaption (such as Phenylpropanoid biosynthesis, Riboflavin metabolism, Isoflavonoid biosynthesis, Plant-pathogen interaction and Tropane, piperidine and pyridine alkaloid biosynthesis) and reproduction (Homologous recombination). Whereas, the NSGs of the three tetraploids significantly enriched in dozens of biologic processes and pathways without clear relationships with evolution. Moreover, the PSGs of B. carinata were found specifically enriched in lipid biosynthesis and metabolism which possibly contributed to the domestication of B. carinata as an oil crop. Our data suggest that selective modes affected the gene feature of Brassica tetraploids, and PSGs contributed in not only the evolution but also the domestication of Brassica tetraploids.

Phylogenetic Analyses and Transcriptional Survey Reveal the Characteristics, Evolution, and Expression Profile of NBS-Type Resistance Genes in Papaya

Carica papaya maintains an abnormally small but complete NLR family while showing weak disease resistance. To better understand their origin, evolution, and biological function, we identified 59 NLR genes via a customized RGAugury and investigated their characteristics, evolutionary history, and expression profiles based on the improved papaya genome and large-scale RNA-seq data. The results indicated that duplication is a major evolutionary force driving the formation of the papaya NLR family. Synteny analyses of papaya and other angiosperms showed that both insertion and inheritance-derived NLRs are present in papaya. Transcriptome-based expression and network analyses revealed that NLRs are actively involved in biotic stress responses. For example, a papaya-specific inserted TNL was up-regulated strongly by the fungal infection. Both transcriptome and qRT-PCR analyses confirmed the expression divergence of an RNL and an RCNL, a pair of tandem duplication genes involved in different co-expression modules. Furthermore, we observed an inserted gene cluster composed of five duplicated CNLs, showing dosage effects and functional differentiation of disease-resistance genes during evolution. This research will enhance our knowledge of the special NLR family in papaya, which may serve as a model plant for disease-resistance genetic studies.

Neofunctionalization of tandem duplicate genes encoding putative β-L-arabinofuranosidases in Arabidopsis.

Tandem duplication, one of the major types of duplication, provides the raw material for the evolution of divergent functions. In this study, we identified 1 pair of tandem duplicate genes (AT5G12950 and AT5G12960) in Arabidopsis (Arabidopsis thaliana) that originated within the last 16 million years after the split of Arabidopsis from the Capsella-Boechera ancestor. We systematically used bioinformatic tools to redefine their putative biochemical function as β-L-arabinofuranosidases that release L-Arabinose from the β-L-Araf-containing molecules in Arabidopsis. Comprehensive transcriptomic and proteomic analyses using various datasets showed divergent expression patterns among tissues between the 2 duplicate genes. We further collected phenotypic data from 2 types of measurements to indicate that AT5G12950 and AT5G12960 have different roles resulting in divergent phenotypic effects. Overall, AT5G12950 and AT5G12960 represent putative β-L-arabinofuranosidase encoding genes in Arabidopsis. After duplication, 1 duplicate copy developed diverged biological functions and contributed to a different phenotypic evolution in Arabidopsis.

Dirigent gene editing of gossypol enantiomers for toxicity-depleted cotton seeds

Axial chirality of biaryls can generate varied bioactivities. Gossypol is a binaphthyl compound made by cotton plants. Of its two axially chiral isomers, (-)-gossypol is the bioactive form in mammals and has antispermatogenic activity, and its accumulation in cotton seeds poses health concerns. Here we identified two extracellular dirigent proteins (DIRs) from Gossypium hirsutum, GhDIR5 and GhDIR6, which impart the hemigossypol oxidative coupling into (-)- and (+)-gossypol, respectively. To reduce cotton seed toxicity, we disrupted GhDIR5 by genome editing, which eliminated (-)-gossypol but had no effects on other phytoalexins, including (+)-gossypol, that provide pest resistance. Reciprocal mutagenesis identified three residues responsible for enantioselectivity. The (-)-gossypol-forming DIRs emerged later than their enantiocomplementary counterparts, from tandem gene duplications that occurred shortly after the cotton genus diverged. Our study offers insight into how plants control enantiomeric ratios and how to selectively modify the chemical spectra of cotton plants and thereby improve crop quality.

Genome-wide analysis of ELOVL family genes in the Pacific oyster genome and functional characterization of CgELOVL2/5 and genomic variations in synthesis of polyunsaturated fatty acids

Long-chain polyunsaturated fatty acids (LC-PUFAs) are essential for cellular structure and function. The elongation of very long chain fatty acids protein (ELOVL) is an enzyme involved in very long-chain fatty acids biosynthesis. In the present study, we performed a genome-wide analysis of the ELOVL family genes in the Pacific oyster Crassostrea gigas and characterized the functions of a CgELOVL2/5 gene and the genomic variations in the synthesis of PUFAs. A total of 9 CgELOVL genes were identified in four chromosomes of C. gigas, and could be grouped into five subfamilies. The CgELOVL genes were all involved in oyster larval development and four involved in response to external stress. Two tandem-duplicated genes showed nearly identical expression profiles and appear essential for the early embryonic development of oysters. Two genes (CgELOVL2/5, CgELOVL1/7) were differentially expressed in oysters with high and low PUFA contents, and CgELOVL2/5 could extend the fatty acids with C18:3n-3, C18:2n-6, and C20:5n-3 as substrates. In the CgELOVL2/5 promoter region, three SNPs and one indel were associated with the PUFA contents and had substantial impacts on the transcriptional expression of CgELOVL2/5. These genomic variations may lead to the PUFA content difference among oysters to some degree and showed potential in breeding applications. This study could contribute to research on the genetic dissection of oyster PUFA metabolism and provide clues for the molecular-based selection of high-PUFA Pacific oysters.

Genome-Wide Mining of the Tandem Duplicated Type III Polyketide Synthases and Their Expression, Structure Analysis of Senna tora

Senna tora is one of the homologous crops used as a medicinal food containing an abundance of anthraquinones. Type III polyketide synthases (PKSs) are key enzymes that catalyze polyketide formation; in particular, the chalcone synthase-like (CHS-L) genes are involved in anthraquinone production. Tandem duplication is a fundamental mechanism for gene family expansion. However, the analysis of the tandem duplicated genes (TDGs) and the identification and characterization of PKSs have not been reported for S. tora. Herein, we identified 3087 TDGs in the S. tora genome; the synonymous substitution rates (Ks) analysis indicated that the TDGs had recently undergone duplication. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the type III PKSs were the most enriched TDGs involved in the biosynthesis of the secondary metabolite pathways, as evidenced by 14 tandem duplicated CHS-L genes. Subsequently, we identified 30 type III PKSs with complete sequences in the S. tora genome. Based on the phylogenetic analysis, the type III PKSs were classified into three groups. The protein conserved motifs and key active residues showed similar patterns in the same group. The transcriptome analysis showed that the chalcone synthase (CHS) genes were more highly expressed in the leaves than in the seeds in S. tora. The transcriptome and qRT-PCR analysis showed that the CHS-L genes had a higher expression in the seeds than in other tissues, particularly seven tandem duplicated CHS-L2/3/5/6/9/10/13 genes. The key active-site residues and three-dimensional models of the CHS-L2/3/5/6/9/10/13 proteins showed slight variation. These results indicated that the rich anthraquinones in S. tora seeds might be ascribed to the PKSs' expansion from tandem duplication, and the seven key CHS-L2/3/5/6/9/10/13 genes provide candidate genes for further research. Our study provides an important basis for further research on the regulation of anthraquinones' biosynthesis in S. tora.

Genome-Wide Identification and Expression Profiling of Glutathione S-Transferase Gene Family in Foxtail Millet (Setaria italica L.)

Glutathione S-transferases (GSTs) are a critical superfamily of multifunctional enzymes in plants. As a ligand or binding protein, GSTs regulate plant growth and development and detoxification. Foxtail millet (Setaria italica (L.) P. Beauv) could respond to abiotic stresses through a highly complex multi-gene regulatory network in which the GST family is also involved. However, GST genes have been scarcely studied in foxtail millet. Genome-wide identification and expression characteristics analysis of the foxtail millet GST gene family were conducted by biological information technology. The results showed that 73 GST genes (SiGSTs) were identified in the foxtail millet genome and were divided into seven classes. The chromosome localization results showed uneven distribution of GSTs on the seven chromosomes. There were 30 tandem duplication gene pairs belonging to 11 clusters. Only one pair of SiGSTU1 and SiGSTU23 were identified as fragment duplication genes. A total of ten conserved motifs were identified in the GST family of foxtail millet. The gene structure of SiGSTs is relatively conservative, but the number and length of exons of each gene are still different. The cis-acting elements in the promoter region of 73 SiGST genes showed that 94.5% of SiGST genes possessed defense and stress-responsive elements. The expression profiles of 37 SiGST genes covering 21 tissues suggested that most SiGST genes were expressed in multiple organs and were highly expressed in roots and leaves. By qPCR analysis, we found that 21 SiGST genes were responsive to abiotic stresses and abscisic acid (ABA). Taken together, this study provides a theoretical basis for identifying foxtail millet GST family information and improving their responses to different stresses.

Upregulation of tandem duplicated BoFLC1 genes is associated with the non-flowering trait in Brassica oleracea var. capitata.

Tandem duplicated BoFLC1 genes(BoFLC1a and BoFLC1b), which were identified as the candidate causal genes for the non-flowering trait in the cabbage mutant 'nfc', were upregulated during winter in 'nfc'. The non-flowering natural cabbage mutant 'nfc' was discovered from the breeding line 'T15' with normal flowering characteristics. In this study, we investigated the molecular basis underlying the non-flowering trait of 'nfc'. First, 'nfc' was induced to flower using the grafting floral induction method, and three F2 populations were generated. The flowering phenotype of each F2 population was widely distributed with non-flowering individuals appearing in two populations. QTL-seq analysis detected a genomic region associated with flowering date at approximately 51Mb on chromosome 9 in two of the three F2 populations. Subsequent validation and fine mapping of the candidate genomic region using QTL analysis identified the quantitative trait loci (QTL) at 50,177,696-51,474,818bp on chromosome 9 covering 241 genes. Additionally, RNA-seq analysis in leaves and shoot apices of 'nfc' and 'T15' plants identified 19 and 15 differentially expressed genes related to floweringtime, respectively. Based on these results, we identified tandem duplicated BoFLC1 genes, which are homologs of thefloral repressor FLOWERING LOCUS C, as the candidate genes responsible for the non-flowering trait of 'nfc'. We designated the tandem duplicated BoFLC1 genes as BoFLC1a and BoFLC1b. Expression analysis revealed that the expression levels of BoFLC1a and BoFLC1b were downregulated during winter in 'T15' but were upregulated and maintained during winter in 'nfc'. Additionally, the expression level of the floral integrator BoFT was upregulated in the spring in 'T15' but hardly upregulated in 'nfc'. These results suggest that the upregulated levels of BoFLC1a and BoFLC1b contributed to the non-flowering trait of 'nfc'.

The high-quality genome of lotus reveals tandem duplicate genes involved in stress response and secondary metabolites biosynthesis.

The high-quality genome of lotus reveals tandem duplicate genes involved in stress response and secondary metabolites biosynthesis.

Genome-wide identification of papain-like cysteine protease family genes in cultivated peanut (Arachis hypogaea L.) and functional characterization of AhRD21B in response to chilling stress

Papain-like cysteine proteases (PLCPs) play essential roles in plant growth, senescence, and stress responses. However, researches on the PLCP genes are limited in cultivated peanut (Arachis hypogaea L.). Here, we identified and analyzed AhPLCP genes in cultivated peanut using bioinformatic approaches and verified the involvement of AhRD21B in response to chilling tolerance using an experimental test. We identified 52 AhPLCP genes in peanut genome and categorized them into eight subfamilies. Six tandem duplicated genes were identified in cultivated peanut. The promoters of AhPLCP genes contained many stress responsive elements. AhPLCP genes were generally expressed constitutively, and many AhPLCP genes were cold-responsive. Among these genes, heterologous expression of AhRD21B significantly increased the tolerance of transgenic Arabidopsis to cold compared with the wild-type plants. In addation, AhRD21B can regulate cold-related genes that participate in enhancing chilling tolerance in transgenic Arabidopsis. These findings provide valuable insights into the functions of AhPLCP genes in plant stress tolerance.

Genome-Wide Identification and Characterization of R2R3-MYB Provide Insight into Anthocyanin Biosynthesis Regulation Mechanism of Ananas comosus var. bracteatus.

The R2R3-MYB proteins comprise the largest class of MYB transcription factors, which play an essential role in regulating anthocyanin synthesis in various plant species. Ananas comosus var. bracteatus is an important colorful anthocyanins-rich garden plant. The spatio-temporal accumulation of anthocyanins in chimeric leaves, bracts, flowers, and peels makes it an important plant with a long ornamental period and highly improves its commercial value. We conducted a comprehensive bioinformatic analysis of the R2R3-MYB gene family based on genome data from A. comosus var. bracteatus. Phylogenetic analysis, gene structure and motif analysis, gene duplication, collinearity, and promoter analysis were used to analyze the characteristics of this gene family. In this work, a total of 99 R2R3-MYB genes were identified and classified into 33 subfamilies according to phylogenetic analysis, and most of them were localized in the nucleus. We found these genes were mapped to 25 chromosomes. Gene structure and protein motifs were conserved among AbR2R3-MYB genes, especially within the same subfamily. Collinearity analysis revealed four pairs of tandem duplicated genes and 32 segmental duplicates in AbR2R3-MYB genes, indicating that segmental duplication contributed to the amplification of the AbR2R3-MYB gene family. A total of 273 ABRE responsiveness, 66 TCA elements, 97 CGTCA motifs, and TGACG motifs were the main cis elements in the promoter region under response to ABA, SA, and MEJA. These results revealed the potential function of AbR2R3-MYB genes in response to hormone stress. Ten R2R3-MYBs were found to have high homology to MYB proteins reported to be involved in anthocyanin biosynthesis from other plants. RT-qPCR results revealed the 10 AbR2R3-MYB genes showed tissue-specific expression patterns, six of them expressed the highest in the flower, two genes in the bract, and two genes in the leaf. These results suggested that these genes may be the candidates that regulate anthocyanin biosynthesis of A. comosus var. bracteatus in the flower, leaf, and bract, respectively. In addition, the expressions of these 10 AbR2R3-MYB genes were differentially induced by ABA, MEJA, and SA, implying that these genes may play crucial roles in hormone-induced anthocyanin biosynthesis. Our study provided a comprehensive and systematic analysis of AbR2R3-MYB genes and identified the AbR2R3-MYB genes regulating the spatial-temporal anthocyanin biosynthesis in A. comosus var. bracteatus, which would be valuable for further study on the anthocyanin regulation mechanism of A. comosus var. bracteatus.

Apiaceae FNS I originated from F3H through tandem gene duplication.

Flavonoids are specialized metabolites with numerous biological functions in stress response and reproduction of plants. Flavones are one subgroup that is produced by the flavone synthase (FNS). Two distinct enzyme families evolved that can catalyze the biosynthesis of flavones. While the membrane-bound FNS II is widely distributed in seed plants, one lineage of soluble FNS I appeared to be unique to Apiaceae species. We show through phylogenetic and comparative genomic analyses that Apiaceae FNS I evolved through tandem gene duplication of flavanone 3-hydroxylase (F3H) followed by neofunctionalization. Currently available datasets suggest that this event happened within the Apiaceae in a common ancestor of Daucus carota and Apium graveolens. The results also support previous findings that FNS I in the Apiaceae evolved independent of FNS I in other plant species. We validated a long standing hypothesis about the evolution of Apiaceae FNS I and predicted the phylogenetic position of this event. Our results explain how an Apiaceae-specific FNS I lineage evolved and confirm independence from other FNS I lineages reported in non-Apiaceae species.

Application of Cloning-Free Genome Engineering to Escherichia coli.

The propagation of foreign DNA in Escherichia coli is central to molecular biology. Recent advances have dramatically expanded the ability to engineer (bacterial) cells; however, most of these techniques remain time-consuming. The aim of the present work was to explore the possibility to use the cloning-free genome editing (CFGE) approach, proposed by Döhlemann and coworkers (2016), for E. coli genetics, and to deepen the knowledge about the homologous recombination mechanism. The E. coli auxotrophic mutant strains FB182 (hisF892) and FB181 (hisI903) were transformed with the circularized wild-type E. coli (i) hisF gene and hisF gene fragments of decreasing length, and (ii) hisIE gene, respectively. His+ clones were selected based on their ability to grow in the absence of histidine, and their hisF/hisIE gene sequences were characterized. CFGE method allowed the recombination of wild-type his genes (or fragments of them) within the mutated chromosomal copy, with a different recombination frequency based on the fragment length, and the generation of clones with a variable number of in tandem his genes copies. Data obtained pave the way to further evolutionary studies concerning the homologous recombination mechanism and the fate of in tandem duplicated genes.

Genome-wide identification glutathione-S-transferase gene superfamily in Daphnia pulex and its transcriptional response to nanoplastics

Glutathione S-transferases (GSTs) are key multifunctional phase II detoxification enzymes involved in the regulation of growth, development, and stress responses. However, the knowledge of GSTs in the model invertebrate organism Daphnia pulex at the genomic level remains limited. In the present study, 35 GST genes were identified in D. pulex (Dp-GST), belonging to eight subfamilies, with the sigma, mu, and delta/epsilon subfamilies constituting approximately 29 %, 20 %, and 20 % of the GST superfamily, respectively. Chromosome tandem duplication of genes within the same subfamily was observed, which may be the main force driving GST expansion in D. pulex. DpGST genes showed different expression patterns in response to nanoplastic exposure for 96 h and 21 days. Some homologous GST genes in D. pulex showed similar expression patterns in response to nanoplastic exposure, likely owing to their unique motifs. For example, motif 9 is found in all delta/epsilon GST genes, whereas motifs 1, 2, 3, 5, and 7 are highly conserved in sigma GST genes. The characterization of D. pulex GSTs extending the knowledge of GST-mediated environmental contaminants, especially nanoplastics.

The ubiquitous position effect, synergistic effect of recent generated tandem duplicated genes in grapevine, and their co-response and overactivity to biotic stress

The ubiquitous position effect, synergistic effect of recent generated tandem duplicated genes in grapevine, and their co-response and overactivity to biotic stress

Genome-wide characterization of the inositol transporters gene family in Populus and functional characterization of PtINT1b in response to salt stress

Inositol transporters (INTs) can mediate the transmembrane transport of inositol, and play crucial roles in plant growth, development and stress resistance. However, the INT gene family in Populus has not been reported. Herein, nine INT genes were identified in the Populus trichocarpa genome and divided into three clades. Tandem duplication and whole-genome duplication events could induce the expansion of PtINT gene family. It was worth noting that PtINT1c* and 1d* formed by twice tandem gene duplication events of PtINT1b, but both had undergone partial structural loss during evolution. PtINT2_p1* and PtINT2_p2* might be originated from one INT2 gene by stop codon- and start codon-gain variants. Different members of PtINTs were localized to the plasma membrane or vacuolar membrane. PtINTs had diversified tissue expression profiles, and many members were significantly induced or suppressed after salt and drought treatments. PtINT1b was induced by drought and salinity stresses, and encoded a vacuolar inositol transporter. Overexpression of PtINT1b rendered the transgenic Arabidopsis plants more resistant to salt stress. In conclusion, this study provides valuable clues for future research on the function of PtINTs, and PtINT1b was identified as a candidate gene for genetic engineering to enhance salinity tolerance in plants.

Genome-wide identification, characterization, evolution, and expression pattern analyses of the typical thioredoxin gene family in wheat (Triticum aestivum L.).

Typical thioredoxin (TRX) plays an important role in maintaining redox balance in plants. However, the typical TRX genes in wheat still need to be comprehensively and deeply studied. In this research, a total of 48 typical TaTRX genes belonging to eight subtypes were identified via a genome-wide search in wheat, and the gene structures, protein conserved motifs, and protein 3D structures of the same subtype were very similar. Evolutionary analysis showed that there are two pairs of tandem duplication genes and 14 clusters of segmental duplication genes in typical TaTRX family members; TaTRX15, TaTRX36, and TaTRX42 had positive selection compared with the orthologs of their ancestral species; rice and maize have 11 and 13 orthologous typical TRXs with wheat, respectively. Gene Ontology (GO) analysis indicated that typical TaTRXs were involved in maintaining redox homeostasis in wheat cells. Estimation of ROS content, determination of antioxidant enzyme activity, and gene expression analysis in a line overexpressing one typical TaTRX confirmed that TRX plays an important role in maintaining redox balance in wheat. A predictive analysis of cis-acting elements in the promoter region showed that typical TaTRXs were extensively involved in various hormone metabolism and response processes to stress. The results predicted using public databases or verified using RT-qPCR show that typical TaTRXs were able to respond to biotic and abiotic stresses, and their expression in wheat was spatiotemporal. A total of 16 wheat proteins belonging to four different families interacting with typical TaTRXs were predicted. The above comprehensive analysis of typical TaTRX genes can enrich our understanding of this gene family in wheat and provide valuable insights for further gene function research.

Genome-wide identification, characterization, and expression profile of NBS-LRR gene family in sweet orange (Citrus sinensis)

BackgroundThe NBS-LRR (nucleotide-binding site-leucine-rich repeat gene) gene family, known as the plant R (resistance) gene family with the most members, plays a significant role in plant resistance to various external adversity stresses. The NBS-LRR gene family has been researched in many plant species. Citrus is one of the most vital global cash crops, the number one fruit group, and the third most traded agricultural product world wild. However, as one of the largest citrus species, a comprehensive study of the NBS-LRR gene family has not been reported on sweet oranges. MethodsIn this study, NBS-LRR genes were identified from the Citrus sinensis genome (v3.0), with a comprehensive analysis of this gene family performed, including phylogenetic analysis, gene structure, cis-acting element of a promoter, and chromosomal localization, among others. The expression pattern of NBS-LRR genes was analyzed when sweet orange fruits were infected by Penicillium digitatum, employing experimental data from our research group. It first reported the expression patterns of NBS-LRR genes under abiotic stresses, using three transcript data from NCBI (National Center for Biotechnology Information). ResultsIn this study, 111 NBS-LRR genes were identified in the C. sinensis genome (v3.0) and classified into seven subfamilies according to their N-terminal and C-terminal domains. The phylogenetic tree results indicate that genes containing only the NBS structural domain are more ancient in the sweet orange NBS-LRR gene family. The chromosome localization results showed that 111 NBS-LRR genes were distributed unevenly on nine chromosomes, with the most genes distributed on chromosome 1. In addition, we identified a total of 18 tandem duplication gene pairs in the sweet orange NBS-LRR gene family, and based on the Ka/Ks ratio, all of the tandem duplication genes underwent purifying selection. Transcriptome data analysis showed a significant number of NBS-LRR genes expressed under biotic and abiotic stresses, and some reached significantly different levels of expression. It indicates that the NBS-LRR gene family is vital in resistance to biotic and abiotic stresses in sweet oranges. ConclusionOur study provides the first comprehensive framework on the NBS-LRR family of genes, which provides a basis for further in-depth studies on the biological functions of NBS-LRR in growth, development, and response to abiotic stresses in sweet orange.

Genome-wide identification and functional analysis of the TCP gene family in rye (Secale cereale L.)

TEOSINTE BRANCHED1/CYCLOIDEA/PCF (TCP) proteins are plant-specific transcription factors that play significant roles in plant growth, development, and stress response. Rye is a high-value crop with strong resistance to adverse environments. However, the functions of TCP proteins in rye are rarely reported. Based on a genome-wide analysis, the present study identified 26 TCP genes (ScTCPs) in rye. Mapping showed an uneven distribution of the ScTCP genes on the seven rye chromosomes and detected three pairs of tandem duplication genes. Phylogenetic analysis divided these genes into PCF (Proliferrating Cell Factors), CIN (CINCINNATA), and CYC (CYCLOIDEA)/TB1 (Teosinte Branched1) classes, which showed the highest homology between rye and wheat genes. Analysis of miRNA targeting sites indicated that five ScTCP genes were identified as potential targets of miRNA319. Promoter cis-acting elements analysis indicated that ScTCPs were regulated by light signals. Further analysis of the gene expression patterns and functional annotations suggested the role of a few ScTCPs in grain development and stress response. In addition, two TB1 homologous genes (ScTCP9 and ScTCP10) were identified in the ScTCP family. Synteny analysis showed that TB1 orthologous gene pairs existed before the ancestral divergence. Finally, the yeast two-hybrid assay and luciferase complementation imaging assay proved that ScTCP9, localized in the nucleus, interacts with ScFT (Flowering locus T), indicating their role in regulating flowering time. Taken together, this comprehensive study of ScTCPs provides important information for further research on gene function and crop improvement.