Hiding in plain sight: Genome-wide recombination and a dynamic accessory genome drive diversity in Fusarium oxysporum f.sp. ciceris.

Understanding the origins of variation in agricultural pathogens is of fundamental interest and practical importance, especially for diseases that threaten food security. Fusarium oxysporum is among the most important of soil-borne pathogens, with a global distribution and an extensive host range. The pathogen is considered to be asexual, with horizontal transfer of chromosomes providing an analog of assortment by meiotic recombination. Here, we challenge those assumptions based on the results of population genomic analyses, describing the pathogen's diversity and inferring its origins and functional consequences in the context of a single, long-standing agricultural system. We identify simultaneously low nucleotide distance among strains, and unexpectedly high levels of genetic and genomic variability. We determine that these features arise from a combination of genome-scale recombination, best explained by widespread sexual reproduction, and presence-absence variation consistent with chromosomal rearrangement. Pangenome analyses document an accessory genome more than twice the size of the core genome, with contrasting evolutionary dynamics. The core genome is stable, with low diversity and high genetic differentiation across geographic space, while the accessory genome is paradoxically more diverse and unstable but with lower genetic differentiation and hallmarks of contemporary gene flow at local scales. We suggest a model in which episodic sexual reproduction generates haplotypes that are selected and then maintained through clone-like dynamics, followed by contemporary genomic rearrangements that reassort the accessory genome among sympatric strains. Taken together, these processes contribute unique genome content, including reassortment of virulence determinants that may explain observed variation in pathogenic potential.

Genome Context Viewer (GCV) version 2: enhanced visual exploration of multiple annotated genomes.

The Genome Context Viewer is a web application for identifying, aligning, and visualizing genomic regions based on their micro and macrosyntenic structures. By using functional elements such as gene annotations as the unit of search and comparison, the Genome Context Viewer can compute and display relationships between regions across many assemblies from federated data sources in real-time, enabling users to rapidly explore multiple annotated genomes and identify divergence and structural events that can help provide insight into evolutionary mechanisms associated with functional consequences. In this work, we introduce version 2 of the Genome Context Viewer and highlight new features that enhance usability, performance, and ease of deployment.

Genome assembly of Medicago truncatula accession SA27063 provides insight into spring black stem and leaf spot disease resistance.

Medicago truncatula, model legume and alfalfa relative, has served as an essential resource for advancing our understanding of legume physiology, functional genetics, and crop improvement traits. Necrotrophic fungus, Ascochyta medicaginicola, the causal agent of spring black stem (SBS) and leaf spot is a devasting foliar disease of alfalfa affecting stand survival, yield, and forage quality. Host resistance to SBS disease is poorly understood, and control methods rely on cultural practices. Resistance has been observed in M. truncatula accession SA27063 (HM078) with two recessively inherited quantitative-trait loci (QTL), rnpm1 and rnpm2, previously reported. To shed light on host resistance, we carried out a de novo genome assembly of HM078. The genome, referred to as MtHM078 v1.0, is comprised of 23 contigs totaling 481.19 Mbp. Notably, this assembly contains a substantial amount of novel centromere-related repeat sequences due to deep long-read sequencing. Genome annotation resulted in 98.4% of BUSCO fabales proteins being complete. The assembly enabled sequence-level analysis of rnpm1 and rnpm2 for gene content, synteny, and structural variation between SBS-resistant accession SA27063 (HM078) and SBS-susceptible accession A17 (HM101). Fourteen candidate genes were identified, and some have been implicated in resistance to necrotrophic fungi. Especially interesting candidates include loss-of-function events in HM078 because they fit the inverse gene-for-gene model, where resistance is recessively inherited. In rnpm1, these include a loss-of-function in a disease resistance gene due to a premature stop codon, and a 10.85 kbp retrotransposon-like insertion disrupting a ubiquitin conjugating E2. In rnpm2, we identified a frameshift mutation causing a loss-of-function in a glycosidase, as well as a missense and frameshift mutation altering an F-box family protein. This study generated a high-quality genome of HM078 and has identified promising candidates, that once validated, could be further studied in alfalfa to enhance disease resistance.

Legacy genetics of Arachis cardenasii in the peanut crop shows the profound benefits of international seed exchange.

The narrow genetics of most crops is a fundamental vulnerability to food security. This makes wild crop relatives a strategic resource of genetic diversity that can be used for crop improvement and adaptation to new agricultural challenges. Here, we uncover the contribution of one wild species accession, Arachis cardenasii GKP 10017, to the peanut crop (Arachis hypogaea) that was initiated by complex hybridizations in the 1960s and propagated by international seed exchange. However, until this study, the global scale of the dispersal of genetic contributions from this wild accession had been obscured by the multiple germplasm transfers, breeding cycles, and unrecorded genetic mixing between lineages that had occurred over the years. By genetic analysis and pedigree research, we identified A. cardenasii-enhanced, disease-resistant cultivars in Africa, Asia, Oceania, and the Americas. These cultivars provide widespread improved food security and environmental and economic benefits. This study emphasizes the importance of wild species and collaborative networks of international expertise for crop improvement. However, it also highlights the consequences of the implementation of a patchwork of restrictive national laws and sea changes in attitudes regarding germplasm that followed in the wake of the Convention on Biological Diversity. Today, the botanical collections and multiple seed exchanges which enable benefits such as those revealed by this study are drastically reduced. The research reported here underscores the vital importance of ready access to germplasm in ensuring long-term world food security.

The genome of a wild Medicago species provides insights into the tolerant mechanisms of legume forage to environmental stress.

BACKGROUND: Medicago ruthenica, a wild and perennial legume forage widely distributed in semi-arid grasslands, is distinguished by its outstanding tolerance to environmental stress. It is a close relative of commonly cultivated forage of alfalfa (Medicago sativa). The high tolerance of M. ruthenica to environmental stress makes this species a valuable genetic resource for understanding and improving traits associated with tolerance to harsh environments. RESULTS: We sequenced and assembled genome of M. ruthenica using an integrated approach, including PacBio, Illumina, 10×Genomics, and Hi-C. The assembled genome was 904.13 Mb with scaffold N50 of 99.39 Mb, and 50,162 protein-coding genes were annotated. Comparative genomics and transcriptomic analyses were used to elucidate mechanisms underlying its tolerance to environmental stress. The expanded FHY3/FAR1 family was identified to be involved in tolerance of M. ruthenica to drought stress. Many genes involved in tolerance to abiotic stress were retained in M. ruthenica compared to other cultivated Medicago species. Hundreds of candidate genes associated with drought tolerance were identified by analyzing variations in single nucleotide polymorphism using accessions of M. ruthenica with varying tolerance to drought. Transcriptomic data demonstrated the involvements of genes related to transcriptional regulation, stress response, and metabolic regulation in tolerance of M. ruthenica. CONCLUSIONS: We present a high-quality genome assembly and identification of drought-related genes in the wild species of M. ruthenica, providing a valuable resource for genomic studies on perennial legume forages.

The genome of cowpea (Vigna unguiculata [L.] Walp.).

Cowpea (Vigna unguiculata [L.] Walp.) is a major crop for worldwide food and nutritional security, especially in sub-Saharan Africa, that is resilient to hot and drought-prone environments. An assembly of the single-haplotype inbred genome of cowpea IT97K-499-35 was developed by exploiting the synergies between single-molecule real-time sequencing, optical and genetic mapping, and an assembly reconciliation algorithm. A total of 519 Mb is included in the assembled sequences. Nearly half of the assembled sequence is composed of repetitive elements, which are enriched within recombination-poor pericentromeric regions. A comparative analysis of these elements suggests that genome size differences between Vigna species are mainly attributable to changes in the amount of Gypsy retrotransposons. Conversely, genes are more abundant in more distal, high-recombination regions of the chromosomes; there appears to be more duplication of genes within the NBS-LRR and the SAUR-like auxin superfamilies compared with other warm-season legumes that have been sequenced. A surprising outcome is the identification of an inversion of 4.2 Mb among landraces and cultivars, which includes a gene that has been associated in other plants with interactions with the parasitic weed Striga gesnerioides. The genome sequence facilitated the identification of a putative syntelog for multiple organ gigantism in legumes. A revised numbering system has been adopted for cowpea chromosomes based on synteny with common bean (Phaseolus vulgaris). An estimate of nuclear genome size of 640.6 Mbp based on cytometry is presented.

A genome-wide association and meta-analysis reveal regions associated with seed size in cowpea [Vigna unguiculata (L.) Walp].

KEY MESSAGE: This paper combined GWAS, meta-analysis and sequence homology comparison with common bean to identify regions associated with seed size variation in domesticated cowpea. Seed size is an important trait for yield and commercial value in dry-grain cowpea. Seed size varies widely among different cowpea accessions, and the genetic basis of such variation is not yet well understood. To better decipher the genetic basis of seed size, a genome-wide association study (GWAS) and meta-analysis were conducted on a panel of 368 cowpea diverse accessions from 51 countries. Four traits, including seed weight, length, width and density were evaluated across three locations. Using 51,128 single nucleotide polymorphisms covering the cowpea genome, 17 loci were identified for these traits. One locus was common to weight, width and length, suggesting pleiotropy. By integrating synteny-based analysis with common bean, six candidate genes (Vigun05g036000, Vigun05g039600, Vigun05g204200, Vigun08g217000, Vigun11g187000, and Vigun11g191300) which are implicated in multiple functional categories related to seed size such as endosperm development, embryo development, and cell elongation were identified. These results suggest that a combination of GWAS meta-analysis with synteny comparison in a related plant is an efficient approach to identify candidate gene (s) for complex traits in cowpea. The identified loci and candidate genes provide useful information for improving cowpea varieties and for molecular investigation of seed size.

Combating a Global Threat to a Clonal Crop: Banana Black Sigatoka Pathogen Pseudocercospora fijiensis (Synonym Mycosphaerella fijiensis) Genomes Reveal Clues for Disease Control.

Black Sigatoka or black leaf streak disease, caused by the Dothideomycete fungus Pseudocercospora fijiensis (previously: Mycosphaerella fijiensis), is the most significant foliar disease of banana worldwide. Due to the lack of effective host resistance, management of this disease requires frequent fungicide applications, which greatly increase the economic and environmental costs to produce banana. Weekly applications in most banana plantations lead to rapid evolution of fungicide-resistant strains within populations causing disease-control failures throughout the world. Given its extremely high economic importance, two strains of P. fijiensis were sequenced and assembled with the aid of a new genetic linkage map. The 74-Mb genome of P. fijiensis is massively expanded by LTR retrotransposons, making it the largest genome within the Dothideomycetes. Melting-curve assays suggest that the genomes of two closely related members of the Sigatoka disease complex, P. eumusae and P. musae, also are expanded. Electrophoretic karyotyping and analyses of molecular markers in P. fijiensis field populations showed chromosome-length polymorphisms and high genetic diversity. Genetic differentiation was also detected using neutral markers, suggesting strong selection with limited gene flow at the studied geographic scale. Frequencies of fungicide resistance in fungicide-treated plantations were much higher than those in untreated wild-type P. fijiensis populations. A homologue of the Cladosporium fulvum Avr4 effector, PfAvr4, was identified in the P. fijiensis genome. Infiltration of the purified PfAVR4 protein into leaves of the resistant banana variety Calcutta 4 resulted in a hypersensitive-like response. This result suggests that Calcutta 4 could carry an unknown resistance gene recognizing PfAVR4. Besides adding to our understanding of the overall Dothideomycete genome structures, the P. fijiensis genome will aid in developing fungicide treatment schedules to combat this pathogen and in improving the efficiency of banana breeding programs.

Correction to: The transcriptome landscape of early maize meiosis.

CORRECTION: Following publication of the original article [1], the authors reported that the number of genes overlaying the bar graph in Fig. 3A were incorrectly counted and inserted (i.e. including a title tile, and in inverse order). The corrected numbers are below and match with the listed genes supplied in Additional File: Table S2.

A transcriptome-SNP-derived linkage map of Apios americana (potato bean) provides insights about genome re-organization and synteny conservation in the phaseoloid legumes.

KEY MESSAGE: We report a linkage map for Apios americana and describe synteny with selected warm-season legumes. A translocation event in common bean and soybean is confirmed against Apios and Vigna species. Apios (Apios americana; "apios"), a tuberous perennial legume in the Phaseoleae tribe, was widely used as a food by Native Americans. Work in the last 40 years has led to several improved breeding lines. Aspects of the pollination biology (complex floral structure and tripping mechanism) have made controlled crosses difficult, and the previous reports indicated that the plant is likely primarily an outcrosser. We used a pseudo-testcross strategy to construct a genetic map specific to the maternal parent. The map was built using single-nucleotide polymorphism markers identified by comparing the expressed sequences of individuals in the mapping population against a de novo maternal reference transcriptome assembly. The apios map consists of 11 linkage groups and 1121 recombinationally distinct loci, covering ~ 938.6 cM. By sequencing the transcriptomes of all potential pollen parents, we were able to identify the probable pollen donors and to discover new aspects of the pollination biology in apios. No selfing was observed, but multiple pollen parents were seen within individual pods. Comparisons with genome sequences in other species in the Phaseoleae showed extended synteny for most apios linkage groups. This synteny supports the robustness of the map, and also sheds light on the history of the Phaseoleae, as apios is relatively early diverging in this tribe. We detected a translocation event that separates apios and two Vigna species from Phaseolus vulgaris and Glycine max. This apios mapping work provides a general protocol for sequencing-based construction of high-density linkage maps in outcrossing species with heterogeneous pollen parents.

Legume information system (LegumeInfo.org): a key component of a set of federated data resources for the legume family.

Legume Information System (LIS), at http://legumeinfo.org, is a genomic data portal (GDP) for the legume family. LIS provides access to genetic and genomic information for major crop and model legumes. With more than two-dozen domesticated legume species, there are numerous specialists working on particular species, and also numerous GDPs for these species. LIS has been redesigned in the last three years both to better integrate data sets across the crop and model legumes, and to better accommodate specialized GDPs that serve particular legume species. To integrate data sets, LIS provides genome and map viewers, holds synteny mappings among all sequenced legume species and provides a set of gene families to allow traversal among orthologous and paralogous sequences across the legumes. To better accommodate other specialized GDPs, LIS uses open-source GMOD components where possible, and advocates use of common data templates, formats, schemas and interfaces so that data collected by one legume research community are accessible across all legume GDPs, through similar interfaces and using common APIs. This federated model for the legumes is managed as part of the 'Legume Federation' project (accessible via http://legumefederation.org), which can be thought of as an umbrella project encompassing LIS and other legume GDPs.

Exploring structural variation and gene family architecture with De Novo assemblies of 15 Medicago genomes.

BACKGROUND: Previous studies exploring sequence variation in the model legume, Medicago truncatula, relied on mapping short reads to a single reference. However, read-mapping approaches are inadequate to examine large, diverse gene families or to probe variation in repeat-rich or highly divergent genome regions. De novo sequencing and assembly of M. truncatula genomes enables near-comprehensive discovery of structural variants (SVs), analysis of rapidly evolving gene families, and ultimately, construction of a pan-genome. RESULTS: Genome-wide synteny based on 15 de novo M. truncatula assemblies effectively detected different types of SVs indicating that as much as 22% of the genome is involved in large structural changes, altogether affecting 28% of gene models. A total of 63 million base pairs (Mbp) of novel sequence was discovered, expanding the reference genome space for Medicago by 16%. Pan-genome analysis revealed that 42% (180 Mbp) of genomic sequences is missing in one or more accession, while examination of de novo annotated genes identified 67% (50,700) of all ortholog groups as dispensable - estimates comparable to recent studies in rice, maize and soybean. Rapidly evolving gene families typically associated with biotic interactions and stress response were found to be enriched in the accession-specific gene pool. The nucleotide-binding site leucine-rich repeat (NBS-LRR) family, in particular, harbors the highest level of nucleotide diversity, large effect single nucleotide change, protein diversity, and presence/absence variation. However, the leucine-rich repeat (LRR) and heat shock gene families are disproportionately affected by large effect single nucleotide changes and even higher levels of copy number variation. CONCLUSIONS: Analysis of multiple M. truncatula genomes illustrates the value of de novo assemblies to discover and describe structural variation, something that is often under-estimated when using read-mapping approaches. Comparisons among the de novo assemblies also indicate that different large gene families differ in the architecture of their structural variation.

The Medicago genome provides insight into the evolution of rhizobial symbioses.

Legumes (Fabaceae or Leguminosae) are unique among cultivated plants for their ability to carry out endosymbiotic nitrogen fixation with rhizobial bacteria, a process that takes place in a specialized structure known as the nodule. Legumes belong to one of the two main groups of eurosids, the Fabidae, which includes most species capable of endosymbiotic nitrogen fixation. Legumes comprise several evolutionary lineages derived from a common ancestor 60 million years ago (Myr ago). Papilionoids are the largest clade, dating nearly to the origin of legumes and containing most cultivated species. Medicago truncatula is a long-established model for the study of legume biology. Here we describe the draft sequence of the M. truncatula euchromatin based on a recently completed BAC assembly supplemented with Illumina shotgun sequence, together capturing ∼94% of all M. truncatula genes. A whole-genome duplication (WGD) approximately 58 Myr ago had a major role in shaping the M. truncatula genome and thereby contributed to the evolution of endosymbiotic nitrogen fixation. Subsequent to the WGD, the M. truncatula genome experienced higher levels of rearrangement than two other sequenced legumes, Glycine max and Lotus japonicus. M. truncatula is a close relative of alfalfa (Medicago sativa), a widely cultivated crop with limited genomics tools and complex autotetraploid genetics. As such, the M. truncatula genome sequence provides significant opportunities to expand alfalfa's genomic toolbox.

The plant-specific protein FEHLSTART controls male meiotic entry, initializing meiotic synchronization in Arabidopsis.

Meiosis marks the transition from the sporophyte to the gametophyte generation in the life cycle of flowering plants, and creates genetic variations through homologous recombination. In most flowering plants, meiosis is highly synchronized within each anther, which is significant for efficient fertilization. To date, little is known about the molecular mechanisms of entry into meiosis and exit from it, and only a few genes in Arabidopsis have been characterized with a role in regulating meiotic progression. In this study, we report the functional characterization of a plant-specific basic helix-loop-helix (bHLH) protein, FEHLSTART (FST), a defect in which leads to premature meiotic entry and asynchronous meiosis, and results in decreased seed yield. Investigation of the time course of meiosis showed that the onset of leptotene, the first stage of prophase I, frequently occurred earlier in fst-1 than in the wild type. Asynchronous meiosis followed, which could manifest in the disruption of regular spindle structures and symmetric cell divisions in fst-1 mutants during the meiosis I/II transition. In accordance with frequently accelerated meiotic entry, whole-transcriptome analysis of fst-1 anthers undergoing meiosis revealed that 19 circadian rhythm genes were affected and 47 pollen-related genes were prematurely expressed at a higher level. Taken together, we propose that FST is required for normal meiotic entry and the establishment of meiotic synchrony.

Multiple post-domestication origins of kabuli chickpea through allelic variation in a diversification-associated transcription factor.

Chickpea (Cicer arietinum) is among the founder crops domesticated in the Fertile Crescent. One of two major forms of chickpea, the so-called kabuli type, has white flowers and light-colored seed coats, properties not known to exist in the wild progenitor. The origin of the kabuli form has been enigmatic. We genotyped a collection of wild and cultivated chickpea genotypes with 538 single nucleotide polymorphisms (SNPs) and examined patterns of molecular diversity relative to geographical sources and market types. In addition, we examined sequence and expression variation in candidate anthocyanin biosynthetic pathway genes. A reduction in genetic diversity and extensive genetic admixture distinguish cultivated chickpea from its wild progenitor species. Among germplasm, the kabuli form is polyphyletic. We identified a basic helix-loop-helix (bHLH) transcription factor at chickpea's B locus that conditions flower and seed colors, orthologous to Mendel's A gene of garden pea, whose loss of function is associated invariantly with the kabuli type of chickpea. From the polyphyletic distribution of the kabuli form in germplasm, an absence of nested variation within the bHLH gene and invariant association of loss of function of bHLH among the kabuli type, we conclude that the kabuli form arose multiple times during the phase of phenotypic diversification after initial domestication of cultivated chickpea.

Genome, epigenome and RNA sequences of monozygotic twins discordant for multiple sclerosis.

Monozygotic or 'identical' twins have been widely studied to dissect the relative contributions of genetics and environment in human diseases. In multiple sclerosis (MS), an autoimmune demyelinating disease and common cause of neurodegeneration and disability in young adults, disease discordance in monozygotic twins has been interpreted to indicate environmental importance in its pathogenesis. However, genetic and epigenetic differences between monozygotic twins have been described, challenging the accepted experimental model in disambiguating the effects of nature and nurture. Here we report the genome sequences of one MS-discordant monozygotic twin pair, and messenger RNA transcriptome and epigenome sequences of CD4(+) lymphocytes from three MS-discordant, monozygotic twin pairs. No reproducible differences were detected between co-twins among approximately 3.6 million single nucleotide polymorphisms (SNPs) or approximately 0.2 million insertion-deletion polymorphisms. Nor were any reproducible differences observed between siblings of the three twin pairs in HLA haplotypes, confirmed MS-susceptibility SNPs, copy number variations, mRNA and genomic SNP and insertion-deletion genotypes, or the expression of approximately 19,000 genes in CD4(+) T cells. Only 2 to 176 differences in the methylation of approximately 2 million CpG dinucleotides were detected between siblings of the three twin pairs, in contrast to approximately 800 methylation differences between T cells of unrelated individuals and several thousand differences between tissues or between normal and cancerous tissues. In the first systematic effort to estimate sequence variation among monozygotic co-twins, we did not find evidence for genetic, epigenetic or transcriptome differences that explained disease discordance. These are the first, to our knowledge, female, twin and autoimmune disease individual genome sequences reported.

Population genomics of the facultatively mutualistic bacteria Sinorhizobium meliloti and S. medicae.

The symbiosis between rhizobial bacteria and legume plants has served as a model for investigating the genetics of nitrogen fixation and the evolution of facultative mutualism. We used deep sequence coverage (>100×) to characterize genomic diversity at the nucleotide level among 12 Sinorhizobium medicae and 32 S. meliloti strains. Although these species are closely related and share host plants, based on the ratio of shared polymorphisms to fixed differences we found that horizontal gene transfer (HGT) between these species was confined almost exclusively to plasmid genes. Three multi-genic regions that show the strongest evidence of HGT harbor genes directly involved in establishing or maintaining the mutualism with host plants. In both species, nucleotide diversity is 1.5-2.5 times greater on the plasmids than chromosomes. Interestingly, nucleotide diversity in S. meliloti but not S. medicae is highly structured along the chromosome - with mean diversity (θ(π)) on one half of the chromosome five times greater than mean diversity on the other half. Based on the ratio of plasmid to chromosome diversity, this appears to be due to severely reduced diversity on the chromosome half with less diversity, which is consistent with extensive hitchhiking along with a selective sweep. Frequency-spectrum based tests identified 82 genes with a signature of adaptive evolution in one species or another but none of the genes were identified in both species. Based upon available functional information, several genes identified as targets of selection are likely to alter the symbiosis with the host plant, making them attractive targets for further functional characterization.

Comprehensive characterization and RNA-Seq profiling of the HD-Zip transcription factor family in soybean (Glycine max) during dehydration and salt stress.

BACKGROUND: The homeodomain leucine zipper (HD-Zip) transcription factor family is one of the largest plant specific superfamilies, and includes genes with roles in modulation of plant growth and response to environmental stresses. Many HD-Zip genes are characterized in Arabidopsis (Arabidopsis thaliana), and members of the family are being investigated for abiotic stress responses in rice (Oryza sativa), maize (Zea mays), poplar (Populus trichocarpa) and cucumber (Cucmis sativus). Findings in these species suggest HD-Zip genes as high priority candidates for crop improvement. RESULTS: In this study we have identified members of the HD-Zip gene family in soybean cv. 'Williams 82', and characterized their expression under dehydration and salt stress. Homology searches with BLASTP and Hidden Markov Model guided sequence alignments identified 101 HD-Zip genes in the soybean genome. Phylogeny reconstruction coupled with domain and gene structure analyses using soybean, Arabidopsis, rice, grape (Vitis vinifera), and Medicago truncatula homologues enabled placement of these sequences into four previously described subfamilies. Of the 101 HD-Zip genes identified in soybean, 88 exist as whole-genome duplication-derived gene pairs, indicating high retention of these genes following polyploidy in Glycine ~13 Mya. The HD-Zip genes exhibit ubiquitous expression patterns across 24 conditions that include 17 tissues of soybean. An RNA-Seq experiment performed to study differential gene expression at 0, 1, 6 and 12 hr soybean roots under dehydration and salt stress identified 20 differentially expressed (DE) genes. Several of these DE genes are orthologs of genes previously reported to play a role under abiotic stress, implying conservation of HD-Zip gene functions across species. Screening of HD-Zip promoters identified transcription factor binding sites that are overrepresented in the DE genes under both dehydration and salt stress, providing further support for the role of HD-Zip genes in abiotic stress responses. CONCLUSIONS: We provide a thorough description of soybean HD-Zip genes, and identify potential candidates with probable roles in dehydration and salt stress. Expression profiles generated for all soybean genes, under dehydration and salt stress, at four time points, will serve as an important resource for the soybean research community, and will aid in understanding plant responses to abiotic stress.

The transcriptome landscape of early maize meiosis.

BACKGROUND: A major step in the higher plant life cycle is the decision to leave the mitotic cell cycle and begin the progression through the meiotic cell cycle that leads to the formation of gametes. The molecular mechanisms that regulate this transition and early meiosis remain largely unknown. To gain insight into gene expression features during the initiation of meiotic recombination, we profiled early prophase I meiocytes from maize (Zea mays) using capillary collection to isolate meiocytes, followed by RNA-seq. RESULTS: We detected ~2,000 genes as preferentially expressed during early meiotic prophase, most of them uncharacterized. Functional analysis uncovered the importance of several cellular processes in early meiosis. Processes significantly enriched in isolated meiocytes included proteolysis, protein targeting, chromatin modification and the regulation of redox homeostasis. The most significantly up-regulated processes in meiocytes were processes involved in carbohydrate metabolism. Consistent with this, many mitochondrial genes were up-regulated in meiocytes, including nuclear- and mitochondrial-encoded genes. The data were validated with real-time PCR and in situ hybridization and also used to generate a candidate maize homologue list of known meiotic genes from Arabidopsis. CONCLUSIONS: Taken together, we present a high-resolution analysis of the transcriptome landscape in early meiosis of an important crop plant, providing support for choosing genes for detailed characterization of recombination initiation and regulation of early meiosis. Our data also reveal an important connection between meiotic processes and altered/increased energy production.