Development of Genome-Wide Intron Length Polymorphism (ILP) Markers in Tea Plant (Camellia sinensis) and Related Applications for Genetics Research.

The market value of tea is largely dependent on the tea species and cultivar. Therefore, it is important to develop efficient molecular markers covering the entire tea genome that can be used for the identification of tea varieties, marker-assisted breeding, and mapping important quantitative trait loci for beneficial traits. In this study, genome-wide molecular markers based on intron length polymorphism (ILP) were developed for tea trees. A total of 479, 1393, and 1342 tea ILP markers were identified using the PCR method in silico from the 'Shuchazao' scaffold genome, the chromosome-level genome of 'Longjing 43', and the ancient tea DASZ chromosome-level genome, respectively. A total of 230 tea ILP markers were used to amplify six tea tree species. Among these, 213 pairs of primers successfully characterize products in all six species, with 112 primer pairs exhibiting polymorphism. The polymorphism rate of primer pairs increased with the improvement in reference genome assembly quality level. The cross-species transferability analysis of 35 primer pairs of tea ILP markers showed an average amplification rate of 85.17% through 11 species in 6 families, with high transferability in Camellia reticulata and tobacco. We also used 40 pairs of tea ILP primers to evaluate the genetic diversity and population structure of C. tetracocca with 176 plants from Puan County, Guizhou Province, China. These genome-wide markers will be a valuable resource for genetic diversity analysis, marker-assisted breeding, and variety identification in tea, providing important information for the tea industry.

FKBP5 intron 7 methylation is associated with higher anxiety proneness and smaller right thalamus volume in adolescents.

Dysregulation of stress response systems may mediate the detrimental effects of childhood trauma (CT) on mental health. FKBP5 regulates glucocorticoid receptor sensitivity and exerts pleiotropic effects on intracellular signaling, neurobiology and behavior. We investigated whether CT, alone and in combination with rs1360780 genotype, is associated with altered FKBP5 methylation and whether CT-associated methylation profiles are associated with anxiety proneness (AP) and structural brain volumes. Ninety-four adolescents completed the Childhood Trauma Questionnaire, and a composite AP score was generated from the Childhood Anxiety Sensitivity Index and the State-Trait Anxiety Inventory-Trait measure. Mean methylation values for 12 regulatory regions and 25 individual CpG sites were determined using high-accuracy measurement via targeted bisulfite sequencing. FKBP5 rs1360780 genotype and structural MRI data were available for a subset of participants (n = 71 and n = 75, respectively). Regression models revealed an inverse association between methylation of three intron 7 CpG sites (35558438, 35558566 and 35558710) and right thalamus volume. CpG35558438 methylation was positively associated with AP scores. Our data indicate that an intron 7 methylation profile, consistent with lower FKBP5 expression and elevated high sensitivity glucocorticoid receptor levels, is associated with higher AP and smaller right thalamus volume. Research into the mechanisms underlying the intron 7 methylation-thalamus volume relationship, and whether it confers increased risk for long-term psychopathology by altering the regulatory threshold of stress responding, is required.

Genome-wide identification and development of miniature inverted-repeat transposable elements and intron length polymorphic markers in tea plant (Camellia sinensis).

Marker-assisted breeding and tagging of important quantitative trait loci for beneficial traits are two important strategies for the genetic improvement of plants. However, the scarcity of diverse and informative genetic markers covering the entire tea genome limits our ability to achieve such goals. In the present study, we used a comparative genomic approach to mine the tea genomes of Camellia sinensis var. assamica (CSA) and C. sinensis var. sinensis (CSS) to identify the markers to differentiate tea genotypes. In our study, 43 and 60 Camellia sinensis miniature inverted-repeat transposable element (CsMITE) families were identified in these two sequenced tea genomes, with 23,170 and 37,958 putative CsMITE sequences, respectively. In addition, we identified 4912 non-redundant, Camellia sinensis intron length polymorphic (CsILP) markers, 85.8% of which were shared by both the CSS and CSA genomes. To validate, a subset of randomly chosen 10 CsMITE markers and 15 CsILP markers were tested and found to be polymorphic among the 36 highly diverse tea genotypes. These genome-wide markers, which were identified for the first time in tea plants, will be a valuable resource for genetic diversity analysis as well as marker-assisted breeding of tea genotypes for quality improvement.

Intron retention is a stress response in sensor genes and is restored by Japanese herbal medicines: A basis for future clinical applications.

Intron retention (IR) is a regulatory mechanism that can retard protein production by acting at the level of mRNA processing. We recently demonstrated that IR occurs at the pre-symptomatic state during the aging process of a mouse model of aging, providing a promising biomarker for that state, and can be restored to the normal state by juzentaihoto (JTT), a Japanese herbal medicine (Kampo) (Okada et al. 2021). Here we characterized the genes that accumulate retained introns, examined the biological significance of increased IR in these genes for the host, and determined whether drugs other than JTT can have this effect. By analyzing RNA-sequencing data generated from the hippocampus of the 19-week-old SAMP8 mouse, a model for studying age-related depression and Alzheimer's disease, we showed that genes with increased IR are generally involved in multiple metabolic pathways and have pivotal roles in sensing homeostasis. We thus propose that IR is a stress response and works to fine-tune the expression of many downstream target genes, leading to lower levels of their translation under stress conditions. Interestingly, Kampo medicines, as well as other organic compounds, restored splicing of a specific set of retained introns in these sensor genes in accordance with the physiological recovery conditions of the host, which corresponds with the recovery of transcripts represented by differentially expressed genes. Thus, analysis of IR genes may have broad applicability in evaluating the pre-symptomatic state based on the extent of IR of selective sensor genes, opening a promising early diagnosis of any diseases and a strategy for evaluating efficacies of several drugs based on the extent of IR restoration of these sensor genes.

Discrimination of Curcuma species from Asia using intron length polymorphism markers in genes encoding diketide-CoA synthase and curcumin synthase.

Recently, Curcuma rhizome-related foods with claimed health benefits have been used worldwide; however, correct identification and quality assessment have not been conducted. Due to the wide distribution and morphological similarities of Curcuma species, the classification of some species is debated and nomenclature is inconsistent among countries. In this study, to elucidate specific molecular markers of medicinally used Curcuma species in Asia, and to solve the confusion on the reported botanical origin of crude drugs, molecular analysis based on the intron length polymorphism (ILP) in genes encoding diketide-CoA synthase and curcumin synthase and the trnK intron sequences was performed using 59 plant specimens and 42 crude drug samples from 13 Curcuma species, obtained from Asian countries. The ILP patterns of the respective species from both plant specimens and crude drug samples revealed high consistency in C. aromatica, C. zedoaria, C. phaeocaulis, C. aeruginosa, C. wenyujin, and C. zanthorrhiza, but showed intraspecies polymorphism in C. longa, C. kwangsiensis, C. amada, C. mangga and C. comosa. The C. longa specimens and samples were separated into three subgroups which were highly consistent with their geographical origins. Based on the ILP markers and the trnK intron sequences, the botanical origins of "Khamin oi" from Thailand were correctly determined to be C. longa or a hybrid between C. longa and other species, and "Wan narn kum" from Thailand and "Kasturi manjal" from India were correctly determined to be C. zanthorrhiza.

Genetic analysis of Curcuma species from Asia based on intron regions of genes encoding diketide-CoA synthase and curcumin synthase.

Intron length polymorphism (ILP) markers in genes encoding diketide-CoA synthase (DCS) and curcumin synthase (CURS) showed high identification rates in 13 Curcuma species from Asia. However, the sequences of the intron regions have not yet been analyzed. To elucidate the sequence differences in intron regions of the DCS and CURS genes and to search for specific sequences suitable for the identification of Curcuma species, a large number of sequences were determined through subcloning coupled with sequencing analysis of six Curcuma plant specimens belonging to five species that showed distinct ILP patterns. More than 30 sequences of each region from each specimen were grouped into genes DCS1, DCS2, or CURS1-3 and subsequently the sequences of the same genes were compared. Sequences belonging to the same gene showed inter-species similarity, and thus, these intron sequences were less informative within each single-gene region. The determined sequences from each specimen showed 3-5 kinds of sequence lengths in DCS intron I region, and 5-7 kinds of sequence lengths in CURS intron region. The length of determined sequences and the fragment number in each intron region were different among species, or specimens in C. longa, which were in accordance with the fragment lengths and numbers in their corresponding ILP patterns.

A dynamic intron retention program regulates the expression of several hundred genes during pollen meiosis.

KEY MESSAGE: Intron retention is a stage-specific mechanism of functional attenuation of a subset of co-regulated, functionally related genes during early stages of pollen development. To improve our understanding of the gene regulatory mechanisms that drive developmental processes, we performed a genome-wide study of alternative splicing and isoform switching during five key stages of pollen development in field mustard, Brassica rapa. Surprisingly, for several hundred genes (12.3% of the genes analysed), isoform switching results in stage-specific expression of intron-retaining transcripts at the meiotic stage of pollen development. In such cases, we report temporally regulated switching between expression of a canonical, translatable isoform and an intron-retaining transcript that is predicted to produce a truncated and presumably inactive protein. The results suggest a new pervasive mechanism underlying modulation of protein levels in a plant developmental program. The effect is not based on gene expression induction but on the type of transcript produced. We conclude that intron retention is a stage-specific mechanism of functional attenuation of a subset of co-regulated, functionally related genes during meiosis, especially genes related to ribosome biogenesis, mRNA transport and nuclear envelope architecture. We also propose that stage-specific expression of a non-functional isoform of Brassica rapa BrSDG8, a non-redundant member of histone methyltransferase gene family, linked to alternative splicing regulation, may contribute to the intron retention observed.

Multiomic analysis elucidates Complex I deficiency caused by a deep intronic variant in NDUFB10.

The diagnosis of Mendelian disorders following uninformative exome and genome sequencing remains a challenging and often unmet need. Following uninformative exome and genome sequencing of a family quartet including two siblings with suspected mitochondrial disorder, RNA sequencing (RNAseq) was pursued in one sibling. Long-read amplicon sequencing was used to determine and quantify transcript structure. Immunoblotting studies and quantitative proteomics were performed to demonstrate functional impact. Differential expression analysis of RNAseq data identified significantly decreased expression of the mitochondrial OXPHOS Complex I subunit NDUFB10 associated with a cryptic exon in intron 1 of NDUFB10, that included an in-frame stop codon. The cryptic exon contained a rare intronic variant that was homozygous in both affected siblings. Immunoblot and quantitative proteomic analysis of fibroblasts revealed decreased abundance of Complex I subunits, providing evidence of isolated Complex I deficiency. Through multiomic analysis we present data implicating a deep intronic variant in NDUFB10 as the cause of mitochondrial disease in two individuals, providing further support of the gene-disease association. This study highlights the importance of transcriptomic and proteomic analyses as complementary diagnostic tools in patients undergoing genome-wide diagnostic evaluation.

Molecular characterization and antiparasitic activity analysis of a novel piscidin 5-like type 4 from Larimichthys crocea.

Cryptocaryon irritans is an obligate parasitic ciliate protozoan that can infect various commercially important mariculture teleosts and cause high lethality and economic loss, especially Larimichthys crocea. Current methods of controlling or preventing this parasite with chemicals or antibiotics are widely considered to be environmentally harmful. The antiparasitic activity of some antimicrobial peptides (AMPs) attracted extensive attention of scholars. In the study, a novel piscidin 5-like type 4 (termed Lc-P5L4) excavated from comparative transcriptome of C. irritans - immuned L. crocea was identified and characterized. Sequence analysis shows the full-length cDNA of Lc-P5L4 is 539 bp containing an open reading frame (ORF) of 198 bp which encodes a peptide of 65 amino acid residues. The genome consists of three exons and two introns which exist in its ORF, and all the exon-intron boundaries are in accordance with classical GT-AG rule (GT/intron/AG). Multiple alignments indicate the signal peptides share highly conserved identity, while mature peptides are more diverse. Phylogenetic analysis displays Lc-P5L4 clusters together with other members of piscidin 5-like family. Next, quantitative Real-time PCR (qRT-PCR) detection found C. irritans infection could upregulate Lc-P5L4 expression level in all tested tissues significantly, it appeared earliest upregulation in the theronts infection stage in the head kidney; the expression contents reached to maximum level in the intestine, gill and muscle during trophonts falling off stage; while it was just upregulated during secondary bacterial infection stage in the liver and spleen. The data showed Lc-P5L4 upregulation time points were in accordance with different infection stages. With recombinant Lc-P5L4 (rLc-P5L4) obtained through Escherichia coli system, in vitro assay showed rLc-P5L4 could cause cilia deactivation, cell bodiesclumping and sticking to each other, then cell membrane rupture and contents leakage. The data illustrated Lc-P5L4 played critical roles in the immune defense against C. irritans infection, and provided another proof that piscidins exhibit multiple anti- C. irritans features.

LILRB1 Intron 1 Has a Polymorphic Regulatory Region That Enhances Transcription in NK Cells and Recruits YY1.

LILRB1 is a highly polymorphic receptor expressed by subsets of innate and adaptive immune cells associated with viral and autoimmune diseases and targeted by pathogens for immune evasion. LILRB1 expression on human NK cells is variegated, and the frequency of LILRB1+ cells differs among people. However, little is known about the processes and factors mediating LILRB1 transcription in NK cells. LILRB1 gene expression in lymphoid and myeloid cells arises from two distinct promoters that are separated by the first exon and intron. In this study, we identified a polymorphic 3-kb region within LILRB1 intron 1 that is epigenetically marked as an active enhancer in human lymphoid cells and not monocytes. This region possesses multiple YY1 sites, and complexes of the promoter/enhancer combination were isolated using anti-YY1 in chromatin immunoprecipitation-loop. CRISPR-mediated deletion of the 3-kb region lowers LILRB1 expression in human NKL cells. Together, these results indicate the enhancer in intron 1 binds YY1 and suggest YY1 provides a scaffold function enabling enhancer function in regulating LILRB1 gene transcription in human NK cells.

Genome-wide identification and expression analysis of the superoxide dismutase (SOD) gene family in Salvia miltiorrhiza.

Adverse environmental conditions, such as salinity, cold, drought, heavy metals, and pathogens affect the yield and quality of Salvia miltiorrhiza, a well-known medicinal plant used for the treatment of cardiovascular and cerebrovascular diseases. Superoxide dismutase (SOD), a key enzyme of antioxidant system in plants, plays a vital role in protecting plants against various biotic and abiotic stresses via scavenging the reactive oxygen species produced by organisms. However, little is known about the SOD gene family in S. miltiorrhiza. In this study, eight SOD genes, including three Cu/Zn-SODs, two Fe-SODs and three Mn-SODs, were identified in the S. miltiorrhiza genome. Their gene structures, promoters, protein features, phylogenetic relationships, and expression profiles were comprehensively investigated. Gene structure analysis implied that most SmSODs have different introns/exons distrbution patterns. Many cis-elements related to different stress responses or plant hormones were found in the promoter of each SmSOD. Expression profile analysis indicated that SmSODs exhibited diverse responses to cold, salt, drought, heavy metal, and plant hormones. Additionally, 31 types of TFs regulating SmSODs were predicted and analyzed. These findings provided valuable information for further researches on the functions and applications of SmSODs in S. miltiorrhiza growth and adaptation to stress.

Exploiting MATE efflux proteins to improve flavonoid accumulation in Camellia sinensis in silico.

Flavonoids in tea plant are the important bioactive compounds for both human health and taste quality. Multidrug and Toxic compound Extrusion (MATE) proteins could improve flavonoid accumulations by transporting and sequestering the flavonoid in vacuoles. We identified 41 putative MATE genes in tea plants. The similar intron-exon structures of tea MATEs clustered within the same gene clade. The correlation analysis of tea flavonoid and transcriptome data showed that TEA006173 might be involve in the tea flavonoid accumulation. The RT-PCR results confirmed that TEA006173 showed high expression in the young leaf tissues. Tertiary structure prediction has shown that TEA006173 contained the 12 helices with three active pockets, comprising 13 critical residues. The present study provided the structural variations and expression patterns of tea MATEs and it would be helpful for taste and nutrient quality improvement in tea plant.

Evolutionary Conservation and Expression Patterns of Neutral/Alkaline Invertases in Solanum.

The invertase gene family in plants is composed of two subfamilies of enzymes, namely, acid- and neutral/alkaline invertases (cytosolic invertase, CIN). Both can irreversibly cleave sucrose into fructose and glucose, which are thought to play key roles in carbon metabolism and plant growth. CINs are widely found in plants, but little is reported about this family. In this paper, a comparative genomic approach was used to analyze the CIN gene family in Solanum, including Solanumtuberosum, Solanumlycopersicum, Solanumpennellii, Solanumpimpinellifolium, and Solanummelongena. A total of 40 CINs were identified in five Solanum plants, and sequence features, phylogenetic relationships, motif compositions, gene structure, collinear relationship, and expression profile were further analyzed. Sequence analysis revealed a remarkable conservation of CINs in sequence length, gene number, and molecular weight. The previously verified four amino acid residues (D188, E414, Arg430, and Ser547) were also observed in 39 out of 40 CINs in our study, showing to be deeply conserved. The CIN gene family could be distinguished into groups α and ß, and α is further subdivided into subgroups α1 and α2 in our phylogenetic tree. More remarkably, each species has an average of four CINs in the α and ß groups. Marked interspecies conservation and collinearity of CINs were also further revealed by chromosome mapping. Exon-intron configuration and conserved motifs were consistent in each of these α and ß groups on the basis of in silico analysis. Expression analysis indicated that CINs were constitutively expressed and share similar expression profiles in all tested samples from S. tuberosum and S.lycopersicum. In addition, in CIN genes of the tomato and potato in response to abiotic and biotic stresses, phytohormones also performed. Overall, CINs in Solanum were encoded by a small and highly conserved gene family, possibly reflecting structural and functional conservation in Solanum. These results lay the foundation for further expounding the functional characterization of CIN genes and are also significant for understanding the evolutionary profiling of the CIN gene family in Solanum.

Genome-wide investigation of superoxide dismutase (SOD) gene family and their regulatory miRNAs reveal the involvement in abiotic stress and hormone response in tea plant (Camellia sinensis).

Superoxide dismutases (SODs), as a family of metalloenzymes related to the removal of reactive oxygen species (ROS), have not previously been investigated at genome-wide level in tea plant. In this study, 10 CsSOD genes were identified in tea plant genome, including 7 Cu/Zn-SODs (CSDs), 2 Fe-SODs (FSDs) and one Mn-SOD (MSD), and phylogenetically classified in three subgroups, respectively. Physico-chemical characteristic, conserved motifs and potential protein interaction analyses about CsSOD proteins were carried out. Exon-intron structures and codon usage bias about CsSOD genes were also examined. Exon-intron structures analysis revealed that different CsSOD genes contained various number of introns. On the basis of the prediction of regulatory miRNAs of CsSODs, a modification 5' RNA ligase-mediated (RLM)-RACE was performed and validated that csn-miR398a-3p-1 directly cleaves CsCSD4. By prediction of cis-acting elements, the expression patterns of 10 CsSOD genes and their regulatory miRNAs were detected under cold, drought, exogenous methyl jasmonate (MeJA) and gibberellin (GA3) treatments. The results showed that most of CsSODs except for CsFSD2 were induced under cold stress and CsCSDs may play primary roles under drought stress; exogenous GA3 and MeJA could also stimulated/inhibited distinct CsSODs at different stages. In addition, we found that csn-miR398a-3p-1 negatively regulated the expression of CsCSD4 may be a crucial regulatory mechanism under cold stress. This study provides a certain basis for the studies about stress resistance in tea plants, even provide insight into comprehending the classification, evolution, diverse functions and influencing factors of expression patterns for CsSOD genes.

Molecular Cloning and Functional Identification of a Squalene Synthase Encoding Gene from Alfalfa (Medicago sativa L.).

The quality of alfalfa, a main forage legume worldwide, is of great importance for the dairy industry and is affected by the content of triterpene saponins. These natural terpenoid products of triterpene aglycones are catalyzed by squalene synthase (SQS), a highly conserved enzyme present in eukaryotes. However, there is scare information on alfalfa SQS. Here, an open reading frame (ORF) of SQS was cloned from alfalfa. Sequence analysis showed MsSQS had the same exon/intron composition and shared high homology with its orthologs. Bioinformatic analysis revealed the deduced MsSQS had two transmembrane domains. When transiently expressed, GFP-MsSQS fusion protein was localized on the plasma membrane of onion epidermal cells. Removal of the C-terminal transmembrane domain of MsSQS improved solubility in Escherichia coli. MsSQS was preferably expressed in roots, followed by leaves and stems. MeJA treatment induced MsSQS expression and increased the content of total saponins. Overexpression of MsSQS in alfalfa led to the accumulation of total saponins, suggesting a correlation between MsSQS expression level with saponins content. Therefore, MsSQS is a canonical squalene synthase and contributes to saponin synthesis in alfalfa. This study provides a key candidate gene for genetic manipulation of the synthesis of triterpene saponins, which impact both plant and animal health.

Disruption of an enhancer associated with addictive behaviour within the cannabinoid receptor-1 gene suggests a possible role in alcohol intake, cannabinoid response and anxiety-related behaviour.

The cannabinoid-1 receptor (CB1) plays a critical role in a number of biological processes including nutrient intake, addiction and anxiety-related behaviour. Numerous studies have shown that expression of the gene encoding CB1 (CNR1) is highly dynamic with changes in the tissue specific expression of CNR1 associated with brain homeostasis and disease progression. However, little is known of the mechanisms regulating this dynamic expression. To gain a better understanding of the genomic mechanisms modulating the expression of CNR1 in health and disease we characterised the role of a highly conserved regulatory sequence (ECR1) in CNR1 intron 2 that contained a polymorphism in linkage disequilibrium with disease associated SNPs. We used CRISPR/CAS9 technology to disrupt ECR1 within the mouse genome. Disruption of ECR1 significantly reduced CNR1 expression in the hippocampus but not in the hypothalamus. These mice also displayed an altered sex-specific anxiety-related behavioural profile (open field test), reduced ethanol intake and a reduced hypothermic response following CB1 agonism. However, no significant changes in feeding patterns were detected. These data suggest that, whilst not all of the expression of CNR1 is modulated by ECR1, this highly conserved enhancer is required for appropriate physiological responses to a number of stimuli. The combination of comparative genomics and CRISPR/CAS9 disruption used in our study to determine the functional effects of genetic and epigenetic changes on the activity of tissue-specific regulatory elements at the CNR1 locus represent an important first step in gaining a mechanistic understanding of cannabinoid regulatory pharmacogenetics.

The complete organelle genomes of Physochlaina orientalis: Insights into short sequence repeats across seed plant mitochondrial genomes.

Short repeats (SR) play an important role in shaping seed plant mitochondrial genomes (mtDNAs). However, their origin, distribution, and relationships across the different plant lineages remain unresolved. We focus on the angiosperm family Solanaceae that shows great variation in repeat content and extend the study to a wide diversity of seed plants. We determined the complete nucleotide sequences of the organellar genomes of the medicinal plant Physochlaina orientalis (Solanaceae), member of the tribe Hyoscyameae. To understand the evolution of the P. orientalis mtDNA we made comparisons with those of five other Solanaceae. P. orientalis mtDNA presents the largest mitogenome (∼685 kb in size) among the Solanaceae and has an unprecedented 8-copy repeat family of ∼8.2 kb in length and a great number of SR arranged in tandem-like structures. We found that the SR in the Solanaceae share a common origin, but these only expanded in members of the tribe Hyoscyameae. We discuss a mechanism that could explain SR formation and expansion in P. orientalis and Hyoscyamus niger. Finally, the great increase in plant mitochondrial data allowed us to systematically extend our repeat analysis to a total of 136 seed plants to characterize and analyze for the first time families of SR among seed plant mtDNAs.

Genome-Wide Analysis of Glycoside Hydrolase Family 1 ß-glucosidase Genes in Brassica rapa and Their Potential Role in Pollen Development.

Glycoside hydrolase family 1 (GH1) ß-glucosidases (BGLUs) are encoded by a large number of genes, and are involved in many developmental processes and stress responses in plants. Due to their importance in plant growth and development, genome-wide analyses have been conducted in model plants (Arabidopsis and rice) and maize, but not in Brassica species, which are important vegetable crops. In this study, we systematically analyzed B. rapa BGLUs (BrBGLUs), and demonstrated the involvement of several genes in pollen development. Sixty-four BrBGLUs were identified in Brassica databases, which were anchored onto 10 chromosomes, with 10 tandem duplications. Phylogenetic analysis revealed that 64 genes were classified into 10 subgroups, and each subgroup had relatively conserved intron/exon structures. Clustering with Arabidopsis BGLUs (AtBGLUs) facilitated the identification of several important subgroups for flavonoid metabolism, the production of glucosinolates, the regulation of abscisic acid (ABA) levels, and other defense-related compounds. At least six BrBGLUs might be involved in pollen development. The expression of BrBGLU10/AtBGLU20, the analysis of co-expressed genes, and the examination of knocked down Arabidopsis plants strongly suggests that BrBGLU10/AtBGLU20 has an indispensable function in pollen development. The results that are obtained from this study may provide valuable information for the further understanding of ß-glucosidase function and Brassica breeding, for nutraceuticals-rich Brassica crops.

Genome-wide identification and expression analysis of the E2 gene family in potato.

E2 (ubiquitin conjugating enzymes) is an important part of the ubiquitin-proteasome pathway. These enzymes have a significant role to play during plant growth and development, which can response to various stresses. To date, the E2 family has been reported in some high plants, but the genome-wide characterization of this gene family in potato remains unknown. In the present study, 57 putative StUBCs were identified, which were clustered into eight subgroups based on phylogeny. The introns varied in numbers 0 to 9. The highest numbers of introns were 5, which accounted for 31.57%. The analysis of gene duplication showed that 22 StUBC genes were involved in 13 segmental duplication events, while no tandem duplication was found in StUBC genes. According to gene ontology analysis (GO), StUBC family major function is protein binding and ion binding. The RNA sequencing data revealed that 15 StUBC genes were highly expressed in different organs and tubers. 27 StUBC genes were up-regulated under 50 µM ABA treatments. Moreover, the RNA-seq data and qRT-PCR analysis indicated that 17 StUBC genes responded to heat stress. 8 StUBC genes responded to salt stress according to qRT-PCR analysis, and StUBC2, StUBC12, StUBC30 and StUBC13 were predominant expression. The result of this research could provide valuable information to insight into potato E2 family and establish a foundation for further to elucidate function of E2 genes.

Exon junction complex (EJC) core genes play multiple developmental roles in Physalis floridana.

KEY MESSAGE: Molecular and functional characterization of four gene families of the Physalis exon junction complex (EJC) core improved our understanding of the evolution and function of EJC core genes in plants. The exon junction complex (EJC) plays significant roles in posttranscriptional regulation of genes in eukaryotes. However, its developmental roles in plants are poorly known. We characterized four EJC core genes from Physalis floridana that were named PFMAGO, PFY14, PFeIF4AIII and PFBTZ. They shared a similar phylogenetic topology and were expressed in all examined organs. PFMAGO, PFY14 and PFeIF4AIII were localized in both the nucleus and cytoplasm while PFBTZ was mainly localized in the cytoplasm. No protein homodimerization was observed, but they could form heterodimers excluding the PFY14-PFBTZ heterodimerization. Virus-induced gene silencing (VIGS) of PFMAGO or PFY14 aborted pollen development and resulted in low plant survival due to a leaf-blight-like phenotype in the shoot apex. Carpel functionality was also impaired in the PFY14 knockdowns, whereas pollen maturation was uniquely affected in PFBTZ-VIGS plants. Once PFeIF4AIII was strongly downregulated, plant survival was reduced via a decomposing root collar after flowering and Chinese lantern morphology was distorted. The expression of Physalis orthologous genes in the DYT1-TDF1-AMS-bHLH91 regulatory cascade that is associated with pollen maturation was significantly downregulated in PFMAGO-, PFY14- and PFBTZ-VIGS flowers. Intron-retention in the transcripts of P. floridana dysfunctional tapetum1 (PFDYT1) occurred in these mutated flowers. Additionally, the expression level of WRKY genes in defense-related pathways in the shoot apex of PFMAGO- or PFY14-VIGS plants and in the root collar of PFeIF4AIII-VIGS plants was significantly downregulated. Taken together, the Physalis EJC core genes play multiple roles including a conserved role in male fertility and newly discovered roles in Chinese lantern development, carpel functionality and defense-related processes. These data increase our understanding of the evolution and functions of EJC core genes in plants.