Comparative linkage maps of the rice and maize genomes

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

Pan-genome inversion index reveals evolutionary insights into the subpopulation structure of Asian rice

Understanding and exploiting genetic diversity is a key factor for the productive and stable production of rice. Here, we utilize 73 high-quality genomes that encompass the subpopulation structure of Asian rice (Oryza sativa), plus the genomes of two wild relatives (O. rufipogon and O. punctata), to build a pan-genome inversion index of 1769 non-redundant inversions that span an average of ~29% of the O. sativa cv. Nipponbare reference genome sequence. Using this index, we estimate an inversion rate of ~700 inversions per million years in Asian rice, which is 16 to 50 times higher than previously estimated for plants. Detailed analyses of these inversions show evidence of their effects on gene expression, recombination rate, and linkage disequilibrium. Our study uncovers the prevalence and scale of large inversions (≥100 bp) across the pan-genome of Asian rice and hints at their largely unexplored role in functional biology and crop performance.

Time series canopy phenotyping enables the identification of genetic variants controlling dynamic phenotypes in soybean

Advances in plant phenotyping technologies are dramatically reducing the marginal costs of collecting multiple phenotypic measurements across several time points. Yet, most current approaches and best statistical practices implemented to link genetic and phenotypic variation in plants have been developed in an era of single-time-point data. Here, we used time-series phenotypic data collected with an unmanned aircraft system for a large panel of soybean (Glycine max (L.) Merr.) varieties to identify previously uncharacterized loci. Specifically, we focused on the dissection of canopy coverage (CC) variation from this rich data set. We also inferred the speed of canopy closure, an additional dimension of CC, from the time-series data, as it may represent an important trait for weed control. Genome-wide association studies (GWASs) identified 35 loci exhibiting dynamic associations with CC across developmental stages. The time-series data enabled the identification of 10 known flowering time and plant height quantitative trait loci (QTLs) detected in previous studies of adult plants and the identification of novel QTLs influencing CC. These novel QTLs were disproportionately likely to act earlier in development, which may explain why they were missed in previous single-time-point studies. Moreover, this time-series data set contributed to the high accuracy of the GWASs, which we evaluated by permutation tests, as evidenced by the repeated identification of loci across multiple time points. Two novel loci showed evidence of adaptive selection during domestication, with different genotypes/haplotypes favored in different geographic regions. In summary, the time-series data, with soybean CC as an example, improved the accuracy and statistical power to dissect the genetic basis of traits and offered a promising opportunity for crop breeding with quantitative growth curves.

Variation in leaf transcriptome responses to elevated ozone corresponds with physiological sensitivity to ozone across maize inbred lines

We examine the impact of sustained elevated ozone concentration on the leaf transcriptome of 5 diverse maize inbred genotypes, which vary in physiological sensitivity to ozone (B73, Mo17, Hp301, C123, and NC338), using long reads to assemble transcripts and short reads to quantify expression of these transcripts. More than 99% of the long reads, 99% of the assembled transcripts, and 97% of the short reads map to both B73 and Mo17 reference genomes. Approximately 95% of the genes with assembled transcripts belong to known B73-Mo17 syntenic loci and 94% of genes with assembled transcripts are present in all temperate lines in the nested association mapping pan-genome. While there is limited evidence for alternative splicing in response to ozone stress, there is a difference in the magnitude of differential expression among the 5 genotypes. The transcriptional response to sustained ozone stress in the ozone resistant B73 genotype (151 genes) was modest, while more than 3,300 genes were significantly differentially expressed in the more sensitive NC338 genotype. There is the potential for tandem duplication in 30% of genes with assembled transcripts, but there is no obvious association between potential tandem duplication and differential expression. Genes with a common response across the 5 genotypes (83 genes) were associated with photosynthesis, in particular photosystem I. The functional annotation of genes not differentially expressed in B73 but responsive in the other 4 genotypes (789) identifies reactive oxygen species. This suggests that B73 has a different response to long-term ozone exposure than the other 4 genotypes. The relative magnitude of the genotypic response to ozone, and the enrichment analyses are consistent regardless of whether aligning short reads to: long read assembled transcripts; the B73 reference; the Mo17 reference. We find that prolonged ozone exposure directly impacts the photosynthetic machinery of the leaf.

Linkage mapping, comparative genome analysis, and QTL detection for growth in a non-model teleost, the meagre Argyrosomus regius, using ddRAD sequencing

Meagre (Argyrosomus regius), is a benthopelagic species rapidly emerging in aquaculture, due to its low food to biomass conversion rate, good fillet yield and ease of production. Tracing a species genomic background along with describing the genetic basis of important traits can greatly influence both conservation strategies and production perspectives. In this study, we employed ddRAD sequencing of 266 fish from six F1 meagre families, to construct a high-density genetic map comprising 4529 polymorphic SNP markers. The QTL mapping analysis provided a genomic appreciation for the weight trait identifying a statistically significant QTL on linkage group 15 (LG15). The comparative genomics analysis with six teleost species revealed an evolutionarily conserved karyotype structure. The synteny observed, verified the already well-known fusion events of the three-spine stickleback genome, reinforced the evidence of reduced evolutionary distance of Sciaenids with the Sparidae family, reflected the evolutionary proximity with Dicentrarchus labrax, traced several putative chromosomal rearrangements and a prominent putative fusion event in meagre’s LG17. This study presents novel elements concerning the genome evolutionary history of a non-model teleost species recently adopted in aquaculture, starts to unravel the genetic basis of the species growth-related traits, and provides a high-density genetic map as a tool that can help to further establish meagre as a valuable resource for research and production.

Molecular sexual determinants in Pistacia genus by KASP assay

BACKGROUND

Pistacia is a genus of dioecious plant species whose trees can take 4-5 years to reach the economically valuable fruit-bearing stage. The fruits have great importance as raw material in the food, healthcare, and baking industries. For that reason, the identification of individual plants in the early juvenile period for the pollination and positioning of trees is crucial for growers. The objective of this study is to develop markers for each Pistacia species that can help in screening the sex of plant seedlings before they reach the reproductive stage, without waiting for morphological characteristics to appear.

METHODS AND RESULTS

Within this context, by using the power of the kompetitive allele-specific PCR (KASP) assay technology as a marker screening system, we successfully discriminated seven out of eight Pistacia species: P. atlantica, P. integerrima, P. khinjuk, P. mutica, P. terebinthus, P. vera, and P. lentiscus. We used a high-throughput DNA sequence read archive (SRA) to assemble a reference genome that was employed in our studies as a de novo bioinformatics method. Four genomic regions from SRA and three single-nucleotide polymorphism (SNP) positions from Kafkas et al. BMC Genomics 16:98, 2015) were selected and sequenced with collected plant material from predominantly the Antepfıstıgı Research Institute Collection Garden, and eight species were aligned intraspecifically for SNP mining. In total, 12 SNP markers were converted to KASP markers, and 5 of them (SNP-PIS-133396, SNP-PIS-167992, P-ATL-91951-565, P-INT-91951-256, P-KHI-91951-115) showed clear allelic discrimination between male and female plants. SNP-PIS-167992 and P-ATL-91951-565 were identified as the best marker assays because they showed allelic frequency differences for all individuals and for both homozygous and heterozygous characters. These markers could be the most comprehensive ones for the whole genus because they showed discriminative power for several species.

CONCLUSIONS

This study is the first one to use the KASP assay for sex discrimination in Pistacia species, and it can be regarded as a precursor study for sex discrimination by KASP for plants in general.

On the Trail of the German Purity Law: Distinguishing the Metabolic Signatures of Wheat, Corn and Rice in Beer

Here, we report a non-targeted analytical approach to investigate the influence of different starch sources on the metabolic signature in the final beer product. An extensive sample set of commercial beers brewed with barley, wheat, corn and/or rice were analyzed by both direct infusion Fourier transform ion cyclotron mass spectrometry (DI-FTICR MS, 400 samples) and UPLC-ToF-MS (100 samples). By its unrivaled mass resolution and accuracy, DI-FTICR-MS was able to uncover the compositional space of both polar and non-polar metabolites that can be traced back to the use of different starch sources. Reversed phase UPLC-ToF-MS was used to access information about molecular structures (MS²-fragmentation spectra) and isomeric separation, with a focus on less polar compounds. Both analytical approaches were able to achieve a clear statistical differentiation (OPLS-DA) of beer samples and reveal metabolic profiles according to the starch source. A mass difference network analysis, applied to the exact marker masses resolved by FTICR, showed a network of potential secondary metabolites specific to wheat, corn and rice. By MS²-similarity networks, database and literature search, we were able to identify metabolites and compound classes significant for the use of the different starch sources. Those were also found in the corresponding brewing raw materials, confirming the potential of our approach for quality control and monitoring. Our results also include the identification of the aspartic acid-conjugate of N-β-D-glucopyranosyl-indole-3-acetic acid as a potential marker for the use of rice in the brewing industry regarding quality control and food inspection purposes.

Genome-wide identification of ZmSnRK2 genes and functional analysis of ZmSnRK2.10 in ABA signaling pathway in maize (Zea mays L)

BACKGROUND

Phytohormone abscisic acid (ABA) is involved in the regulation of a wide range of biological processes. In Arabidopsis, it has been well-known that SnRK2s are the central components of the ABA signaling pathway that control the balance between plant growth and stress response, but the functions of ZmSnRK2 in maize are rarely reported. Therefore, the study of ZmSnRK2 is of great importance to understand the ABA signaling pathways in maize.

RESULTS

In this study, 14 ZmSnRK2 genes were identified in the latest version of maize genome database. Phylogenetic analysis revealed that ZmSnRK2s are divided into three subclasses based on their diversity of C-terminal domains. The exon-intron structures, phylogenetic, synteny and collinearity analysis indicated that SnRK2s, especially the subclass III of SnRK2, are evolutionally conserved in maize, rice and Arabidopsis. Subcellular localization showed that ZmSnRK2 proteins are localized in the nucleus and cytoplasm. The RNA-Seq datasets and qRT-PCR analysis showed that ZmSnRK2 genes exhibit spatial and temporal expression patterns during the growth and development of different maize tissues, and the transcript levels of some ZmSnRK2 genes in kernel are significantly induced by ABA and sucrose treatment. In addition, we found that ZmSnRK2.10, which belongs to subclass III, is highly expressed in kernel and activated by ABA. Overexpression of ZmSnRK2.10 partially rescued the ABA-insensitive phenotype of snrk2.2/2.3 double and snrk2.2/2.3/2.6 triple mutants and led to delaying plant flowering in Arabidopsis.

CONCLUSION

The SnRK2 gene family exhibits a high evolutionary conservation and has expanded with whole-genome duplication events in plants. The ZmSnRK2s expanded in maize with whole-genome and segmental duplication, not tandem duplication. The expression pattern analysis of ZmSnRK2s in maize offers important information to study their functions. Study of the functions of ZmSnRK.10 in Arabidopsis suggests that the ABA-dependent members of SnRK2s are evolutionarily conserved in plants. Our study elucidated the structure and evolution of SnRK2 genes in plants and provided a basis for the functional study of ZmSnRK2s protein in maize.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

Frequency, Origins, and Evolutionary Role of Chromosomal Inversions in Plants

Chromosomal inversions have the potential to play an important role in evolution by reducing recombination between favorable combinations of alleles. Until recently, however, most evidence for their likely importance derived from dipteran flies, whose giant larval salivary chromosomes aided early cytogenetic studies. The widespread application of new genomic technologies has revealed that inversions are ubiquitous across much of the plant and animal kingdoms. Here we review the rapidly accumulating literature on inversions in the plant kingdom and discuss what we have learned about their establishment and likely evolutionary role. We show that inversions are prevalent across a wide range of plant groups. We find that inversions are often associated with locally favored traits, as well as with traits that contribute to assortative mating, suggesting that they may be key to adaptation and speciation in the face of gene flow. We also discuss the role of inversions in sex chromosome formation, and explore possible parallels with inversion establishment on autosomes. The identification of inversion origins, as well as the causal variants within them, will advance our understanding of chromosomal evolution in plants.

Mapping and Identifying a Candidate Gene Plr4, a Recessive Gene Regulating Purple Leaf in Rice, by Using Bulked Segregant and Transcriptome Analysis with Next-Generation Sequencing

The anthocyanin biosynthesis of rice is a major concern due to the potential nutritional value. Purple appears in various organs and tissues of rice such as pericarp, flower organs, leaves, leaf sheaths, internodes, ligules, apex, and stigma. At present, there are many studies on the color of rice pericarp, but the gene and mechanism of other organs such as leaves are still unclear, and the gene regulatory network of specific organ coloring has not been systematically understood. In this study, genetic analysis demonstrated that the purple leaf traits of rice were regulated by a recessive gene. The green leaf cultivar Y58S and purple leaf cultivar XianHongB were used to construct the mapping population. A set of near isogenicline (NIL) (BC₃F₁) was bred via crossing and back-crossing. The generations of BC₃F₂ appeared to separate four phenotypes, pl1, pl2, pl3, and pl4, due to the occurrence of a purple color in different organs. We constructed three bulked segregant analysis (BSA) pools (pl1-pl2, pl1-pl3, and pl1-pl4) by using the separated generations of BC₃F₅ and mapped the purple leaf gene plr4 to the vicinity of 27.9-31.1 Mb on chromosome 4. Subsequently, transcriptome sequencing (RNA-Seq) for pl3 and pl2 was used to analyze the differentially expressed genes in the localization interval, where 12 unigenes exhibited differential expression in which two genes (Os04g0577800, Os04g0616400) were downregulated. The two downregulated genes (Os04g0577800 and Os04g0616400) are possible candidate genes because of the recessive genetic characteristics of the purple leaf genes. These results will facilitate the cloning of plr4 and illustrate the molecular mechanisms of the anthocyanin synthesis pathway.

A study of CCD8 genes/proteins in seven monocots and eight dicots

In plants, the enzyme CCD8 (carotenoid cleavage dioxygenase 8) is involved in the synthesis of an important hormone, strigolactone, and therefore, plays an important role in controlling growth and development. Using cDNA and protein sequence derived from the gene ZmCCD8 from maize, we identified putative orthologs of the gene encoding CCD8 in six other monocots and eight dicots; the sequence similarity ranged from 52–75.9% at the gene level and 60.9–93.7% at the protein level. The average length of the gene was ~3.3 kb (range: 2.08 to 3.98 kb), although the number of introns within the genes differed (4 or 5 in dicots and 3 or 4 in monocots, except in T. urartu with 6 introns). Several cis-acting regulatory elements were identified in the promoters of CCD8 genes, which are known to respond to biotic and abiotic stresses. The N-terminal end (up to ~70 amino acids) of CCD8 proteins was highly variable due to insertions, deletions and mismatches. The variation in genes and proteins were particularly conspicuous in T. urartu and Ae. tauschii among the monocots and A. thaliana and P. persica among the dicots. In CCD8 proteins, 12 motifs were also identified, of which 6 were novel; 4 of these novel motifs occurred in all the 15 species. The 3D structures of proteins had the characteristic features of the related enzyme apocarotenoid oxygenase (ACO) of Synechocystis (a representative of cyanobacteria). The results of qRT-PCR in wheat revealed that under phosphorous (P)-starved condition (relative to expression under optimum P used as control), the expression of TaCCD8 genes increased ~37 fold in root tissue of the cultivar C306 and ~33 fold in shoot tissue of the cultivar HUW468 (the two cultivars differed in their P-use efficiency). This suggested that expression of TaCCD8 genes is genotype-dependent and tissue-specific and is regulated under different levels of P supply.

Identification of Soybean Genes Related to Soybean Seed Protein Content Based on Quantitative Trait Loci Collinearity Analysis

Increasing the protein content of soybean seeds through a higher ratio of glycinin is important for soybean breeding and food processing; therefore, the integration of different quantitative trait loci (QTLs) is of great significance. In this study, we investigated the collinearity of seed protein QTLs. We identified 192 collinear protein QTLs that formed six hotspot regions. The two most important regions had seed protein 36-10 and seed protein 36-20 as hub nodes. We used a chromosome segment substitution line (CSSL) population for QTL validation and identified six CSSL materials with collinear QTLs. Five materials with higher protein and glycinin contents in comparison to the recurrent parent were analyzed. A total of 13 candidate genes related to seed protein from the QTL hotspot intervals were detected, 8 of which had high expression in mature soybean seeds. These results offer a new analysis method for molecular-assisted selection (MAS) and improvement of soybean product quality.

Metabolic GWAS-based dissection of genetic bases underlying the diversity of plant metabolism

Plants have served as sources providing humans with metabolites for food and nutrition, biomaterials for living, and treatment for pain and disease. Plants produce a huge array of metabolites, with an immense diversity at both the population and individual levels. Dissection of the genetic bases for metabolic diversity has attracted increasing research attention. The concept of genome-wide association study (GWAS) was extended to studies on the diversity of plant metabolome that benefitted from the development of mass-spectrometry-based analytical systems and genome sequencing technologies. Metabolic genome-wide association study (mGWAS) is one of the most powerful tools for global identification of genetic determinants for diversity of plant metabolism. Recently, mGWAS has been performed for various species with continuous improvements, providing deeper insights into the genetic bases of metabolic diversity. In this review, we discuss fully the achievements to date and remaining challenges that are associated with both mGWAS and mGWAS-based multi-dimensional analysis. We begin with a summary of GWAS and its development based on statistical methods and populations. As variation in targeted traits is essential for GWAS, we review metabolic diversity and its rise at both the population and individual levels. Subsequently, the application of mGWAS for plants and its corresponding achievements are fully discussed. We address the current knowledge on mGWAS-based multi-dimensional analysis and emerging insights into the diversity of metabolism.

Exploring the Diversity of Plant Metabolism

Plants produce a huge array of metabolites, far more than those produced by most other organisms. Unraveling this diversity and its underlying genetic variation has attracted increasing research attention. Post-genomic profiling platforms have enabled the marriage and mining of the enormous amount of phenotypic and genetic diversity. We review here achievements to date and challenges remaining that are associated with plant metabolic research using multi-omic strategies. We focus mainly on strategies adopted in investigating the diversity of plant metabolism and its underlying features. Recent advances in linking metabotypes with phenotypic and genotypic traits are also discussed. Taken together, we conclude that exploring the diversity of metabolism could provide new insights into plant evolution and domestication.

Effect of Co-segregating Markers on High-Density Genetic Maps and Prediction of Map Expansion Using Machine Learning Algorithms

Advances in sequencing and genotyping methods have enable cost-effective production of high throughput single nucleotide polymorphism (SNP) markers, making them the choice for linkage mapping. As a result, many laboratories have developed high-throughput SNP assays and built high-density genetic maps. However, the number of markers may, by orders of magnitude, exceed the resolution of recombination for a given population size so that only a minority of markers can accurately be ordered. Another issue attached to the so-called 'large p, small n' problem is that high-density genetic maps inevitably result in many markers clustering at the same position (co-segregating markers). While there are a number of related papers, none have addressed the impact of co-segregating markers on genetic maps. In the present study, we investigated the effects of co-segregating markers on high-density genetic map length and marker order using empirical data from two populations of wheat, Mohawk × Cocorit (durum wheat) and Norstar × Cappelle Desprez (bread wheat). The maps of both populations consisted of 85% co-segregating markers. Our study clearly showed that excess of co-segregating markers can lead to map expansion, but has little effect on markers order. To estimate the inflation factor (IF), we generated a total of 24,473 linkage maps (8,203 maps for Mohawk × Cocorit and 16,270 maps for Norstar × Cappelle Desprez). Using seven machine learning algorithms, we were able to predict with an accuracy of 0.7 the map expansion due to the proportion of co-segregating markers. For example in Mohawk × Cocorit, with 10 and 80% co-segregating markers the length of the map inflated by 4.5 and 16.6%, respectively. Similarly, the map of Norstar × Cappelle Desprez expanded by 3.8 and 11.7% with 10 and 80% co-segregating markers. With the increasing number of markers on SNP-chips, the proportion of co-segregating markers in high-density maps will continue to increase making map expansion unavoidable. Therefore, we suggest developers improve linkage mapping algorithms for efficient analysis of high-throughput data. This study outlines a practical strategy to estimate the IF due to the proportion of co-segregating markers and outlines a method to scale the length of the map accordingly.

Quantitative trait loci from identification to exploitation for crop improvement

Advancement in the field of genetics and genomics after the discovery of Mendel's laws of inheritance has led to map the genes controlling qualitative and quantitative traits in crop plant species. Mapping of genomic regions controlling the variation of quantitatively inherited traits has become routine after the advent of different types of molecular markers. Recently, the next generation sequencing methods have accelerated the research on QTL analysis. These efforts have led to the identification of more closely linked molecular markers with gene/QTLs and also identified markers even within gene/QTL controlling the trait of interest. Efforts have also been made towards cloning gene/QTLs or identification of potential candidate genes responsible for a trait. Further new concepts like crop QTLome and QTL prioritization have accelerated precise application of QTLs for genetic improvement of complex traits. In the past years, efforts have also been made in exploitation of a number of QTL for improving grain yield or other agronomic traits in various crops through markers assisted selection leading to cultivation of these improved varieties at farmers' field. In present article, we reviewed QTLs from their identification to exploitation in plant breeding programs and also reviewed that how improved cultivars developed through introgression of QTLs have improved the yield productivity in many crops.

Comparative genomics of Eucalyptus and Corymbia reveals low rates of genome structural rearrangement

BACKGROUND

Previous studies suggest genome structure is largely conserved between Eucalyptus species. However, it is unknown if this conservation extends to more divergent eucalypt taxa. We performed comparative genomics between the eucalypt genera Eucalyptus and Corymbia. Our results will facilitate transfer of genomic information between these important taxa and provide further insights into the rate of structural change in tree genomes.

RESULTS

We constructed three high density linkage maps for two Corymbia species (Corymbia citriodora subsp. variegata and Corymbia torelliana) which were used to compare genome structure between both species and Eucalyptus grandis. Genome structure was highly conserved between the Corymbia species. However, the comparison of Corymbia and E. grandis suggests large (from 1-13 MB) intra-chromosomal rearrangements have occurred on seven of the 11 chromosomes. Most rearrangements were supported through comparisons of the three independent Corymbia maps to the E. grandis genome sequence, and to other independently constructed Eucalyptus linkage maps.

CONCLUSIONS

These are the first large scale chromosomal rearrangements discovered between eucalypts. Nonetheless, in the general context of plants, the genomic structure of the two genera was remarkably conserved; adding to a growing body of evidence that conservation of genome structure is common amongst woody angiosperms.

Comparative and parallel genome-wide association studies for metabolic and agronomic traits in cereals

The plant metabolome is characterized by extensive diversity and is often regarded as a bridge between genome and phenome. Here we report metabolic and phenotypic genome-wide studies (mGWAS and pGWAS) in rice grain that, in addition to previous metabolic GWAS in rice leaf and maize kernel, show both distinct and overlapping aspects of genetic control of metabolism within and between species. We identify new candidate genes potentially influencing important metabolic and/or morphological traits. We show that the differential genetic architecture of rice metabolism between different tissues is in part determined by tissue specific expression. Using parallel mGWAS and pGWAS we identify new candidate genes potentially responsible for variation in traits such as grain colour and size, and provide evidence of metabotype-phenotype linkage. Our study demonstrates a powerful strategy for interactive functional genomics and metabolomics in plants, especially the cloning of minor QTLs for complex phenotypic traits.

Whole Genome Duplications Shaped the Receptor Tyrosine Kinase Repertoire of Jawed Vertebrates

The receptor tyrosine kinase (RTK) gene family, involved primarily in cell growth and differentiation, comprises proteins with a common enzymatic tyrosine kinase intracellular domain adjacent to a transmembrane region. The amino-terminal portion of RTKs is extracellular and made of different domains, the combination of which characterizes each of the 20 RTK subfamilies among mammals. We analyzed a total of 7,376 RTK sequences among 143 vertebrate species to provide here the first comprehensive census of the jawed vertebrate repertoire. We ascertained the 58 genes previously described in the human and mouse genomes and established their phylogenetic relationships. We also identified five additional RTKs amounting to a total of 63 genes in jawed vertebrates. We found that the vertebrate RTK gene family has been shaped by the two successive rounds of whole genome duplications (WGD) called 1R and 2R (1R/2R) that occurred at the base of the vertebrates. In addition, the Vegfr and Ephrin receptor subfamilies were expanded by single gene duplications. In teleost fish, 23 additional RTK genes have been retained after another expansion through the fish-specific third round (3R) of WGD. Several lineage-specific gene losses were observed. For instance, birds have lost three RTKs, and different genes are missing in several fish sublineages. The RTK gene family presents an unusual high gene retention rate from the vertebrate WGDs (58.75% after 1R/2R, 64.4% after 3R), resulting in an expansion that might be correlated with the evolution of complexity of vertebrate cellular communication and intracellular signaling.

Paralog analyses reveal gene duplication events and genes under positive selection in Ixodes scapularis and other ixodid ticks

Background

Hard ticks (family Ixodidae) are obligatory hematophagous ectoparasites of worldwide medical and veterinary importance. The haploid genomes of multiple species of ixodid ticks exceed 1 Gbp, prompting questions regarding gene, segmental and whole genome duplication in this phyletic group. The availability of the genome assembly for the black legged tick, Ixodes scapularis, and transcriptome datasets for multiple species of ticks offers an opportunity to assess the contribution of gene duplication to the genome. Here we developed a bioinformatics pipeline to identify and analyze duplicated genes (paralogs) using gene models from the prostriate tick, I. scapularis IscaW1.1 annotation and expressed sequence tags (ESTs) from I. scapularis and the metastriate ticks, Rhipicephalus microplus (southern cattle tick), R. appendiculatus (brown ear tick) and Amblyomma variegatum (tropical bont tick).

Results

Approximately 1-2 % of I. scapularis gene models and 2-14 % of ESTs from the four species represent duplicated genes. The ratio of non-synonymous to synonymous nucleotide substitution rates suggests ~ 25 % of duplicated genes are under positive selection pressure in each species. Analyses of synonymous substitution rates provide evidence for two duplication events in I. scapularis and R. microplus involving several hundred genes. Conservative molecular clock estimates based on synonymous substitution rates for species of Anopheles mosquitoes and the fruit fly, Drosophila melanogaster, suggest these events occurred within the last 50 MYA. Mapping of paralogs to the I. scapularis genome assembly supports tandem, or possibly segmental duplication events.

Conclusions

The present study marks the first genome-level analyses of gene duplication for the Ixodidae and provides insights into mechanisms shaping genome evolution in this group. At least two duplication events involving hundreds of genes may have occurred independently in the lineages prostriata and metastriata, with tandem and segmental duplication the most likely mechanisms for paralog generation. Duplicated genes under positive selection pressure may be linked to emerging functions in the tick and represent important candidates for further study.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2350-2) contains supplementary material, which is available to authorized users.

Genome-Wide Identification, Phylogenetic and Expression Analyses of the Ubiquitin-Conjugating Enzyme Gene Family in Maize

BACKGROUND

Ubiquitination is a post-translation modification where ubiquitin is attached to a substrate. Ubiquitin-conjugating enzymes (E2s) play a major role in the ubiquitin transfer pathway, as well as a variety of functions in plant biological processes. To date, no genome-wide characterization of this gene family has been conducted in maize (Zea mays).

METHODOLOGY/PRINCIPAL FINDINGS

In the present study, a total of 75 putative ZmUBC genes have been identified and located in the maize genome. Phylogenetic analysis revealed that ZmUBC proteins could be divided into 15 subfamilies, which include 13 ubiquitin-conjugating enzymes (ZmE2s) and two independent ubiquitin-conjugating enzyme variant (UEV) groups. The predicted ZmUBC genes were distributed across 10 chromosomes at different densities. In addition, analysis of exon-intron junctions and sequence motifs in each candidate gene has revealed high levels of conservation within and between phylogenetic groups. Tissue expression analysis indicated that most ZmUBC genes were expressed in at least one of the tissues, indicating that these are involved in various physiological and developmental processes in maize. Moreover, expression profile analyses of ZmUBC genes under different stress treatments (4°C, 20% PEG6000, and 200 mM NaCl) and various expression patterns indicated that these may play crucial roles in the response of plants to stress.

CONCLUSIONS

Genome-wide identification, chromosome organization, gene structure, evolutionary and expression analyses of ZmUBC genes have facilitated in the characterization of this gene family, as well as determined its potential involvement in growth, development, and stress responses. This study provides valuable information for better understanding the classification and putative functions of the UBC-encoding genes of maize.

Genome Alignment Spanning Major Poaceae Lineages Reveals Heterogeneous Evolutionary Rates and Alters Inferred Dates for Key Evolutionary Events

Multiple comparisons among genomes can clarify their evolution, speciation, and functional innovations. To date, the genome sequences of eight grasses representing the most economically important Poaceae (grass) clades have been published, and their genomic-level comparison is an essential foundation for evolutionary, functional, and translational research. Using a formal and conservative approach, we aligned these genomes. Direct comparison of paralogous gene pairs all duplicated simultaneously reveal striking variation in evolutionary rates among whole genomes, with nucleotide substitution slowest in rice and up to 48% faster in other grasses, adding a new dimension to the value of rice as a grass model. We reconstructed ancestral genome contents for major evolutionary nodes, potentially contributing to understanding the divergence and speciation of grasses. Recent fossil evidence suggests revisions of the estimated dates of key evolutionary events, implying that the pan-grass polyploidization occurred ∼96 million years ago and could not be related to the Cretaceous-Tertiary mass extinction as previously inferred. Adjusted dating to reflect both updated fossil evidence and lineage-specific evolutionary rates suggested that maize subgenome divergence and maize-sorghum divergence were virtually simultaneous, a coincidence that would be explained if polyploidization directly contributed to speciation. This work lays a solid foundation for Poaceae translational genomics.

Comparative mapping combined with homology-based cloning of the rice genome reveals candidate genes for grain zinc and iron concentration in maize

BACKGROUND

Grain zinc and iron concentration is a complex trait that is controlled by quantitative trait loci (QTL) and is important for maintaining body health. Despite the substantial effort that has been put into identifying QTL for grain zinc and iron concentration, the integration of independent QTL is useful for understanding the genetic foundation of traits. The number of QTL for grain zinc and iron concentration is relatively low in a single species. Therefore, combined analysis of different genomes may help overcome this challenge.

RESULTS

As a continuation of our work on maize, meta-analysis of QTL for grain zinc and iron concentration in rice was performed to identify meta-QTL (MQTL). Based on MQTL in rice and maize, comparative mapping combined with homology-based cloning was performed to identify candidate genes for grain zinc and iron concentration in maize. In total, 22 MQTL in rice, 4 syntenic MQTL-related regions, and 3 MQTL-containing candidate genes in maize (ortho-mMQTL) were detected. Two maize orthologs of rice, GRMZM2G366919 and GRMZM2G178190, were characterized as natural resistance-associated macrophage protein (NRAMP) genes and considered to be candidate genes. Phylogenetic analysis of NRAMP genes among maize, rice, and Arabidopsis thaliana further demonstrated that they are likely responsible for the natural variation of maize grain zinc and iron concentration.

CONCLUSIONS

Syntenic MQTL-related regions and ortho-mMQTL are prime areas for future investigation as well as for marker-assisted selection breeding programs. Furthermore, the combined method using the rice genome that was used in this study can shed light on other species and help direct future quantitative trait research. In conclusion, these results help elucidate the molecular mechanism that underlies grain zinc and iron concentration in maize.

A powerful tool for genome analysis in maize: development and evaluation of the high density 600 k SNP genotyping array

Background

High density genotyping data are indispensable for genomic analyses of complex traits in animal and crop species. Maize is one of the most important crop plants worldwide, however a high density SNP genotyping array for analysis of its large and highly dynamic genome was not available so far.

Results

We developed a high density maize SNP array composed of 616,201 variants (SNPs and small indels). Initially, 57 M variants were discovered by sequencing 30 representative temperate maize lines and then stringently filtered for sequence quality scores and predicted conversion performance on the array resulting in the selection of 1.2 M polymorphic variants assayed on two screening arrays. To identify high-confidence variants, 285 DNA samples from a broad genetic diversity panel of worldwide maize lines including the samples used for sequencing, important founder lines for European maize breeding, hybrids, and proprietary samples with European, US, semi-tropical, and tropical origin were used for experimental validation. We selected 616 k variants according to their performance during validation, support of genotype calls through sequencing data, and physical distribution for further analysis and for the design of the commercially available Affymetrix® Axiom® Maize Genotyping Array. This array is composed of 609,442 SNPs and 6,759 indels. Among these are 116,224 variants in coding regions and 45,655 SNPs of the Illumina® MaizeSNP50 BeadChip for study comparison. In a subset of 45,974 variants, apart from the target SNP additional off-target variants are detected, which show only a minor bias towards intermediate allele frequencies. We performed principal coordinate and admixture analyses to determine the ability of the array to detect and resolve population structure and investigated the extent of LD within a worldwide validation panel.

Conclusions

The high density Affymetrix® Axiom® Maize Genotyping Array is optimized for European and American temperate maize and was developed based on a diverse sample panel by applying stringent quality filter criteria to ensure its suitability for a broad range of applications. With 600 k variants it is the largest currently publically available genotyping array in crop species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-823) contains supplementary material, which is available to authorized users.

Diversification of the plant-specific hybrid glycine-rich protein (HyGRP) genes in cereals

Plant-specific hybrid proline- or glycine-rich proteins (HyP/GRPs) are involved in diverse gene functions including plant development and responses to biotic and abiotic stresses. The quantitative trait locus, qLTG3-1, enhances seed germination in rice under low-temperature conditions and encodes a member with a glycine-rich motif of the HyP/GRP family. The function of this gene may be related to the weakening of tissue covering the embryo during seed germination. In the present study, the diversification of the HyP/GRP gene family was elucidated in rice based on phylogenetic relationships and gene expression levels. At least 21 members of the HyP/GRP family have been identified in the rice genome and clustered in five regions on four chromosomes by tandem and chromosomal duplications. Of these, OsHyPRP05 (qLTG3-1) and its paralogous gene, OsHyPRP21, had a glycine-rich motif. Furthermore, orthologous genes with a glycine-rich motif and the HyP/GRP gene family were detected in four genome-sequenced monocots: 12 in barley, 10 in Brachypodium, 20 in maize, and 28 in sorghum, using a BLAST search of qLTG3-1 as the query. All members of the HyP/GRP family in these five species were classified into seven main groups, which were clustered together in these species. These results suggested that the HyP/GRP gene family was formed in the ancestral genome before the divergence of these species. The collinearity of chromosomal regions around qLTG3-1 and its orthologous genes were conserved among rice, Brachypodium, sorghum, and maize, indicating that qLTG3-1 and orthologous genes conserve gene function during seed germination.

Comparative genome analysis between Agrostis stolonifera and members of the Pooideae subfamily, including Brachypodium distachyon

Creeping bentgrass (Agrostis stolonifera, allotetraploid 2n = 4x = 28) is one of the major cool-season turfgrasses. It is widely used on golf courses due to its tolerance to low mowing and aggressive growth habit. In this study, we investigated genome relationships of creeping bentgrass relative to the Triticeae (a consensus map of Triticum aestivum, T. tauschii, Hordeum vulgare, and H. spontaneum), oat, rice, and ryegrass maps using a common set of 229 EST-RFLP markers. The genome comparisons based on the RFLP markers revealed large-scale chromosomal rearrangements on different numbers of linkage groups (LGs) of creeping bentgrass relative to the Triticeae (3 LGs), oat (4 LGs), and rice (8 LGs). However, we detected no chromosomal rearrangement between creeping bentgrass and ryegrass, suggesting that these recently domesticated species might be closely related, despite their memberships to different Pooideae tribes. In addition, the genome of creeping bentgrass was compared with the complete genome sequence of Brachypodium distachyon in Pooideae subfamily using both sequences of the above-mentioned mapped EST-RFLP markers and sequences of 8,470 publicly available A. stolonifera ESTs (AgEST). We discovered large-scale chromosomal rearrangements on six LGs of creeping bentgrass relative to B. distachyon. Also, a total of 24 syntenic blocks based on 678 orthologus loci were identified between these two grass species. The EST orthologs can be utilized in further comparative mapping of Pooideae species. These results will be useful for genetic improvement of Agrostis species and will provide a better understanding of evolution within Pooideae species.

Karyotype stability in the genus Phaseolus evidenced by the comparative mapping of the wild species Phaseolus microcarpus

The genus Phaseolus L. (Fabaceae) is monophyletic and comprises approximately 75 species distributed into two principal clades. The five cultivated species, including the common bean (Phaseolus vulgaris), were placed in clade B. Clade A comprises only wild species, with more limited distribution. In the present work, bacterial artificial chromosomes (BACs) previously mapped in common bean (2n = 22) were used as probes in fluorescent in situ hybridization (FISH) in this comparative study of Phaseolus microcarpus (2n = 22), a species from clade A. We also analyzed the chromomycin A3 (CMA)/4',6-diamidino-2-phenylindole (DAPI) banding pattern and the localization of rDNA and telomeric DNA sites. The single 45S rDNA site from P. microcarpus was mapped to chromosome 6, showing conservation to the P. vulgaris homeolog. Of the two 5S rDNA sites identified in both species, only the site on chromosome 10 appeared conserved. In spite of the phylogenetic distance between the two species, all of the single-copy BACs demonstrated conservation of synteny. However, four collinearity breaks were observed, probably caused by para- and pericentric inversions. Some variation in the repetitive fraction of the genome was also observed. Thus, a broader analysis of the genus confirms that few, rare inversions seem to represent the main karyotype changes during the evolution of this genus.

Characterization of starch phosphorylases in barley grains

BACKGROUND

Starch is synthesized in both leaves and storage tissues of plants. The role of starch syntheses and branching enzymes is well understood; however, the role of starch phosphorylase is not clear.

RESULTS

A gene encoding Pho1 from barley was characterized and starch phosphorylases from both developing and germinating grain were characterized and purified. Two activities were detected: one with a molecular mass of 110 kDa and the other of 95 kDa. It was demonstrated through the use of antisera that the 110 kDa activity was located in the amyloplast and could correspond to the polypeptide encoded by the Pho1 gene cloned. The 95 kDa activity was localized to the cytoplasm, most strongly expressed in germinating grain, and was classified as a Pho2-type sequence. Using RNAi technology to reduce the content of Pho1 in the grain to less than 30% of wild type did not lead to any visible phenotype, and no dramatic alterations in the structure of the starch were observed.

CONCLUSION

Two starch phosphorylase activities were identified and characterized in barley grains, and shown to be present during starch synthesis. However, their role in starch synthesis still remains to be elucidated.

Genome conservation among three legume genera detected with DNA markers

A set of 219 DNA clones derived from mungbean (Vigna radiata), cowpea (V. unguiculata), common bean (Phaseolus vulgaris), and soybean (Glycine max) were used to generate comparative linkage maps among mungbean, common bean, and soybean. The maps allowed an assessment of linkage conservation and collinearity among the three genomes. Mungbean and common bean, both of the subtribe Phaseolinae, exhibited a high degree of linkage conservation and preservation of marker order. Most linkage groups of mungbean consisted of only one or two linkage blocks from common bean (and vice versa). The situation was significantly different with soybean, a member of the subtribe Glycininae. Mungbean and common bean linkage groups were generally mosaics of short soybean linkage blocks, each only a few centimorgans in length. These results suggest that it would be fruitful to join maps of mungbean and common bean, while knowledge of conserved genomic blocks would be useful in increasing marker density in specific genomic regions for all three genera. These comparative maps may also contribute to enhanced understanding of legume evolution.

Comparison off genetic distance and order off DNA markers in five populations of rice

A group of about 300 evenly distributed DNA markers from a high density RFLP linkage map of rice constructed using an F2 population derived from a japonica variety, Nipponbare, and an indica variety, Kasalath, were used to evaluate gene order and genetic distance in four other rice mapping populations. The purpose of this study was to determine the degree to which information gained from the high density linkage map could be applied to other mapping populations, particularly with regard to its utility in bridging quantitative traits and molecular and physical mapping information. The mapping populations consisted of two F2 populations derived from Dao Ren Qiao/Fl-1084 and Kinandangputi/Fl-1007, recombinant inbred lines from Asominori/IR24, and a backcross population from Sasanishiki/Habataki//Sasanishiki. All DNA markers commonly mapped in the four populations showed the same linkage groups as in the Nipponbare/Kasalath linkage map with conserved linkage order. The genetic distance between markers among the different populations did not vary to a significant level in any of the 12 chromosomes. The differences in some markers could be attributed to the size of the population used in the construction of the linkage maps. Furthermore, the conservation of linkage order found in the distal region of chromosomes 11 and 12 was also confirmed in the RFLP maps based on the four populations of rice. These suggest that any major genetic information from the Nipponbare/Kasalath map can be expected to be approximately the same in other crosses or populations. This high density RFLP linkage map, which is being utilized in constructing a physical map of rice, can be very useful in interpreting genome structure with great accuracy in other populations. Key words : linkage map, japonica, indica, gene order, genetic distance.

Molecular-genetic maps for group 1 chromosomes of Triticeae species and their relation to chromosomes in rice and oat

Group 1 chromosomes of the Triticeae tribe have been studied extensively because many important genes have been assigned to them. In this paper, chromosome 1 linkage maps of Triticum aestivum, T. tauschii, and T. monococcum are compared with existing barley and rye maps to develop a consensus map for Triticeae species and thus facilitate the mapping of agronomic genes in this tribe. The consensus map that was developed consists of 14 agronomically important genes, 17 DNA markers that were derived from known-function clones, and 76 DNA markers derived from anonymous clones. There are 12 inconsistencies in the order of markers among seven wheat, four barley, and two rye maps. A comparison of the Triticeae group 1 chromosome consensus map with linkage maps of homoeologous chromosomes in rice indicates that the linkage maps for the long arm and the proximal portion of the short arm of group 1 chromosomes are conserved among these species. Similarly, gene order is conserved between Triticeae chromosome 1 and its homoeologous chromosome in oat. The location of the centromere in rice and oat chromosomes is estimated from its position in homoeologous group 1 chromosomes of Triticeae.

Molecular mapping of wheat. Homoeologous group 3

A prerequisite for molecular level genetic studies and breeding in wheat is a molecular marker map detailing its similarities with those of other grass species in the Gramineae family. We have constructed restriction fragment length polymorphism maps of the A-, B-, and D-genome chromosomes of homoeologous group 3 of hexaploid wheat (Triticum aestivum L. em. Thell) using 114 F7-8 lines from a synthetic x bread wheat cross. The map consists of 58 markers spanning 230 cM on chromosome 3A, 62 markers spanning 260 cM on 3B, and 40 markers spanning 171 cM on 3D. Thirteen libraries of genomic or cDNA clones from wheat, barley, and T. tauschii, the wheat D genome donor, are represented, facilitating the alignment and comparison of these maps with maps of other grass species. Twenty-four clones reveal homoeoloci on two of the three genomes and the associated linkages are largely comparable across genomes. A consensus sequence of orthologous loci in grass species genomes is assembled from this map and from existing maps of the chromosome-3 homoeologs in barley (Hordeum spp.), T. tauschii, and rice (Oryza spp.). It illustrates the close homoeology among the four species and the partial homoeology of wheat chromosome 3 with oat (Avena spp.) chromosome C. Two orthologous red grain color genes, R3 and R1, are mapped on chromosome arms 3BL and 3DL.

A composite map of expressed sequences in maize

A maize genetic map based mainly on expressed sequences has been constructed. The map incorporates data from four segregating populations. Three recombinant inbred line populations were derived from the nonreciprocal crosses between three inbred lines. A map derived from an independent F2 progeny from one of the crosses was also used. With a total of 521 genotyped individuals, accuracy in gene order is expected. Five sources of markers were used: (i) 109 loci corresponding to 69 genes of known function, (ii) 39 loci controlling protein position shifts revealed by two-dimensional electrophoresis, (iii) 8 isozyme loci, (iv) 17 loci corresponding to 14 sequenced cDNAs for which no homology was found in gene banks, and (v) 102 loci corresponding to 81 anonymous probes. As many loci were common to all maps, we tested heterogeneity between recombination fractions. The comparison of recombination fractions revealed: (i) a good correspondence between the maps derived from the same cross, (ii) few significant differences in interval distances, and (iii) global differences, which can reach 20% of the total map length. A composite map of 275 loci covering 1765 cM has been constructed. Key words : Zea mays L., RFLP, genetic map, molecular markers, proteins.

Molecular mapping of wheat. Homoeologous group 2

A molecular-marker map of bread wheat having many markers in common with other grasses in the Gramineae family is a prerequisite for molecular level genetic studies and breeding in this crop species. We have constructed restriction fragment length polymorphism maps of the A-, B-, and D-genome chromosomes of homoeologous group 2 of hexaploid wheat (Triticum aestivum L. em. Thell) using 114 F7 lines from a synthetic x bread wheat cross and clones from 11 libraries. Chromosomes 2A, 2B, and 2D comprise 57, 60, and 56 markers and each spans about 200 cM. Comparisons between chromosomes are facilitated by 26 sets of homoeoloci. Genes mapped include a heterologous abscisic acid responsive locus cloned as pBS128, the epidermal waxiness inhibitor W21, and two presumed leaf rust and stem rust resistance genes. Anomalies suggesting ancestral rearrangements in chromosome 2B are pointed out and features of wheat group 2 chromosomes that are common to barley (Hordeum vulgare L.), rice (Oryza spp.), and T. tauschii are discussed.

Molecular genetic maps of the group 6 chromosomes of hexaploid wheat (Triticum aestivum L. em. Thell.)

Restriction fragment length polymorphism (RFLP) maps of chromosomes 6A, 6B, and 6D of hexaploid wheat (Triticum aestivum L. em. Thell.) have been produced. They were constructed using a population of F7-8 recombinant inbred lines derived from a synthetic wheat x bread wheat cross. The maps consist of 74 markers assigned to map positions at a LOD >= 3 (29 markers assigned to 6A, 24 to 6B, and 21 to 6D) and 2 markers assigned to 6D ordered at a LOD of 2.7. Another 78 markers were assigned to intervals on the maps. The maps of 6A, 6B, and 6D span 178, 132, and 206 cM, respectively. Twenty-one clones detected orthologous loci in two homoeologues and 3 detected an orthologous locus in each chromosome. Orthologous loci are located at intervals of from 1.5 to 26 cM throughout 70% of the length of the linkage maps. Within this portion of the maps, colinearity (homosequentiality) among the three homoeologues is strongly indicated. The remainder of the linkage maps consists of three segments ranging in length from 47 to 60 cM. Colinearity among these chromosomes and other Triticeae homoeologous group 6 chromosomes is indicated and a consensus RFLP map derived from maps of the homoeologous group 6 chromosomes of hexaploid wheat, tetraploid wheat, Triticum tauschii, and barley is presented. Key words : RFLP, wheat, linkage maps, molecular markers.

Epigenetic neofunctionalisation and regulatory gene evolution in grasses

During plant evolution, genome duplication and subsequent selection acting on new gene pairs has frequently resulted in partition of gene functions, or acquisition of new functions. This 'sub- and neofunctionalisation' (subF and neoF) is held to have driven the expansion of key gene classes. One such gene class in maize (Zea mays) includes a pair of Polycomb group (PcG) protein genes that, unlike their single Arabidopsis (Arabidopsis thaliana) counterpart, are both parentally imprinted with only the maternal alleles being expressed in the seed endosperm. Surprisingly, this imprinting is regulated by different mechanisms in the two genes, resulting in different phasing of parent-specific expression. In this opinion article we propose that recruitment of different imprinting systems constitutes 'epigenetic neoF', and has enhanced maternal control over seed development, with a potential impact on the evolution of the large and persistent endosperms of cereal grains.

An EST-derived SNP and SSR genetic linkage map of cassava (Manihot esculenta Crantz)

Cassava (Manihot esculenta Crantz) is one of the most important food security crops in the tropics and increasingly being adopted for agro-industrial processing. Genetic improvement of cassava can be enhanced through marker-assisted breeding. For this, appropriate genomic tools are required to dissect the genetic architecture of economically important traits. Here, a genome-wide SNP-based genetic map of cassava anchored in SSRs is presented. An outbreeder full-sib (F1) family was genotyped on two independent SNP assay platforms: an array of 1,536 SNPs on Illumina's GoldenGate platform was used to genotype a first batch of 60 F1. Of the 1,358 successfully converted SNPs, 600 which were polymorphic in at least one of the parents and was subsequently converted to KBiosciences' KASPar assay platform for genotyping 70 additional F1. High-precision genotyping of 163 informative SSRs using capillary electrophoresis was also carried out. Linkage analysis resulted in a final linkage map of 1,837 centi-Morgans (cM) containing 568 markers (434 SNPs and 134 SSRs) distributed across 19 linkage groups. The average distance between adjacent markers was 3.4 cM. About 94.2% of the mapped SNPs and SSRs have also been localized on scaffolds of version 4.1 assembly of the cassava draft genome sequence. This more saturated genetic linkage map of cassava that combines SSR and SNP markers should find several applications in the improvement of cassava including aligning scaffolds of the cassava genome sequence, genetic analyses of important agro-morphological traits, studying the linkage disequilibrium landscape and comparative genomics.

Genome-enabled insights into legume biology

Legumes are the third-largest family of angiosperms, the second-most-important crop family, and a key source of biological nitrogen in agriculture. Recently, the genome sequences of Glycine max (soybean), Medicago truncatula, and Lotus japonicus were substantially completed. Comparisons among legume genomes reveal a key role for duplication, especially a whole-genome duplication event approximately 58 Mya that is shared by most agriculturally important legumes. A second and more recent genome duplication occurred only in the lineage leading to soybean. Outcomes of genome duplication, including gene fractionation and sub- and neofunctionalization, have played key roles in shaping legume genomes and in the evolution of legume-specific traits. Analysis of legume genome sequences also enables the discovery of legume-specific gene families and provides a framework for genome-wide association mapping that will target phenotypes of special importance in legumes. Translating genomic resources from sequenced species to less studied but still important "orphan" legumes will enhance prospects for world food production.

Cloning and characterization of a maize SnRK2 protein kinase gene confers enhanced salt tolerance in transgenic Arabidopsis

SnRK2 (sucrose non-fermenting 1-related protein kinases 2) represents a unique family of protein kinase in regulating signaling transduction in plants. Although the regulatory mechanisms of SnRK2 have been well demonstrated in Arabidopsis thaliana, their functions in maize are still unknown. In our study, we cloned an SnRK2 gene from maize, ZmSAPK8, which encoded a putative homolog of the rice SAPK8 protein. ZmSAPK8 had two copies in the maize genome and harbored eight introns in its coding region. We demonstrated that ZmSAPK8 expressed differentially in various organs of maize plants and was up-regulated by high-salinity and drought treatment. A green fluorescent protein (GFP)-tagged ZmSAPK8 showed subcellular localization in the cell membrane, cytoplasm and nucleus. In vitro kinase assays indicated that ZmSAPK8 preferred Mn(2+) to Mg(2+) as cofactor for phosphorylation, and Ser-182 and Thr-183 in activation loop was important for its activity. Heterologous overexpression of ZmSAPK8 in Arabidopsis could significantly strengthen tolerance to salt stress. Under salt treatment, ZmSAPK8-overexpressed transgenic plants exhibited higher germination rate and proline content, low electrolyte leakage and higher survival rate than wild type. Further analysis indicated that transgenic plants showed increased transcription of the stress-related genes, RD29A, RD29B, RAB18, ABI1, DREB2A and P5CS1, under high-salinity conditions. The results demonstrated that ZmSAPK8 was involved in diverse stress signal transduction. Moreover, no obvious adverse effects on growth and development in the ZmSAPK8-overexpressed transgenic plants implied that ZmSAPK8 was potentially useful in transgenic breeding to improve salt tolerance in crops.

Syntenic relationships between cucumber (Cucumis sativus L.) and melon (C. melo L.) chromosomes as revealed by comparative genetic mapping

BACKGROUND

Cucumber, Cucumis sativus L. (2n = 2 × = 14) and melon, C. melo L. (2n = 2 × = 24) are two important vegetable species in the genus Cucumis (family Cucurbitaceae). Both species have an Asian origin that diverged approximately nine million years ago. Cucumber is believed to have evolved from melon through chromosome fusion, but the details of this process are largely unknown. In this study, comparative genetic mapping between cucumber and melon was conducted to examine syntenic relationships of their chromosomes.

RESULTS

Using two melon mapping populations, 154 and 127 cucumber SSR markers were added onto previously reported F(2)- and RIL-based genetic maps, respectively. A consensus melon linkage map was developed through map integration, which contained 401 co-dominant markers in 12 linkage groups including 199 markers derived from the cucumber genome. Syntenic relationships between melon and cucumber chromosomes were inferred based on associations between markers on the consensus melon map and cucumber draft genome scaffolds. It was determined that cucumber Chromosome 7 was syntenic to melon Chromosome I. Cucumber Chromosomes 2 and 6 each contained genomic regions that were syntenic with melon chromosomes III+V+XI and III+VIII+XI, respectively. Likewise, cucumber Chromosomes 1, 3, 4, and 5 each was syntenic with genomic regions of two melon chromosomes previously designated as II+XII, IV+VI, VII+VIII, and IX+X, respectively. However, the marker orders in several syntenic blocks on these consensus linkage maps were not co-linear suggesting that more complicated structural changes beyond simple chromosome fusion events have occurred during the evolution of cucumber.

CONCLUSIONS

Comparative mapping conducted herein supported the hypothesis that cucumber chromosomes may be the result of chromosome fusion from a 24-chromosome progenitor species. Except for a possible inversion, cucumber Chromosome 7 has largely remained intact in the past nine million years since its divergence from melon. Meanwhile, many structural changes may have occurred during the evolution of the remaining six cucumber chromosomes. Further characterization of the genomic nature of Cucumis species closely related to cucumber and melon might provide a better understanding of the evolutionary history leading to modern cucumber.

QTL mapping in white spruce: gene maps and genomic regions underlying adaptive traits across pedigrees, years and environments

BACKGROUND

The genomic architecture of bud phenology and height growth remains poorly known in most forest trees. In non model species, QTL studies have shown limited application because most often QTL data could not be validated from one experiment to another. The aim of our study was to overcome this limitation by basing QTL detection on the construction of genetic maps highly-enriched in gene markers, and by assessing QTLs across pedigrees, years, and environments.

RESULTS

Four saturated individual linkage maps representing two unrelated mapping populations of 260 and 500 clonally replicated progeny were assembled from 471 to 570 markers, including from 283 to 451 gene SNPs obtained using a multiplexed genotyping assay. Thence, a composite linkage map was assembled with 836 gene markers.For individual linkage maps, a total of 33 distinct quantitative trait loci (QTLs) were observed for bud flush, 52 for bud set, and 52 for height growth. For the composite map, the corresponding numbers of QTL clusters were 11, 13, and 10. About 20% of QTLs were replicated between the two mapping populations and nearly 50% revealed spatial and/or temporal stability. Three to four occurrences of overlapping QTLs between characters were noted, indicating regions with potential pleiotropic effects. Moreover, some of the genes involved in the QTLs were also underlined by recent genome scans or expression profile studies.Overall, the proportion of phenotypic variance explained by each QTL ranged from 3.0 to 16.4% for bud flush, from 2.7 to 22.2% for bud set, and from 2.5 to 10.5% for height growth. Up to 70% of the total character variance could be accounted for by QTLs for bud flush or bud set, and up to 59% for height growth.

CONCLUSIONS

This study provides a basic understanding of the genomic architecture related to bud flush, bud set, and height growth in a conifer species, and a useful indicator to compare with Angiosperms. It will serve as a basic reference to functional and association genetic studies of adaptation and growth in Picea taxa. The putative QTNs identified will be tested for associations in natural populations, with potential applications in molecular breeding and gene conservation programs. QTLs mapping consistently across years and environments could also be the most important targets for breeding, because they represent genomic regions that may be least affected by G × E interactions.

Genome-wide identification, classification and analysis of heat shock transcription factor family in maize

Background

Heat shock response in eukaryotes is transcriptionally regulated by conserved heat shock transcription factors (Hsfs). Hsf genes are represented by a large multigene family in plants and investigation of the Hsf gene family will serve to elucidate the mechanisms by which plants respond to stress. In recent years, reports of genome-wide structural and evolutionary analysis of the entire Hsf gene family have been generated in two model plant systems, Arabidopsis and rice. Maize, an important cereal crop, has represented a model plant for genetics and evolutionary research. Although some Hsf genes have been characterized in maize, analysis of the entire Hsf gene family were not completed following Maize (B73) Genome Sequencing Project.

Results

A genome-wide analysis was carried out in the present study to identify all Hsfs maize genes. Due to the availability of complete maize genome sequences, 25 nonredundant Hsf genes, named ZmHsfs were identified. Chromosomal location, protein domain and motif organization of ZmHsfs were analyzed in maize genome. The phylogenetic relationships, gene duplications and expression profiles of ZmHsf genes were also presented in this study. Twenty-five ZmHsfs were classified into three major classes (class A, B, and C) according to their structural characteristics and phylogenetic comparisons, and class A was further subdivided into 10 subclasses. Moreover, phylogenetic analysis indicated that the orthologs from the three species (maize, Arabidopsis and rice) were distributed in all three classes, it also revealed diverse Hsf gene family expression patterns in classes and subclasses. Chromosomal/segmental duplications played a key role in Hsf gene family expansion in maize by investigation of gene duplication events. Furthermore, the transcripts of 25 ZmHsf genes were detected in the leaves by heat shock using quantitative real-time PCR. The result demonstrated that ZmHsf genes exhibit different expression levels in heat stress treatment.

Conclusions

Overall, data obtained from our investigation contributes to a better understanding of the complexity of the maize Hsf gene family and provides the first step towards directing future experimentation designed to perform systematic analysis of the functions of the Hsf gene family.

Following tetraploidy in maize, a short deletion mechanism removed genes preferentially from one of the two homologs

Following genome duplication and selfish DNA expansion, maize used a heretofore unknown mechanism to shed redundant genes and functionless DNA with bias toward one of the parental genomes.

Previous work in Arabidopsis showed that after an ancient tetraploidy event, genes were preferentially removed from one of the two homeologs, a process known as fractionation. The mechanism of fractionation is unknown. We sought to determine whether such preferential, or biased, fractionation exists in maize and, if so, whether a specific mechanism could be implicated in this process. We studied the process of fractionation using two recently sequenced grass species: sorghum and maize. The maize lineage has experienced a tetraploidy since its divergence from sorghum approximately 12 million years ago, and fragments of many knocked-out genes retain enough sequence similarity to be easily identifiable. Using sorghum exons as the query sequence, we studied the fate of both orthologous genes in maize following the maize tetraploidy. We show that genes are predominantly lost, not relocated, and that single-gene loss by deletion is the rule. Based on comparisons with orthologous sorghum and rice genes, we also infer that the sequences present before the deletion events were flanked by short direct repeats, a signature of intra-chromosomal recombination. Evidence of this deletion mechanism is found 2.3 times more frequently on one of the maize homeologs, consistent with earlier observations of biased fractionation. The over-fractionated homeolog is also a greater than 3-fold better target for transposon removal, but does not have an observably higher synonymous base substitution rate, nor could we find differentially placed methylation domains. We conclude that fractionation is indeed biased in maize and that intra-chromosomal or possibly a similar illegitimate recombination is the primary mechanism by which fractionation occurs. The mechanism of intra-chromosomal recombination explains the observed bias in both gene and transposon loss in the maize lineage. The existence of fractionation bias demonstrates that the frequency of deletion is modulated. Among the evolutionary benefits of this deletion/fractionation mechanism is bulk DNA removal and the generation of novel combinations of regulatory sequences and coding regions.

All genomes can accumulate dispensable DNA in the form of duplications of individual genes or even partial or whole genome duplications. Genomes also can accumulate selfish DNA elements. Duplication events specifically are often followed by extensive gene loss. The maize genome is particularly extreme, having become tetraploid 10 million years ago and played host to massive transposon amplifications. We compared the genome of sorghum (which is homologous to the pre-tetraploid maize genome) with the two identifiable parental genomes retained in maize. The two maize genomes differ greatly: one of the parental genomes has lost 2.3 times more genes than the other, and the selfish DNA regions between genes were even more frequently lost, suggesting maize can distinguish between the parental genomes present in the original tetraploid. We show that genes are actually lost, not simply relocated. Deletions were rarely longer than a single gene, and occurred between repeated DNA sequences, suggesting mis-recombination as a mechanism of gene removal. We hypothesize an epigenetic mechanism of genome distinction to account for the selective loss. To the extent that the rate of base substitutions tracks time, we neither support nor refute claims of maize allotetraploidy. Finally, we explain why it makes sense that purifying selection in mammals does not operate at all like the gene and genome deletion program we describe here.

On the tetraploid origin of the maize genome

Data from cytological and genetic mapping studies suggest that maize arose as a tetraploid. Two previous studies investigating the most likely mode of maize origin arrived at different conclusions. Gaut and Doebley proposed a segmental allotetraploid origin of the maize genome and estimated that the two maize progenitors diverged at 20.5 million years ago (mya). In a similar study, using larger data set, Brendel and colleagues suggested a single genome duplication at 16 mya. One of the key components of such analyses is to examine sequence divergence among strictly orthologous genes. In order to identify such genes, Lai and colleagues sequenced five duplicated chromosomal regions from the maize genome and the orthologous counterparts from the sorghum genome. They also identified the orthologous regions in rice. Using positional information of genetic components, they identified 11 orthologous genes across the two duplicated regions of maize, and the sorghum and rice regions. Swigonova et al. analyzed the 11 orthologues, and showed that all five maize chromosomal regions duplicated at the same time, supporting a tetraploid origin of maize, and that the two maize progenitors diverged from each other at about the same time as each of them diverged from sorghum, about 11.9 mya.