Transcriptome variations in hybrids of wild emmer wheat (Triticum turgidum ssp. dicoccoides)

BACKGROUND

Wild emmer wheat is a great candidate to revitalize domesticated wheat genetic diversity. Recent years have seen intensive investigation into the evolution and domestication of wild emmer wheat, including whole-genome DNA and transcriptome sequencing. However, the impact of intraspecific hybridization on the transcriptome of wild emmer wheat has been poorly studied. In this study, we assessed changes in methylation patterns and transcriptomic variations in two accessions of wild emmer wheat collected from two marginal populations, Mt. Hermon and Mt. Amasa, and in their stable F4 hybrid.

RESULTS

Methylation-Sensitive Amplified Polymorphism (MSAP) detected significant cytosine demethylation in F4 hybrids vs. parental lines, suggesting potential transcriptome variation. After a detailed analysis, we examined nine RNA-Seq samples, which included three biological replicates from the F4 hybrid and its parental lines. RNA-Seq databases contained approximately 200 million reads, with each library consisting of 15 to 25 million reads. There are a total of 62,490 well-annotated genes in these databases, with 6,602 genes showing differential expression between F4 hybrid and parental lines Mt. Hermon and Mt. Amasa. The differentially expressed genes were classified into four main categories based on their expression patterns. Gene ontology (GO) analysis revealed that differentially expressed genes are associated with DNA/RNA metabolism, photosynthesis, stress response, phosphorylation and developmental processes.

CONCLUSION

This study highlights the significant transcriptomic changes resulting from intraspecific hybridization within natural plant populations, which might aid the nascent hybrid in adapting to various environmental conditions.

OUP accepted manuscript

Interspecific hybridization is often seen as a genomic stress that may lead to new gene expression patterns and deregulation of transposable elements (TEs). The understanding of expression changes in hybrids compared with parental species is essential to disentangle their putative role in speciation processes. However, to date we ignore the detailed mechanisms involved in genomic deregulation in hybrids. We studied the ovarian transcriptome and epigenome of the Drosophila buzzatii and Drosophila koepferae species together with their F₁ hybrid females. We found a trend toward underexpression of genes and TE families in hybrids. The epigenome in hybrids was highly similar to the parental epigenomes and showed intermediate histone enrichments between parental species in most cases. Differential gene expression in hybrids was often associated only with changes in H3K4me3 enrichments, whereas differential TE family expression in hybrids may be associated with changes in H3K4me3, H3K9me3, or H3K27me3 enrichments. We identified specific genes and TE families, which their differential expression in comparison with the parental species was explained by their differential chromatin mark combination enrichment. Finally, cis–trans compensatory regulation could also contribute in some way to the hybrid deregulation. This work provides the first study of histone content in Drosophila interspecific hybrids and their effect on gene and TE expression deregulation.

Transcriptome and DNA methylome analyses provide insight into the heterosis in flag leaf of inter-subspecific hybrid rice

Flag leaf heterosis of inter-subspecific hybrid rice is suggested to be related to leaf area, gene expression pattern and allele-specific expression, putatively related to DNA methylation differences between the hybrid and its parents. Inter-subspecific hybrid rice combinations of indica × japonica have great potential to broaden genetic diversity and enhance the heterosis. The genetic and epigenetic molecular mechanism of its heterosis is not completely understood. Here, the dissection of gene expression and epigenetic regulation of an elite inter-subspecific hybrid rice were reported. In the hybrid, plant height, flag leaf area and Pn showed significant heterosis at the heading stage. Chloroplast-related differentially expressed genes (DEGs) and 530 allele-specific expression genes in hybrid were identified. Analysis of the genome-wide distribution of DNA methylation (5-methylcytosine, 5mC) and its association with transcription showed that there were variant DNA methylation maps and that the regulation of gene expression levels was negatively regulated by DNA methylation in the inter-subspecific hybrid rice. Differentially methylated DEGs were significantly enriched in photosynthetic functions. Moreover, distinct 5mC sequence contexts and distinct functional elements (promoter/gene body) may have different influences on heterosis related genes. The data identified heterosis related molecular mechanisms in inter-subspecific hybrid rice and suggested that epigenetic changes could extensively influence the flag leaf gene expression and heterosis.

Inheritance of gene expression throughout fruit development in chili pepper

Gene expression is the primary molecular phenotype and can be estimated in specific organs or tissues at particular times. Here we analyzed genome-wide inheritance of gene expression in fruits of chili pepper (Capsicum annuum L.) in reciprocal crosses between a domesticated and a wild accession, estimating this parameter during fruit development. We defined a general hierarchical schema to classify gene expression inheritance which can be employed for any quantitative trait. We found that inheritance of gene expression is affected by both, the time of fruit development as well as the direction of the cross, and propose that such variations could be common in many developmental processes. We conclude that classification of inheritance patterns is important to have a better understanding of the mechanisms underlying gene expression regulation, and demonstrate that sets of genes with specific inheritance pattern at particular times of fruit development are enriched in different biological processes, molecular functions and cell components. All curated data and functions for analysis and visualization are publicly available as an R package.

Imprints of independent allopolyploid formations on patterns of gene expression in two sibling yarrow species (Achillea, Asteraceae)

BACKGROUND

Polyploid species often originate recurrently. While this is well known, there is little information on the extent to which distinct allotetraploid species formed from the same parent species differ in gene expression. The tetraploid yarrow species Achillea alpina and A. wilsoniana arose independently from allopolyploidization between diploid A. acuminata and A. asiatica. The genetics and geography of these origins are clear from previous studies, providing a solid basis for comparing gene expression patterns of sibling allopolyploid species that arose independently.

RESULTS

We conducted comparative RNA-sequencing analyses on the two Achillea tetraploid species and their diploid progenitors to evaluate: 1) species-specific gene expression and coexpression across the four species; 2) patterns of inheritance of parental gene expression; 3) parental contributions to gene expression in the allotetraploid species, and homeolog expression bias. Diploid A. asiatica showed a higher contribution than diploid A. acuminata to the transcriptomes of both tetraploids and also greater homeolog bias in these transcriptomes, possibly reflecting a maternal effect. Comparing expressed genes in the two allotetraploids, we found expression of ca. 30% genes were species-specific in each, which were most enriched for GO terms pertaining to "defense response". Despite species-specific and differentially expressed genes between the two allotetraploids, they display similar transcriptome changes in comparison to their diploid progenitors.

CONCLUSION

Two independently originated Achillea allotetraploid species exhibited difference in gene expression, some of which must be related to differential adaptation during their post-speciation evolution. On the other hand, they showed similar expression profiles when compared to their progenitors. This similarity might be expected when pairs of merged diploid genomes in tetraploids are similar, as is the case in these two particular allotetraploids.

Rearing environment affects the genetic architecture and plasticity of DNA methylation in Chinook salmon

Genetic architecture and phenotypic plasticity are important considerations when studying trait variation within and among populations. Since environmental change can induce shifts in the genetic architecture and plasticity of traits, it is important to consider both genetic and environmental sources of phenotypic variation. While there is overwhelming evidence for environmental effects on phenotype, the underlying mechanisms are less clear. Variation in DNA methylation is a potential mechanism mediating environmental effects on phenotype due to its sensitivity to environmental stimuli, transgenerational inheritance, and influences on transcription. To characterize the effect of environment on methylation, we created two 6 × 6 (North Carolina II) Chinook salmon breeding crosses and reared the offspring in two environments: uniform hatchery tanks and seminatural stream channels. We sampled the fish twice during development, at the alevin (larval) and fry (juvenile) stages. We measured DNA methylation at 13 genes using a PCR-based bisulfite sequencing protocol. The genetic architecture of DNA methylation differed between rearing environments, with greater additive and nonadditive genetic variance in hatchery fish and greater maternal effects in seminatural channel fish, though gene-specific variation was evident. We observed plasticity in methylation across all assayed genes, as well as gene-specific effects at two genes in alevin and six genes in fry, indicating developmental stage-specific effects of rearing environment on methylation. Characterizing genetic and environmental influences on methylation is critical for future studies on DNA methylation as a potential mechanism for acclimation and adaptation.

Maternal effects in gene expression of interspecific coral hybrids

Maternal effects have been well documented for offspring morphology and life history traits in plants and terrestrial animals, yet little is known about maternal effects in corals. Further, few studies have explored maternal effects in gene expression. In a previous study, F1 interspecific hybrid and purebred larvae of the coral species Acropora tenuis and Acropora loripes were settled and exposed to ambient or elevated temperature and pCO₂ conditions for 7 months. At this stage, the hybrid coral recruits from both ocean conditions exhibited strong maternal effects in several fitness traits. We conducted RNA-sequencing on these corals and showed that gene expression of the hybrid Acropora also exhibited clear maternal effects. Only 40 genes were differentially expressed between hybrids and their maternal progenitor. In contrast, ~2000 differentially expressed genes were observed between hybrids and their paternal progenitors, and between the reciprocal F1 hybrids. These results indicate that maternal effects in coral gene expression can be long-lasting. Unlike findings from most short-term stress experiments in corals, no genes were differentially expressed in the hybrid nor purebred offspring after seven months of exposure to elevated temperature and pCO₂ conditions.

DNA Methylation Profiles Suggest Intergenerational Transfer of Maternal Effects

The view of maternal effects (nongenetic maternal environmental influence on offspring phenotype) has changed from one of distracting complications in evolutionary genetics to an important evolutionary mechanism for improving offspring fitness. Recent studies have shown that maternal effects act as an adaptive mechanism to prepare offspring for stressful environments. Although research into the magnitude of maternal effects is abundant, the molecular mechanisms of maternal influences on offspring phenotypic variation are not fully understood. Despite recent work identifying DNA methylation as a potential mechanism of nongenetic inheritance, currently proposed links between DNA methylation and parental effects are indirect and primarily involve genomic imprinting. We combined a factorial breeding design and gene-targeted sequencing methods to assess inheritance of methylation during early life stages at 14 genes involved in growth, development, metabolism, stress response, and immune function of Chinook salmon (Oncorhynchus tshawytscha). We found little evidence for additive or nonadditive genetic effects acting on methylation levels during early development; however, we detected significant maternal effects. Consistent with conventional maternal effect data, maternal effects on methylation declined through development and were replaced with nonadditive effects when offspring began exogenous feeding. We mapped methylation at individual CpG sites across the selected candidate genes to test for variation in site-specific methylation profiles and found significant maternal effects at selected CpG sites that also declined with development stage. While intergenerational inheritance of methylated DNA is controversial, we show that CpG-specific methylation may function as an underlying molecular mechanism for maternal effects, with important implications for offspring fitness.

Some Like It Hot: Maternal-Switching With Climate Change Modifies Formation of Invasive Spartina Hybrids

Climate change can induce temporary, spatial or behavioral changes in species, so that only some species can adapt to the new climatic conditions. In the case of invasive species, it is expected that they will be promoted in a context of global change, given their high tolerance to environmental factors and phenotypic plasticity. Once in the invaded range, these species can hybridize with native species thus introducing their genotype in the native biota. However, the effects that climate change will have on this process of invasion by hybridization remain unclear. We evaluated the historical establishment of the reciprocal hybrids between the native Spartina maritima and the invasive S. densiflora in the Gulf of Cadiz (SW Iberian Peninsula) and we related it to climatic changes during the period 1955-2017. Our results showed that, according to their dating based on their rate of lateral expansion rates, the establishment of S. maritima × densiflora and S. densiflora × maritima in the Gulf of Cadiz has occurred in the last two centuries and has been related to changes in air temperature and rainfall during the flowering periods of their parental species, with antagonist impacts on both hybrids. Thus, the hybrid S. densiflora × maritima has been established in years with mild ends of spring and beginning of summer when the flowering of S. maritima lengthened and its pollen production was higher, and it coincided with the beginning of the flowering period of S. densiflora. Moreover, the establishment of this hybrid was related to higher spring/summer rainfalls, probably due to the reduction in salinity in middle marshes. However, the hybrid S. maritima × densiflora, was established mainly in warmer spring/summers in which the proportion of pollen:ovule of S. maritima was reduced favoring its pollination by S. densiflora. As a consequence of the promotion of S. maritima × densiflora with climate change, the native and endangered species S. maritima would be threatened, as both taxa share the same habitat and the hybrid shows a remarkably higher competitive potential.

Asymmetric expression patterns reveal a strong maternal effect and dosage compensation in polyploid hybrid fish

Background

Hybridization and polyploidization are regarded as the major driving forces in plant speciation, diversification, and ecological adaptation. Our knowledge regarding the mechanisms of duplicated-gene regulation following genomic merging or doubling is primarily derived from plants and is sparse for vertebrates.

Results

We successfully obtained an F1 generation (including allodiploid hybrids and triploid hybrids) from female Megalobrama amblycephala Yih (BSB, 2n = 48) × male Xenocypri davidi Bleeker (YB, 2n = 48). The duplicated-gene expression patterns of the two types of hybrids were explored using RNA-Seq data. In total, 5.44 × 10⁸ (69.32 GB) clean reads and 499,631 assembled unigenes were obtained from the testis transcriptomes. The sequence similarity analysis of 4265 orthologs revealed that the merged genomes were dominantly expressed in different ploidy hybrids. The differentially expressed genes in the two types of hybrids were asymmetric compared with those in both parents. Furthermore, the genome-wide expression level dominance (ELD) was biased toward the maternal BSB genome in both the allodiploid and triploid hybrids. In addition, the dosage-compensation mechanisms that reduced the triploid expression levels to the diploid state were determined in the triploid hybrids.

Conclusions

Our results indicate that divergent genomes undergo strong interactions and domination in allopolyploid offspring. Genomic merger has a greater effect on the gene-expression patterns than genomic doubling. The various expression mechanisms (including maternal effect and dosage compensation) in different ploidy hybrids suggest that the initial genomic merger and doubling play important roles in polyploidy adaptation and evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4883-7) contains supplementary material, which is available to authorized users.

An Arabidopsis introgression zone studied at high spatio-temporal resolution: interglacial and multiple genetic contact exemplified using whole nuclear and plastid genomes

Background

Gene flow between species, across ploidal levels, and even between evolutionary lineages is a common phenomenon in the genus Arabidopsis. However, apart from two genetically fully stabilized allotetraploid species that have been investigated in detail, the extent and temporal dynamics of hybridization are not well understood. An introgression zone, with tetraploid A. arenosa introgressing into A. lyrata subsp. petraea in the Eastern Austrian Forealps and subsequent expansion towards pannonical lowlands, was described previously based on morphological observations as well as molecular data using microsatellite and plastid DNA markers. Here we investigate the spatio-temporal context of this suture zone, making use of the potential of next-generation sequencing and whole-genome data. By utilizing a combination of nuclear and plastid genomic data, the extent, direction and temporal dynamics of gene flow are elucidated in detail and Late Pleistocene evolutionary processes are resolved.

Results

Analysis of nuclear genomic data significantly recognizes the clinal structure of the introgression zone, but also reveals that hybridization and introgression is more common and substantial than previously thought. Also tetraploid A. lyrata and A. arenosa subsp. borbasii from outside the previously defined suture zone show genomic signals of past introgression. A. lyrata is shown to serve usually as the maternal parent in these hybridizations, but one exception is identified from plastome-based phylogenetic reconstruction. Using plastid phylogenomics with secondary time calibration, the origin of A. lyrata and A. arenosa lineages is pre-dating the last three glaciation complexes (approx. 550,000 years ago). Hybridization and introgression followed during the last two glacial-interglacial periods (since approx. 300,000 years ago) with later secondary contact at the northern and southern border of the introgression zone during the Holocene.

Conclusions

Footprints of adaptive introgression in the Northeastern Forealps are older than expected and predate the Last Glaciation Maximum. This correlates well with high genetic diversity found within areas that served as refuge area multiple times. Our data also provide some first hints that early introgressed and presumably preadapted populations account for successful and rapid postglacial re-colonization and range expansion.

Electronic supplementary material

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

Genetic architecture of gene transcription in two Atlantic salmon (Salmo salar) populations

Gene expression regulation has an important role in short-term acclimation and long-term adaptation to changing environments. However, the genetic architecture of gene expression has received much less attention than that of traditional phenotypic traits. In this study, we used a 5 × 5 full-factorial breeding design within each of two Atlantic salmon (Salmo salar) populations to characterize the genetic architecture of gene transcription. The two populations (LaHave and Sebago) are being used for reintroduction efforts into Lake Ontario, Canada. We used high-throughput quantitative real-time PCR to measure gene transcription levels for 22 genes in muscle tissue of Atlantic salmon fry. We tested for population differences in gene transcription and partitioned the transcription variance into additive genetic, non-additive genetic and maternal effects within each population. Interestingly, average additive genetic effects for gene transcription were smaller than those reported for traditional phenotypic traits in salmonids, suggesting that the evolutionary potential of gene transcription is lower than that of traditional traits. Contrary to expectations for early life stage traits, maternal effects were small. In general, the LaHave population had higher additive genetic effects for gene transcription than the Sebago population had, indicating that the LaHave fish have a higher adaptive potential to respond to the novel selection pressures associated with reintroduction into a novel environment. This study highlights not only the profound variation in gene transcription possible among salmonid populations but also the among-population variation in the underlying genetic architecture of such traits.

Field Guide to Plant Model Systems

For the past several decades, advances in plant development, physiology, cell biology, and genetics have relied heavily on the model (or reference) plant Arabidopsis thaliana. Arabidopsis resembles other plants, including crop plants, in many but by no means all respects. Study of Arabidopsis alone provides little information on the evolutionary history of plants, evolutionary differences between species, plants that survive in different environments, or plants that access nutrients and photosynthesize differently. Empowered by the availability of large-scale sequencing and new technologies for investigating gene function, many new plant models are being proposed and studied.

Genomic Imprinting in the Endosperm Is Systematically Perturbed in Abortive Hybrid Tomato Seeds

Hybrid seed failure represents an important postzygotic barrier to interbreeding among species of wild tomatoes (Solanum section Lycopersicon) and other flowering plants. We studied genome-wide changes associated with hybrid seed abortion in the closely related Solanum peruvianum and S. chilense where hybrid crosses yield high proportions of inviable seeds due to endosperm failure and arrested embryo development. Based on differences of seed size in reciprocal hybrid crosses and developmental evidence implicating endosperm failure, we hypothesized that perturbed genomic imprinting is involved in this strong postzygotic barrier. Consequently, we surveyed the transcriptomes of developing endosperms from intra- and inter-specific crosses using tissues isolated by laser-assisted microdissection. We implemented a novel approach to estimate parent-of-origin–specific expression using both homozygous and heterozygous nucleotide differences between parental individuals and identified candidate imprinted genes. Importantly, we uncovered systematic shifts of “normal” (intraspecific) maternal:paternal transcript proportions in hybrid endosperms; the average maternal proportion of gene expression increased in both crossing directions but was stronger with S. peruvianum in the maternal role. These genome-wide shifts almost entirely eliminated paternally expressed imprinted genes in S. peruvianum hybrid endosperm but also affected maternally expressed imprinted genes and all other assessed genes. These profound, systematic changes in parental expression proportions suggest that core processes of transcriptional regulation are functionally compromised in hybrid endosperm and contribute to hybrid seed failure.

Recent mating-system evolution in Eichhornia is accompanied by cis-regulatory divergence

The evolution of predominant self-fertilization from cross-fertilization in plants is accompanied by diverse changes to morphology, ecology and genetics, some of which likely result from regulatory changes in gene expression. We examined changes in gene expression during early stages in the transition to selfing in populations of animal-pollinated Eichhornia paniculata with contrasting mating patterns. We crossed plants from outcrossing and selfing populations and tested for the presence of allele-specific expression (ASE) in floral buds and leaf tissue of F1 offspring, indicative of cis-regulatory changes. We identified 1365 genes exhibiting ASE in floral buds and leaf tissue. These genes preferentially expressed alleles from outcrossing parents. Moreover, we found evidence that genes exhibiting ASE had a greater nonsynonymous diversity compared to synonymous diversity in the selfing parents. Our results suggest that the transition from outcrossing to high rates of self-fertilization may have the potential to shape the cis-regulatory genomic landscape of angiosperm species, but that the changes in ASE may be moderate, particularly during the early stages of this transition.