Association analysis of production traits of Japanese quail (Coturnix japonica) using restriction-site associated DNA sequencing.

This study was designed to perform an association analysis and identify SNP markers associated with production traits of Japanese quail using restriction-site-associated DNA sequencing. Weekly body weight data from 805 quail were collected from hatching to 16 weeks of age. A total number of 3990 eggs obtained from 399 female quail were used to assess egg quality traits. Egg-related traits were measured at the beginning of egg production (first stage) and at 12 weeks of age (second stage). Five eggs were analyzed at each stage. Traits, such as egg weight, egg length and short axes, eggshell strength and weight, egg equator thickness, yolk weight, diameter, and colour, albumen weight, age of first egg, total number of laid eggs, and egg production rate, were assessed. A total of 383 SNPs and 1151 associations as well as 734 SNPs and 1442 associations were identified in relation to quail production traits using general linear model (GLM) and mixed linear model (MLM) approaches, respectively. The GLM-identified SNPs were located on chromosomes 1-13, 15, 17-20, 24, 26-28, and Z, underlying phenotypic traits, except for egg and albumen weight at the first stage and yolk yellowness at the second stage. The MLM-identified SNPs were positioned on defined chromosomes associated with phenotypic traits except for the egg long axis at the second stage of egg production. Finally, 35 speculated genes were identified as candidate genes for the targeted traits based on their nearest positions. Our findings provide a deeper understanding and allow a more precise genetic improvement of production traits of Galliformes, particularly in Japanese quail.

A genome-wide association and fine-mapping study of white rust resistance in hexaploid chrysanthemum cultivars with a wild diploid reference genome.

White rust caused by Puccinia horiana is one of the most serious diseases of chrysanthemum (Chrysanthemum × morifolium). In this study, we report the DNA markers associated with resistance against P. horiana via a simple approach using the genome of a wild diploid relative, Chrysanthemum seticuspe. First, we identified the important region of the genome in the resistant cultivar "Ariesu" via a genome-wide association study. Simplex single nucleotide polymorphism (SNP) markers mined from ddRAD-Seq were used in a biparental population originating from crosses between resistant "Ariesu" and susceptible "Yellow Queen". The C. seticuspe genome was used as a reference. For the fine mapping of P. horiana resistance locus 2 (Phr2), a comparative whole genome sequencing study was conducted. Although the genome sequences of chrysanthemum cultivars assembled via the short-read approach were fragmented, reliable genome alignments were reconstructed by mapping onto the chromosome level of the C. seticuspe pseudomolecule. Base variants were then identified by comparing the assembled genome sequences of resistant "Ariesu" and susceptible "Yellow Queen". Consequently, SNP markers that were closer to Phr2 compared with ddRAD-Seq markers were obtained. These SNP markers co-segregated with resistance in F1 progenies originating from resistant "Ariesu" and showed robust transferability for detecting Phr2-conferring resistance among chrysanthemum genetic resources. The wild C. seticuspe pseudomolecule, a de facto monoploid genome used for ddRAD-Seq analysis and assembled genome sequence comparison, demonstrated this method's utility as a model for developing DNA markers in hexaploid chrysanthemum cultivars.

The complete sequence of the chloroplast genome of Chrysanthemum rupestre, a diploid disciform capitula species of Chrysanthemum.

In this study, we analyzed the complete sequence of the chloroplast genome of Chrysanthemum rupestre Matsum. et Koidz., 1910, a diploid disciform capitula species of Chrysanthemum endemic to Japan, formerly classified as Ajania rupestris (Matsum. & Koidz.) Muldashev, Bot. Zhurn. (Moscow & Leningrad), 1983. The chloroplast genome of C. rupestre has a typical conserved quadripartite structure of 151,061 bp in length, comprising a large single-copy region (82,846 bp), a small single-copy region (18,301 bp), and a pair of inverted repeat regions (each 24,957 bp). Phylogenetic analysis indicated that C. rupestre clustered with other Chrysanthemum species, including another former Ajania species, Chrysanthemum pacificum Nakai, 1928. However, Ajania variifolia (C.C.Chang) Tzvelev, 1961, which is a synonym of Phaeostigma variifolium (C.C.Chang) Muldashev, 1981, was placed outside the Chrysanthemum clade, thereby implying that the former genus Ajania includes heterogeneous species.

Intragenic complementation at the Lotus japonicus CELLULOSE SYNTHASE-LIKE D1 locus rescues root hair defects.

Root hair cells form the primary interface of plants with the soil environment, playing key roles in nutrient uptake and plant defense. In legumes, they are typically the first cells to become infected by nitrogen-fixing soil bacteria during root nodule symbiosis. Here, we report a role for the CELLULOSE SYNTHASE-LIKE D1 (CSLD1) gene in root hair development in the legume species Lotus japonicus. CSLD1 belongs to the cellulose synthase protein family that includes cellulose synthases and cellulose synthase-like proteins, the latter thought to be involved in the biosynthesis of hemicellulose. We describe 11 Ljcsld1 mutant alleles that impose either short (Ljcsld1-1) or variable (Ljcsld1-2 to 11) root hair length phenotypes. Examination of Ljcsld1-1 and one variable-length root hair mutant, Ljcsld1-6, revealed increased root hair cell wall thickness, which in Ljcsld1-1 was significantly more pronounced and also associated with a strong defect in root nodule symbiosis. Lotus japonicus plants heterozygous for Ljcsld1-1 exhibited intermediate root hair lengths, suggesting incomplete dominance. Intragenic complementation was observed between alleles with mutations in different CSLD1 domains, suggesting CSLD1 function is modular and that the protein may operate as a homodimer or multimer during root hair development.

A chromosome-level genome sequence of Chrysanthemum seticuspe, a model species for hexaploid cultivated chrysanthemum.

Chrysanthemums are one of the most industrially important cut flowers worldwide. However, their segmental allopolyploidy and self-incompatibility have prevented the application of genetic analysis and modern breeding strategies. We thus developed a model strain, Gojo-0 (Chrysanthemum seticuspe), which is a diploid and self-compatible pure line. Here, we present the 3.05 Gb chromosome-level reference genome sequence, which covered 97% of the C. seticuspe genome. The genome contained more than 80% interspersed repeats, of which retrotransposons accounted for 72%. We identified recent segmental duplication and retrotransposon expansion in C. seticuspe, contributing to arelatively large genome size. Furthermore, we identified a retrotransposon family, SbdRT, which was enriched in gene-dense genome regions and had experienced a very recent transposition burst. We also demonstrated that the chromosome-level genome sequence facilitates positional cloning in C. seticuspe. The genome sequence obtained here can greatly contribute as a reference for chrysanthemum in front-line breeding including genome editing.

Quantitative trait loci for growth-related traits in Japanese quail (Coturnix japonica) using restriction-site associated DNA sequencing.

This study aimed to identify quantitative trait loci (QTLs) for growth-related traits by constructing a genetic linkage map based on single nucleotide polymorphism (SNP) markers in Japanese quail. A QTL mapping population of 277 F2 birds was obtained from an intercross between a male of a large-sized strain and three females of a normal-sized strain. Body weight (BW) was measured weekly from hatching to 16 weeks of age. Non-linear regression growth models of Weibull, Logistic, Gompertz, Richards, and Brody were analyzed, and growth curve parameters of Richards was selected as the best model to describe the quail growth curve of the F2 birds. Restriction-site associated DNA sequencing developed 125 SNP markers that were informative between their parental strains. The SNP markers were distributed on 16 linkage groups that spanned 795.9 centiMorgan (cM) with an average marker interval of 7.3 cM. QTL analysis of phenotypic traits revealed four main-effect QTLs. Detected QTLs were located on chromosomes 1 and 3 and were associated with BW from 4 to 16 weeks of age and asymptotic weight of Richards model at genome-wide significant at 1% or 5% level. No QTL was detected for BW from 0 to 3 weeks of age. This is the first report identified QTLs for asymptotic weight of the Richards parameter in Japanese quail. These results highlight that the combination of QTL studies and the RAD-seq method will aid future breeding programs identify genes underlying the QTL and the application of marker-assisted selection in the poultry industry, particularly the Japanese quail.

Mapping of Quantitative Trait Loci Controlling Egg-Quality and -Production Traits in Japanese Quail (Coturnix japonica) Using Restriction-Site Associated DNA Sequencing.

This research was conducted to identify quantitative trait loci (QTL) associated with egg-related traits by constructing a genetic linkage map based on single nucleotide polymorphism (SNP) markers using restriction-site associated DNA sequencing (RAD-seq) in Japanese quail. A total of 138 F2 females were produced by full-sib mating of F1 birds derived from an intercross between a male of the large-sized strain with three females of the normal-sized strain. Eggs were investigated at two different stages: the beginning stage of egg-laying and at 12 weeks of age (second stage). Five eggs were analyzed for egg weight, lengths of the long and short axes, egg shell strength and weight, yolk weight and diameter, albumen weight, egg equator thickness, and yolk color (L*, a*, and b* values) at each stage. Moreover, the age at first egg, the cumulative number of eggs laid, and egg production rate were recorded. RAD-seq developed 118 SNP markers and mapped them to 13 linkage groups using the Map Manager QTX b20 software. Markers were spanned on 776.1 cM with an average spacing of 7.4 cM. Nine QTL were identified on chromosomes 2, 4, 6, 10, 12, and Z using the simple interval mapping method in the R/qtl package. The QTL detected affected 10 egg traits of egg weight, lengths of the long and short axes of egg, egg shell strength, yolk diameter and weight, albumen weight, and egg shell weight at the beginning stage, yellowness of the yolk color at the second stage, and age at first egg. This is the first report to perform a quail QTL analysis of egg-related traits using RAD-seq. These results highlight the effectiveness of RAD-seq associated with targeted QTL and the application of marker-assisted selection in the poultry industry, particularly in the Japanese quail.

RAD-seq-Based High-Density Linkage Map Construction and QTL Mapping of Biomass-Related Traits in Sorghum using the Japanese Landrace Takakibi NOG.

Sorghum [Sorghum bicolor (L.) Moench] grown locally by Japanese farmers is generically termed Takakibi, although its genetic diversity compared with geographically distant varieties or even within Takakibi lines remains unclear. To explore the genomic diversity and genetic traits controlling biomass and other physiological traits in Takakibi, we focused on a landrace, NOG, in this study. Admixture analysis of 460 sorghum accessions revealed that NOG belonged to the subgroup that represented Asian sorghums, and it was only distantly related to American/African accessions including BTx623. In an attempt to dissect major traits related to biomass, we generated a recombinant inbred line (RIL) from a cross between BTx623 and NOG, and we constructed a high-density linkage map based on 3,710 single-nucleotide polymorphisms obtained by restriction-site-associated DNA sequencing of 213 RIL individuals. Consequently, 13 fine quantitative trait loci (QTLs) were detected on chromosomes 2, 3, 6, 7, 8 and 9, which included five QTLs for days to heading, three for plant height (PH) and total shoot fresh weight and two for Brix. Furthermore, we identified two dominant loci for PH as being identical to the previously reported dw1 and dw3. Together, these results corroborate the diversified genome of Japanese Takakibi, while the RIL population and high-density linkage map generated in this study will be useful for dissecting other important traits in sorghum.

A pure line derived from a self-compatible Chrysanthemum seticuspe mutant as a model strain in the genus Chrysanthemum.

Asteraceae is the largest family of angiosperms, comprising approximately 24,000 species. Molecular genetic studies of Asteraceae are essential for understanding plant diversity. Chrysanthemum morifolium is the most industrially important ornamental species in Asteraceae. Most cultivars of C. morifolium are autohexaploid and self-incompatible. These properties are major obstacles to the genetic analysis and modern breeding of C. morifolium. Furthermore, high genome heterogeneity complicates molecular biological analyses. In this study, we developed a model strain in the genus Chrysanthemum. C. seticuspe is a diploid species with a similar flowering property and morphology to C. morifolium and can be subjected to Agrobacterium-mediated transformation. We isolated a natural self-compatible mutant of C. seticuspe and established a pure line through repeated selfing and selection. The resultant strain, named Gojo-0, was favorable for genetic analyses, including isolation of natural and induced mutants, and facilitated molecular biological analysis, including whole genome sequencing, owing to the simplicity and homogeneity of its genome. Interspecific hybridization with Chrysanthemum species was possible, enabling molecular genetic analysis of natural interspecific variations. The accumulation of research results and resources using Gojo-0 as a platform is expected to promote molecular genetic studies on the genus Chrysanthemum and the genetic improvement of chrysanthemum cultivars.

De novo whole-genome assembly in Chrysanthemum seticuspe, a model species of Chrysanthemums, and its application to genetic and gene discovery analysis.

Cultivated chrysanthemum (Chrysanthemum morifolium Ramat.) is one of the most economically important ornamental crops grown worldwide. It has a complex hexaploid genome (2n = 6x = 54) and large genome size. The diploid Chrysanthemum seticuspe is often used as a model of cultivated chrysanthemum, since the two species are closely related. To expand our knowledge of the cultivated chrysanthemum, we here performed de novo whole-genome assembly in C. seticuspe using the Illumina sequencing platform. XMRS10, a C. seticuspe accession developed by five generations of self-crossing from a self-compatible strain, AEV2, was used for genome sequencing. The 2.72 Gb of assembled sequences (CSE_r1.0), consisting of 354,212 scaffolds, covered 89.0% of the 3.06 Gb C. seticuspe genome estimated by k-mer analysis. The N50 length of scaffolds was 44,741 bp. For protein-encoding genes, 71,057 annotated genes were deduced (CSE_r1.1_cds). Next, based on the assembled genome sequences, we performed linkage map construction, gene discovery and comparative analyses for C. seticuspe and cultivated chrysanthemum. The generated C. seticuspe linkage map revealed skewed regions in segregation on the AEV2 genome. In gene discovery analysis, candidate flowering-related genes were newly found in CSE_r1.1_cds. Moreover, single nucleotide polymorphism identification and annotation on the C. × morifolium genome showed that the C. seticuspe genome was applicable to genetic analysis in cultivated chrysanthemums. The genome sequences assembled herein are expected to contribute to future chrysanthemum studies. In addition, our approach demonstrated the usefulness of short-read genome assembly and the importance of choosing an appropriate next genome sequencing technology based on the purpose of the post-genome analysis.

MITE insertion-dependent expression of CitRKD1 with a RWP-RK domain regulates somatic embryogenesis in citrus nucellar tissues.

BACKGROUND: Somatic embryogenesis in nucellar tissues is widely recognized to induce polyembryony in major citrus varieties such as sweet oranges, satsuma mandarins and lemons. This capability for apomixis is attractive in agricultural production systems using hybrid seeds, and many studies have been performed to elucidate the molecular mechanisms of various types of apomixis. To identify the gene responsible for somatic embryogenesis in citrus, a custom oligo-DNA microarray including predicted genes in the citrus polyembryonic locus was used to compare the expression profiles in reproductive tissues between monoembryonic and polyembryonic varieties. The full length of CitRKD1, which was identified as a candidate gene responsible for citrus somatic embryogenesis, was isolated from satsuma mandarin and its molecular function was investigated using transgenic 'Hamlin' sweet orange by antisense-overexpression. RESULTS: The candidate gene CitRKD1, predominantly transcribed in reproductive tissues of polyembryonic varieties, is a member of the plant RWP-RK domain-containing protein. CitRKD1 of satsuma mandarin comprised two alleles (CitRKD1-mg1 and CitRKD1-mg2) at the polyembryonic locus controlling embryonic type (mono/polyembryony) that were structurally divided into two types with or without a miniature inverted-repeat transposable element (MITE)-like insertion in the upstream region. CitRKD1-mg2 with the MITE insertion was the predominant transcript in flowers and young fruits where somatic embryogenesis of nucellar cells occurred. Loss of CitRKD1 function by antisense-overexpression abolished somatic embryogenesis in transgenic sweet orange and the transgenic T1 plants were confirmed to derive from zygotic embryos produced by self-pollination by DNA diagnosis. Genotyping PCR analysis of 95 citrus traditional and breeding varieties revealed that the CitRKD1 allele with the MITE insertion (polyembryonic allele) was dominant and major citrus varieties with the polyembryonic allele produced polyembryonic seeds. CONCLUSION: CitRKD1 at the polyembryonic locus plays a principal role in regulating citrus somatic embryogenesis. CitRKD1 comprised multiple alleles that were divided into two types, polyembryonic alleles with a MITE insertion in the upstream region and monoembryonic alleles without it. CitRKD1 was transcribed in reproductive tissues of polyembryonic varieties with the polyembryonic allele. The MITE insertion in the upstream region of CitRKD1 might be involved in regulating the transcription of CitRKD1.

Draft Genome Sequence of Mn(II)-Oxidizing Pseudomonas resinovorans Strain MO-1.

Pseudomonas resinovorans strain MO-1, which possesses a high ability to oxidize Mn(II), has been isolated from oligotrophic pond sediment. The draft genome sequence consists of 6,252,942 bp and has a G+C content of 63.4%. Strain MO-1 has 5,694 coding sequences, including 13 putative Mn(II) oxidation genes.

Impairment of Lhca4, a subunit of LHCI, causes high accumulation of chlorophyll and the stay-green phenotype in rice.

Chlorophyll is an essential molecule for acquiring light energy during photosynthesis. Mutations that result in chlorophyll retention during leaf senescence are called 'stay-green' mutants. One of the several types of stay-green mutants, Type E, accumulates high levels of chlorophyll in the pre-senescent leaves, resulting in delayed yellowing. We isolated delayed yellowing1-1 (dye1-1), a rice mutant whose yellowing is delayed in the field. dye1-1 accumulated more chlorophyll than the wild-type in the pre-senescent and senescent leaves, but did not retain leaf functionality in the 'senescent green leaves', suggesting that dye1-1 is a Type E stay-green mutant. Positional cloning revealed that DYE1 encodes Lhca4, a subunit of the light-harvesting complex I (LHCI). In dye1-1, amino acid substitution occurs at the location of a highly conserved amino acid residue involved in pigment binding; indeed, a severely impaired structure of the PSI-LHCI super-complex in dye1-1 was observed in a blue native PAGE analysis. Nevertheless, the biomass and carbon assimilation rate of dye1-1 were comparable to those in the wild-type. Interestingly, Lhcb1, a trimeric LHCII protein, was highly accumulated in dye1-1, in the chlorophyll-protein complexes. The high accumulation of LHCII in the LHCI mutant dye1 suggests a novel functional interaction between LHCI and LHCII.

The Non-Mendelian Green Cotyledon Gene in Soybean Encodes a Small Subunit of Photosystem II.

Chlorophyll degradation plays important roles in leaf senescence including regulation of degradation of chlorophyll-binding proteins. Although most genes encoding enzymes of the chlorophyll degradation pathway have been identified, the regulation of their activity has not been fully understood. Green cotyledon mutants in legume are stay-green mutants, in which chlorophyll degradation is impaired during leaf senescence and seed maturation. Among them, the soybean (Glycine max) green cotyledon gene cytG is unique because it is maternally inherited. To isolate cytG, we extensively sequenced the soybean chloroplast genome, and detected a 5-bp insertion causing a frame-shift in psbM, which encodes one of the small subunits of photosystem II. Mutant tobacco plants (Nicotiana tabacum) with a disrupted psbM generated using a chloroplast transformation technique had green senescent leaves, confirming that cytG encodes PsbM. The phenotype of cytG was very similar to that of mutant of chlorophyll b reductase catalyzing the first step of chlorophyll b degradation. In fact, chlorophyll b-degrading activity in dark-grown cytG and psbM-knockout seedlings was significantly lower than that of wild-type plants. Our results suggest that PsbM is a unique protein linking photosynthesis in presenescent leaves with chlorophyll degradation during leaf senescence and seed maturation. Additionally, we discuss the origin of cytG, which may have been selected during domestication of soybean.

Overexpression of a citrus basic helix-loop-helix transcription factor (CubHLH1), which is homologous to Arabidopsis activation-tagged bri1 suppressor 1 interacting factor genes, modulates carotenoid metabolism in transgenic tomato.

To explore the transcription factors associated with carotenoid metabolism in citrus fruit, one transcription factor (CubHLH1) was selected through microarray screening in Satsuma mandarin (Citrus unshiu Marc.) fruit, which was treated with exogenous ethylene or gibberellin (GA), accelerating or retarding carotenoid accumulation in peel, respectively. The amino acid sequence of CubHLH1 has homology to Arabidopsis activation-tagged bri1 suppressor 1 (ATBS1) interacting factor (AIF), which is functionally characterized as a negative regulator of the brassinolide (BR) signalling pathway. Yeast two-hybrid analysis revealed that protein for CubHLH1 could interact with Arabidopsis and tomato ATBS1. Overexpression of CubHLH1 caused a dwarf phenotype in transgenic tomato (Solanum lycopersicum L.), suggesting that CubHLH1 has a similar function to Arabidopsis AIF. In the transgenic tomato fruit at ripening stage, the lycopene content was reduced along with the changes in carotenoid biosynthetic gene expression. The abscisic acid (ABA) content of all the transgenic tomato fruit was higher than that of the wild type. These results implied that CubHLH1 is considered to have a similar function to Arabidopsis AIFs and might be directly involved in carotenoid metabolism in mature citrus fruit.

A green-cotyledon/stay-green mutant exemplifies the ancient whole-genome duplications in soybean.

The recent whole-genome sequencing of soybean (Glycine max) revealed that soybean experienced whole-genome duplications 59 million and 13 million years ago, and it has an octoploid-like genome in spite of its diploid nature. We analyzed a natural green-cotyledon mutant line, Tenshin-daiseitou. The physiological analysis revealed that Tenshin-daiseitou shows a non-functional stay-green phenotype in senescent leaves, which is similar to that of the mutant of Mendel's green-cotyledon gene I, the ortholog of SGR in pea. The identification of gene mutations and genetic segregation analysis suggested that defects in GmSGR1 and GmSGR2 were responsible for the green-cotyledon/stay-green phenotype of Tenshin-daiseitou, which was confirmed by RNA interference (RNAi) transgenic soybean experiments using GmSGR genes. The characterized green-cotyledon double mutant d1d2 was found to have the same mutations, suggesting that GmSGR1 and GmSGR2 are D1 and D2. Among the examined d1d2 strains, the d1d2 strain K144a showed a lower Chl a/b ratio in mature seeds than other strains but not in senescent leaves, suggesting a seed-specific genetic factor of the Chl composition in K144a. Analysis of the soybean genome sequence revealed four genomic regions with microsynteny to the Arabidopsis SGR1 region, which included the GmSGR1 and GmSGR2 regions. The other two regions contained GmSGR3a/GmSGR3b and GmSGR4, respectively, which might be pseudogenes or genes with a function that is unrelated to Chl degradation during seed maturation and leaf senescence. These GmSGR genes were thought to be produced by the two whole-genome duplications, and they provide a good example of such whole-genome duplication events in the evolution of the soybean genome.

Increased levels of IAA are required for system 2 ethylene synthesis causing fruit softening in peach (Prunus persica L. Batsch).

The fruit of melting-flesh peach (Prunus persica L. Batsch) cultivars produce high levels of ethylene caused by high expression of PpACS1 (an isogene of 1-aminocyclopropane-1-carboxylic acid synthase), resulting in rapid fruit softening at the late-ripening stage. In contrast, the fruit of stony hard peach cultivars do not soften and produce little ethylene due to low expression of PpACS1. To elucidate the mechanism for suppressing PpACS1 expression in stony hard peaches, a microarray analysis was performed. Several genes that displayed similar expression patterns as PpACS1 were identified and shown to be indole-3-acetic acid (IAA)-inducible genes (Aux/IAA, SAUR). That is, expression of IAA-inducible genes increased at the late-ripening stage in melting flesh peaches; however, these transcripts were low in mature fruit of stony hard peaches. The IAA concentration increased suddenly just before harvest time in melting flesh peaches exactly coinciding with system 2 ethylene production. In contrast, the IAA concentration did not increase in stony hard peaches. Application of 1-naphthalene acetic acid, a synthetic auxin, to stony hard peaches induced a high level of PpACS1 expression, a large amount of ethylene production and softening. Application of an anti-auxin, α-(phenylethyl-2-one)-IAA, to melting flesh peaches reduced levels of PpACS1 expression and ethylene production. These observations indicate that suppression of PpACS1 expression at the late-ripening stage of stony hard peach may result from a low level of IAA and that a high concentration of IAA is required to generate a large amount of system 2 ethylene in peaches.

Characterization of genomic sequence showing strong association with polyembryony among diverse Citrus species and cultivars, and its synteny with Vitis and Populus.

Polyembryony, in which multiple somatic nucellar cell-derived embryos develop in addition to the zygotic embryo in a seed, is common in the genus Citrus. Previous genetic studies indicated polyembryony is mainly determined by a single locus, but the underlying molecular mechanism is still unclear. As a step towards identification and characterization of the gene or genes responsible for nucellar embryogenesis in Citrus, haplotype-specific physical maps around the polyembryony locus were constructed. By sequencing three BAC clones aligned on the polyembryony haplotype, a single contiguous draft sequence consisting of 380 kb containing 70 predicted open reading frames (ORFs) was reconstructed. Single nucleotide polymorphism genotypes detected in the sequenced genomic region showed strong association with embryo type in Citrus, indicating a common polyembryony locus is shared among widely diverse Citrus cultivars and species. The arrangement of the predicted ORFs in the characterized genomic region showed high collinearity to the genomic sequence of chromosome 4 of Vitis vinifera and linkage group VI of Populus trichocarpa, suggesting that the syntenic relationship among these species is conserved even though V. vinifera and P. trichocarpa are non-apomictic species. This is the first study to characterize in detail the genomic structure of an apomixis locus determining adventitious embryony.