V-primer: software for the efficient design of genome-wide InDel and SNP markers from multi-sample variant call format (VCF) genotyping data.

DNA markers are indispensable tools in genetics and genomics research as well as in crop breeding, particularly for marker-assisted selection. Recent advances in next-generation sequencing technology have made it easier to obtain genome sequences for various crop species, enabling the large-scale identification of DNA polymorphisms among varieties, which in turn has made DNA marker design more accessible. However, existing primer design software is not suitable for designing many types of genome-wide DNA markers from next-generation sequencing data. Here, we describe the development of V-primer, high-throughput software for designing insertion/deletion, cleaved amplified polymorphic sequence, and single-nucleotide polymorphism (SNP) markers. We validated the applicability of these markers in different crops. In addition, we performed multiplex PCR targeted amplicon sequencing using SNP markers designed with V-primer. Our results demonstrate that V-primer facilitates the efficient and accurate design of primers and is thus a useful tool for genetics, genomics, and crop breeding. V-primer is freely available at https://github.com/ncod3/vprimer.

Genome Analysis Revives a Forgotten Hybrid Crop Edo-dokoro in the Genus Dioscorea.

A rhizomatous Dioscorea crop 'Edo-dokoro' was described in old records of Japan, but its botanical identity has not been characterized. We found that Edo-dokoro is still produced by four farmers in Tohoku-machi of the Aomori prefecture, Japan. The rhizomes of Edo-dokoro are a delicacy to the local people and are sold in the markets. Morphological characters of Edo-dokoro suggest its hybrid origin between the two species, Dioscorea tokoro and Dioscorea tenuipes. Genome analysis revealed that Edo-dokoro likely originated by hybridization of a male D. tokoro to a female D. tenuipes, followed by a backcross with a male plant of D. tokoro. Edo-dokoro is a typical minor crop possibly maintained for more than 300 years but now almost forgotten by the public. We hypothesize that there are many such uncharacterized genetic heritages passed over generations by small-scale farmers that await serious scientific investigation for future use and improvement by using modern genomics information.

High-performance pipeline for MutMap and QTL-seq.

Summary: Bulked segregant analysis implemented in MutMap and QTL-seq is a powerful and efficient method to identify loci contributing to important phenotypic traits. However, the previous pipelines were not user-friendly to install and run. Here, we describe new pipelines for MutMap and QTL-seq. These updated pipelines are approximately 5-8 times faster than the previous pipeline, are easier for novice users to use, and can be easily installed through bioconda with all dependencies. Availability: The new pipelines of MutMap and QTL-seq are written in Python and can be installed via bioconda. The source code and manuals are available online (MutMap: https://github.com/YuSugihara/MutMap, QTL-seq: https://github.com/YuSugihara/QTL-seq).

Genome analyses reveal the hybrid origin of the staple crop white Guinea yam (Dioscorea rotundata).

White Guinea yam (Dioscorea rotundata) is an important staple tuber crop in West Africa. However, its origin remains unclear. In this study, we resequenced 336 accessions of white Guinea yam and compared them with the sequences of wild Dioscorea species using an improved reference genome sequence of D. rotundata In contrast to a previous study suggesting that D. rotundata originated from a subgroup of Dioscorea praehensilis, our results suggest a hybrid origin of white Guinea yam from crosses between the wild rainforest species D. praehensilis and the savannah-adapted species Dioscorea abyssinica We identified a greater genomic contribution from D. abyssinica in the sex chromosome of Guinea yam and extensive introgression around the SWEETIE gene. Our findings point to a complex domestication scenario for Guinea yam and highlight the importance of wild species as gene donors for improving this crop through molecular breeding.

Identification of candidate flowering and sex genes in white Guinea yam (D. rotundata Poir.) by SuperSAGE transcriptome profiling.

Dioecy (distinct male and female individuals) and scarce to non-flowering are common features of cultivated yam (Dioscorea spp.). However, the molecular mechanisms underlying flowering and sex determination in Dioscorea are largely unknown. We conducted SuperSAGE transcriptome profiling of male, female and monoecious individuals to identify flowering and sex-related genes in white Guinea yam (D. rotundata), generating 20,236 unique tags. Of these, 13,901 were represented by a minimum of 10 tags. A total 88 tags were significantly differentially expressed in male, female and monoecious plants, of which 18 corresponded to genes previously implicated in flower development and sex determination in multiple plant species. We validated the SuperSAGE data with quantitative real-time PCR (qRT-PCR)-based analysis of the expression of three candidate genes. We further investigated the flowering patterns of 1938 D. rotundata accessions representing diverse geographical origins over two consecutive years. Over 85% of accessions were either male or non-flowering, less than 15% were female, while monoecious plants were rare. Intensity of flowering varied between male and female plants, with the former flowering more abundantly than the latter. Candidate genes identified in this study can be targeted for further validation and to induce regular flowering in poor to non-flowering cultivars. Findings of the study provide important inputs for further studies aiming to overcome the challenge of flowering in yams and to improve efficiency of yam breeding.

Genome analysis of the foxtail millet pathogen Sclerospora graminicola reveals the complex effector repertoire of graminicolous downy mildews.

BACKGROUND: Downy mildew, caused by the oomycete pathogen Sclerospora graminicola, is an economically important disease of Gramineae crops including foxtail millet (Setaria italica). Plants infected with S. graminicola are generally stunted and often undergo a transformation of flower organs into leaves (phyllody or witches' broom), resulting in serious yield loss. To establish the molecular basis of downy mildew disease in foxtail millet, we carried out whole-genome sequencing and an RNA-seq analysis of S. graminicola. RESULTS: Sequence reads were generated from S. graminicola using an Illumina sequencing platform and assembled de novo into a draft genome sequence comprising approximately 360 Mbp. Of this sequence, 73% comprised repetitive elements, and a total of 16,736 genes were predicted from the RNA-seq data. The predicted genes included those encoding effector-like proteins with high sequence similarity to those previously identified in other oomycete pathogens. Genes encoding jacalin-like lectin-domain-containing secreted proteins were enriched in S. graminicola compared to other oomycetes. Of a total of 1220 genes encoding putative secreted proteins, 91 significantly changed their expression levels during the infection of plant tissues compared to the sporangia and zoospore stages of the S. graminicola lifecycle. CONCLUSIONS: We established the draft genome sequence of a downy mildew pathogen that infects Gramineae plants. Based on this sequence and our transcriptome analysis, we generated a catalog of in planta-induced candidate effector genes, providing a solid foundation from which to identify the effectors causing phyllody.

Genome sequencing of the staple food crop white Guinea yam enables the development of a molecular marker for sex determination.

BACKGROUND: Root and tuber crops are a major food source in tropical Africa. Among these crops are several species in the monocotyledonous genus Dioscorea collectively known as yam, a staple tuber crop that contributes enormously to the subsistence and socio-cultural lives of millions of people, principally in West and Central Africa. Yam cultivation is constrained by several factors, and yam can be considered a neglected "orphan" crop that would benefit from crop improvement efforts. However, the lack of genetic and genomic tools has impeded the improvement of this staple crop. RESULTS: To accelerate marker-assisted breeding of yam, we performed genome analysis of white Guinea yam (Dioscorea rotundata) and assembled a 594-Mb genome, 76.4% of which was distributed among 21 linkage groups. In total, we predicted 26,198 genes. Phylogenetic analyses with 2381 conserved genes revealed that Dioscorea is a unique lineage of monocotyledons distinct from the Poales (rice), Arecales (palm), and Zingiberales (banana). The entire Dioscorea genus is characterized by the occurrence of separate male and female plants (dioecy), a feature that has limited efficient yam breeding. To infer the genetics of sex determination, we performed whole-genome resequencing of bulked segregants (quantitative trait locus sequencing [QTL-seq]) in F1 progeny segregating for male and female plants and identified a genomic region associated with female heterogametic (male = ZZ, female = ZW) sex determination. We further delineated the W locus and used it to develop a molecular marker for sex identification of Guinea yam plants at the seedling stage. CONCLUSIONS: Guinea yam belongs to a unique and highly differentiated clade of monocotyledons. The genome analyses and sex-linked marker development performed in this study should greatly accelerate marker-assisted breeding of Guinea yam. In addition, our QTL-seq approach can be utilized in genetic studies of other outcrossing crops and organisms with highly heterozygous genomes. Genomic analysis of orphan crops such as yam promotes efforts to improve food security and the sustainability of tropical agriculture.

Lentinula edodes Genome Survey and Postharvest Transcriptome Analysis.

Lentinula edodes is a popular, cultivated edible and medicinal mushroom. Lentinula edodes is susceptible to postharvest problems, such as gill browning, fruiting body softening, and lentinan degradation. We constructed a de novo assembly draft genome sequence and performed gene prediction for Lentinula edodesDe novo assembly was carried out using short reads from paired-end and mate-paired libraries and by using long reads by PacBio, resulting in a contig number of 1,951 and an N50 of 1 Mb. Furthermore, we predicted genes by Augustus using transcriptome sequencing (RNA-seq) data from the whole life cycle of Lentinula edodes, resulting in 12,959 predicted genes. This analysis revealed that Lentinula edodes lacks lignin peroxidase. To reveal genes involved in the loss of quality of Lentinula edodes postharvest fruiting bodies, transcriptome analysis was carried out using serial analysis of gene expression (SuperSAGE). This analysis revealed that many cell wall-related enzymes are upregulated after harvest, such as ß-1,3-1,6-glucan-degrading enzymes in glycoside hydrolase (GH) families GH5, GH16, GH30, GH55, and GH128, and thaumatin-like proteins. In addition, we found that several chitin-related genes are upregulated, such as putative chitinases in GH family 18, exochitinases in GH20, and a putative chitosanase in GH family 75. The results suggest that cell wall-degrading enzymes synergistically cooperate for rapid fruiting body autolysis. Many putative transcription factor genes were upregulated postharvest, such as genes containing high-mobility-group (HMG) domains and zinc finger domains. Several cell death-related proteins were also upregulated postharvest.IMPORTANCE Our data collectively suggest that there is a rapid fruiting body autolysis system in Lentinula edodes The genes for the loss of postharvest quality newly found in this research will be targets for the future breeding of strains that keep fresh longer than present strains. De novoLentinula edodes genome assembly data will be used for the construction of a complete Lentinula edodes chromosome map for future breeding.

Screening of rice mutants with improved saccharification efficiency results in the identification of CONSTITUTIVE PHOTOMORPHOGENIC 1 and GOLD HULL AND INTERNODE 1.

MAIN CONCLUSION: The screening of rice mutants with improved cellulose to glucose saccharification efficiency (SE) identifies reduced xylan and/or ferulic acid, and a qualitative change of lignin to impact SE. To ensure the availability of sustainable energy, considerable effort is underway to utilize lignocellulosic plant biomass as feedstock for the production of biofuels. However, the high cost of degrading plant cell wall components to fermentable sugars (saccharification) has been problematic. One way to overcome this barrier is to develop plants possessing cell walls that are amenable to saccharification. In this study, we aimed to identify new molecular factors that influence saccharification efficiency (SE) in rice. By screening 22 rice mutants, we identified two lines, 122 and 108, with improved SE. Reduced xylan and ferulic acid within the cell wall of line 122 were probable reasons of improved SE. Line 108 showed reduced levels of thioglycolic-released lignin; however, the amount of Klason lignin was comparable to the wild-type, indicating that structural changes had occurred in the 108 lignin polymer which resulted in improved SE. Positional cloning revealed that the genes responsible for improved SE in 122 and 108 were rice CONSTITUTIVE PHOTOMORPHOGENIC 1 (OsCOP1) and GOLD HULL AND INTERNODE 1 (GH1), respectively, which have not been previously reported to influence SE. The screening of mutants for improved SE is an efficient approach to identify novel genes that affect SE, which is relevant in the development of crops as biofuel sources.

Draft genome sequence of bitter gourd (Momordica charantia), a vegetable and medicinal plant in tropical and subtropical regions.

Bitter gourd (Momordica charantia) is an important vegetable and medicinal plant in tropical and subtropical regions globally. In this study, the draft genome sequence of a monoecious bitter gourd inbred line, OHB3-1, was analyzed. Through Illumina sequencing and de novo assembly, scaffolds of 285.5 Mb in length were generated, corresponding to ∼84% of the estimated genome size of bitter gourd (339 Mb). In this draft genome sequence, 45,859 protein-coding gene loci were identified, and transposable elements accounted for 15.3% of the whole genome. According to synteny mapping and phylogenetic analysis of conserved genes, bitter gourd was more related to watermelon (Citrullus lanatus) than to cucumber (Cucumis sativus) or melon (C. melo). Using RAD-seq analysis, 1507 marker loci were genotyped in an F2 progeny of two bitter gourd lines, resulting in an improved linkage map, comprising 11 linkage groups. By anchoring RAD tag markers, 255 scaffolds were assigned to the linkage map. Comparative analysis of genome sequences and predicted genes determined that putative trypsin-inhibitor and ribosome-inactivating genes were distinctive in the bitter gourd genome. These genes could characterize the bitter gourd as a medicinal plant.

Host specialization of the blast fungus Magnaporthe oryzae is associated with dynamic gain and loss of genes linked to transposable elements.

BACKGROUND: Magnaporthe oryzae (anamorph Pyricularia oryzae) is the causal agent of blast disease of Poaceae crops and their wild relatives. To understand the genetic mechanisms that drive host specialization of M. oryzae, we carried out whole genome resequencing of four M. oryzae isolates from rice (Oryza sativa), one from foxtail millet (Setaria italica), three from wild foxtail millet S. viridis, and one isolate each from finger millet (Eleusine coracana), wheat (Triticum aestivum) and oat (Avena sativa), in addition to an isolate of a sister species M. grisea, that infects the wild grass Digitaria sanguinalis. RESULTS: Whole genome sequence comparison confirmed that M. oryzae Oryza and Setaria isolates form a monophyletic and close to another monophyletic group consisting of isolates from Triticum and Avena. This supports previous phylogenetic analysis based on a small number of genes and molecular markers. When comparing the host specific subgroups, 1.2-3.5 % of genes showed presence/absence polymorphisms and 0-6.5 % showed an excess of non-synonymous substitutions. Most of these genes encoded proteins whose functional domains are present in multiple copies in each genome. Therefore, the deleterious effects of these mutations could potentially be compensated by functional redundancy. Unlike the accumulation of nonsynonymous nucleotide substitutions, gene loss appeared to be independent of divergence time. Interestingly, the loss and gain of genes in pathogens from the Oryza and Setaria infecting lineages occurred more frequently when compared to those infecting Triticum and Avena even though the genetic distance between Oryza and Setaria lineages was smaller than that between Triticum and Avena lineages. In addition, genes showing gain/loss and nucleotide polymorphisms are linked to transposable elements highlighting the relationship between genome position and gene evolution in this pathogen species. CONCLUSION: Our comparative genomics analyses of host-specific M. oryzae isolates revealed gain and loss of genes as a major evolutionary mechanism driving specialization to Oryza and Setaria. Transposable elements appear to facilitate gene evolution possibly by enhancing chromosomal rearrangements and other forms of genetic variation.

Changes in the Proteome of Xylem Sap in Brassica oleracea in Response to Fusarium oxysporum Stress.

Fusarium oxysporum f.sp. conlutinans (Foc) is a serious root-invading and xylem-colonizing fungus that causes yellowing in Brassica oleracea. To comprehensively understand the interaction between F. oxysporum and B. oleracea, composition of the xylem sap proteome of the non-infected and Foc-infected plants was investigated in both resistant and susceptible cultivars using liquid chromatography-tandem mass spectrometry (LC-MS/MS) after in-solution digestion of xylem sap proteins. Whole genome sequencing of Foc was carried out and generated a predicted Foc protein database. The predicted Foc protein database was then combined with the public B. oleracea and B. rapa protein databases downloaded from Uniprot and used for protein identification. About 200 plant proteins were identified in the xylem sap of susceptible and resistant plants. Comparison between the non-infected and Foc-infected samples revealed that Foc infection causes changes to the protein composition in B. oleracea xylem sap where repressed proteins accounted for a greater proportion than those of induced in both the susceptible and resistant reactions. The analysis on the proteins with concentration change > = 2-fold indicated a large portion of up- and down-regulated proteins were those acting on carbohydrates. Proteins with leucine-rich repeats and legume lectin domains were mainly induced in both resistant and susceptible system, so was the case of thaumatins. Twenty-five Foc proteins were identified in the infected xylem sap and 10 of them were cysteine-containing secreted small proteins that are good candidates for virulence and/or avirulence effectors. The findings of differential response of protein contents in the xylem sap between the non-infected and Foc-infected samples as well as the Foc candidate effectors secreted in xylem provide valuable insights into B. oleracea-Foc interactions.

A chloroplast-localized protein LESION AND LAMINA BENDING affects defence and growth responses in rice.

Understanding how plants allocate their resources to growth or defence is of long-term importance to the development of new and improved varieties of different crops. Using molecular genetics, plant physiology, hormone analysis and Next-Generation Sequencing (NGS)-based transcript profiling, we have isolated and characterized the rice (Oryza sativa) LESION AND LAMINA BENDING (LLB) gene that encodes a chloroplast-targeted putative leucine carboxyl methyltransferase. Loss of LLB function results in reduced growth and yield, hypersensitive response (HR)-like lesions, accumulation of the antimicrobial compounds momilactones and phytocassanes, and constitutive expression of pathogenesis-related genes. Consistent with these defence-associated responses, llb shows enhanced resistance to rice blast (Magnaporthe oryzae) and bacterial blight (Xanthomonas oryzae pv. oryzae). The lesion and resistance phenotypes are likely to be caused by the over-accumulation of jasmonates (JAs) in the llb mutant including the JA precursor 12-oxo-phytodienoic acid. Additionally, llb shows an increased lamina inclination and enhanced early seedling growth due to elevated brassinosteroid (BR) synthesis and/or signalling. These findings show that LLB functions in the chloroplast to either directly or indirectly repress both JA- and BR-mediated responses, revealing a possible mechanism for controlling how plants allocate resources for defence and growth.

Grouping of multicopper oxidases in Lentinula edodes by sequence similarities and expression patterns.

The edible white rot fungus Lentinula edodes possesses a variety of lignin degrading enzymes such as manganese peroxidases and laccases. Laccases belong to the multicopper oxidases, which have a wide range of catalytic activities including polyphenol degradation and synthesis, lignin degradation, and melanin formation. The exact number of laccases in L. edodes is unknown, as are their complete properties and biological functions. We analyzed the draft genome sequence of L. edodes D703PP-9 and identified 13 multicopper oxidase-encoding genes; 11 laccases in sensu stricto, of which three are new, and two ferroxidases. lcc8, a laccase previously reported in L. edodes, was not identified in D703PP-9 genome. Phylogenetic analysis showed that the 13 multicopper oxidases can be classified into laccase sensu stricto subfamily 1, laccase sensu stricto subfamily 2 and ferroxidases. From sequence similarities and expression patterns, laccase sensu stricto subfamily 1 can be divided into two subgroups. Laccase sensu stricto subfamily 1 group A members are mainly secreted from mycelia, while laccase sensu stricto subfamily 1 group B members are expressed mainly in fruiting bodies during growth or after harvesting but are lowly expressed in mycelia. Laccase sensu stricto subfamily 2 members are mainly expressed in mycelia, and two ferroxidases are mainly expressed in the fruiting body during growth or after harvesting, and are expressed at very low levels in mycelium. Our data suggests that L. edodes laccases in same group share expression patterns and would have common biological functions.

A cytochrome P450, OsDSS1, is involved in growth and drought stress responses in rice (Oryza sativa L.).

Cytochrome P450s are among the largest protein coding gene families in plant genomes. However, majority of the genes remain uncharacterized. Here, we report the characterization of dss1, a rice mutant showing dwarfism and reduced grain size. The dss1 phenotype is caused by a non-synonymous point mutation we identified in DSS1, which is member of a P450 gene cluster located on rice chromosome 3 and corresponds to the previously reported CYP96B4/SD37 gene. Phenotypes of several dwarf mutants characterized in rice are associated with defects in the biosynthesis or perception of the phytohormones gibberellins (GAs) and brassinosteroids (BRs). However, both GA and BR failed to rescue the dss1 phenotype. Hormone profiling revealed the accumulation of abscisic acid (ABA) and ABA metabolites, as well as significant reductions in GA19 and GA53 levels, precursors of the bioactive GA1, in the mutant. The dss1 contents of cytokinin and auxins were not significantly different from wild-type plants. Consistent with the accumulation of ABA and metabolites, germination and early growth was delayed in dss1, which also exhibited an enhanced tolerance to drought. Additionally, expressions of members of the DSS1/CYP96B gene cluster were regulated by drought stress and exogenous ABA. RNA-seq-based transcriptome profiling revealed, among others, that cell wall-related genes and genes involved in lipid metabolism were up- and down-regulated in dss1, respectively. Taken together, these findings suggest that DSS1 mediates growth and stress responses in rice by fine-tuning GA-to-ABA balance, and might as well play a role in lipid metabolism.

The rice (Oryza sativa L.) LESION MIMIC RESEMBLING, which encodes an AAA-type ATPase, is implicated in defense response.

Lesion mimic mutants (LMMs) provide a useful tool to study defense-related programmed cell death (PCD) in plants. Although a number of LMMs have been identified in multiple species, most of the candidate genes are yet to be isolated. Here, we report the identification and characterization of a novel rice (Oryza sativa L.) lesion mimic resembling (lmr) mutant, and cloning of the corresponding LMR gene. The LMR locus was initially delineated to 1.2 Mb region on chromosome 6, which was further narrowed down to 155-kb using insertions/deletions (INDELs) and cleavage amplified polymorphic sequence markers developed in this study. We sequenced the open reading frames predicted within the candidate genomic region, and identified a G-A base substitution causing a premature translation termination in a gene that encodes an ATPase associated with various cellular activities type (AAA-type) protein. RNA interference transgenic lines with reduced LMR transcripts exhibited the lesion mimic phenotype similar to that of lmr plants. Furthermore, expression of the wild-type LMR in the mutant background complemented the lesion phenotype, confirming that the mutation identified in LMR is responsible for the mutant phenotype. The pathogenesis-related (PR) genes PBZ1 and PR1 were induced in lmr, which also showed enhanced resistance to rice blast (Magnaporthe oryzae) and bacterial blight (Xanthomonas oryzae pv. oryzae), suggesting LMR is a negative regulator of cell death in rice. The identification of lmr and cloning of the corresponding LMR gene provide an additional resource for the study of PCD in plants.

The Draft Genome of Hop (Humulus lupulus), an Essence for Brewing.

The female flower of hop (Humulus lupulus var. lupulus) is an essential ingredient that gives characteristic aroma, bitterness and durability/stability to beer. However, the molecular genetic basis for identifying DNA markers in hop for breeding and to study its domestication has been poorly established. Here, we provide draft genomes for two hop cultivars [cv. Saazer (SZ) and cv. Shinshu Wase (SW)] and a Japanese wild hop [H. lupulus var. cordifolius; also known as Karahanasou (KR)]. Sequencing and de novo assembly of genomic DNA from heterozygous SW plants generated scaffolds with a total size of 2.05 Gb, corresponding to approximately 80% of the estimated genome size of hop (2.57 Gb). The scaffolds contained 41,228 putative protein-encoding genes. The genome sequences for SZ and KR were constructed by aligning their short sequence reads to the SW reference genome and then replacing the nucleotides at single nucleotide polymorphism (SNP) sites. De novo RNA sequencing (RNA-Seq) analysis of SW revealed the developmental regulation of genes involved in specialized metabolic processes that impact taste and flavor in beer. Application of a novel bioinformatics tool, phylogenetic comparative RNA-Seq (PCP-Seq), which is based on read depth of genomic DNAs and RNAs, enabled the identification of genes related to the biosynthesis of aromas and flavors that are enriched in SW compared to KR. Our results not only suggest the significance of historical human selection process for enhancing aroma and bitterness biosyntheses in hop cultivars, but also serve as crucial information for breeding varieties with high quality and yield.

Genome sequence comparison reveals a candidate gene involved in male-hermaphrodite differentiation in papaya (Carica papaya) trees.

The sex type of papaya (Carica papaya) is determined by the pair of sex chromosomes (XX, female; XY, male; and XY(h), hermaphrodite), in which there is a non-recombining genomic region in the Y and Y(h) chromosomes. This region is presumed to be involved in determination of males and hermaphrodites; it is designated as the male-specific region in the Y chromosome (MSY) and the hermaphrodite-specific region in the Y(h) chromosome (HSY). Here, we identified the genes determining male and hermaphrodite sex types by comparing MSY and HSY genomic sequences. In the MSY and HSY genomic regions, we identified 14,528 nucleotide substitutions and 965 short indels with a large gap and two highly diverged regions. In the predicted genes expressed in flower buds, we found no nucleotide differences leading to amino acid changes between the MSY and HSY. However, we found an HSY-specific transposon insertion in a gene (SVP like) showing a similarity to the Short Vegetative Phase (SVP) gene. Study of SVP-like transcripts revealed that the MSY allele encoded an intact protein, while the HSY allele encoded a truncated protein. Our findings demonstrated that the SVP-like gene is a candidate gene for male-hermaphrodite determination in papaya.