Identification of Quantitative Trait Loci and Candidate Genes Controlling Seed Dormancy in Eggplant (Solanum melongena L.).

Seed dormancy is a life adaptation trait exhibited by plants in response to environmental changes during their growth and development. The dormancy of commercial seeds is the key factor affecting seed quality. Eggplant seed dormancy is controlled by quantitative trait loci (QTLs), but reliable QTLs related to eggplant dormancy are still lacking. In this study, F2 populations obtained through the hybridization of paternally inbred lines with significant differences in dormancy were used to detect regulatory sites of dormancy in eggplant seeds. Three QTLs (dr1.1, dr2.1, and dr6.1) related to seed dormancy were detected on three chromosomes of eggplant using the QTL-Seq technique. By combining nonsynonymous sites within the candidate regions and gene functional annotation analysis, nine candidate genes were selected from three QTL candidate regions. According to the germination results on the eighth day, the male parent was not dormant, but the female parent was dormant. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to verify the expression of nine candidate genes, and the Smechr0201082 gene showed roughly the same trend as that in the phenotypic data. We proposed Smechr0201082 as the potential key gene involved in regulating the dormancy of eggplant seeds. The results of seed experiments with different concentrations of gibberellin A3 (GA3) showed that, within a certain range, the higher the gibberellin concentration, the earlier the emergence and the higher the germination rate. However, higher concentrations of GA3 may have potential effects on eggplant seedlings. We suggest the use of GA3 at a concentration of 200-250 mg·L-1 to treat dormant seeds. This study provides a foundation for the further exploration of genes related to the regulation of seed dormancy and the elucidation of the molecular mechanism of eggplant seed dormancy and germination.

Construction of SNP fingerprints and genetic diversity analysis of radish (Raphanus sativus L.).

Radish (Raphanus sativus L.) is a vegetable crop with economic value and ecological significance in the genus Radish, family Brassicaceae. In recent years, developed countries have attached great importance to the collection and conservation of radish germplasm resources and their research and utilization, but the lack of population genetic information and molecular markers has hindered the development of the genetic breeding of radish. In this study, we integrated the radish genomic data published in databases for the development of single-nucleotide polymorphism (SNP) markers, and obtained a dataset of 308 high-quality SNPs under strict selection criteria. With the support of Kompetitive Allele-Specific PCR (KASP) technology, we screened a set of 32 candidate core SNP marker sets to analyse the genetic diversity of the collected 356 radish varieties. The results showed that the mean values of polymorphism information content (PIC), minor allele frequency (MAF), gene diversity and heterozygosity of the 32 candidate core SNP markers were 0.32, 0.30, 0.40 and 0.25, respectively. Population structural analysis, principal component analysis and genetic evolutionary tree analysis indicated that the 356 radish materials were best classified into two taxa, and that the two taxa of the material were closely genetically exchanged. Finally, on the basis of 32 candidate core SNP markers we calculated 15 core markers using a computer algorithm to construct a fingerprint map of 356 radish varieties. Furthermore, we constructed a core germplasm population consisting of 71 radish materials using 32 candidate core markers. In this study, we developed SNP markers for radish cultivar identification and genetic diversity analysis, and constructed DNA fingerprints, providing a basis for the identification of radish germplasm resources and molecular marker-assisted breeding as well as genetic research.

Metabolic and Transcriptomic Analyses Reveal Different Metabolite Biosynthesis Profiles between Purple and Green Pak Choi.

Pak choi is one of the most important leafy vegetables planted in East Asia and provides essential nutrients for the human body. Purple pak choi differs mainly in leaf colour but exhibits distinct nutritional profiles from green pak choi. In this study, we performed metabolic and transcriptomic analyses to uncover the mechanisms underlying the differences in metabolite biosynthesis profiles between the two pak choi varieties. Metabolite profiling revealed significant differences in the levels of metabolites, mainly amino acids and their derivatives and flavonoids. Furthermore, 34 flavonoids significantly differed between green and purple pak choi leaves, and cyanidin and its derivative anthocyanins were abundant in purple pak choi. In addition, we found that the structural genes CHS, DFR, ANS, and UGT75C1, as well as the transcription factor MYB2, play a major role in anthocyanin synthesis. These results provide insight into the molecular mechanisms underlying leaf pigmentation in pak choi and offer a platform for assessing related varieties.

Metabolomic and Transcriptomic Analyses Reveal the Effects of Grafting on Nutritional Properties in Eggplant.

Grafting has a significant impact on the botany properties, commercial character, disease resistance, and productivity of eggplants. However, the mechanism of phenotypic modulation on grafted eggplants is rarely reported. In this study, a widely cultivated eggplant (Solanum. melongena cv. 'Zheqie No.10') was selected as the scion and grafted, respectively, onto four rootstocks of TOR (S. torvum), Sa (S. aculeatissimum), SS (S. sisymbriifolium), and Sm64R (S. melongena cv. 'Qiezhen No. 64R') for phenotypic screening. Physiological and biochemical analysis showed the rootstock Sm64R could improve the fruit quality with the increasing of fruit size, yield, and the contents of total soluble solid, phenolic acid, total amino acid, total sugar, and vitamin C. To further investigate the improvement of fruit quality on Sm64R, a transcriptome and a metabolome between the Sm64R-grafted eggplant and self-grafted eggplant were performed. Significant differences in metabolites, such as phenolic acids, lipids, nucleotides and derivatives, alkaloids, terpenoids, and amino acids, were observed. Differential metabolites and differentially expressed genes were found to be abundant in three core pathways of nutritional qualities, including biosynthesis of phenylpropanoids, phospholipids, and nucleotide metabolism. Thus, this study may provide a novel insight into the effects of grafting on the fruit quality in eggplant.

Understanding the nutraceutical diversity through a comparative analysis of the taproot metabolomes of different edible radish types via UHPLC-Q-TOF-MS.

Radishes are root vegetables that are rich in bioactive compounds and provide numerous health benefits, but the overall metabolic profiles of radish taproots and the metabolic differences among different edible types are not fully understood. In this research, we used UHPLC-Q-TOF-MS to identify the metabolites in cooked, processed and fruit radishes of ten varieties. In total, 264 metabolites belonging to 18 categories were detected. A multivariate analysis revealed that the metabolite composition differed among the three radish groups, and a comparative analysis showed that the significantly differentially accumulated metabolites were mainly amino acids and derivatives, lipids, flavonoids, hydroxycinnamate derivatives and carbohydrates. The accumulation of metabolites, particularly flavonoids, was greater in fruit radishes than in cooked and processed radishes. This work provides novel insights into the radish metabolomic profiles for assessment of the nutritional value of different edible radish types for humans.

Characterization of repetitive sequences in Dendrobium officinale and comparative chromosomal structures in Dendrobium species using FISH.

Tandem repeats are one of the most conserved features in the eukaryote genomes. Dendrobium is the third largest genus in family Orchidaceae compromising over 1,200 species. However, the organization of repetitive sequences in Dendrobium species remains unclear. In this study, we performed the identification and characterization of the tandem repeats in D. officinale genome using graph-based clustering and Fluorescence in situ hybridization (FISH). Six major clusters including five satellite DNAs (DofSat1-5) and one 5S rDNA repeat (Dof5S) were identified as tandem repeats. The tandem organization of DofSat5 was verified by PCR amplification and southern blotting. The chromosomal locations of the repetitive DNAs in D. officinale were investigated by FISH using the tandem repeats and oligos probes. The results showed that each of the DofSat5, 5S and 45S rDNA had one pair of strong signals on D. officinale chromosomes. The distribution of repetitive DNAs along chromosomes was also investigated based on genomic in situ hybridization (GISH) among four Dendrobium species. The results suggested complex chromosomal fusion/segmentation and rearrangements during the evolution of Dendrobium species. In conclusion, the present study provides new landmarks for unequival differentiation of the Dendrobium chromosomes and facilitate the understanding the chromosome evolution in Dendrobium speceis.

Shotgun metagenomic sequencing reveals skin microbial variability from different facial sites.

Biogeography (body site) is known to be one of the main factors influencing the composition of the skin microbial community. However, site-associated microbial variability at a fine-scale level was not well-characterized since there was a lack of high-resolution recognition of facial microbiota across kingdoms by shotgun metagenomic sequencing. To investigate the explicit microbial variance in the human face, 822 shotgun metagenomic sequencing data from Han Chinese recently published by our group, in combination with 97 North American samples from NIH Human Microbiome Project (HMP), were reassessed. Metagenomic profiling of bacteria, fungi, and bacteriophages, as well as enriched function modules from three facial sites (forehead, cheek, and the back of the nose), was analyzed. The results revealed that skin microbial features were more alike in the forehead and cheek while varied from the back of the nose in terms of taxonomy and functionality. Analysis based on biogeographic theories suggested that neutral drift with niche selection from the host could possibly give rise to the variations. Of note, the abundance of porphyrin-producing species, i.e., Cutibacterium acnes, Cutibacterium avidum, Cutibacterium granulosum, and Cutibacterium namnetense, was all the highest in the back of the nose compared with the forehead/cheek, which was consistent with the highest porphyrin level on the nose in our population. Sequentially, the site-associated microbiome variance was confirmed in American populations; however, it was not entirely consistent. Furthermore, our data revealed correlation patterns between Propionibacterium acnes bacteriophages with genus Cutibacterium at different facial sites in both populations; however, C. acnes exhibited a distinct correlation with P. acnes bacteriophages in Americans/Chinese. Taken together, in this study, we explored the fine-scale facial site-associated changes in the skin microbiome and provided insight into the ecological processes underlying facial microbial variations.

The New Variation in the Promoter Region of FLOWERING LOCUS T Is Involved in Flowering in Brassica rapa.

Flowering time is an important agronomic trait in Brassica rapa and has a wide range of variation. The change from vegetative to reproductive development is a major transition period, especially in flowering vegetable crops. In this study, two non-heading Chinese cabbage varieties with significantly different flowering times, Pak-choi (B. rapa var. communis Tesn et Lee) and Caitai (B. rapa var. tsaitai Hort.), were used to construct segregated F2 populations. The bulk-segregant approach coupled with whole genome re-sequencing was used for QTL sequencing (QTL-seq) analysis to map flowering time traits. The candidate genes controlling flowering time in B. rapa were predicted by homologous gene alignment and function annotation. The major-effect QTL ft7.1 was detected on chromosome A07 of B. rapa, and the FT family gene BrFT was predicted as the candidate gene. Moreover, a new promoter regional difference of 1577 bp was revealed by analyzing the sequence of the BrFT gene. The promoter region activity analysis and divergent gene expression levels indicated that the difference in the promoter region may contribute to different flowering times. These findings provide insights into the mechanisms underlying the flowering time in Brassica and the candidate genes regulating flowering in production.

Species-Specific Gene Expansion of the Cellulose synthase Gene Superfamily in the Orchidaceae Family and Functional Divergence of Mannan Synthesis-Related Genes in Dendrobium officinale.

Plant Cellulose synthase genes constitute a supergene family that includes the Cellulose synthase (CesA) family and nine Cellulose synthase-like (Csl) families, the members of which are widely involved in the biosynthesis of cellulose and hemicellulose. However, little is known about the Cellulose synthase superfamily in the family Orchidaceae, one of the largest families of angiosperms. In the present study, we identified and systematically analyzed the CesA/Csl family members in three fully sequenced Orchidaceae species, i.e., Dendrobium officinale, Phalaenopsis equestris, and Apostasia shenzhenica. A total of 125 Cellulose synthase superfamily genes were identified in the three orchid species and classified into one CesA family and six Csl families: CslA, CslC, CslD, CslE, CslG, and CslH according to phylogenetic analysis involving nine representative plant species. We found species-specific expansion of certain gene families, such as the CslAs in D. officinale (19 members). The CesA/Csl families exhibited sequence divergence and conservation in terms of gene structure, phylogeny, and deduced protein sequence, indicating multiple origins via different evolutionary processes. The distribution of the DofCesA/DofCsl genes was investigated, and 14 tandemly duplicated genes were detected, implying that the expansion of DofCesA/DofCsl genes may have originated via gene duplication. Furthermore, the expression profiles of the DofCesA/DofCsl genes were investigated using transcriptome sequencing and quantitative Real-time PCR (qRT-PCR) analysis, which revealed functional divergence in different tissues and during different developmental stages of D. officinale. Three DofCesAs were highly expressed in the flower, whereas DofCslD and DofCslC family genes exhibited low expression levels in all tissues and at all developmental stages. The 19 DofCslAs were differentially expressed in the D. officinale stems at different developmental stages, among which six DofCslAs were expressed at low levels or not at all. Notably, two DofCslAs (DofCslA14 and DofCslA15) showed significantly high expression in the stems of D. officinale, indicating a vital role in mannan synthesis. These results indicate the functional redundancy and specialization of DofCslAs with respect to polysaccharide accumulation. In conclusion, our results provide insights into the evolution, structure, and expression patterns of CesA/Csl genes and provide a foundation for further gene functional analysis in Orchidaceae and other plant species.

Comparative transcriptome analysis reveals distinct responsive biological processes in radish genotypes contrasting for Plasmodiophora brassicae interaction.

Plasmodiophora brassicae is a protozoan pathogen that causes clubroot disease, which is one of the most destructive diseases for Brassica crops, including radish. However, little is known about the molecular mechanism of clubroot resistance in radish. In this study, we performed a comparative transcriptome analysis between resistant and susceptible radish inoculated with P. brassicae. More differentially expressed genes (DEGs) were identified at 28 days after inoculation (DAI) compared to 7 DAI in both genotypes. Gene ontology (GO) and KEGG enrichment indicated that stress/defense response, secondary metabolic biosynthesis, hormone metabolic process, and cell periphery are directly involved in the defense response process. Further analysis of the transcriptome revealed that effector-triggered immunity (ETI) plays key roles in the defense response. The plant hormones jasmonic acid (JA), ethylene (ET), and abscisic acid (ABA) related genes are activated in clubroot defense in the resistant line. Auxin (AUX) hormone related genes are activated in the developing galls of susceptible radish. Our study provides a global transcriptional overview for clubroot development for insights into the P. brassicae defense mechanisms in radish.

Phylogenetic relationships of sucrose transporters (SUTs) in plants and genome-wide characterization of SUT genes in Orchidaceae reveal roles in floral organ development.

Sucrose is the primary form of photosynthetically produced carbohydrates transported long distance in many plant species and substantially affects plant growth, development and physiology. Sucrose transporters (SUTs or SUCs) are a group of membrane proteins that play vital roles in mediating sucrose allocation within cells and at the whole-plant level. In this study, we investigated the relationships among SUTs in 24 representative plant species and performed an analysis of SUT genes in three sequenced Orchidaceae species: Dendrobium officinale, Phalaenopsis equestris, and Apostasia shenzhenica. All the SUTs from the 24 plant species were classified into three groups and five subgroups, subgroups A, B1, B2.1, B2.2, and C, based on their evolutionary relationships. A total of 22 SUT genes were identified among Orchidaceae species, among which D. officinale had 8 genes (DoSUT01-08), P. equestris had eight genes (PeqSUT01-08) and A. shenzhenica had 6 genes (AsSUT01-06). For the 22 OrchidaceaeSUTs, subgroups A, B2.2 and C contained three genes, whereas the SUT genes were found to have significantly expanded in the monocot-specific subgroup B2.1, which contained 12 genes. To understand sucrose partitioning and the functions of sucrose transporters in Orchidaceae species, we analyzed the water-soluble sugar content and performed RNA sequencing of different tissues of D. officinale, including leaves, stems, flowers and roots. The results showed that although the total content of water-soluble polysaccharides was highest in the stems of D. officinale, the sucrose content was highest in the flowers. Moreover, gene expression analysis showed that most of the DoSUTs were expressed in the flowers, among which DoSUT01,DoSUT07 and DoSUT06 had significantly increased expression levels. These results indicated that stems are used as the main storage sinks for photosynthetically produced sugar in D. officinale and that DoSUTs mainly function in the cellular machinery and development of floral organs. Our findings provide valuable information on sucrose partitioning and the evolution and functions of SUT genes in Orchidaceae and other species.

A high-quality chromosome-level genome assembly reveals genetics for important traits in eggplant.

Eggplant (Solanum melongena L.) is an economically important vegetable crop in the Solanaceae family, with extensive diversity among landraces and close relatives. Here, we report a high-quality reference genome for the eggplant inbred line HQ-1315 (S. melongena-HQ) using a combination of Illumina, Nanopore and 10X genomics sequencing technologies and Hi-C technology for genome assembly. The assembled genome has a total size of ~1.17 Gb and 12 chromosomes, with a contig N50 of 5.26 Mb, consisting of 36,582 protein-coding genes. Repetitive sequences comprise 70.09% (811.14 Mb) of the eggplant genome, most of which are long terminal repeat (LTR) retrotransposons (65.80%), followed by long interspersed nuclear elements (LINEs, 1.54%) and DNA transposons (0.85%). The S. melongena-HQ eggplant genome carries a total of 563 accession-specific gene families containing 1009 genes. In total, 73 expanded gene families (892 genes) and 34 contraction gene families (114 genes) were functionally annotated. Comparative analysis of different eggplant genomes identified three types of variations, including single-nucleotide polymorphisms (SNPs), insertions/deletions (indels) and structural variants (SVs). Asymmetric SV accumulation was found in potential regulatory regions of protein-coding genes among the different eggplant genomes. Furthermore, we performed QTL-seq for eggplant fruit length using the S. melongena-HQ reference genome and detected a QTL interval of 71.29-78.26 Mb on chromosome E03. The gene Smechr0301963, which belongs to the SUN gene family, is predicted to be a key candidate gene for eggplant fruit length regulation. Moreover, we anchored a total of 210 linkage markers associated with 71 traits to the eggplant chromosomes and finally obtained 26 QTL hotspots. The eggplant HQ-1315 genome assembly can be accessed at http://eggplant-hq.cn. In conclusion, the eggplant genome presented herein provides a global view of genomic divergence at the whole-genome level and powerful tools for the identification of candidate genes for important traits in eggplant.

Genome-Wide Identification and Functional Characterization of the Heat Shock Factor Family in Eggplant (Solanum melongena L.) under Abiotic Stress Conditions.

Plant heat shock factors (Hsfs) play crucial roles in various environmental stress responses. Eggplant (Solanum melongena L.) is an agronomically important and thermophilic vegetable grown worldwide. Although the functions of Hsfs under environmental stress conditions have been characterized in the model plant Arabidopsis thaliana and tomato, their roles in responding to various stresses remain unclear in eggplant. Therefore, we characterized the eggplant SmeHsf family and surveyed expression profiles mediated by the SmeHsfs under various stress conditions. Here, using reported Hsfs from other species as queries to search SmeHsfs in the eggplant genome and confirming the typical conserved domains, we identified 20 SmeHsf genes. The SmeHsfs were further classified into 14 subgroups on the basis of their structure. Additionally, quantitative real-time PCR revealed that SmeHsfs responded to four stresses-cold, heat, salinity and drought-which indicated that SmeHsfs play crucial roles in improving tolerance to various abiotic stresses. The expression pattern of SmeHsfA6b exhibited the most immediate response to the various environmental stresses, except drought. The genome-wide identification and abiotic stress-responsive expression pattern analysis provide clues for further analysis of the roles and regulatory mechanism of SmeHsfs under environmental stresses.

Construction of a SNP-Based Genetic Map Using SLAF-Seq and QTL Analysis of Morphological Traits in Eggplant.

Eggplant (Solanum melongena; 2n = 24) is an economically important fruit crop of the family Solanaceae that was domesticated in India and Southeast Asia. Construction of a high-resolution genetic map and map-based gene mining in eggplant have lagged behind other crops within the family such as tomato and potato. In this study, we conducted high-throughput single nucleotide polymorphism (SNP) discovery in the eggplant genome using specific length amplified fragment (SLAF) sequencing and constructed a high-density genetic map for the quantitative trait locus (QTL) analysis of multiple traits. An interspecific F2 population of 121 individuals was developed from the cross between cultivated eggplant "1836" and the wild relative S. linnaeanum "1809." Genomic DNA extracted from parental lines and the F2 population was subjected to high-throughput SLAF sequencing. A total of 111.74 Gb of data and 487.53 million pair-end reads were generated. A high-resolution genetic map containing 2,122 SNP markers and 12 linkage groups was developed for eggplant, which spanned 1530.75 cM, with an average distance of 0.72 cM between adjacent markers. A total of 19 QTLs were detected for stem height and fruit and leaf morphology traits of eggplant, explaining 4.08-55.23% of the phenotypic variance. These QTLs were distributed on nine linkage groups (LGs), but not on LG2, 4, and 9. The number of SNPs ranged from 2 to 11 within each QTL, and the genetic interval varied from 0.15 to 10.53 cM. Overall, the results establish a foundation for the fine mapping of complex QTLs, candidate gene identification, and marker-assisted selection of favorable alleles in eggplant breeding.

An irregularly striped rind mutant reveals new insight into the function of PG1ß in cucumber (Cucumis sativus L.).

KEY MESSAGE: Via bulked segregant analysis sequencing combined with linkage mapping, the ist gene responsible for the irregularly striped rind mutation was delimited to a 144-kb region in cucumber. Sequencing and expression analysis identified Csa1G005490 as the candidate gene. The rind appearance of cucumber is one of the most important commercial quality traits. Usually, an immature cucumber fruit has a uniform rind that varies from green to yellow to white among different cultivated varieties. In the present paper, we isolated a novel fruit appearance cucumber mutant, ist, that has an irregularly striped rind pattern. The mutant displayed green irregular stripes on a yellow-green background at the immature fruit stage. Genetic analysis revealed that a single recessive gene, ist, is responsible for this mutation. A BSA (bulked segregant analysis) sequencing approach combined with genetic mapping delimited the ist locus to an interval with a length of 144 kb, and 21 predicted genes were annotated in the region. Based on mutation site screening and expression analysis, two single-nucleotide polymorphisms within the candidate gene, Csa1G005490, were identified as constituting the mutation. Csa1G005490 encodes a polygalacturonase-1 noncatalytic subunit beta protein (PG1ß) known to be involved in fruit softening. The expression of Csa1G005490 was significantly lower in the ist mutant than in the wild type. Transcriptome analysis identified 1796 differentially expressed genes (DEGs) between the ist mutant and wild type. Gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that these DEGs were enriched mostly in photosynthesis and chlorophyll metabolism pathways. Decreased expression patterns of several chlorophyll synthesis genes in the mutant suggest that ist plays a key role in chlorophyll biosynthesis. These results will provide new insight into the molecular mechanism underlying rind appearance polymorphisms in cucumber.

Bioinformatics Analysis of the Lipoxygenase Gene Family in Radish (Raphanus sativus) and Functional Characterization in Response to Abiotic and Biotic Stresses.

Lipoxygenases (LOXs) are non-heme iron-containing dioxygenases involved in many developmental and stress-responsive processes in plants. However, little is known about the radish LOX gene family members and their functions in response to biotic and abiotic stresses. In this study, we completed a genome-wide analysis and expression profiling of RsLOX genes under abiotic and biotic stress conditions. We identified 11 RsLOX genes, which encoded conserved domains, and classified them in 9-LOX and 13-LOX categories according to their phylogenetic relationships. The characteristic structural features of 9-LOX and 13-LOX genes and the encoded protein domains as well as their evolution are presented herein. A qRT-PCR analysis of RsLOX expression levels in the roots under simulated drought, salinity, heat, and cold stresses, as well as in response to a Plasmodiophora brassicae infection, revealed three tandem-clustered RsLOX genes that are involved in responses to various environmental stresses via the jasmonic acid pathway. Our findings provide insights into the evolution and potential biological roles of RsLOXs related to the adaptation of radish to stress conditions.

Patterns of human social contact and contact with animals in Shanghai, China.

East Asia is as a principal hotspot for emerging zoonotic infections. Understanding the likely pathways for their emergence and spread requires knowledge on human-human and human-animal contacts, but such studies are rare. We used self-completed and interviewer-completed contact diaries to quantify patterns of these contacts for 965 individuals in 2017/2018 in a high-income densely-populated area of China, Shanghai City. Interviewer-completed diaries recorded more social contacts (19.3 vs. 18.0) and longer social contact duration (35.0 vs. 29.1 hours) than self-reporting. Strong age-assortativity was observed in all age groups especially among young participants (aged 7-20) and middle aged participants (25-55 years). 17.7% of participants reported touching animals (15.3% (pets), 0.0% (poultry) and 0.1% (livestock)). Human-human contact was very frequent but contact with animals (especially poultry) was rare although associated with frequent human-human contact. Hence, this densely populated area is more likely to act as an accelerator for human-human spread but less likely to be at the source of a zoonosis outbreak. We also propose that telephone interview at the end of reporting day is a potential improvement of the design of future contact surveys.

Transcriptome-Wide Identification and Characterization of Circular RNAs in Leaves of Chinese Cabbage (Brassica rapa L. ssp. pekinensis) in Response to Calcium Deficiency-Induced Tip-burn.

Circular RNA (circRNA) is a newly discovered non-coding RNA, which play significant roles in the function and transcriptional regulation of microRNA. To date, in Chinese cabbage, the functional characteristic of circRNAs in response to calcium deficiency-induced tip-burn have not been reported. In this study, 730 circRNAs were isolated from Chinese cabbage leaves, of which 23 and 22 were differentially expressed in different calcium deficiency stages compared with the control. Forty-six host genes of the differentially expressed circRNAs were identified, and one circRNA was found to act as miRNAs sponges. Based on the functional analysis of host genes and target mRNAs of the corresponding miRNAs, the identified circRNAs might participated in response to stimulus, electron carrier activity, ATPase activity, cell wall metabolism, transcription factors and plant hormone signal transduction. ABF2, a positive regulator of the abiotic stress response in the abscisic acid (ABA) pathway, may play a role in calcium deficiency tolerance through a circRNA regulatory pathway. Correspondingly, the concentration of ABA is also increased during the Ca2+ deficiency stress. Our results suggest that circRNAs participate in a broad range of biological processes and physiological functions in the response to calcium deficiency-induced tip-burn and provide a basis for further studies of the biological roles that circRNAs play in the plant stress response.

Comparative Transcriptome Analysis in Eggplant Reveals Selection Trends during Eggplant Domestication.

Eggplant (Solanum melongena L.) is an economically and nutritionally important fruit crop of the Solanaceae family, which was domesticated in India and southern China. However, the genome regions subjected to selective sweeps in eggplant remain unknown. In the present study, we performed comparative transcriptome analysis of cultivated and wild eggplant species with emphasis on the selection pattern during domestication. In total, 44,073 (S. sisymbriifolium) to 58,677 (S. melongena cultivar S58) unigenes were generated for the six eggplant accessions with total lengths of 36.6-46 Mb. The orthologous genes were assessed using the ratio of nonsynonymous (K a) to synonymous (K s) nucleotide substitutions to characterize selective patterns during eggplant domestication. We identified 19 genes under positive selection across the phylogeny that were classified into four groups. The gene (OG12205) under positive selection was possibly associated with fruit-related traits in eggplant, which may have resulted from human manipulation. Eight positive selected genes were potentially involved in stress tolerance or disease resistance, suggesting that environmental changes and biotic stresses were important selective pressures in eggplant domestication. Taken together, our results shed light on the effects of artificial and natural selection on the transcriptomes of eggplant and its wild relatives. Identification of the selected genes will facilitate the understanding of genetic architecture of domesticated-related traits and provide resources for resistant breeding in eggplant.

Investigation of evolutionary and expressional relationships in the function of the leucine-rich repeat receptor-like protein kinase gene family (LRR-RLK) in the radish (Raphanus sativus L.).

The leucine-rich repeat receptor-like protein kinase (LRR-RLK) plays an important role in plant development and disease defence. Although genome-wide studies of LRR-RLKs have been performed in several species, a comprehensive analysis, including evolutionary, structural and expressional analyses and their relationships to function, has not been carried out in the radish (Raphanus sativus L.). In this study, we identified 292 LRR-RLK genes in the R. sativus genome and classified them into 23 subgroups. The subgroups containing genes involved in defence were more likely to evolve from tandem duplication rather than whole genome triplication (WGT), had lower expression profiles and were expressed in fewer tissues than the subgroups related to development. Gene structures and conserved domains did not differ in the defence-related or development-related subgroups, but they were distinct in each subgroup. This study sheds light on the evolutionary and expressional relationships with the functions of R. sativus LRR-RLKs and provides an integrated framework for additional investigation into these functions.