BackgroundThe honey bee is an important model system for increasing understanding of molecular and neural mechanisms underlying social behaviors relevant to the agricultural industry and basic science. The western honey bee, Apis mellifera, has served as a model species, and its genome sequence has been published. In contrast, the genome of the Asian honey bee, Apis cerana, has not yet been sequenced. A. cerana has been raised in Asian countries for thousands of years and has brought considerable economic benefits to the apicultural industry. A cerana has divergent biological traits compared to A. mellifera and it has played a key role in maintaining biodiversity in eastern and southern Asia. Here we report the first whole genome sequence of A. cerana.ResultsUsing de novo assembly methods, we produced a 238 Mbp draft of the A. cerana genome and generated 10,651 genes. A.cerana-specific genes were analyzed to better understand the novel characteristics of this honey bee species. Seventy-two percent of the A. cerana-specific genes had more than one GO term, and 1,696 enzymes were categorized into 125 pathways. Genes involved in chemoreception and immunity were carefully identified and compared to those from other sequenced insect models. These included 10 gustatory receptors, 119 odorant receptors, 10 ionotropic receptors, and 160 immune-related genes.ConclusionsThis first report of the whole genome sequence of A. cerana provides resources for comparative sociogenomics, especially in the field of social insect communication. These important tools will contribute to a better understanding of the complex behaviors and natural biology of the Asian honey bee and to anticipate its future evolutionary trajectory.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-16-1) contains supplementary material, which is available to authorized users.

Although polyploidy is common in plants and some animals, mechanisms for functional divergence between homoeologous genes are poorly understood. MYB2 gene promotes cotton fibre development and is functionally homologous to Arabidopsis GLABROUS1 (GL1) in trichome formation. The most widely cultivated cotton is an allotetraploid (Gossypium hirsutum, AADD) that contains GhMYB2A and GhMYB2D homoeologs. Here we show that GhMYB2D mRNA accumulates more than GhMYB2A during fibre initiation and is targeted by miR828 and miR858, generating trans-acting siRNAs (ta-siRNAs) in the TAS4 family. Overexpressing GhMYB2A but not GhMYB2D complements the gl1 phenotype. Mutating the miR828-binding site or replacing its downstream sequence in GhMYB2D abolishes ta-siRNA production and restores trichome development in gl1 mutants. Moreover, disrupting Dicer-like protein 4 or RDR6, the biogenesis genes for ta-siRNAs, in the gl1 GhMYB2D overexpressors restores trichome development. These data support a unique role for microRNAs in functional divergence between target homoeologous genes that are important for evolution and selection of morphological traits.

BackgroundThe morphogenesis of single-celled cotton fiber includes extreme elongation and staged cell wall differentiation. Designing strategies for improving cotton fiber for textiles and other uses relies on uncovering the related regulatory mechanisms. In this research we compared the transcriptomes and metabolomes of two Gossypium genotypes, Gossypium barbadense cv Phytogen 800 and G. hirsutum cv Deltapine 90. When grown in parallel, the two types of fiber developed similarly except for prolonged fiber elongation in the G. barbadense cultivar. The data were collected from isolated fibers between 10 to 28 days post anthesis (DPA) representing: primary wall synthesis to support elongation; transitional cell wall remodeling; and secondary wall cellulose synthesis, which was accompanied by continuing elongation only in G. barbadense fiber.ResultsOf 206 identified fiber metabolites, 205 were held in common between the two genotypes. Approximately 38,000 transcripts were expressed in the fiber of each genotype, and these were mapped to the reference set and interpreted by homology to known genes. The developmental changes in the transcriptomes and the metabolomes were compared within and across genotypes with several novel implications. Transitional cell wall remodeling is a distinct stable developmental stage lasting at least four days (18 to 21 DPA). Expression of selected cell wall related transcripts was similar between genotypes, but cellulose synthase gene expression patterns were more complex than expected. Lignification was transcriptionally repressed in both genotypes. Oxidative stress was lower in the fiber of G. barbadense cv Phytogen 800 as compared to G. hirsutum cv Deltapine 90. Correspondingly, the G. barbadense cultivar had enhanced capacity for management of reactive oxygen species during its prolonged elongation period, as indicated by a 138-fold increase in ascorbate concentration at 28 DPA.ConclusionsThe parallel data on deep-sequencing transcriptomics and non-targeted metabolomics for two genotypes of single-celled cotton fiber showed that a discrete developmental stage of transitional cell wall remodeling occurs before secondary wall cellulose synthesis begins. The data showed how lignification can be transcriptionally repressed during secondary cell wall synthesis, and they implicated enhanced capacity to manage reactive oxygen species through the ascorbate-glutathione cycle as a positive contributor to fiber length.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1708-9) contains supplementary material, which is available to authorized users.

Summary• Flowering time is an important adaptive trait and varies among Arabidopsis thaliana and its related species, including allopolyploids that are formed between A. thaliana and Arabidopsis arenosa. FLOWERING LOCUS C (FLC) inhibits early flowering in A. thaliana. A previous study showed that late-flowering A. arenosa contained two or more FLC alleles that were differentially expressed in Arabidopsis allotetraploids, but the genomic organization and evolution of FLC locus were unknown.• Comparative sequence and evolutionary analyses were performed in FLC-containing genomic regions in A. thaliana, A. arenosa and Arabidopsis lyrata, and expression of FLC loci and alleles was examined in Arabidopsis allopolyploids.• The FLC locus was tandemly duplicated in A. lyrata and triplicated in A. arenosa, and the tandem duplication event occurred after divergence from A. thaliana. Although FLC duplicates were highly conserved, their upstream sequences rapidly diverged. The third FLC copy in A. arenosa acquired a new splicing site through a point mutation in the intron and generated the new exon followed by an early stop codon, resulting in a novel MADS box gene.• Flowering time variation in Arabidopsis allopolyploids is probably related to the expression diversity and ⁄ or copy number of multiple FLC loci. Moreover, exonization of intronic sequence is a mechanism for the origin of new genes.

Single nucleotide polymorphisms (SNPs) play important roles as molecular markers in plant genomics and breeding studies. Although onion (Allium cepa L.) is an important crop globally, relatively few molecular marker resources have been reported due to its large genome and high heterozygosity. Genotyping-by-sequencing (GBS) offers a greater degree of complexity reduction followed by concurrent SNP discovery and genotyping for species with complex genomes. In this study, GBS was employed for SNP mining in onion, which currently lacks a reference genome. A segregating F2 population, derived from a cross between ‘NW-001’ and ‘NW-002,’ as well as multiple parental lines were used for GBS analysis. A total of 56.15 Gbp of raw sequence data were generated and 1,851,428 SNPs were identified from the de novo assembled contigs. Stringent filtering resulted in 10,091 high-fidelity SNP markers. Robust SNPs that satisfied the segregation ratio criteria and with even distribution in the mapping population were used to construct an onion genetic map. The final map contained eight linkage groups and spanned a genetic length of 1,383 centiMorgans (cM), with an average marker interval of 8.08 cM. These robust SNPs were further analyzed using the high-throughput Fluidigm platform for marker validation. This is the first study in onion to develop genome-wide SNPs using GBS. The resulting SNP markers and developed linkage map will be valuable tools for genetic mapping of important agronomic traits and marker-assisted selection in onion breeding programs.