Extensive allelic variation in gene expression in populus F1 hybrids

Identification and Functional Prediction of CircRNAs in Leaves of F1 Hybrid Poplars with Different Growth Potential and Their Parents

Circular RNAs (CircRNAs) regulate plant growth and development; however, their role in poplar heterosis is unclear. We identified 3722 circRNAs in poplar leaves, most of which were intergenic (57.2%) and exonic (40.2%). The expression of circRNAs in F1 hybrids with high growth potential was higher than that in those with low growth potential. Non-additive expression of circRNAs and single-parent expression of circRNAs (SPE-circRNAs) might regulate poplar heterosis through microRNA sponging and protein translation, respectively. DECs among F1 hybrids with different growth potentials might regulate the growth potential of poplar via microRNA sponging. Correlation analysis between circRNA expression and its parent gene expression showed that SPE-M circRNA (circRNAs expressed by male parent only) might regulate poplar heterosis by inhibiting parent gene expression, while other circRNAs might regulate poplar heterosis by enhancing parent gene expression. Weighted correlation network analysis of gene/circRNA expression showed that circRNAs mainly regulate poplar heterosis via carbohydrate metabolism, amino acid metabolism, energy metabolism, and material transport. In addition, we identified seven circRNAs that positively or negatively regulate poplar heterosis. Thus, non-additively expressed circRNAs and SPE circRNAs are involved in regulating poplar heterosis, and DECs among F1 hybrids with different growth potentials were involved in regulating poplar growth potential.

Nonadditive gene expression is correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow (Salix spp.)

Many studies have highlighted the complex and diverse basis for heterosis in inbred crops. Despite the lack of a consensus model, it is vital that we turn our attention to understanding heterosis in undomesticated, heterozygous, and polyploid species, such as willow (Salix spp.). Shrub willow is a dedicated energy crop bred to be fast-growing and high yielding on marginal land without competing with food crops. A trend in willow breeding is the consistent pattern of heterosis in triploids produced from crosses between diploid and tetraploid species. Here, we test whether differentially expressed genes are associated with heterosis in triploid families derived from diploid Salix purpurea, diploid Salix viminalis, and tetraploid Salix miyabeana parents. Three biological replicates of shoot tips from all family progeny and parents were collected after 12 weeks in the greenhouse and RNA extracted for RNA-Seq analysis. This study provides evidence that nonadditive patterns of gene expression are correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow. Expression-level dominance was most correlated with heterosis for biomass yield traits and was highly enriched for processes involved in starch and sucrose metabolism. In addition, there was a global dosage effect of parent alleles in triploid hybrids, with expression proportional to copy number variation. Importantly, differentially expressed genes between family parents were most predictive of heterosis for both field and greenhouse collected traits. Altogether, these data will be used to progress models of heterosis to complement the growing genomic resources available for the improvement of heterozygous perennial bioenergy crops.

Homology-Based Interactions between Small RNAs and Their Targets Control Dominance Hierarchy of Male Determinant Alleles of Self-Incompatibility in Arabidopsis lyrata

Self-incompatibility (SI) is conserved among members of the Brassicaceae plant family. This trait is controlled epigenetically by the dominance hierarchy of the male determinant alleles. We previously demonstrated that a single small RNA (sRNA) gene is sufficient to control the linear dominance hierarchy in Brassica rapa and proposed a model in which a homology-based interaction between sRNAs and target sites controls the complicated dominance hierarchy of male SI determinants. In Arabidopsis halleri, male dominance hierarchy is reported to have arisen from multiple networks of sRNA target gains and losses. Despite these findings, it remains unknown whether the molecular mechanism underlying the dominance hierarchy is conserved among Brassicaceae. Here, we identified sRNAs and their target sites that can explain the linear dominance hierarchy of Arabidopsis lyrata, a species closely related to A. halleri. We tested the model that we established in Brassica to explain the linear dominance hierarchy in A. lyrata. Our results suggest that the dominance hierarchy of A. lyrata is also controlled by a homology-based interaction between sRNAs and their targets.

Identification of Imprinted Genes Based on Homology: An Example of Fragaria vesca

Genomic imprinting has drawn increasing attention in plant biology in recent years. At present, hundreds of imprinted genes have been identified in various plants, and some of them have been reported to be evolutionarily conserved in plant species. In this research, 17 candidate genes in Fragaria vesca were obtained based on the homologous imprinted genes in Arabidopsis thaliana and other species. We further constructed reciprocal crosses of diploid strawberry (F. vesca) using the varieties 10-41 and 18-86 as the parents to investigate the conservation of these imprinted genes. Potentially informative single nucleotide polymorphisms (SNPs) were used as molecular markers of two parents obtained from candidate imprinted genes which have been cloned and sequenced. Meanwhile, we analyzed the SNP site variation ratios and parent-of-origin expression patterns of candidate imprinted genes at 10 days after pollination (DAP) endosperm and embryo for the hybrids of reciprocal cross, respectively. A total of five maternally expressed genes (MEGs), i.e., FvARI8, FvKHDP-2, FvDRIP2, FvBRO1, and FvLTP3, were identified in the endosperm, which did not show imprinting in the embryo. Finally, tissues expression analysis indicated that the five imprinted genes excluding FvDRIP2 mainly expressed in the endosperm. This is the first report on imprinted genes of Fragaria, and we provide a simple and rapid method based on homologous conservation to screen imprinted genes. The present study will provide a basis for further study of function and mechanism of genomic imprinting in F. vesca.

The Effects of Gene Duplication Modes on the Evolution of Regulatory Divergence in Wild and Cultivated Soybean

Regulatory changes include divergence in both cis-elements and trans-factors, which play roles in organismal evolution. Whole genome duplications (WGD) followed by diploidization are a recurrent feature in the evolutionary history of angiosperms. Prior studies have shown that duplicated genes have different evolutionary fates due to variable selection constraints and results in genomic compositions with hallmarks of paleopolyploidy. The recent sequential WGDs and post-WGD evolution in the common ancestor of cultivated soybean (Glycine max) and wild soybean (Glycine soja), together with other models of gene duplication, have resulted in a highly duplicated genome. In this study, we investigated the transcriptional changes in G. soja and G. max. We identified a sizable proportion of interspecific differentially expressed genes (DEGs) and found parental expression level dominance of G. max in their F1 hybrids. By classifying genes into different regulatory divergence types, we found the trans-regulatory changes played a predominant role in transcriptional divergence between wild and cultivated soybean. The same gene ontology (GO) and protein family (Pfam) terms were found to be over-represented in DEGs and genes of cis-only between JY47 and GS, suggesting the substantial contribution of cis-regulatory divergences to the evolution of wild and cultivated soybeans. By further dissecting genes into five different duplication modes, we found genes in different duplication modes tend to accumulate different types of regulatory differences. A relatively higher proportion of cis-only regulatory divergences was detected in singleton, dispersed, proximal, and tandem duplicates than WGD duplicates and genome-wide level, which is in line with the prediction of gene balance hypothesis for the differential fates of duplicated genes post-WGD. The numbers of cis-only and trans-only regulated genes were similar for singletons, whereas there were more genes of trans-only than cis-only in the rest duplication types, especially in WGD in which there were two times more trans-only genes than that in cis-only type. Tandem duplicates showed the highest proportion of trans-only genes probably due to some special features of this class. In summary, our results demonstrate that genes in different duplication modes have different fates in transcriptional evolution underpinned by cis- or trans-regulatory divergences in soybean and likely in other paleopolyploid higher organisms.

Genome-wide analysis of spatiotemporal allele-specific expression in F1 hybrids of meat- and egg-type chickens

In diploid organisms, each gene locus is composed of two parental alleles, which would interact with each other for determining the phenotypic variation. Better understanding of the allele-specific expression (ASE) in farm animals is much important to explore the genetic basis underlying economically important traits, which have been poorly understood yet. In this study, genome-wide analysis was applied to explore the spatiotemporal pattern of ASE in the F1 hybrids of chicken. First, meat- and egg-type chickens were selected for producing a full-sib F1 hybrid population (n = 57). Then, genome resequencing of two parents and 38 offspring were performed and liver and breast muscle samples (n = 38) were subjected to strand-specific RNA sequencing (ssRNA-seq) for ASE detection at 1, 28, and 56 days of age, respectively. The results accurately identified a total of 465 informative genes that could be distinguished with respect to their parental origins. There were 0.4% - 4.1% of informative genes showing ASE, and 57 of them were found across different tissues and time points. Besides, most ASE genes in chickens were tissue-specific, and no matter what the time-point pattern of one ASE gene, the same parental allele of this gene almost showed consistently higher or lower expression across all time points in the same type tissue. In conclusion, this study indicated that most of ASE genes were tissue-specific and time-dependent.

Cas9-mediated mutagenesis of potato starch-branching enzymes generates a range of tuber starch phenotypes

We investigated whether Cas9-mediated mutagenesis of starch-branching enzymes (SBEs) in tetraploid potatoes could generate tuber starches with a range of distinct properties. Constructs containing the Cas9 gene and sgRNAs targeting SBE1, SBE2 or both genes were introduced by Agrobacterium-mediated transformation or by PEG-mediated delivery into protoplasts. Outcomes included lines with mutations in all or only some of the homoeoalleles of SBE genes and lines in which homoeoalleles carried several different mutations. DNA delivery into protoplasts resulted in mutants with no detectable Cas9 gene, suggesting the absence of foreign DNA. Selected mutants with starch granule abnormalities had reductions in tuber SBE1 and/or SBE2 protein that were broadly in line with expectations from genotype analysis. Strong reduction in both SBE isoforms created an extreme starch phenotype, as reported previously for low-SBE potato tubers. HPLC-SEC and 1 H NMR revealed a decrease in short amylopectin chains, an increase in long chains and a large reduction in branching frequency relative to wild-type starch. Mutants with strong reductions in SBE2 protein alone had near-normal amylopectin chain-length distributions and only small reductions in branching frequency. However, starch granule initiation was enormously increased: cells contained many granules of <4 μm and granules with multiple hila. Thus, large reductions in both SBEs reduce amylopectin branching during granule growth, whereas reduction in SBE2 alone primarily affects numbers of starch granule initiations. Our results demonstrate that Cas9-mediated mutagenesis of SBE genes has the potential to generate new, potentially valuable starch properties without integration of foreign DNA into the genome.

Transcriptional regulation of anthocyanin biosynthesis in a high-anthocyanin resynthesized Brassica napus cultivar

BACKGROUND

Anthocyanins are plant secondary metabolites with key roles in attracting insect pollinators and protecting against biotic and abiotic stresses. They have potential health-promoting effects as part of the human diet. Anthocyanin biosynthesis has been elucidated in many species, enabling the development of anthocyanin-enriched fruits, vegetables, and grains; however, few studies have investigated Brassica napus anthocyanin biosynthesis.

RESULTS

We developed a high-anthocyanin resynthesized B. napus line, Rs035, by crossing anthocyanin-rich B. rapa (A genome) and B. oleracea (C genome) lines, followed by chromosome doubling. We identified and characterized 73 and 58 anthocyanin biosynthesis genes in silico in the A and C genomes, respectively; these genes showed syntenic relationships with 41 genes in Arabidopsis thaliana and B. napus. Among the syntenic genes, twelve biosynthetic and six regulatory genes showed transgressively higher expression in Rs035, and eight structural genes and one regulatory gene showed additive expression. We identified three early-, four late-biosynthesis pathways, three transcriptional regulator genes, and one transporter as putative candidates enhancing anthocyanin accumulation in Rs035. Principal component analysis and Pearson's correlation coefficients corroborated the contribution of these genes to anthocyanin accumulation.

CONCLUSIONS

Our study lays the foundation for producing high-anthocyanin B. napus cultivars. The resynthesized lines and the differentially expressed genes we have identified could be used to transfer the anthocyanin traits to other commercial rapeseed lines using molecular and conventional breeding.

Meiotic chromosome stability of a newly formed allohexaploid wheat is facilitated by selection under abiotic stress as a spandrel

Polyploidy is a prominent route to speciation in plants; however, this entails resolving the challenges of meiotic instability facing abrupt doubling of chromosome complement. This issue remains poorly understood. We subjected progenies of a synthetic hexaploid wheat, analogous to natural common wheat, but exhibiting extensive meiotic chromosome instability, to heat or salt stress. We selected stress-tolerant cohorts and generated their progenies under normal condition. We conducted fluorescent in situ hybridization/genomic in situ hybridization-based meiotic/mitotic analysis, RNA-Seq and virus-induced gene silencing (VIGS)-mediated assay of meiosis candidate genes. We show that heritability of stress tolerance concurred with increased euploidy frequency due to enhanced meiosis stability. We identified a set of candidate meiosis genes with altered expression in the stress-tolerant plants vs control, but the expression was similar to that of common wheat (cv Chinese Spring, CS). We demonstrate VIGS-mediated downregulation of individual candidate meiosis genes in CS is sufficient to confer an unstable meiosis phenotype mimicking the synthetic wheat. Our results suggest that heritable regulatory changes of preexisting meiosis genes may be hitchhiked as a spandrel of stress tolerance, which significantly improves meiosis stability in the synthetic wheat. Our findings implicate a plausible scenario that the meiosis machinery in hexaploid wheat may have already started to evolve at its onset stage.

Transcriptomic analysis of transgressive segregants revealed the central role of photosynthetic capacity and efficiency in biomass accumulation in sugarcane

Sugarcane is among the most efficient crops in converting solar energy into chemical energy. However, due to its complex genome structure and inheritance, the genetic and molecular basis of biomass yield in sugarcane is still largely unknown. We created an F2 segregating population by crossing S. officinarum and S. spontaneum and evaluated the biomass yield of the F2 individuals. The F2 individuals exhibited clear transgressive segregation in biomass yield. We sequenced transcriptomes of source and sink tissues from 12 selected extreme segregants to explore the molecular basis of high biomass yield for future breeding of high-yielding energy canes. Among the 103,664 assembled unigenes, 10,115 and 728 showed significant differential expression patterns between the two extreme segregating groups in the top visible dewlap leaf and the 9^th culm internode, respectively. The most enriched functional categories were photosynthesis and fermentation in the high-biomass and the low-biomass groups, respectively. Our results revealed that high-biomass yield was mainly determined by assimilation of carbon in source tissues. The high-level expression of fermentative genes in the low-biomass group was likely induced by their low-energy status. Group-specific expression alleles which can be applied in the development of new high-yielding energy cane varieties via molecular breeding were identified.

Adiposity QTL Adip20 decomposes into at least four loci when dissected using congenic strains

An average mouse in midlife weighs between 25 and 30 g, with about a gram of tissue in the largest adipose depot (gonadal), and the weight of this depot differs between inbred strains. Specifically, C57BL/6ByJ mice have heavier gonadal depots on average than do 129P3/J mice. To understand the genetic contributions to this trait, we mapped several quantitative trait loci (QTLs) for gonadal depot weight in an F2 intercross population. Our goal here was to fine-map one of these QTLs, Adip20 (formerly Adip5), on mouse chromosome 9. To that end, we analyzed the weight of the gonadal adipose depot from newly created congenic strains. Results from the sequential comparison method indicated at least four rather than one QTL; two of the QTLs were less than 0.5 Mb apart, with opposing directions of allelic effect. Different types of evidence (missense and regulatory genetic variation, human adiposity/body mass index orthologues, and differential gene expression) implicated numerous candidate genes from the four QTL regions. These results highlight the value of mouse congenic strains and the value of this sequential method to dissect challenging genetic architecture.

Segmental allotetraploidy generates extensive homoeologous expression rewiring and phenotypic diversity at the population level in rice

Allopolyploidization, that is, concomitant merging and doubling of two or more divergent genomes in a common nucleus/cytoplasm, is known to instantly alter genomewide transcriptome dynamics, a phenomenon referred to as "transcriptomic shock." However, the immediate effects of transcriptomic alteration in generating phenotypic diversity at the population level remain underinvestigated. Here, we employed the MassARRAY-based Sequenom platform to assess and compare orthologous, allelic and homoeologous gene expression status in two tissues (leaf and root) of a set of randomly chosen individuals from populations of parental rice subspecies (indica and japonica), in vitro "hybrids" (parental mixes), reciprocal F1 hybrids and reciprocal tetraploids at the 5th-selfed generation (S5). We show that hybridization and whole genome duplication (WGD) have opposing effects on allelic and homoeologous expression in the F1 hybrids and tetraploids, respectively. Whereas hybridization exerts strong attenuating effects on allelic expression differences in diploid hybrids, WGD augments the intrinsic parental differences and generates extensive and variable homoeolog content which triggers diversification in expression patterning among the tetraploid plants. Coupled with the vast phenotypic diversity observed among the tetraploid individuals, our results provide experimental evidence in support of the notion that allopolyploidy catalyses rapid phenotypic diversification in higher plants. Our data further suggest that largely stochastic homoeolog content reshuffling rather than alteration in total expression level may be an important feature of evolution in young segmental allopolyploids, which underlies rapid expression diversity at the population level.

Evolutionary Dynamics of Regulatory Changes Underlying Gene Expression Divergence among Saccharomyces Species

Heritable changes in gene expression are important contributors to phenotypic differences within and between species and are caused by mutations in cis-regulatory elements and trans-regulatory factors. Although previous work has suggested that cis-regulatory differences preferentially accumulate with time, technical restrictions to closely related species and limited comparisons have made this observation difficult to test. To address this problem, we used allele-specific RNA-seq data from Saccharomyces species and hybrids to expand both the evolutionary timescale and number of species in which the evolution of regulatory divergence has been investigated. We find that as sequence divergence increases, cis-regulatory differences do indeed become the dominant type of regulatory difference between species, ultimately becoming a better predictor of expression divergence than trans-regulatory divergence. When both cis- and trans-regulatory differences accumulate for the same gene, they more often have effects in opposite directions than in the same direction, indicating widespread compensatory changes underlying the evolution of gene expression. The frequency of compensatory changes within and between species and the magnitude of effect for the underlying cis- and trans-regulatory differences suggests that compensatory changes accumulate primarily due to selection against divergence in gene expression as a result of weak stabilizing selection on gene expression levels. These results show that cis-regulatory differences and compensatory changes in regulation play increasingly important roles in the evolution of gene expression as time increases.

Genetic variation of growth dynamics in maize (Zea mays L.) revealed through automated non-invasive phenotyping

Hitherto, most quantitative trait loci of maize growth and biomass yield have been identified for a single time point, usually the final harvest stage. Through this approach cumulative effects are detected, without considering genetic factors causing phase-specific differences in growth rates. To assess the genetics of growth dynamics, we employed automated non-invasive phenotyping to monitor the plant sizes of 252 diverse maize inbred lines at 11 different developmental time points; 50 k SNP array genotype data were used for genome-wide association mapping and genomic selection. The heritability of biomass was estimated to be over 71%, and the average prediction accuracy amounted to 0.39. Using the individual time point data, 12 main effect marker-trait associations (MTAs) and six pairs of epistatic interactions were detected that displayed different patterns of expression at various developmental time points. A subset of them also showed significant effects on relative growth rates in different intervals. The detected MTAs jointly explained up to 12% of the total phenotypic variation, decreasing with developmental progression. Using non-parametric functional mapping and multivariate mapping approaches, four additional marker loci affecting growth dynamics were detected. Our results demonstrate that plant biomass accumulation is a complex trait governed by many small effect loci, most of which act at certain restricted developmental phases. This highlights the need for investigation of stage-specific growth affecting genes to elucidate important processes operating at different developmental phases.

Phenotypic Variations in the Foliar Chemical Profile of Persea americana Mill. cv. Hass

The Hass avocado tree Persea americana cv. Hass was derived from a single hybrid tree of P. americana var. drymifolia and P. americana var. guatemalensis, and it is propagated clonally by grafting. This cultivar is the most widely planted in the world but its profile of secondary metabolites has been studied rarely despite of its importance in plant protection. We illustrate the variability of the volatilome of mature leaves by describing the average chemical composition and the phenotypic variability found in 70 trees. Contrary to the uniformity expected in the Hass cultivar, high variability coefficients were found for most of the 36 detected foliar volatile compounds; furthermore we found six chemotypes grouping the foliar phenotypes of the sampled trees using hierarchical cluster analysis. About 48% of trees were grouped in one chemotype; five chemotypes grouped the remaining trees. The compounds that determined these chemotypes were: estragole, α-farnesene, β-caryophyllene, germacrene D, α-cubebene and eugenol. This striking variation in a cultivar propagated clonally is discussed in terms of somatic mutation.

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.

Gene regulation and speciation in house mice

One approach to understanding the process of speciation is to characterize the genetic architecture of post-zygotic isolation. As gene regulation requires interactions between loci, negative epistatic interactions between divergent regulatory elements might underlie hybrid incompatibilities and contribute to reproductive isolation. Here, we take advantage of a cross between house mouse subspecies, where hybrid dysfunction is largely unidirectional, to test several key predictions about regulatory divergence and reproductive isolation. Regulatory divergence between Mus musculus musculus and M. m. domesticus was characterized by studying allele-specific expression in fertile hybrid males using mRNA-sequencing of whole testes. We found extensive regulatory divergence between M. m. musculus and M. m. domesticus, largely attributable to cis-regulatory changes. When both cis and trans changes occurred, they were observed in opposition much more often than expected under a neutral model, providing strong evidence of widespread compensatory evolution. We also found evidence for lineage-specific positive selection on a subset of genes related to transcriptional regulation. Comparisons of fertile and sterile hybrid males identified a set of genes that were uniquely misexpressed in sterile individuals. Lastly, we discovered a nonrandom association between these genes and genes showing evidence of compensatory evolution, consistent with the idea that regulatory interactions might contribute to Dobzhansky-Muller incompatibilities and be important in speciation.

Distinct Regulatory Changes Underlying Differential Expression of TEOSINTE BRANCHED1-CYCLOIDEA-PROLIFERATING CELL FACTOR Genes Associated with Petal Variations in Zygomorphic Flowers of Petrocosmea spp. of the Family Gesneriaceae

CYCLOIDEA (CYC)-like genes, belonging to the plant-specific TCP transcription factor family that is named after TEOSINTE BRANCHED1 (TB1) from maize (Zea mays), CYC from Antirrhinum majus, and the PROLIFERATING CELL FACTORS (PCF) from rice (Oryza sativa), have conserved dorsal identity function in patterning floral zygomorphy mainly through specific expression in dorsal petals of a flower. Their expression changes are usually related to morphological diversity of zygomorphic flowers. However, it is still a challenge to elucidate the molecular mechanism underlying their expression differentiation. It is also unknown whether CINCINNATA (CIN)-like TCP genes, locally controlling cell growth and proliferation, are involved in the evolution of floral zygomorphy. To address these questions, we selected two closely related species, i.e. Petrocosmea glabristoma and Petrocosmea sinensis, with distinct petal morphology to conduct expression, hybridization, mutant, and allele-specific expression analyses. The results show that the size change of the dorsal petals between the two species is mainly mediated by the expression differentiation of CYC1C and CYC1D, while the shape variation of all petals is related to the expression change of CIN1. In reciprocal F1 hybrids, the expression of CYC1C, CYC1D, and CIN1 conforms to an additive inheritance mode, consistent with the petal phenotypes of hybrids. Through allele-specific expression analyses, we find that the expression differentiation of these TCP genes is underlain by distinctly different types of regulatory changes. We suggest that highly redundant paralogs with identical expression patterns and interspecific expression differentiation may be controlled by remarkably different regulatory pathways because natural selection may favor different regulatory modifications rather than coding sequence changes of key developmental genes in generating morphological diversity.

Transcriptomic changes following synthesis of a Populus full-sib diploid and allotriploid population with different heterozygosities driven by three types of 2n female gamete

Diploid gametes are usually applied to produce triploids of Populus [originating from first-division restitution (FDR), second-division restitution (SDR), and postmeiotic restitution (PMR) 2n eggs]. Three types of 2n gametes transmitted different parental heterozygosities in Populus. Failed spindle formation and no chromosomal separation to opposite poles during meiosis I mean that FDR 2n gametes carry nonsister chromatids that are potentially heterozygous. By contrast, SDR 2n gametes result from failed sister chromatid separation in meiosis II, and therefore, they carry sister chromatid that are potentially homozygous. Completely homozygous 2n gametes can arise from the PMR mechanism. The alteration of gene expression resulting from allopolyploidization is a prominent feature in plants. We compared gene expression in the full-sib progeny of three allotriploid Populus populations (triploid-F, triploid-S, and triploid-P) with that in its parent species, and their full-sib diploid F1 hybrid. Genome-wide expression level dominance was biased toward the maternal in the diploid F1 hybrid and three allotriploid populations, whereas our data indicated important, but different, effects of the transmission of different heterozygosity by 2n female gametes in the expression patterns of allopolyploids. Because of the higher level of heterozygosity, the triploids had higher rates of non-additive and transgressive expression patterns in the triploid-F than in triploid-S and triploid-P. Compared with diploid F1, about 30-fold more genes (251) were differently expressed in the triploid-F than in the triploid-S (9) and triploid-P (8), respectively. These findings indicate that hybridization and polyploidization have immediate and distinct effects on the large-scale patterns of gene expression, and different effects on the transmission of heterozygosity by three 2n female gametes.

Heterosis in early maize ear inflorescence development: a genome-wide transcription analysis for two maize inbred lines and their hybrid

Heterosis, or hybrid vigor, contributes to superior agronomic performance of hybrids compared to their inbred parents. Despite its importance, little is known about the genetic and molecular basis of heterosis. Early maize ear inflorescences formation affects grain yield, and are thus an excellent model for molecular mechanisms involved in heterosis. To determine the parental contributions and their regulation during maize ear-development-genesis, we analyzed genome-wide digital gene expression profiles in two maize elite inbred lines (B73 and Mo17) and their F₁ hybrid using deep sequencing technology. Our analysis revealed 17,128 genes expressed in these three genotypes and 22,789 genes expressed collectively in the present study. Approximately 38% of the genes were differentially expressed in early maize ear inflorescences from heterotic cross, including many transcription factor genes and some presence/absence variations (PAVs) genes, and exhibited multiple modes of gene action. These different genes showing differential expression patterns were mainly enriched in five cellular component categories (organelle, cell, cell part, organelle part and macromolecular complex), five molecular function categories (structural molecule activity, binding, transporter activity, nucleic acid binding transcription factor activity and catalytic activity), and eight biological process categories (cellular process, metabolic process, biological regulation, regulation of biological process, establishment of localization, cellular component organization or biogenesis, response to stimulus and localization). Additionally, a significant number of genes were expressed in only one inbred line or absent in both inbred lines. Comparison of the differences of modes of gene action between previous studies and the present study revealed only a small number of different genes had the same modes of gene action in both maize seedlings and ear inflorescences. This might be an indication that in different tissues or developmental stages, different global expression patterns prevail, which might nevertheless be related to heterosis. Our results support the hypotheses that multiple molecular mechanisms (dominance and overdominance modes) contribute to heterosis.

Analysis of the leaf methylomes of parents and their hybrids provides new insight into hybrid vigor in Populus deltoides

Background

Plants with heterosis/hybrid vigor perform better than their parents in many traits. However, the biological mechanisms underlying heterosis remain unclear. To investigate the significance of DNA methylation to heterosis, a comprehensive analysis of whole-genome DNA methylome profiles of Populus deltoides cl.'55/65' and '10/17' parental lines and their intraspecific F1 hybrids lines was performed using methylated DNA immunoprecipitation (MeDIP) and high-throughput sequencing.

Results

Here, a total of 486.27 million reads were mapped to the reference genome of Populus trichocarpa, with an average unique mapping rate of 57.8%. The parents with similar genetic background had distinct DNA methylation levels. F1 hybrids with hybrid vigor possessed non-additive DNA methylation level (their levels were higher than mid-parent values). The DNA methylation levels in promoter and repetitive sequences and transposable element of better-parent F1 hybrids and parents and lower-parent F1 hybrids were different. Compared with the maternal parent, better-parent F1 hybrids had fewer hypermethylated genes and more hypomethylated ones. Compared with the paternal parent and lower-parent L1, better-parent F1 hybrids had more hypermethylated genes and fewer hypomethylated ones. The differentially methylated genes between better-parent F1 hybrids, the parents and lower-parent F1 hybrids were enriched in the categories metabolic processes, response to stress, binding, and catalytic activity, development, and involved in hormone biosynthesis, signaling pathway.

Conclusions

The methylation patterns of the parents both partially and dynamically passed onto their hybrids, and F1 hybrids has a non-additive mathylation level. A multidimensional process is involved in the formation of heterosis.

Pollen transcriptome analysis of Solanum tuberosum (2n = 4x = 48), S. demissum (2n = 6x = 72), and their reciprocal F1 hybrids

Pollen mRNAs were different in reciprocal F 1 hybrids, which were probably caused by a cytoplasm-nuclear chromosomal genes interaction. We have found reciprocal differences in crossability between F1 hybrids of Solanum tuberosum (T) and a Mexican wild potato species S. demissum (D). When the reciprocal hybrids were crossed as pollen parents with S. demissum, a significantly higher berry-setting rate was obtained in TD compared with DT. In this study, we performed a whole-genome transcript analysis of the pollen mRNA using a high-throughput sequencer. We obtained 12.6 billion bases that were aligned into 13,020 transcripts with 9,366 loci. All possible genetic modes were observed between the parents and their progeny, where over-dominance and under-recessive types were relatively frequent (15.7 and 15.3 %, respectively). We found that 59.1 % of transcripts were more abundant in TD and over fourfold higher transcription levels were found in 66 TD transcripts and three DT transcripts. A higher proportion of over-dominance and a lower proportion of under-recessive transcription types were also observed in TD. The percentage contributions of multiple transcripts at the same locus varied greatly and were transcribed differently between species. In the new allelic combinations created by hybridization, approximately three-fourth of the transcripts had intermediate percentage contributions between the parents but no differential transcription patterns were apparent between the reciprocal hybrids. A broad spectrum of functionally different nuclear genes was over-represented in TD pollen, some of which were directly related to pollen behavior. Since TD and DT pollen had the same composition of nuclear genes, a cytoplasm-nuclear chromosomal genes interaction is suggested for the cause of transcriptional and phenotypic differences between reciprocal hybrids.

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

Hybridization, a common process in nature, can give rise to a vast reservoir of allelic variants. Combination of these allelic variants may result in novel patterns of gene action and is thought to contribute to heterosis. In this study, we analyzed genome-wide allele-specific gene expression (ASGE) in the super-hybrid rice variety Xieyou9308 using RNA sequencing technology (RNA-Seq). We identified 9325 reliable single nucleotide polymorphisms (SNPs) distributed throughout the genome. Nearly 68% of the identified polymorphisms were CT and GA SNPs between R9308 and Xieqingzao B, suggesting the existence of DNA methylation, a heritable epigenetic mark, in the parents and their F1 hybrid. Of 2793 identified transcripts with consistent allelic biases, only 480 (17%) showed significant allelic biases during tillering and/or heading stages, implying that trans effects may mediate most transcriptional differences in hybrid offspring. Approximately 67% and 62% of the 480 transcripts showed R9308 allelic expression biases at tillering and heading stages, respectively. Transcripts with higher levels of gene expression in R9308 also exhibited R9308 allelic biases in the hybrid. In addition, 125 transcripts were identified with significant allelic expression biases at both stages, of which 74% showed R9308 allelic expression biases. R9308 alleles may tend to preserve their characteristic states of activity in the hybrid and may play important roles in hybrid vigor at both stages. The allelic expression of 355 transcripts was highly stage-specific, with divergent allelic expression patterns observed at different developmental stages. Many transcripts associated with stress resistance were differently regulated in the F1 hybrid. The results of this study may provide valuable insights into molecular mechanisms of heterosis.

Genome anchored QTLs for biomass productivity in hybrid Populus grown under contrasting environments

Traits related to biomass production were analyzed for the presence of quantitative trait loci (QTLs) in a Populus trichocarpa × P. deltoides F(2) population. A genetic linkage map composed of 841 SSR, AFLP, and RAPD markers and phenotypic data from 310 progeny were used to identify genomic regions harboring biomass QTLs. Twelve intervals were identified, of which BM-1, BM-2, and BM-7 were identified in all three years for both height and diameter. One putative QTL, BM-7, and one suggestive QTL exhibited significant evidence of over-dominance in all three years for both traits. Conversely, QTLs BM-4 and BM-6 exhibited evidence of under-dominance in both environments for height and diameter. Seven of the nine QTLs were successfully anchored, and QTL peak positions were estimated for each one on the P. trichocarpa genome assembly using flanking SSR markers with known physical positions. Of the 3,031 genes located in genome-anchored QTL intervals, 1,892 had PFAM annotations. Of these, 1,313, representing 255 unique annotations, had at least one duplicate copy in a QTL interval identified on a separate scaffold. This observation suggests that some QTLs identified in this study may have shared the same ancestral sequence prior to the salicoid genome duplication in Populus.

A 34K SNP genotyping array for Populus trichocarpa: design, application to the study of natural populations and transferability to other Populus species

Genetic mapping of quantitative traits requires genotypic data for large numbers of markers in many individuals. For such studies, the use of large single nucleotide polymorphism (SNP) genotyping arrays still offers the most cost-effective solution. Herein we report on the design and performance of a SNP genotyping array for Populus trichocarpa (black cottonwood). This genotyping array was designed with SNPs pre-ascertained in 34 wild accessions covering most of the species latitudinal range. We adopted a candidate gene approach to the array design that resulted in the selection of 34 131 SNPs, the majority of which are located in, or within 2 kb of, 3543 candidate genes. A subset of the SNPs on the array (539) was selected based on patterns of variation among the SNP discovery accessions. We show that more than 95% of the loci produce high quality genotypes and that the genotyping error rate for these is likely below 2%. We demonstrate that even among small numbers of samples (n = 10) from local populations over 84% of loci are polymorphic. We also tested the applicability of the array to other species in the genus and found that the number of polymorphic loci decreases rapidly with genetic distance, with the largest numbers detected in other species in section Tacamahaca. Finally, we provide evidence for the utility of the array to address evolutionary questions such as intraspecific studies of genetic differentiation, species assignment and the detection of natural hybrids.

RNA-seq analysis of allele-specific expression, hybrid effects, and regulatory divergence in hybrids compared with their parents from natural populations

Hybridization is a prominent process among natural plant populations that can result in phenotypic novelty, heterosis, and changes in gene expression. The effects of intraspecific hybridization on F1 hybrid gene expression were investigated using parents from divergent, natural populations of Cirsium arvense, an invasive Compositae weed. Using an RNA-seq approach, the expression of 68,746 unigenes was quantified in parents and hybrids. The expression levels of 51% of transcripts differed between parents, a majority of which had less than 1.25× fold-changes. More unigenes had higher expression in the invasive parent (P1) than the noninvasive parent (P2). Of those that were divergently expressed between parents, 10% showed additive and 81% showed nonadditive (transgressive or dominant) modes of gene action in the hybrids. A majority of the dominant cases had P2-like expression patterns in the hybrids. Comparisons of allele-specific expression also enabled a survey of cis- and trans-regulatory effects. Cis- and trans-regulatory divergence was found at 70% and 68% of 62,281 informative single-nucleotide polymorphism sites, respectively. Of the 17% of sites exhibiting both cis- and trans-effects, a majority (70%) had antagonistic regulatory interactions (cis x trans); trans-divergence tended to drive higher expression of the P1 allele, whereas cis-divergence tended to increase P2 transcript abundance. Trans-effects correlated more highly than cis with parental expression divergence and accounted for a greater proportion of the regulatory divergence at sites with additive compared with nonadditive inheritance patterns. This study explores the nature of, and types of mechanisms underlying, expression changes that occur in upon intraspecific hybridization in natural populations.

The inheritance pattern of 24 nt siRNA clusters in arabidopsis hybrids is influenced by proximity to transposable elements

Hybrids often display increased size and growth, and thus are widely cultivated in agriculture and horticulture. Recent discoveries demonstrating the important regulatory roles of small RNAs have greatly improved our understanding of many basic biological questions, and could illuminate the molecular basis for the enhanced growth and size of hybrid plants. We profiled small RNAs by deep sequencing to characterize the inheritance patterns of small RNA levels in reciprocal hybrids of two Arabidopsis thaliana accessions, Columbia and Landsberg erecta. We find 24-nt siRNAs predominate among those small RNAs that are differentially expressed between the parents. Following hybridization, the transposable element (TE)-derived siRNAs are often inherited in an additive manner, whereas siRNAs associated with protein-coding genes are often down-regulated in hybrids to the levels observed for the parent with lower relative siRNA levels. Among the protein-coding genes that exhibit this pattern, genes that function in pathogen defense, abiotic stress tolerance, and secondary metabolism are significantly enriched. Small RNA clusters from protein-coding genes where a TE is present within one kilobase show a different predominant inheritance pattern (additive) from those that do not (low-parent dominance). Thus, down-regulation in the form of low-parent dominance is likely the default pattern of inheritance for genic siRNA, and a different inheritance mechanism for TE siRNA is suggested.

Complementation contributes to transcriptome complexity in maize (Zea mays L.) hybrids relative to their inbred parents

Typically, F₁-hybrids are more vigorous than their homozygous, genetically distinct parents, a phenomenon known as heterosis. In the present study, the transcriptomes of the reciprocal maize (Zea mays L.) hybrids B73×Mo17 and Mo17×B73 and their parental inbred lines B73 and Mo17 were surveyed in primary roots, early in the developmental manifestation of heterotic root traits. The application of statistical methods and a suitable experimental design established that 34,233 (i.e., 86%) of all high-confidence maize genes were expressed in at least one genotype. Nearly 70% of all expressed genes were differentially expressed between the two parents and 42%–55% of expressed genes were differentially expressed between one of the parents and one of the hybrids. In both hybrids, ∼10% of expressed genes exhibited nonadditive gene expression. Consistent with the dominance model (i.e., complementation) for heterosis, 1124 genes that were expressed in the hybrids were expressed in only one of the two parents. For 65 genes, it could be shown that this was a consequence of complementation of genomic presence/absence variation. For dozens of other genes, alleles from the inactive inbred were activated in the hybrid, presumably via interactions with regulatory factors from the active inbred. As a consequence of these types of complementation, both hybrids expressed more genes than did either parental inbred. Finally, in hybrids, ∼14% of expressed genes exhibited allele-specific expression (ASE) levels that differed significantly from the parental-inbred expression ratios, providing further evidence for interactions of regulatory factors from one parental genome with target genes from the other parental genome.

Genome-wide analysis of cis-regulatory divergence between species in the Arabidopsis genus

Cis-regulatory DNA has been suspected to play a preeminent role in adaptive evolution, but understanding the role of cis-regulatory mutations in gene expression divergence first requires an accurate analysis of the functional differences associated with these regions. We analyzed allele-specific expression (ASE) in leaf and floral tissues of F1 interspecific hybrids generated between the two closely related species Arabidopsis thaliana and A. lyrata with a whole-genome SNP (single nucleotide polymorphism) tiling array. We observed 2,205 genes showing ASE pattern in at least one tissue. Nearly 90% of genes displaying ASE preferentially expressed the allele of A. lyrata. Genome-wide comparison of sequence divergence revealed that genes displaying ASE had a higher ratio of nonsynonymous to synonymous substitutions in coding regions. We further observe that the epigenetic landscape of histone methylation in A. thaliana genome associate with ASE. The asymmetry in the direction of allele-specific expression suggests interspecific differences in the efficiency of gene silencing in F1 hybrids.

Extreme multiple reticulate origins of the Pteris cadieri complex (Pteridaceae)

The Pteris cadieri complex displays extensive morphological variation and seems to have originated through hybridization. However, the members of this complex reproduce by apogamy, which usually limits genetic variation. To evaluate the hypotheses of hybrid origins, the pattern of evolution in this species complex is reconstructed. Multiple methodologies were used. Diploids, triploids, and tetraploids were identified by chromosome counts and flow cytometry. Nuclear DNA markers (cytosolic phosphoglucose isomerase gene, PgiC) were used, together with chloroplast DNA markers (atpB-rbcL spacer and rbcL gene) to infer the biparental and maternal lineages of the Pteris cadieri complex. The three cpDNA haplotype groups and five PgiC alleles found in this study indicate that the evolution of the Pteris cadieri complex has been extremely reticulate. Up to 11 taxa belonging to eight morphs were identified. By comparing genetic variation in the Pteris cadieri in two independent areas, Hainan and Taiwan, we inferred that hybridization has occurred independently in different areas. Furthermore, we found evidence for phenological divergence (evergreen and deciduous) within Taiwan. We propose that the Pteris cadieri complex originated from different genetic lineages through multiple hybridizations in different geographical areas, leading to its present morphological diversity.

The inheritance pattern of 24 nt siRNA clusters in arabidopsis hybrids is influenced by proximity to transposable elements

Hybrids often display increased size and growth, and thus are widely cultivated in agriculture and horticulture. Recent discoveries demonstrating the important regulatory roles of small RNAs have greatly improved our understanding of many basic biological questions, and could illuminate the molecular basis for the enhanced growth and size of hybrid plants. We profiled small RNAs by deep sequencing to characterize the inheritance patterns of small RNA levels in reciprocal hybrids of two Arabidopsis thaliana accessions, Columbia and Landsberg erecta. We find 24-nt siRNAs predominate among those small RNAs that are differentially expressed between the parents. Following hybridization, the transposable element (TE)-derived siRNAs are often inherited in an additive manner, whereas siRNAs associated with protein-coding genes are often down-regulated in hybrids to the levels observed for the parent with lower relative siRNA levels. Among the protein-coding genes that exhibit this pattern, genes that function in pathogen defense, abiotic stress tolerance, and secondary metabolism are significantly enriched. Small RNA clusters from protein-coding genes where a TE is present within one kilobase show a different predominant inheritance pattern (additive) from those that do not (low-parent dominance). Thus, down-regulation in the form of low-parent dominance is likely the default pattern of inheritance for genic siRNA, and a different inheritance mechanism for TE siRNA is suggested.

Isolation and identification of an AP2/ERF factor that binds an allelic cis-element of rice gene LRK6

Allelic expression of the rice yield-related gene, leucine-rich receptor-like kinase 6 (LRK6), in the hybrid of 93-11 (Oryza sativa L. subsp. Indica var. 93-11) and Nipponbare (O. sativa L. subsp. Japonica var. Nipponbare) is determined by allelic promoter cis-elements. Using deletion analysis of the LRK6 promoter, we identified two distinct regions that might contribute to LRK6 expression. Sequence alignment revealed differences in these LRK6 promoter regions in 93-11 and Nipponbare. One of the segments, named differential sequence of LRK6 promoter 2 (DSLP2), contains potential transcription factor binding sites. Using a yeast one-hybrid assay, we isolated an ethylene-responsive factor (ERF) protein that binds to DSLP2. Sequence analysis and a GCC-box assay showed that the ERF gene, O. sativa ERF 3 (OsERF3), which belongs to ERF subfamily class II, has a conserved ERF domain and an ERF-associated amphiphilic repression repressor motif. We used an in vivo mutation assay to identify a new motif (5′-TAA(A)GT-3′) located in DSLP2, which interacts with OsERF3. These results suggest that OsERF3, an AP2 (APETALA 2 Gene)/ERF transcription factor, binds the LRK6 promoter at this new motif, which might cause differential expression of LRK6 in the 93-11/Nipponbare hybrid.

Heterosis

Heterosis refers to the phenomenon that progeny of diverse varieties of a species or crosses between species exhibit greater biomass, speed of development, and fertility than both parents. Various models have been posited to explain heterosis, including dominance, overdominance, and pseudo-overdominance. In this Perspective, we consider that it might be useful to the field to abandon these terms that by their nature constrain data interpretation and instead attempt a progression to a quantitative genetic framework involving interactions in hierarchical networks. While we do not provide a comprehensive model to explain the phenomenology of heterosis, we provide the details of what needs to be explained and a direction of pursuit that we feel should be fruitful.

Transcriptomic changes following recent natural hybridization and allopolyploidy in the salt marsh species Spartina x townsendii and Spartina anglica (Poaceae)

Allopolyploidy results from two events: the merger of divergent genomes and genome duplication. Both events have important functional consequences for the evolution and adaptation of newly formed allopolyploid species. In spite of the significant progress made in recent years, few studies have decoupled the effects of hybridization from genome duplication in the observed patterns of expression changes accompanying allopolyploidy in natural conditions. We used Agilent rice oligomicroarrays to explore gene expression changes following allopolyploidy in Spartina that includes a classic example of recent allopolyploid speciation: S. anglica formed during the 19th century following genome duplication of the hybrid S. x townsendii. Our data indicate important, but different, effects of hybridization and genome duplication in the expression patterns of the hybrid and allopolyploid. Deviation from parental additivity was most important following hybridization and was accompanied by maternal expression dominance, although transgressively expressed genes were also encountered. Maternal dominance was attenuated following genome duplication in S. anglica, but this species exhibits an increased number of transgressively overexpressed genes. These results reflect the decoupled effects of the 'genomic shock' following hybridization and genome redundancy on the genetic, epigenetic and regulatory mechanisms characterizing transcriptomic evolution in allopolyploids.

The application of restriction landmark genome scanning method for surveillance of non-mendelian inheritance in f(1) hybrids

We analyzed inheritance of DNA methylation in reciprocal F₁ hybrids (subsp. japonica cv. Nipponbare × subsp. indica cv. Kasalath) of rice (Oryza sativa L.) using restriction landmark genome scanning (RLGS), and detected differing RLGS spots between the parents and reciprocal F₁ hybrids. MspI/HpaII restriction sites in the DNA from these different spots were suspected to be heterozygously methylated in the Nipponbare parent. These spots segregated in F₁ plants, but did not segregate in selfed progeny of Nipponbare, showing non-Mendelian inheritance of the methylation status. As a result of RT-PCR and sequencing, a specific allele of the gene nearest to the methylated sites was expressed in reciprocal F₁ plants, showing evidence of biased allelic expression. These results show the applicability of RLGS for scanning of non-Mendelian inheritance of DNA methylation and biased allelic expression.

An interspecific plant hybrid shows novel changes in parental splice forms of genes for splicing factors

Interspecific hybridization plays an important role in plant adaptive evolution and speciation, and the process often results in phenotypic novelty. Hybrids can show changes in genome structure and gene expression compared with their parents including chromosomal rearrangments, changes in cytosine methylation, up- and downregulation of gene expression, and gene silencing. Alternative splicing (AS) is a fundamental aspect of the expression of many genes. However alternative splicing patterns have not been examined in multiple genes in an interspecific plant hybrid compared with its parents. Here we studied alternative splicing patterns in an interspecific Populus hybrid and its parents by assaying 40 genes using reverse transcription PCR. Most of the genes showed identical alternative splicing patterns between the parents and the hybrid. We found new alternative splicing variants present in the hybrid in two SR genes involved in the regulation of splicing and alternative splicing. The novel alternative splicing patterns included changes in donor and acceptor sites to create a new exon in one allele of PtRSZ22 in the hybrid and retention of an intron in both alleles of PtSR34a.1 in the hybrid, with effects on the function of the corresponding truncated proteins, if present. Our results suggest that novel alternative splicing patterns are present in a small percentage of genes in hybrids, but they could make a considerable impact on the expression of some genes. Changes in alternative splicing are likely to be an important component of the genetic changes that occur upon interspecific hybridization.

Asymmetric allele-specific expression in relation to developmental variation and drought stress in barley hybrids

In the present study, we analysed allele-specific expression (ASE) in the selfing species barley to assess the frequency of cis-acting regulatory variation and the effects of genetic background, developmental differences and drought stress on allelic expression levels. We measured ASE ratios in 30 genes putatively involved in stress responses in five hybrids and their reciprocals, namely Hordeum spontaneum 41-1/Alexis (HAl), Hordeum spontaneum 41-1/Arta (HAr), Sloop/WI3408 (SW), Tadmor/Sloop (TS) and Tadmor/WI3408 (TW). In order to detect cis-acting variation related to drought and developmental changes, the barley hybrids were grown under control and water-limited conditions, and leaf tissue was harvested at two developmental stages. The analysis demonstrated that more than half of the genes measured (63%) showed allelic differences in expression of up to 19-fold due to cis-regulatory variation in at least one cross by treatment/stage combination. Drought stress induced changes in allelic expression ratios, indicating differences between drought responsive cis-elements. In addition, ASE differences between developmental stages suggested the presence of cis-acting elements interacting with developmental cues. We were also able to demonstrate that the levels and frequency of allelic imbalance and hence differences in cis-regulatory elements are correlated with the genetic divergence between the parental lines, but may also arise as an adaptation to diverse habitats. Our findings suggest that cis-regulatory variation is a common phenomenon in barley, and may provide a molecular basis of transgression. Differential expression of near-isogenic members of the same gene family could potentially result in hybrid lines out performing their parents in terms of expression level, timing and response to developmental and environmental cues. Identification and targeted manipulation of cis-regulatory elements will assist in breeding improved crops with a better adaptation to changing environments.

Reciprocal silencing, transcriptional bias and functional divergence of homeologs in polyploid cotton (gossypium)

Polyploidy is an important force in the evolution of flowering plants. Genomic merger and doubling induce an extensive array of genomic effects, including immediate and long-term alterations in the expression of duplicate genes ("homeologs"). Here we employed a novel high-resolution, genome-specific, mass-spectrometry technology and a well-established phylogenetic framework to investigate relative expression levels of each homeolog for 63 gene pairs in 24 tissues in naturally occurring allopolyploid cotton (Gossypium L.), a synthetic allopolyploid of the same genomic composition, and models of the diploid progenitor species. Results from a total of 2177 successful expression assays permitted us to determine the extent of expression evolution accompanying genomic merger of divergent diploid parents, genome doubling, and genomic coevolution in a common nucleus subsequent to polyploid formation. We demonstrate that 40% of homeologs are transcriptionally biased in at least one stage of cotton development, that genome merger per se has a large effect on relative expression of homeologs, and that the majority of these alterations are caused by cis-regulatory divergence between the diploid progenitors. We describe the scope of transcriptional subfunctionalization and 15 cases of probable neofunctionalization among 8 tissues. To our knowledge, this study represents the first characterization of transcriptional neofunctionalization in an allopolyploid. These results provide a novel temporal perspective on expression evolution of duplicate genomes and add to our understanding of the importance of polyploidy in plants.

Genomics of growth traits in forest trees

Growth traits in trees are fundamental components of adaptation in a forest ecosystem and of productivity in planted forests. A number of processes determine tree growth, which are controlled by genetic and epigenetic factors that respond dynamically to environmental signals throughout centuries. Advances in genomics have allowed an increased comprehension of the complex mechanisms of tree growth and adaptation. Yet, the application of genomics to improving forest productivity and sustainability still entails capturing a large proportion of the total genetic variation controlling the component traits. Nonetheless, genetics and genomics are unifying disciplines that will serve well to dissect the variables and mechanisms of tree growth and development.

Understanding quantitative genetics in the systems biology era

Biology is now entering the new era of systems biology and exerting a growing influence on the future development of various disciplines within life sciences. In early classical and molecular periods of Biology, the theoretical frames of classical and molecular quantitative genetics have been systematically established, respectively. With the new advent of systems biology, there is occurring a paradigm shift in the field of quantitative genetics. Where and how the quantitative genetics would develop after having undergone its classical and molecular periods? This is a difficult question to answer exactly. In this perspective article, the major effort was made to discuss the possible development of quantitative genetics in the systems biology era, and for which there is a high potentiality to develop towards "systems quantitative genetics". In our opinion, the systems quantitative genetics can be defined as a new discipline to address the generalized genetic laws of bioalleles controlling the heritable phenotypes of complex traits following a new dynamic network model. Other issues from quantitative genetic perspective relating to the genetical genomics, the updates of network model, and the future research prospects were also discussed.

Evidence for a positive regulatory role of strawberry (Fragaria x ananassa) Fa WRKY1 and Arabidopsis At WRKY75 proteins in resistance

Knowledge of the molecular basis of plant resistance to pathogens in species other than Arabidopsis is limited. The function of Fa WRKY1, the first WRKY gene isolated from strawberry (Fragaria x ananassa), an important agronomical fruit crop, has been investigated here. Fa WRKY1 encodes a IIc WRKY transcription factor and is up-regulated in strawberry following Colletotrichum acutatum infection, treatments with elicitors, and wounding. Its Arabidopsis sequence homologue, At WRKY75, has been described as playing a role in regulating phosphate starvation responses. However, using T-DNA insertion mutants, a role for the At WRKY75 and Fa WRKY1 in the activation of basal and R-mediated resistance in Arabidopsis is demonstrated. At wrky75 mutants are more susceptible to virulent and avirulent isolates of Pseudomonas syringae. Overexpression of Fa WRKY1 in At wrky75 mutant and wild type reverts the enhanced susceptible phenotype of the mutant, and even increases resistance to avirulent strains of P. syringae. The resistance phenotype is uncoupled to PATHOGENESIS-RELATED (PR) gene expression, but it is associated with a strong oxidative burst and glutathione-S-transferase (GST) induction. Taken together, these results indicate that At WRKY75 and Fa WRKY1 act as positive regulators of defence during compatible and incompatible interactions in Arabidopsis and, very likely, Fa WRKY1 is an important element mediating defence responses to C. acutatum in strawberry. Moreover, these results provide evidence that Arabidopsis can be a useful model for functional studies in Rosacea species like strawberry.

Cis-regulatory variation is typically polyallelic in Drosophila

Gene expression levels vary heritably, with approximately 25-35% of the loci affecting expression acting in cis. We characterized standing cis-regulatory variation among 16 wild-derived strains of Drosophila melanogaster. Our experiment's robust biological and technical replication enabled precise estimates of variation in allelic expression on a high-throughput SNP genotyping platform. We observed concordant, significant differential allelic expression (DAE) in 7/10 genes queried with multiple SNPs, and every member of a set of eight additional, one-assay genes suggest significant DAE. Four of the high-confidence, multiple-assay genes harbor three or more statistically distinguishable allelic classes, often at intermediate frequency. Numerous intermediate-frequency, detectable regulatory polymorphisms cast doubt on a model in which cis-acting variation is a product of deleterious mutations of large effect. Comparing our data to predictions of population genetics theory using coalescent simulations, we estimate that a typical gene harbors 7-15 cis-regulatory sites (nucleotides) at which a selectively neutral mutation would elicit an observable expression phenotype. If standing cis-regulatory variation is actually slightly deleterious, the true mutational target size is larger.