Gibberellins regulate masculinization through the SpGAI-SpSTM module in dioecious spinach.

The sex of dioecious plants is mainly determined by genetic factors, but it can also be converted by environmental cues such as exogenous phytohormones. Gibberellic acids (GAs) are well-known inducers of flowering and sexual development, yet the pathway of gibberellin-induced sex conversion in dioecious spinach (Spinacia oleracea L.) remains elusive. Based on sex detection before and after GA3 application using T11A and SSR19 molecular markers, we confirmed and elevated the masculinization effect of GA on a single female plant through exogenous applications of GA3, showing complete conversion and functional stamens. Silencing of GIBBERELLIC ACID INSENSITIVE (SpGAI), a single DELLA family protein that is a central GA signaling repressor, results in similar masculinization. We also show that SpGAI can physically interact with the spinach KNOX transcription factor SHOOT MERISTEMLESS (SpSTM), which is a homolog of the flower meristem identity regulator STM in Arabidopsis. The silencing of SpSTM also masculinized female flowers in spinach. Furthermore, SpSTM could directly bind the intron of SpPI to repress SpPI expression in developing female flowers. Overall, our results suggest that GA induces a female masculinization process through the SpGAI-SpSTM-SpPI regulatory module in spinach. These insights may help to clarify the molecular mechanism underlying the sex conversion system in dioecious plants while also elucidating the physiological basis for the generation of unisexual flowers so as to establish dioecy in plants.

Pathogenic microorganism DNA high-throughput genetic sequencing to diagnose peritoneal dialysis-associated peritonitis due to Mycobacterium tuberculosis infection.

Peritoneal dialysis-associated peritonitis is a serious complication of peritoneal dialysis, and the prevention and treatment of this condition are important for improving the long-term survival and quality of life of patients. However, peritoneal dialysis-associated peritonitis due to Mycobacterium tuberculosis infection is relatively rare and not easily diagnosed. Here, we present a case of peritoneal dialysis-associated peritonitis caused by Mycobacterium tuberculosis identified by pathogenic microbial DNA high-throughput genetic sequencing. This case demonstrates that pathogenic microbial DNA high-throughput genetic sequencing could be used to improve the detection rate of pathogenic microorganisms in patients with complex conditions, thereby allowing for earlier initiation of treatment.

Lineage-specific amplification and epigenetic regulation of LTR-retrotransposons contribute to the structure, evolution, and function of Fabaceae species.

BACKGROUND: Long terminal repeat (LTR)-retrotransposons (LTR-RTs) are ubiquitous and make up the majority of nearly all sequenced plant genomes, whereas their pivotal roles in genome evolution, gene expression regulation as well as their epigenetic regulation are still not well understood, especially in a large number of closely related species. RESULTS: Here, we analyzed the abundance and dynamic evolution of LTR-RTs in 54 species from an economically and agronomically important family, Fabaceae, and also selected two representative species for further analysis in expression of associated genes, transcriptional activity and DNA methylation patterns of LTR-RTs. Annotation results revealed highly varied proportions of LTR-RTs in these genomes (5.1%~68.4%) and their correlation with genome size was highly positive, and they were significantly contributed to the variance in genome size through species-specific unique amplifications. Almost all of the intact LTR-RTs were inserted into the genomes 4 Mya (million years ago), and more than 50% of them were inserted in the last 0.5 million years, suggesting that recent amplifications of LTR-RTs were an important force driving genome evolution. In addition, expression levels of genes with intronic, promoter, and downstream LTR-RT insertions of Glycine max and Vigna radiata, two agronomically important crops in Fabaceae, showed that the LTR-RTs located in promoter or downstream regions suppressed associated gene expression. However, the LTR-RTs within introns promoted gene expression or had no contribution to gene expression. Additionally, shorter and younger LTR-RTs maintained higher mobility and transpositional potential. Compared with the transcriptionally silent LTR-RTs, the active elements showed significantly lower DNA methylation levels in all three contexts. The distributions of transcriptionally active and silent LTR-RT methylation varied across different lineages due to the position of LTR-RTs located or potentially epigenetic regulation. CONCLUSION: Lineage-specific amplification patterns were observed and higher methylation level may repress the activity of LTR-RTs, further influence evolution in Fabaceae species. This study offers valuable clues into the evolution, function, transcriptional activity and epigenetic regulation of LTR-RTs in Fabaceae genomes.

Endosomal Sorting Complex Required for Transport: The Molecular Machine for Ensuring Proper Mitosis.

Mitosis is a key step of eukaryote proliferation.Endosomal sorting complex required for transport, a protein complex closely associated with membrane shearing, is involved in endosome maturation, virus budding, and autophagy.The structural and functional abnormalities of the complex are associated with the occurrence and progression of cancer and other diseases.In this paper, we summarized the roles of the endosomal sorting complex required for transport in different stages of mitosis and reviewed the studies about the role of the complex in regulating mitosis in diseases.

Impact of LTR-Retrotransposons on Genome Structure, Evolution, and Function in Curcurbitaceae Species.

Long terminal repeat (LTR)-retrotransposons (LTR-RTs) comprise a major portion of many plant genomes and may exert a profound impact on genome structure, function, and evolution. Although many studies have focused on these elements in an individual species, their dynamics on a family level remains elusive. Here, we investigated the abundance, evolutionary dynamics, and impact on associated genes of LTR-RTs in 16 species in an economically important plant family, Cucurbitaceae. Results showed that full-length LTR-RT numbers and LTR-RT content varied greatly among different species, and they were highly correlated with genome size. Most of the full-length LTR-RTs were amplified after the speciation event, reflecting the ongoing rapid evolution of these genomes. LTR-RTs highly contributed to genome size variation via species-specific distinct proliferations. The Angela and Tekay lineages with a greater evolutionary age were amplified in Trichosanthes anguina, whereas a recent activity burst of Reina and another ancient round of Tekay activity burst were examined in Sechium edule. In addition, Tekay and Retand lineages belonging to the Gypsy superfamily underwent a recent burst in Gynostemma pentaphyllum. Detailed investigation of genes with intronic and promoter LTR-RT insertion showed diverse functions, but the term of metabolism was enriched in most species. Further gene expression analysis in G.pentaphyllum revealed that the LTR-RTs within introns suppress the corresponding gene expression, whereas the LTR-RTs within promoters exert a complex influence on the downstream gene expression, with the main function of promoting gene expression. This study provides novel insights into the organization, evolution, and function of LTR-RTs in Cucurbitaceae genomes.

Development and Validation of a Random Forest Risk Prediction Pneumothorax Model in Percutaneous Transthoracic Needle Biopsy.

BACKGROUND Computed tomography (CT)-guided percutaneous transthoracic needle biopsy (PTNB) is an effective means for diagnosing various thoracic diseases. Pneumothorax is the most common complication, and when it becomes life-threatening, urgent medical intervention is required. The purpose of this study was to develop and validate a model that can be used to predict postoperative pneumothorax following CT-guided PTNB. MATERIAL AND METHODS We enrolled 245 patients who completed CT-guided PTNB to develop the model. A random forest (RF) model was built using 15 risk factors (15-RFs). The 7 most critical risk factors (7-RFs) were extracted by feature selection and used to build a new model. The independent external validation data contained 97 patients. Logistic regression (LR), support vector machine (SVM), and decision tree (DT) models were also developed using both 15-RFs and 7-RFs, and their performance was compared with the RF models. RESULTS The length of the aerated lung traversed was identified as the most important risk factor for developing pneumothorax, followed by angle of pleural puncture, lesion depth, lesion size, age, procedure time, and sex. The RF model demonstrated better performance in the development and validation datasets when compared with the LR, SVM, and DT based on 15-RFs and 7-RFs. According to DeLong's test for difference in ROC curves, the RF models based on the 15-RFs and 7-RFs achieved similar classification performance (P>0.05). CONCLUSIONS This study demonstrated the feasibility of using the 7-RFs RF model for predicting postoperative pneumothorax before patients undergo CT-guided PTNB.

Centromere protein U enhances the progression of bladder cancer by promoting mitochondrial ribosomal protein s28 expression.

Bladder cancer is one of the most common types of cancer. Most gene mutations related to bladder cancer are dominantly acquired gene mutations and are not inherited. Previous comparative transcriptome analysis of urinary bladder cancer and control samples has revealed a set of genes that may play a role in tumor progression. Here we set out to investigate further the expression of two candidate genes, centromere protein U (CENPU) and mitochondrial ribosomal protein s28 (MRPS28) to better understand their role in bladder cancer pathogenesis. Our results confirmed that CENPU is up-regulated in human bladder cancer tissues at mRNA and protein levels. Gain-of-function and loss-of-function studies in T24 human urinary bladder cancer cell line revealed a hierarchical relationship between CENPU and MRPS28 in the regulation of cell viability, migration and invasion activity. CENPU expression was also up-regulated in in vivo nude mice xenograft model of bladder cancer and mice overexpressing CENPU had significantly higher tumor volume. In summary, our findings identify CENPU and MRPS28 in the molecular pathogenesis of bladder cancer and suggest that CENPU enhances the progression of bladder cancer by promoting MRPS28 expression.

Nuclear Integrants of Organellar DNA Contribute to Genome Structure and Evolution in Plants.

The transfer of genetic material from the mitochondria and plastid to the nucleus gives rise to nuclear integrants of mitochondrial DNA (NUMTs) and nuclear integrants of plastid DNA (NUPTs). This frequently occurring DNA transfer is ongoing and has important evolutionary implications. In this review, based on previous studies and the analysis of NUMT/NUPT insertions of more than 200 sequenced plant genomes, we analyzed and summarized the general features of NUMTs/NUPTs and highlighted the genetic consequence of organellar DNA insertions. The statistics of organellar DNA integrants among various plant genomes revealed that organellar DNA-derived sequence content is positively correlated with the nuclear genome size. After integration, the nuclear organellar DNA could undergo different fates, including elimination, mutation, rearrangement, fragmentation, and proliferation. The integrated organellar DNAs play important roles in increasing genetic diversity, promoting gene and genome evolution, and are involved in sex chromosome evolution in dioecious plants. The integrating mechanisms, involving non-homologous end joining at double-strand breaks were also discussed.

[Risk factors for cow's milk protein allergy in infants: a multicenter survey].

OBJECTIVE: To investigate the risk factors for cow's milk protein allergy (CMPA) among infants through a multicenter clinical study. METHODS: A total of 1 829 infants, aged 1-12 months, who attended the outpatient service of the pediatric department in six hospitals in Shenzhen, China from June 2016 to May 2017 were enrolled as subjects. A questionnaire survey was performed to screen out suspected cases of CMPA. Food avoidance and oral food challenge tests were used to make a confirmed diagnosis of CMPA CMPA. A multivariate logistic regression analysis was used to investigate the risk factors for CMPA. RESULTS: Among the 1 829 infants, 82 (4.48%) were diagnosed with CMPA. The multivariate logistic regression analysis showed that maternal food allergy (OR=4.91, 95%CI: 2.24-10.76, P<0.05), antibiotic exposure during pregnancy (OR=3.18, 95%CI: 1.32-7.65, P<0.05), and the introduction of complementary food at an age of <4 months (OR=3.55, 95%CI: 1.52-8.27, P<0.05) were risk factors for CMPA, while exclusive breastfeeding (OR=0.21, 95%CI: 0.08-0.58, P<0.05) and the introduction of complementary food at an age of >6 months (OR=0.38, 95%CI: 0.17-0.86, P<0.05) were protective factors. CONCLUSIONS: The introduction of complementary food at an age of <4 months, maternal food allergy, and antibiotic exposure during pregnancy are risk factors for CMPA in infants.

Cytogenetic and genomic organization analyses of chloroplast DNA invasions in the nuclear genome of Asparagus officinalis L. provides signatures of evolutionary complexity and informativity in sex chromosome evolution.

BACKGROUND: The transfer of chloroplast DNA into nuclear genome is a common process in plants. These transfers form nuclear integrants of plastid DNAs (NUPTs), which are thought to be driving forces in genome evolution, including sex chromosome evolution. In this study, NUPTs in the genome of a dioecious plant Asparagus officinalis L. were systematically analyzed, in order to investigate the characteristics of NUPTs in the nuclear genome and the relationship between NUPTs and sex chromosome evolution in this species. RESULTS: A total of 3155 NUPT insertions were detected, and they represented approximated 0.06% of the nuclear genome. About 45% of the NUPTs were organized in clusters. These clusters were derived from various evolutionary events. The Y chromosome contained the highest number and largest proportion of NUPTs, suggesting more accumulation of NUPTs on sex chromosomes. NUPTs were distributed widely in all of the chromosomes, and some regions preferred these insertions. The highest density of NUPTs was found in a 47 kb region in the Y chromosome; more than 75% of this region was occupied by NUPTs. Further cytogenetic and sequence alignment analysis revealed that this region was likely the centromeric region of the sex chromosomes. On the other hand, the male-specific region of the Y chromosome (MSY) and the adjacent regions did not have NUPT insertions. CONCLUSIONS: These results indicated that NUPTs were involved in shaping the genome of A. officinalis through complicated process. NUPTs may play important roles in the centromere shaping of the sex chromosomes of A. officinalis, but were not implicated in MSY formation.

Genome-wide characterization of microsatellites and genetic diversity assessment of spinach in the Chinese germplasm collection.

Spinach is a nutritional leafy green vegetable, and it also serves as a model species for studying sex chromosome evolution. Genetic marker development and genome structure analysis are important in breeding practice and theoretical evolution studies of spinach. In this study, the frequency and distribution of different microsatellites in the recently released draft spinach genome were characterized. A total of 261,002 perfect microsatellites were identified (estimated frequency: ~262.1 loci/Mbp). The most abundant microsatellites were tetranucleotide and trinucleotide, accounting for 33.2% and 27.7% of the total number of microsatellites, respectively. A total of 105 primer pairs were designed and screened, and 34 were polymorphic among the detected spinach cultivars. Combined with seven primer sets developed previously, 41 primer pairs were used to investigate genetic diversity among 43 spinach cultivars in China. The average polymorphism information content value of the 41 markers was 0.43, representing an intermediate level. The spinach cultivars had a low genetic diversity, and no detectable common factors were shared by each group in the UPGMA dendrogram. This study's findings facilitate further investigations on the organization of the microsatellites in spinach genome and provide clues for future breeding applications of spinach in China.

Comparative transcriptome analysis reveals differentially expressed genes associated with sex expression in garden asparagus (Asparagus officinalis).

BACKGROUND: Garden asparagus (Asparagus officinalis) is a highly valuable vegetable crop of commercial and nutritional interest. It is also commonly used to investigate the mechanisms of sex determination and differentiation in plants. However, the sex expression mechanisms in asparagus remain poorly understood. RESULTS: De novo transcriptome sequencing via Illumina paired-end sequencing revealed more than 26 billion bases of high-quality sequence data from male and female asparagus flower buds. A total of 72,626 unigenes with an average length of 979 bp were assembled. In comparative transcriptome analysis, 4876 differentially expressed genes (DEGs) were identified in the possible sex-determining stage of female and male/supermale flower buds. Of these DEGs, 433, including 285 male/supermale-biased and 149 female-biased genes, were annotated as flower related. Of the male/supermale-biased flower-related genes, 102 were probably involved in anther development. In addition, 43 DEGs implicated in hormone response and biosynthesis putatively associated with sex expression and reproduction were discovered. Moreover, 128 transcription factor (TF)-related genes belonging to various families were found to be differentially expressed, and this finding implied the essential roles of TF in sex determination or differentiation in asparagus. Correlation analysis indicated that miRNA-DEG pairs were also implicated in asparagus sexual development. CONCLUSIONS: Our study identified a large number of DEGs involved in the sex expression and reproduction of asparagus, including known genes participating in plant reproduction, plant hormone signaling, TF encoding, and genes with unclear functions. We also found that miRNAs might be involved in the sex differentiation process. Our study could provide a valuable basis for further investigations on the regulatory networks of sex determination and differentiation in asparagus and facilitate further genetic and genomic studies on this dioecious species.

Identification of male-specific AFLP and SCAR markers in the dioecious plant Humulus scandens.

In this study, 17 male-specific amplified fragment length polymorphism (AFLP) markers were identified between male and female Humulus scandens plants. BLAST analysis revealed that 7 of the 17 sex-linked sequences were highly similar to retrotransposons. Two stable male-specific sequence-characterized amplified regions (SCAR) markers were developed. These AFLP and SCAR markers are novel molecular probes that can be used efficiently to identify the genetic gender of H. scandens and may provide a basis for further investigations on the evolution of sex chromosomes.

The application of MutMap in forward genetic studies based on whole-genome sequencing.

Classical forward genetic analysis relies on construction of complicated progeny populations and development of many molecular markers for linkage analysis in genetic mapping, which is both time- and cost-consuming. The recently developed MutMap is a new forward genetic approach based on high-throughput next-generation sequencing technologies. It is more efficient and affordable than traditional methods. Moreover, new extended methods based on MutMap have been developed: MutMap+, which is based on self-crossing; MutMap-Gap, which is used to recognize the causative variations occurring in genome gap regions; QTL-seq, a method similar to MutMap for mapping quantitative trait loci. These methods are free from constructing complicated mapping population, genetic hybridization and linkage information. They have greatly accelerated the identification of genetic elements associated with interested phenotypic variation. Here, we review the basic principles of MutMap, and discuss their future applications in next generation sequencing-based forward genetic mapping and crop improvement.

[Analysis on medication rules of modern traditional Chinese medicines in treating palpitations based on traditional Chinese medicine inheritance support system].

To analyze the prescription and medication rules of Chinese medicines in the treatment of palpitations in the Chinese journal full text database(CNKI) by using traditional Chinese medicine inheritance system, and provide a reference for further research and development of modern traditional Chinese medicines(TCMs) in treatment of palpitations. In order to give better guidance for clinical mediation, prescriptions used for treatment of palpitations in CNKI were collected, and then were input to the TCM inheritance support system for establishing a Chinese medicine prescription database for palpitations. The software's revised mutual information, complex system entropy clustering and other data mining methods were adopted to analyze the prescriptions according to the frequencies of herbs, "four natures", "five flavors" and "meridians" of the high-frequency medicines in the database, identify the core herbs and application characteristics, and analyze the prescription rules and medication experience. Totally, 545 prescriptions used for palpitation were included in this study and involved 247 Chinese herbs. The analysis results showed that the herbs in prescriptions for palpitation mostly had the warm property, and the herbs in heart and spleen meridian accounted for a larger proportion, indicating that the treatment was mainly to nourish heart and strengthen spleen. The top 11 herbs in usage frequency were consistent with the high-frequency medicines in medication patterns of common herbal pairs; therefore, we considered that these 11 herbs were the core herbs; the core herbal combination included Cassia Twig, Licorice, fossil fragments, Ostreae decoction, and evolved into 9 new prescriptions for treating palpitation. Our results objectively presented the prescription and medication rules for treating palpitation and provided extremely effective guidance for the clinical therapy.

A Systematic Review and Meta-Analysis of C-Arm Cone-Beam CT-Guided Percutaneous Transthoracic Needle Biopsy of Lung Nodules.

BACKGROUND: A systematic review and meta-analysis of all available publications was performed to evaluate the diagnostic accuracy of percutaneous transthoracic needle biopsy (PTNB) using a C-Arm Cone-Beam CT (CBCT) system in patients with lung nodules. MATERIAL/METHODS: Thedatabases of PUBMED, OVID, EBSCO, EMBASE, and China National Knowledge Infrastructure (CNKI) were systematically searched for relevant original articles on the diagnostic accuracy of CBCT-guided PTNB for the diagnosis of nodules in the lungs. Diagnostic indices including sensitivity, specificity, positive likelihood ratio (PLR), negative likelihood ratio (NLR), diagnostic odds ratio (DOR) and diagnostic score (DS) were calculated. Moreover,summary receiver operating characteristic curves (SROC) were constructed with Stata (version 13.0), Rev Man (version 5.3), and Meta-disc (version 1.4) software. Other clinical indices such as incidence of complications were also recorded. RESULTS: Eight studies met the inclusion and exclusion criteria for the meta-analysis. The pooled sensitivity, specificity, PLR, NLR, DOR, DS, and SROC with 95% confidence intervals were 0.96 (0.93-0.98), 1.00 (0.91-1.00), 711.15 (9.48-53325.89), 0.04 (0.02-0.07), 16585.29 (284.88-9.7e+05), 9.72 (5.65-13.78), and 0.99 (0.97-0.99), respectively. The incidence of pneumothorax and hemorrhage was 10-29.27% and 1.22-47.25%, respectively. CONCLUSIONS: CBCT-guided PTNB has an acceptable rate of complications and is associated with a reasonable radiation exposure. Moreover, it is a highly accurate and safe technique for the diagnosis of lung nodules and can be recommended to be used in routine clinical practice.

Repetitive sequences and epigenetic modification: inseparable partners play important roles in the evolution of plant sex chromosomes.

MAIN CONCLUSION: The present review discusses the roles of repetitive sequences played in plant sex chromosome evolution, and highlights epigenetic modification as potential mechanism of repetitive sequences involved in sex chromosome evolution. Sex determination in plants is mostly based on sex chromosomes. Classic theory proposes that sex chromosomes evolve from a specific pair of autosomes with emergence of a sex-determining gene(s). Subsequently, the newly formed sex chromosomes stop recombination in a small region around the sex-determining locus, and over time, the non-recombining region expands to almost all parts of the sex chromosomes. Accumulation of repetitive sequences, mostly transposable elements and tandem repeats, is a conspicuous feature of the non-recombining region of the Y chromosome, even in primitive one. Repetitive sequences may play multiple roles in sex chromosome evolution, such as triggering heterochromatization and causing recombination suppression, leading to structural and morphological differentiation of sex chromosomes, and promoting Y chromosome degeneration and X chromosome dosage compensation. In this article, we review the current status of this field, and based on preliminary evidence, we posit that repetitive sequences are involved in sex chromosome evolution probably via epigenetic modification, such as DNA and histone methylation, with small interfering RNAs as the mediator.

Genes encoding Δ(8)-sphingolipid desaturase from various plants: identification, biochemical functions, and evolution.

∆(8)-sphingolipid desaturase catalyzes the C8 desaturation of a long chain base, which is the characteristic structure of various complex sphingolipids. The genes of 20 ∆(8)-sphingolipid desaturases from 12 plants were identified and functionally detected by using Saccharomyces cerevisiae system to elucidate the relationship between the biochemical function and evolution of this enzyme. Results showed that the 20 genes all can encode a functional ∆(8)-sphingolipid desaturase, which catalyzes different ratios of two products, namely, 8(Z) and 8(E)-C18-phytosphingenine. The coded enzymes could be divided into two groups on the basis of biochemical functions: ∆(8)-sphingolipid desaturase with a preference for an E-isomer product and ∆(8)-sphingolipid desaturase with a preference for a Z-isomer product. The conversion rate of the latter was generally lower than that of the former. Phylogenetic analysis revealed that the 20 desaturases could also be clustered into two groups, and this grouping is consistent with that of the biochemical functions. Thus, the biochemical function of ∆(8)-sphingolipid desaturase is correlated with its evolution. The two groups of ∆(8)-sphingolipid desaturases could arise from distinct ancestors in higher plants. However, they might have initially evolved from ∆(8)-sphingolipid desaturases in lower organisms, such as yeasts, which can produce E-isomer products only. Furthermore, almost all of the transgenic yeasts harboring ∆(8)-sphingolipid desaturase genes exhibit an improvement in aluminum tolerance. Our study provided new insights into the biochemical function and evolution of ∆(8)-sphingolipid desaturases in plants.

[The Design of an ATR Probe for Online Monitoring of Biological Process].

The in situ biological process relays on a high performance attenuated total reflection (ATR) probe while ATR technique is a powerful tool for inline analysis. The basic principle and characteristic of the ATR technique was analyzed in this paper. A low cost and high throughout ATR probe was designed which can be inserted into the biological reactor directly with nondestructive and zero delay monitoring. It shows that the online ATR probe has a wide spectrum range, and the collected spectra has a high signal-to-noise ratio, which can be used in in-line quantitative analysis for biological process.

Isolation of differentially expressed sex genes in garden asparagus using suppression subtractive hybridization.

Garden asparagus (Asparagus officinalis L.) is a dioecious species whose male and female flowers are found in separate unisexual individuals. A region called the M-locus, located on a pair of homomorphic sex chromosomes, controls sexual dimorphism in asparagus. To date, no sex determining gene has been isolated from asparagus. To identify more genes involved in flower development in asparagus, subtractive hybridization library of male flowers in asparagus was constructed by suppression subtraction hybridization. A total of 107 expressed sequence tags (ESTs) were identified. BLASTX analysis showed that the library contained several genes that could be related to flower development. The expression patterns of seven selected genes believed to be involved in the development of asparagus male flower were further analyzed by semi-quantitative or real-time reverse-transcription polymerase chain reaction (RT-PCR). Results showed that AOEST4-5, AOEST12-40, and AOEST13-38 were strongly expressed in the male flower stage, whereas no transcript level of AOEST13-38 was detected in the female flower stage. The expression levels of AOEST13-87, AOEST13-92, AOEST13-40, and AOEST18-87 in the male flower stage were also higher than those in the female flower stage, although these transcripts were also expressed in other tissues. The identified genes can provide a strong starting point for further studies on the underlying molecular differences between the male and female flowers of asparagus.