Search

Safflower (Carthamus tinctorius) is a member of the Asteraceae family that is grown in temperate climates as an oil seed crop. Most commercially grown safflower varieties can be sown in late winter or early spring and flower rapidly in the absence of overwintering. There are winter-hardy safflower accessions that can be sown in autumn and survive over-wintering. Here, we show that a winter-hardy safflower possesses a vernalization response, whereby flowering is accelerated by exposing germinating seeds to prolonged cold. The impact of vernalization was quantitative, such that increasing the duration of cold treatment accelerated flowering to a greater extent, until the response was saturated after 2 weeks exposure to low-temperatures. To investigate the molecular-basis of the vernalization-response in safflower, transcriptome activity was compared and contrasted between vernalized versus non-vernalized plants, in both ‘winter hardy’ and ‘spring’ cultivars. These genome-wide expression analyses identified a small set of transcripts that are both differentially expressed following vernalization and that also have different expression levels in the spring versus winter safflowers. Four of these transcripts were quantitatively induced by vernalization in a winter hardy safflower but show high basal levels in spring safflower. Phylogenetic analyses confidently assigned that the nucleotide sequences of the four differentially expressed transcripts are related to FLOWERING LOCUS T (FT), FRUITFUL (FUL), and two genes within the MADS-like clade genes. Gene models were built for each of these sequences by assembling an improved safflower reference genome using PacBio-based long-read sequencing, covering 85% of the genome, with N50 at 594,000 bp in 3000 contigs. Possible evolutionary relationships between the vernalization response of safflower and those of other plants are discussed.

Background The legume cowpea (Vigna unguiculata L.) is extensively grown in sub-Saharan Africa. Cowpea, like many legumes has proved recalcitrant to plant transformation. A rapid transient leaf assay was developed for testing gene expression and editing constructs prior to stable cowpea transformation, to accelerate cowpea and legume crop improvement. Results Attempts to develop a transient protoplast system for cowpea were unsuccessful. Leaflets from plants 3–4 weeks post-germination were age selected to establish a rapid Agrobacterium (Agro) infiltration-mediated transient system for efficacy testing of gene expression and CRISPR/Cas9 gene editing constructs. In planta, Agro-infiltration of leaflets with fluorescent expression constructs, resulted in necrosis. By contrast, Agro-infiltration of detached leaflets with an Arabidopsis (At) ubiquitin3 promoter:ZsGreen construct, followed by culture on solid nutrient medium resulted in fluorescence in over 48% of leaf cells. Expression efficiency was leaf age-dependent. Three cowpea meiosis genes were identified for CRISPR/Cas9 gene-editing, with the forward aim of meiosis-knock out for asexual seed induction in cowpea. Constructs were designed and tested containing candidate gene-specific guide RNAs, expressed using either the cowpea or Arabidopsis U6 promoters with Cas9 expression directed by either the Arabidopsis 40S ribosomal protein or parsley ubiquitin4-2 promoters. Leaflets were infiltrated with test gene-editing constructs and analytical methods developed to identify gene-specific mutations. A construct that produced mutations predicted to induce functional knockout of in the VuSPO11-1 meiosis gene was tested for efficacy in primary transgenic cowpea plants using a previously established stable transformation protocol. Vuspo11-1 mutants were identified, that cytologically phenocopied spo11-1 mutants previously characterized in Arabidopsis, and rice. Importantly, a biallelic male and female sterile mutant was identified in primary transgenics, exhibiting the expected defects in 100% of examined male and female meiocytes. Conclusion The transient, detached cowpea leaf assay, and supporting analytical methods developed, provide a rapid and reproducible means for testing gene expression constructs, and constructs for inducing mutagenesis in genes involved in both vegetative and reproductive developmental programs. The method and tested editing constructs and components have potential application for a range of crop legumes.

Apomixis results in asexual seed formation where progeny are identical to the maternal plant. In ovules of apomictic species of the Hieracium subgenus Pilosella, meiosis of the megaspore mother cell generates four megaspores. Aposporous initial (AI) cells form during meiosis in most ovules. The sexual pathway terminates during functional megaspore (FM) differentiation, when an enlarged AI undergoes mitosis to form an aposporous female gametophyte. Then, the mitotically programmed FM dies along with the three other megaspores by unknown mechanisms. Transcriptomes of laser-dissected AIs, ovule cells, and ovaries from apomicts and AI-deficient mutants were analyzed to understand the pathways involved. The steps leading to AI mitosis and sexual pathway termination were determined using antibodies against arabinogalactan protein epitopes found to mark both sexual and aposporous female gametophyte lineages at inception. At most, four AIs differentiated near developing megaspores. The first expanding AI cell to contact the FM formed a functional AI that underwent mitosis soon after megaspore degeneration. Transcriptome analyses indicated that the enlarged, laser-captured AIs were arrested in the S/G2 phase of the cell cycle and were metabolically active. Further comparisons with AI-deficient mutants showed that AIs were enriched in transcripts encoding homologs of genes involved in, and potentially antagonistic to, known FM specification pathways. We propose that AI and FM cell contact provides cues required for AI mitosis and megaspore degeneration. Specific candidates to further interrogate AI-FM interactions were identified here and include Hieracium arabinogalactan protein family genes.

Molecular knowledge of pathways regulating seed formation in legumes, remains scarce. Thirteen isolated cell-type transcriptomes were developed, spanning temporal events of male and female gametogenesis and seed initiation, to examine pathways involved in cowpea seed formation.In situhybridization confirmed localization ofin silicoidentified cell-specific genes, verifying transcriptome utility. Cowpea andArabidopsisreproductive cells showed some conservation in regulators enabling cell-type expression as some cowpea cell-specific genes promoters and theirArabidopsishomologs directed expression to identical reproductive cell-types in transgenic plants.In silicoanalyses revealed gene expression similarities and differences with genes in pathways regulating reproductive events in other plants. Meiosis-related genes were expressed at mitotic stages of gametogenesis and during sporophytic development in cowpea. Plant hormone pathways showing preferential expression at particular reproductive stages were identified. Expression of epigenetic pathways, resembling those found inArabidopsis,including microRNA mediated gene silencing, RNA directed DNA methylation and histone modification were associated with particular stages of male and female gametophyte development, suggesting roles in gametogenic cell specification and elaboration. Analyses of cell-cycle related gene expression in mature cowpea female gametophytes, indicated that the egg and central cell were arrested at the G1/S and G2/M cell cycle phases, respectively, prior to fertilization. Pre-fertilization female gametophyte arrest was characterized by barely detectable auxin biosynthesis gene expression levels, and elevated expression of genes involved in RNA-mediated gene silencing and histone modification. These transcriptomes provide a useful resource for additional interrogation to support functional analyses for development of higher yielding cowpea and syntenic legume crops.One sentence summaryAnalyses of laser capture derived cell-type transcriptomes spanning meiosis to seed initiation revealed gene expression profiles during cell specification and reproductive development in cowpea.

The legume cowpea (Vigna unguiculata, 2n=2x=22) has significant tolerance to drought and heat stress. Here we analysed and manipulated cowpea centromere-specific histone H3 (CENH3) genes, aiming to establish a centromere-based doubled-haploid method for use in genetic improvement of this dryland crop in future. Cowpea encodes two functional CENH3 variants (CENH3.1 and CENH3.2) and two CENH3 pseudogenes. Phylogenetic analysis suggests that gene duplication of CENH3 occurred independently during the speciation ofV. unguiculataand the relatedV. mungowithout a genome duplication event. Both functional cowpeaCENH3variants are transcribed, and the corresponding proteins are intermingled in subdomains of different types of centromere sequences in a tissue-specific manner together with the outer kinetochore protein CENPC. CENH3.2 is removed from the generative cell of mature pollen, while CENH3.1 persists. Differences between both CENH3 paralogs are restricted to the N-terminus. The complete CRISPR/Cas9-based inactivation ofCENH3.1resulted in delayed vegetative growth and sterility, indicating that this variant is needed for plant development and reproduction. By contrast,CENH3.2knockout individuals did not show obvious defects during vegetative and reproductive development, suggesting that the gene is an early stage of subfunctionalization or pseudogenization.One-sentence summaryThe two paralogous centromeric histones (CENH3) of cowpea contribute unequal to the function of the centromere.

High-throughput genotyping arrays provide a standardized resource for plant breeding communities that are useful for a breadth of applications including high-density genetic mapping, genome-wide association studies (GWAS), genomic selection (GS), complex trait dissection, and studying patterns of genomic diversity among cultivars and wild accessions. We have developed the CottonSNP63K, an Illumina Infinium array containing assays for 45,104 putative intraspecific single nucleotide polymorphism (SNP) markers for use within the cultivated cotton species Gossypium hirsutum L. and 17,954 putative interspecific SNP markers for use with crosses of other cotton species with G. hirsutum. The SNPs on the array were developed from 13 different discovery sets that represent a diverse range of G. hirsutum germplasm and five other species: G. barbadense L., G. tomentosum Nuttal × Seemann, G. mustelinum Miers × Watt, G. armourianum Kearny, and G. longicalyx J.B. Hutchinson and Lee. The array was validated with 1,156 samples to generate cluster positions to facilitate automated analysis of 38,822 polymorphic markers. Two high-density genetic maps containing a total of 22,829 SNPs were generated for two F2 mapping populations, one intraspecific and one interspecific, and 3,533 SNP markers were co-occurring in both maps. The produced intraspecific genetic map is the first saturated map that associates into 26 linkage groups corresponding to the number of cotton chromosomes for a cross between two G. hirsutum lines. The linkage maps were shown to have high levels of collinearity to the JGI G. raimondii Ulbrich reference genome sequence. The CottonSNP63K array, cluster file and associated marker sequences constitute a major new resource for the global cotton research community.

� Background and Aims Apomixis in plants generates clonal progeny with a maternal genotype through asexual seed formation. Hieracium subgenus Pilosella (Asteraceae) contains polyploid, highly heterozygous apomictic and sexual species. Within apomictic Hieracium, dominant genetic loci independently regulate the qualitative developmental components of apomixis. In H. praealtum, LOSS OF APOMEIOSIS (LOA) enables formation of embryo sacs without meiosis and LOSS OF PARTHENOGENESIS (LOP) enables fertilization-independent seed formation. A locus required for fertilization-independent endosperm formation (AutE) has been identified in H. piloselloides. Additional quantitative loci appear to influence the penetrance of the qualitative loci, although the controlling genes remain unknown. This study aimed to develop the first genetic linkage maps for sexual and apomictic Hieracium species using simple sequence repeat (SSR) markers derived from expressed transcripts within the developing ovaries. � Methods RNA from microdissected Hieracium ovule cell types and ovaries was sequenced and SSRs were identified. Two different F1 mapping populations were created to overcome difficulties associated with genome complexity and asexual reproduction. SSR markers were analysed within each mapping population to generate draft linkage maps for apomictic and sexual Hieracium species. � Key Results A collection of 14 684 Hieracium expressed SSR markers were developed and linkage maps were constructed for Hieracium species using a subset of the SSR markers. Both the LOA and LOP loci were successfully assigned to linkage groups; however, AutE could not be mapped using the current populations. Comparisons with lettuce (Lactuca sativa) revealed partial macrosynteny between the two Asteraceae species. � Conclusions A collection of SSR markers and draft linkage maps were developed for two apomictic and one sexual Hieracium species. These maps will support cloning of controlling genes at LOA and LOP loci in Hieracium and should also assist with identification of quantitative loci that affect the expressivity of apomixis. Future work will focus on mapping AutE using alternative populations.

Varietal single nucleotide polymorphisms (SNPs) are the differences within one of the two subgenomes between different tetraploid cotton varieties and have not been practically used in cotton genetics and breeding because they are difficult to identify due to low genetic diversity and very high sequence identity between homeologous genes in cotton. We have used transcriptome and restriction site−associated DNA sequencing to identify varietal SNPs among 18 G. hirsutum varieties based on the rationale that varietal SNPs can be more confidently called when flanked by subgenome-specific SNPs. Using transcriptome data, we successfully identified 37,413 varietal SNPs and, of these, 22,121 did not have an additional varietal SNP within their 20-bp flanking regions so can be used in most SNP genotyping assays. From restriction site−associated DNA sequencing data, we identified an additional 3090 varietal SNPs between two of the varieties. Of the 1583 successful SNP assays achieved using different genotyping platforms, 1363 were verified. Many of the SNPs behaved as dominant markers because of coamplification from homeologous loci, but the number of SNPs acting as codominant markers increased when one or more subgenome-specific SNP(s) were incorporated in their assay primers, giving them greater utility for breeding applications. A G. hirsutum genetic map with 1244 SNP markers was constructed covering 5557.42 centiMorgan and used to map qualitative and quantitative traits. This collection of G. hirsutum varietal SNPs complements existing intra-specific SNPs and provides the cotton community with a valuable marker resource applicable to genetic analyses and breeding programs.

Hieracium praealtum forms seeds asexually by apomixis. During ovule development, sexual reproduction initiates with megaspore mother cell entry into meiosis and formation of a tetrad of haploid megaspores. The sexual pathway ceases when a diploid aposporous initial (AI) cell differentiates, enlarges, and undergoes mitosis, forming an aposporous embryo sac that displaces sexual structures. Embryo and endosperm development in aposporous embryo sacs is fertilization independent. Transcriptional data relating to apomixis initiation in Hieracium spp. ovules is scarce and the functional identity of the AI cell relative to other ovule cell types is unclear. Enlarging AI cells with undivided nuclei, early aposporous embryo sacs containing two to four nuclei, and random groups of sporophytic ovule cells not undergoing these events were collected by laser capture microdissection. Isolated amplified messenger RNA samples were sequenced using the 454 pyrosequencing platform and comparatively analyzed to establish indicative roles of the captured cell types. Transcriptome and protein motif analyses showed that approximately one-half of the assembled contigs identified homologous sequences in Arabidopsis (Arabidopsis thaliana), of which the vast majority were expressed during early Arabidopsis ovule development. The sporophytic ovule cells were enriched in signaling functions. Gene expression indicative of meiosis was notably absent in enlarging AI cells, consistent with subsequent aposporous embryo sac formation without meiosis. The AI cell transcriptome was most similar to the early aposporous embryo sac transcriptome when comparing known functional annotations and both shared expressed genes involved in gametophyte development, suggesting that the enlarging AI cell is already transitioning to an embryo sac program prior to mitotic division.

Background Application of apomixis, or asexual seed formation, in crop breeding would allow rapid fixation of complex traits, economizing improved crop delivery. Identification of apomixis genes is confounded by the polyploid nature, high genome complexity and lack of genomic sequence integration with reproductive tissue transcriptomes in most apomicts. Results A genomic and transcriptomic resource was developed for Hieracium subgenus Pilosella (Asteraceae) which incorporates characterized sexual, apomictic and mutant apomict plants exhibiting reversion to sexual reproduction. Apomicts develop additional female gametogenic cells that suppress the sexual pathway in ovules. Disrupting small RNA pathways in sexual Arabidopsis also induces extra female gametogenic cells; therefore, the resource was used to examine if changes in small RNA pathways correlate with apomixis initiation. An initial characterization of small RNA pathway genes within Hieracium was undertaken, and ovary-expressed ARGONAUTE genes were identified and cloned. Comparisons of whole ovary transcriptomes from mutant apomicts, relative to the parental apomict, revealed that differentially expressed genes were enriched for processes involved in small RNA biogenesis and chromatin silencing. Small RNA profiles within mutant ovaries did not reveal large-scale alterations in composition or length distributions; however, a small number of differentially expressed, putative small RNA targets were identified. Conclusions The established Hieracium resource represents a substantial contribution towards the investigation of early sexual and apomictic female gamete development, and the generation of new candidate genes and markers. Observed changes in small RNA targets and biogenesis pathways within sexual and apomictic ovaries will underlie future functional research into apomixis initiation in Hieracium. Electronic supplementary material The online version of this article (doi:10.1186/s12915-016-0311-0) contains supplementary material, which is available to authorized users.

Motivation: We have recently characterized an instance of alternative splicing that differs from the canonical gene transcript by deletion of a length of sequence not divisible by three, but where translation can be rescued by an alternative start codon. This results in a predicted protein in which the amino terminus differs markedly in sequence from the known protein product(s), as it is translated from an alternative reading frame. Automated pipelines have annotated thousands of splice variants but have overlooked these protein isoforms, leading to them being underrepresented in current databases. Results: Here we describe 1849 human and 733 mouse transcripts that can be transcribed from an alternate ATG. Of these, >80% have not been annotated previously. Those conserved between human and mouse genomes (and hence under likely evolutionary selection) are identified. We provide mass spectroscopy evidence for translation of selected transcripts. Of the described splice variants, only one has previously been studied in detail and converted the encoded protein from an activator of cell-function to a suppressor, demonstrating that these splice variants can result in profound functional change. We investigate the potential functional effects of this splicing using a variety of bioinformatic tools. The 2582 variants we describe are involved in a wide variety of biological processes, and therefore open many new avenues of research. Contact: aude.fahrer@anu.edu.au Supplementary Inforation: Supplementary data are available at Bioinformatics online.

Background In legumes, seed storage proteins are important for the developing seedling and are an important source of protein for humans and animals. Lupinus angustifolius (L.), also known as narrow-leaf lupin (NLL) is a grain legume crop that is gaining recognition as a potential human health food as the grain is high in protein and dietary fibre, gluten-free and low in fat and starch. Results Genes encoding the seed storage proteins of NLL were characterised by sequencing cDNA clones derived from developing seeds. Four families of seed storage proteins were identified and comprised three unique α, seven β, two γ and four δ conglutins. This study added eleven new expressed storage protein genes for the species. A comparison of the deduced amino acid sequences of NLL conglutins with those available for the storage proteins of Lupinus albus (L.), Pisum sativum (L.), Medicago truncatula (L.), Arachis hypogaea (L.) and Glycine max (L.) permitted the analysis of a phylogenetic relationships between proteins and demonstrated, in general, that the strongest conservation occurred within species. In the case of 7S globulin (β conglutins) and 2S sulphur-rich albumin (δ conglutins), the analysis suggests that gene duplication occurred after legume speciation. This contrasted with 11S globulin (α conglutin) and basic 7S (γ conglutin) sequences where some of these sequences appear to have diverged prior to speciation. The most abundant NLL conglutin family was β (56%), followed by α (24%), δ (15%) and γ (6%) and the transcript levels of these genes increased 103 to 106 fold during seed development. We used the 16 NLL conglutin sequences identified here to determine that for individuals specifically allergic to lupin, all seven members of the β conglutin family were potential allergens. Conclusion This study has characterised 16 seed storage protein genes in NLL including 11 newly-identified members. It has helped lay the foundation for efforts to use molecular breeding approaches to improve lupins, for example by reducing allergens or increasing the expression of specific seed storage protein(s) with desirable nutritional properties.

Genomic imprinting causes the expression of an allele depending on its parental origin. In plants, most imprinted genes have been identified in Arabidopsis endosperm, a transient structure consumed by the embryo during seed formation. We identified imprinted genes in rice seed where both the endosperm and embryo are present at seed maturity. RNA was extracted from embryos and endosperm of seeds obtained from reciprocal crosses between two subspecies Nipponbare (Japonica rice) and 93-11 (Indica rice). Sequenced reads from cDNA libraries were aligned to their respective parental genomes using single-nucleotide polymorphisms (SNPs). Reads across SNPs enabled derivation of parental expression bias ratios. A continuum of parental expression bias states was observed. Statistical analyses indicated 262 candidate imprinted loci in the endosperm and three in the embryo (168 genic and 97 non-genic). Fifty-six of the 67 loci investigated were confirmed to be imprinted in the seed. Imprinted loci are not clustered in the rice genome as found in mammals. All of these imprinted loci were expressed in the endosperm, and one of these was also imprinted in the embryo, confirming that in both rice and Arabidopsis imprinted expression is primarily confined to the endosperm. Some rice imprinted genes were also expressed in vegetative tissues, indicating that they have additional roles in plant growth. Comparison of candidate imprinted genes found in rice with imprinted candidate loci obtained from genome-wide surveys of imprinted genes in Arabidopsis to date shows a low degree of conservation, suggesting that imprinting has evolved independently in eudicots and monocots., Author Summary The expression of maternal or paternal alleles in either a preferentially or exclusively uniparental manner, termed imprinting, is prevalent in the transient endosperm of seeds in the model plant Arabidopsis. Cereals form seeds where both the embryo and endosperm are present at seed maturity. They are an important world food source. To date, very few imprinted genes have been identified in cereal seeds. How parental gene expression biases contribute to rice seed development has not yet been studied in detail. The deep resolution of transcript sequencing platforms was used to identify loci expressed in a parentally biased manner in the embryo and endosperm of Indica and Japonica rice at a genome-wide level. We identified 262 candidate imprinted loci expressed in the endosperm, experimentally verified 56 of these, and found novel features pertaining to their expression. Only one gene was found to be imprinted in the rice embryo. Imprinting in Arabidopsis and rice seeds is confined primarily to the endosperm, but the identified loci do not share extensive sequence conservation. Imprinting thus appears to have evolved independently in these plant species.

Low-oxygen (hypoxia) stress associated with natural phenomena such as waterlogging, results in widespread transcriptome changes and a metabolic switch from aerobic respiration to anaerobic fermentation. High-throughput sequencing of small RNA libraries obtained from hypoxia-treated and control root tissue identified a total of 65 unique microRNA (miRNA) sequences from 46 families, and 14 trans-acting small interfering RNA (tasiRNA) from three families. Hypoxia resulted in changes to the abundance of 46 miRNAs from 19 families, and all three tasiRNA families. Chemical inhibition of mitochondrial respiration caused similar changes in expression in a majority of the hypoxia-responsive small RNAs analysed. Our data indicate that miRNAs and tasiRNAs play a role in gene regulation and possibly developmental responses to hypoxia, and that a major signal for these responses is likely to be dependent on mitochondrial function.

Endogenous small RNAs, including microRNAs (miRNAs) and short-interfering RNAs (siRNAs), function as post-transcriptional or transcriptional regulators in plants. miRNA function is essential for normal plant development and therefore is likely to be important in the growth of the rice grain. To investigate the roles of miRNAs in rice grain development, we carried out deep sequencing of the small RNA populations of rice grains at two developmental stages. In a data set of ∼5.5 million sequences, we found representatives of all 20 conserved plant miRNA families. We used an approach based on the presence of miRNA and miRNA* sequences to identify 39 novel, nonconserved rice miRNA families expressed in grains. Cleavage of predicted target mRNAs was confirmed for a number of the new miRNAs. We identified a putative mirtron, indicating that plants may also use spliced introns as a source of miRNAs. We also identified a miRNA-like long hairpin that generates phased 21 nt small RNAs, strongly expressed in developing grains, and show that these small RNAs act in trans to cleave target mRNAs. Comparison of the population of miRNAs and miRNA-like siRNAs in grains to those in other parts of the rice plant reveals that many are expressed in an organ-specific manner.

Low-oxygen (hypoxia) stress associated with natural phenomena such as waterlogging, results in widespread transcriptome changes and a metabolic switch from aerobic respiration to anaerobic fermentation. High-throughput sequencing of small RNA libraries obtained from hypoxia-treated and control root tissue identified a total of 65 unique microRNA (miRNA) sequences from 46 families, and 14 trans-acting small interfering RNA (tasiRNA) from three families. Hypoxia resulted in changes to the abundance of 46 miRNAs from 19 families, and all three tasiRNA families. Chemical inhibition of mitochondrial respiration caused similar changes in expression in a majority of the hypoxia-responsive small RNAs analysed. Our data indicate that miRNAs and tasiRNAs play a role in gene regulation and possibly developmental responses to hypoxia, and that a major signal for these responses is likely to be dependent on mitochondrial function.

Background Cotton germplasm resources contain beneficial alleles that can be exploited to develop germplasm adapted to emerging environmental and climate conditions. Accessions and lines have traditionally been characterized based on phenotypes, but phenotypic profiles are limited by the cost, time, and space required to make visual observations and measurements. With advances in molecular genetic methods, genotypic profiles are increasingly able to identify differences among accessions due to the larger number of genetic markers that can be measured. A combination of both methods would greatly enhance our ability to characterize germplasm resources. Recent efforts have culminated in the identification of sufficient SNP markers to establish high-throughput genotyping systems, such as the CottonSNP63K array, which enables a researcher to efficiently analyze large numbers of SNP markers and obtain highly repeatable results. In the current investigation, we have utilized the SNP array for analyzing genetic diversity primarily among cotton cultivars, making comparisons to SSR-based phylogenetic analyses, and identifying loci associated with seed nutritional traits. Results The SNP markers distinctly separated G. hirsutum from other Gossypium species and distinguished the wild from cultivated types of G. hirsutum. The markers also efficiently discerned differences among cultivars, which was the primary goal when designing the CottonSNP63K array. Population structure within the genus compared favorably with previous results obtained using SSR markers, and an association study identified loci linked to factors that affect cottonseed protein content. Conclusions Our results provide a large genome-wide variation data set for primarily cultivated cotton. Thousands of SNPs in representative cotton genotypes provide an opportunity to finely discriminate among cultivated cotton from around the world. The SNPs will be relevant as dense markers of genome variation for association mapping approaches aimed at correlating molecular polymorphisms with variation in phenotypic traits, as well as for molecular breeding approaches in cotton. Electronic supplementary material The online version of this article (doi:10.1186/s12870-017-0981-y) contains supplementary material, which is available to authorized users.

Background During the early stages of seed development many genes are under dynamic regulation to ensure the proper differentiation and establishment of the tissue that will constitute the mature grain. To investigate how miRNA regulation contributes to this process in barley, a combination of small RNA and mRNA degradome analyses were used to identify miRNAs and their targets. Results Our analysis identified 84 known miRNAs and 7 new miRNAs together with 96 putative miRNA target genes regulated through a slicing mechanism in grain tissues during the first 15 days post anthesis. We also identified many potential miRNAs including several belonging to known miRNA families. Our data gave us evidence for an increase in miRNA-mediated regulation during the transition between pre-storage and storage phases. Potential miRNA targets were found in various signalling pathways including components of four phytohormone pathways (ABA, GA, auxin, ethylene) and the defence response to powdery mildew infection. Among the putative miRNA targets we identified were two essential genes controlling the GA response, a GA3oxidase1 and a homolog of the receptor GID1, and a homolog of the ACC oxidase which catalyses the last step of ethylene biosynthesis. We found that two MLA genes are potentially miRNA regulated, establishing a direct link between miRNAs and the R gene response. Conclusion Our dataset provides a useful source of information on miRNA regulation during the early development of cereal grains and our analysis suggests that miRNAs contribute to the control of development of the cereal grain, notably through the regulation of phytohormone response pathways.