Identification of transcriptome SNPs for assessing allele-specific gene expression in a super-hybrid rice Xieyou9308

Integrated analysis of mRNA and miRNA expression profiling in rice backcrossed progenies (BC2F12) with different plant height

Inter-specific hybridization and backcrossing commonly occur in plants. The use of progeny generated from inter-specific hybridization and backcrossing has been developed as a novel model system to explore gene expression divergence. The present study investigated the analysis of gene expression and miRNA regulation in backcrossed introgression lines constructed from cultivated and wild rice. High-throughput sequencing was used to compare gene and miRNA expression profiles in three progeny lines (L1710, L1817 and L1730), with different plant heights resulting from the backcrossing of introgression lines (BC₂F₁₂) and their parents (O. sativa and O. longistaminata). A total of 25,387 to 26,139 mRNAs and 379 to 419 miRNAs were obtained in these rice lines. More differentially expressed genes and miRNAs were detected in progeny/O. longistaminata comparison groups than in progeny/O. sativa comparison groups. Approximately 80% of the genes and miRNAs showed expression level dominance to O. sativa, indicating that three progeny lines were closer to the recurrent parent, which might be influenced by their parental genome dosage. Approximately 16% to 64% of the differentially expressed miRNAs possessing coherent target genes were predicted, and many of these miRNAs regulated multiple target genes. Most genes were up-regulated in progeny lines compared with their parents, but down-regulated in the higher plant height line in the comparison groups among the three progeny lines. Moreover, certain genes related to cell walls and plant hormones might play crucial roles in the plant height variations of the three progeny lines. Taken together, these results provided valuable information on the molecular mechanisms of hybrid backcrossing and plant height variations based on the gene and miRNA expression levels in the three progeny lines.

Changes in gene expression between a soybean F1 hybrid and its parents are associated with agronomically valuable traits

Soybean [Glycine max (L.) Merr.] genetic diversity is limited because domesticated soybean has undergone multiple genetic bottlenecks. Its progenitor, the wild soybean [Glycine soja Siebold & Zucc], has not undergone the same intense selection and is much more genetically diverse than domesticated soybean. However, the agronomic importance of diversity in wild soybean is unclear, and its weedy nature makes assessment difficult. To address this issue, we chose for study a highly selected, adapted F4-derived progeny of wild soybean, NMS4-44-329. This breeding line is derived from the hybridization between G. max cultivar N7103 and G. soja PI 366122. Agronomic comparisons were made among N7103, NMS4-44-329 and their F1 and F2 progeny in replicated yield trials at two North Carolina locations. Significant F1 mid-parent heterosis was observed at each location for seed yield (189 and 223 kgha-1, P<0.05 and P<0.10, respectively), seed protein content (1.1g/100g, P<0.01) and protein production per hectare (101 and 100 kgha-1, P<0.01 and P<0.06, respectively). Increased yield, seed protein content and protein production per hectare in the hybrids suggested that wild soybean has the potential to improve agronomic traits in applied breeding. Comparisons of differentially-expressed genes in the hybrid vs. parents identified genes associated with N metabolism. Non-additive changes in gene expression in the hybrids relative to the parents could reasonably explain the improved protein levels in the F1 hybrids. Changes in gene expression were influenced by environmental effects; however, allele specific bias in the hybrids were well correlated between environments. We propose that changes in gene expression, both additive and non-additive, and changes in allele specific expression bias may explain agronomic traits, and be valuable tools for plant breeders in the assessment of breeding populations.

Transcriptome profiling and comparison of maize ear heterosis during the spikelet and floret differentiation stages

BACKGROUND

Hybridization is a prominent process in the evolution of crop plants that can give rise to gene expression variation, phenotypic novelty and heterosis. Maize is the most successful crop in utilizing heterosis. The development of hybrid maize ears exhibits strong heterotic vigor and greatly affects maize yield. However, a comprehensive perspective on transcriptional variation and its correlation with heterosis during maize ear development is not available.

RESULTS

Using RNA sequencing technology, we investigated the transcriptome profiles of maize ears in the spikelet and floret differentiation stages of hybrid ZD808 and its parents CL11 and NG5. Our results revealed that 53.9% (21,258) of maize protein-coding genes were transcribed in at least one genotype. In both development stages, significant numbers of genes were differentially expressed between the hybrid and its parents. Gene expression inheritance analysis revealed approximately 80% of genes were expressed additively, which suggested that the complementary effect may play a foundation role in maize ear heterosis. Among non-additively expressed genes, NG5-dominant genes were predominant. Analyses of the allele-specific gene expression in hybrid identified pervasive allelic imbalance and significant preferential expression of NG5 alleles in both developmental stages. The results implied that NG5 may provide beneficial alleles that contribute greatly to heterosis. Further comparison of parental and hybrid allele-specific expression suggested that gene expression variation is largely attributable to cis-regulatory variation in maize. The cis-regulatory variations tend to preserve the allelic expression levels in hybrid and result in additive expression. Comparison between the two development stages revealed that allele-specific expression and cis-/trans-regulatory variations responded differently to developmental cues, which may lead to stage-specific vigor phenotype during maize ear development.

CONCLUSION

Our research suggests that cis-regulated additive expression may fine-tune gene expression level into an optimal status and play a foundation role in maize ear heterosis. Our work provides a comprehensive insight into transcriptional variation and its correlation with heterosis during maize ear development. The knowledge gained from this study presents novel opportunity to improve our maize varieties.

Allele Workbench: transcriptome pipeline and interactive graphics for allele-specific expression

Sequencing the transcriptome can answer various questions such as determining the transcripts expressed in a given species for a specific tissue or condition, evaluating differential expression, discovering variants, and evaluating allele-specific expression. Differential expression evaluates the expression differences between different strains, tissues, and conditions. Allele-specific expression evaluates expression differences between parental alleles. Both differential expression and allele-specific expression have been studied for heterosis (hybrid vigor), where the hybrid has improved performance over the parents for one or more traits. The Allele Workbench software was developed for a heterosis study that evaluated allele-specific expression for a mouse F1 hybrid using libraries from multiple tissues with biological replicates. This software has been made into a distributable package, which includes a pipeline, a Java interface to build the database, and a Java interface for query and display of the results. The required input is a reference genome, annotation file, and one or more RNA-Seq libraries with optional replicates. It evaluates allelic imbalance at the SNP and transcript level and flags transcripts with significant opposite directional allele-specific expression. The Java interface allows the user to view data from libraries, replicates, genes, transcripts, exons, and variants, including queries on allele imbalance for selected libraries. To determine the impact of allele-specific SNPs on protein folding, variants are annotated with their effect (e.g., missense), and the parental protein sequences may be exported for protein folding analysis. The Allele Workbench processing results in transcript files and read counts that can be used as input to the previously published Transcriptome Computational Workbench, which has a new algorithm for determining a trimmed set of gene ontology terms. The software with demo files is available from https://code.google.com/p/allele-workbench. Additionally, all software is ready for immediate use from an Atmosphere Virtual Machine Image available from the iPlant Collaborative (www.iplantcollaborative.org).

Genome-wide transcriptome profiles of rice hybrids and their parents

Heterosis is a widely studied phenomenon in several plant species. However, its genetic basis still remains to be elucidated. In this study, we used RNA-seq data from two rice genotypes and their reciprocal hybrids, and used a combination of transcriptome profiling and allele-specific expression analysis to identify genes that are differentially expressed in the hybrids and their parents or expressed in an allele-specific manner. The differentially expressed genes (DEGs) were identified by a pairwise comparison of the four genotypes. Detailed annotation of DEGs suggested that these genes showed enrichment in some gene ontology categories, and they tend to have tissue-specific expression patterns compared to all genes. A total of 1033 (10.24%) of 10,195 genes with informative single nucleotide polymorphism (SNPs) were identified as ASE genes. These allele-specific expessed (ASE) genes showed a broader expression breadth suggesting that they function in diverse developmental stages. Among 1033 ASE genes, we also identified 45 ASE transcription factors belonging to 17 transcription factor families. These ASE transcription factors may act in trans to regulate gene expression in filial 1 (F₁) hybrids. Our analyses provide a comprehensive transcriptome profile of rice hybrids and their parents, and would be a useful resource for the rice research community.

Single nucleotide polymorphisms generated by genotyping by sequencing to characterize genome-wide diversity, linkage disequilibrium, and selective sweeps in cultivated watermelon

Background

A large single nucleotide polymorphism (SNP) dataset was used to analyze genome-wide diversity in a diverse collection of watermelon cultivars representing globally cultivated, watermelon genetic diversity. The marker density required for conducting successful association mapping depends on the extent of linkage disequilibrium (LD) within a population. Use of genotyping by sequencing reveals large numbers of SNPs that in turn generate opportunities in genome-wide association mapping and marker-assisted selection, even in crops such as watermelon for which few genomic resources are available. In this paper, we used genome-wide genetic diversity to study LD, selective sweeps, and pairwise F_ST distributions among worldwide cultivated watermelons to track signals of domestication.

Results

We examined 183 Citrullus lanatus var. lanatus accessions representing domesticated watermelon and generated a set of 11,485 SNP markers using genotyping by sequencing. With a diverse panel of worldwide cultivated watermelons, we identified a set of 5,254 SNPs with a minor allele frequency of ≥ 0.05, distributed across the genome. All ancestries were traced to Africa and an admixture of various ancestries constituted secondary gene pools across various continents. A sliding window analysis using pairwise F_ST values was used to resolve selective sweeps. We identified strong selection on chromosomes 3 and 9 that might have contributed to the domestication process. Pairwise analysis of adjacent SNPs within a chromosome as well as within a haplotype allowed us to estimate genome-wide LD decay. LD was also detected within individual genes on various chromosomes. Principal component and ancestry analyses were used to account for population structure in a genome-wide association study. We further mapped important genes for soluble solid content using a mixed linear model.

Conclusions

Information concerning the SNP resources, population structure, and LD developed in this study will help in identifying agronomically important candidate genes from the genomic regions underlying selection and for mapping quantitative trait loci using a genome-wide association study in sweet watermelon.

Electronic supplementary material

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

What is crop heterosis: new insights into an old topic

Heterosis (or hybrid vigor) refers to a natural phenomenon whereby hybrid offspring of genetically diverse individuals out-perform their parents in multiple traits including yield, adaptability and resistances to biotic and abiotic stressors. Innovations in technology and research continue to clarify the mechanisms underlying crop heterosis, however the intrinsic relationship between the biological basis of heterosis remain unclear. In this review, we aim to provide insight into the molecular genetic basis of heterosis by presenting recent advances in the 'omics' of heterosis and the role of non-coding regions, particularly in relation to energy-use efficiency. We propose that future research should focus on integrating the expanding datasets from different species and hybrid combinations, to mine key heterotic genes and unravel interactive 'omics' networks associated with heterosis. Improved understanding of heterosis and the biological basis for its manipulation in agriculture should help to streamline its use in enhancing crop productivity.

Novel method for analysis of allele specific expression in triploid Oryzias latipes reveals consistent pattern of allele exclusion

Assessing allele-specific gene expression (ASE) on a large scale continues to be a technically challenging problem. Certain biological phenomena, such as X chromosome inactivation and parental imprinting, affect ASE most drastically by completely shutting down the expression of a whole set of alleles. Other more subtle effects on ASE are likely to be much more complex and dependent on the genetic environment and are perhaps more important to understand since they may be responsible for a significant amount of biological diversity. Tools to assess ASE in a diploid biological system are becoming more reliable. Non-diploid systems are, however, not uncommon. In humans full or partial polyploid states are regularly found in both healthy (meiotic cells, polynucleated cell types) and diseased tissues (trisomies, non-disjunction events, cancerous tissues). In this work we have studied ASE in the medaka fish model system. We have developed a method for determining ASE in polyploid organisms from RNAseq data and we have implemented this method in a software tool set. As a biological model system we have used nuclear transplantation to experimentally produce artificial triploid medaka composed of three different haplomes. We measured ASE in RNA isolated from the livers of two adult, triploid medaka fish that showed a high degree of similarity. The majority of genes examined (82%) shared expression more or less evenly among the three alleles in both triploids. The rest of the genes (18%) displayed a wide range of ASE levels. Interestingly the majority of genes (78%) displayed generally consistent ASE levels in both triploid individuals. A large contingent of these genes had the same allele entirely suppressed in both triploids. When viewed in a chromosomal context, it is revealed that these genes are from large sections of 4 chromosomes and may be indicative of some broad scale suppression of gene expression.