Cytonuclear interactions can favor the evolution of genomic imprinting

Comparative Hox genes expression within the dimorphic annelid Streblospio benedicti reveals patterning variation during development

Hox genes are transcriptional regulators that elicit cell positional identity along the anterior-posterior region of the body plan across different lineages of Metazoan. Comparison of Hox gene expression across distinct species reveals their evolutionary conservation; however, their gains and losses in different lineages can correlate with body plan modifications and morphological novelty. We compare the expression of 11 Hox genes found within Streblospio benedicti, a marine annelid that produces two types of offspring with distinct developmental and morphological features. For these two distinct larval types, we compare Hox gene expression through ontogeny using hybridization chain reaction (HCR) probes for in situ hybridization and RNA-seq data. We find that Hox gene expression patterning for both types is typically similar at equivalent developmental stages. However, some Hox genes have spatial or temporal differences between the larval types that are associated with morphological and life-history differences. This is the first comparison of developmental divergence in Hox gene expression within a single species and these changes reveal how body plan differences may arise in larval evolution.

Comparative Hox genes expression within the dimorphic annelid Streblospio benedicti reveals patterning variation during development

Hox genes are transcriptional regulators that elicit cell positional identity along the anterior-posterior region of the body plan across different lineages of Metazoan. Comparison of Hox gene expression across distinct species reveals their evolutionary conservation, however their gains and losses in different lineages can correlate with body plan modifications and morphological novelty. We compare the expression of eleven Hox genes found within Streblospio benedicti, a marine annelid that produces two types of offspring with distinct developmental and morphological features. For these two distinct larval types, we compare Hox gene expression through ontogeny using HCR (hybridization chain reaction) probes for in-situ hybridization and RNA-seq data. We find that Hox gene expression patterning for both types is typically similar at equivalent developmental stages. However, some Hox genes have spatial or temporal differences between the larval types that are associated with morphological and life-history differences. This is the first comparison of developmental divergence in Hox genes expression within a single species and these changes reveal how body plan differences may arise in larval evolution.

Plant biochemical genetics in the multiomics era

Our understanding of plant biology has been revolutionized by modern genetics and biochemistry. However, biochemical genetics can be traced back to the foundation of Mendelian genetics; indeed, one of Mendel's milestone discoveries of seven characteristics of pea plants later came to be ascribed to a mutation in a starch branching enzyme. Here, we review both current and historical strategies for the elucidation of plant metabolic pathways and the genes that encode their component enzymes and regulators. We use this historical review to discuss a range of classical genetic phenomena including epistasis, canalization, and heterosis as viewed through the lens of contemporary high-throughput data obtained via the array of approaches currently adopted in multiomics studies.

Parent-of-origin effects, allele-specific expression, genomic imprinting and paternal manipulation in social insects

Haplo-diploidy and the relatedness asymmetries it generates mean that social insects are prime candidates for the evolution of genomic imprinting. In single-mating social insect species, some genes may be selected to evolve genomic mechanisms that enhance reproduction by workers when they are inherited from a female. This situation reverses in multiple mating species, where genes inherited from fathers can be under selection to enhance the reproductive success of daughters. Reciprocal crosses between subspecies of honeybees have shown strong parent-of-origin effects on worker reproductive phenotypes, and this could be evidence of such genomic imprinting affecting genes related to worker reproduction. It is also possible that social insect fathers directly affect gene expression in their daughters, for example, by placing small interfering RNA molecules in semen. Gene expression studies have repeatedly found evidence of parent-specific gene expression in social insects, but it is unclear at this time whether this arises from genomic imprinting, paternal manipulation, an artefact of cyto-nuclear interactions, or all of these. This article is part of the theme issue 'How does epigenetics influence the course of evolution?'

Detection of subgenome bias using an anchored syntenic approach in Eleusine coracana (finger millet)

Background

Finger millet (Eleusine coracana 2n = 4x = 36) is a hardy, nutraceutical, climate change tolerant, orphan crop that is consumed throughout eastern Africa and India. Its genome has been sequenced multiple times, but A and B subgenomes could not be separated because no published genome for E. indica existed. The classification of A and B subgenomes is important for understanding the evolution of this crop and provide a means to improve current and future breeding programs.

Results

We produced subgenome calls for 704 syntenic blocks and inferred A or B subgenomic identity for 59,377 genes 81% of the annotated genes. Phylogenetic analysis of a super matrix containing 455 genes shows high support for A and B divergence within the Eleusine genus. Synonymous substitution rates between A and B genes support A and B calls. The repetitive content on highly supported B contigs is higher than that on similar A contigs. Analysis of syntenic singletons showed evidence of biased fractionation showed a pattern of A genome dominance, with 61% A, 37% B and 1% unassigned, and was further supported by the pattern of loss observed among cyto-nuclear interacting genes.

Conclusion

The evidence of individual gene calls within each syntenic block, provides a powerful tool for inference for subgenome classification. Our results show the utility of a draft genome in resolving A and B subgenomes calls, primarily it allows for the proper polarization of A and B syntenic blocks. There have been multiple calls for the use of phylogenetic inference in subgenome classification, our use of synteny is a practical application in a system that has only one parental genome available.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07447-y.

Coevolution in Hybrid Genomes: Nuclear-Encoded Rubisco Small Subunits and Their Plastid-Targeting Translocons Accompanying Sequential Allopolyploidy Events in Triticum

The Triticum/Aegilops complex includes hybrid species resulting from homoploid hybrid speciation and allopolyploid speciation. Sequential allotetra- and allohexaploidy events presumably result in two challenges for the hybrids, which involve 1) cytonuclear stoichiometric disruptions caused by combining two diverged nuclear genomes with the maternal inheritance of the cytoplasmic organellar donor; and 2) incompatibility of chimeric protein complexes with diverged subunits from nuclear and cytoplasmic genomes. Here, we describe coevolution of nuclear rbcS genes encoding the small subunits of Rubisco (ribulose 1,5-bisphosphate carboxylase/oxygenase) and nuclear genes encoding plastid translocons, which mediate recognition and translocation of nuclear-encoded proteins into plastids, in allopolyploid wheat species. We demonstrate that intergenomic paternal-to-maternal gene conversion specifically occurred in the genic region of the homoeologous rbcS3 gene from the D-genome progenitor of wheat (abbreviated as rbcS3D) such that it encodes a maternal-like or B-subgenome-like SSU3D transit peptide in allohexaploid wheat but not in allotetraploid wheat. Divergent and limited interaction between SSU3D and the D-subgenomic TOC90D translocon subunit is implicated to underpin SSU3D targeting into the chloroplast of hexaploid wheat. This implicates early selection favoring individuals harboring optimal maternal-like organellar SSU3D targeting in hexaploid wheat. These data represent a novel dimension of cytonuclear evolution mediated by organellar targeting and transportation of nuclear proteins.

Paternally-biased gene expression follows kin-selected predictions in female honey bee embryos

The Kinship Theory of Genomic Imprinting (KTGI) posits that, in species where females mate with multiple males, there is selection for a male to enhance the reproductive success of his offspring at the expense of other males and his mating partner. Reciprocal crosses between honey bee subspecies show parent-of-origin effects for reproductive traits, suggesting that males modify the expression of genes related to female function in their female offspring. This effect is likely to be greater in the Cape honey bee (Apis mellifera capensis), because a male's daughters have the unique ability to produce female offspring that can develop into reproductive workers or the next queen without mating. We generated reciprocal crosses between Capensis and another subspecies and used RNA-seq to identify transcripts that are over- or underexpressed in the embryos, depending on the parental origin of the gene. As predicted, 21 genes showed expression bias towards the Capensis father's allele in colonies with a Capensis father, with no such bias in the reciprocal cross. A further six genes showed a consistent bias towards expression of the father's allele across all eight colonies examined, regardless of the direction of the cross. Consistent with predictions of the KTGI, six of the 21 genes are associated with female reproduction. No gene consistently showed overexpression of the maternal allele.

The Legacy of Sexual Ancestors in Phenotypic Variability, Gene Expression, and Homoeolog Regulation of Asexual Hybrids and Polyploids

Hybridization and polyploidization are important evolutionary processes whose impacts range from the alteration of gene expression and phenotypic variation to the triggering of asexual reproduction. We investigated fishes of the Cobitis taenia-elongatoides hybrid complex, which allowed us to disentangle the direct effects of both processes, due to the co-occurrence of parental species with their diploid and triploid hybrids. Employing morphological, ecological, and RNAseq approaches, we investigated the molecular determinants of hybrid and polyploid forms. In contrast with other studies, hybridization and polyploidy induced relatively very little transgressivity. Instead, Cobitis hybrids appeared intermediate with a clear effect of genomic dosing when triploids expressed higher similarity to the parent contributing two genome sets. This dosage effect was symmetric in the germline (oocyte gene expression), interestingly though, we observed an overall bias toward C. taenia in somatic tissues and traits. At the level of individual genes, expression-level dominance vastly prevailed over additivity or transgressivity. Also, trans-regulation of gene expression was less efficient in diploid hybrids than in triploids, where the expression modulation of homoeologs derived from the "haploid" parent was stronger than those derived from the "diploid" parent. Our findings suggest that the apparent intermediacy of hybrid phenotypes results from the combination of individual genes with dominant expression rather than from simple additivity. The efficiency of cross-talk between trans-regulatory elements further appears dosage dependent. Important effects of polyploidization may thus stem from changes in relative concentrations of trans-regulatory elements and their binding sites between hybridizing genomes. Links between gene regulation and asexuality are discussed.

Allele-specific expression variation at different ploidy levels in Squalius alburnoides

Allopolyploid plants are long known to be subject to a homoeolog expression bias of varying degree. The same phenomenon was only much later suspected to occur also in animals based on studies of single selected genes in an allopolyploid vertebrate, the Iberian fish Squalius alburnoides. Consequently, this species became a good model for understanding the evolution of gene expression regulation in polyploid vertebrates. Here, we analyzed for the first time genome-wide allele-specific expression data from diploid and triploid hybrids of S. alburnoides and compared homoeolog expression profiles of adult livers and of juveniles. Co-expression of alleles from both parental genomic types was observed for the majority of genes, but with marked homoeolog expression bias, suggesting homoeolog specific reshaping of expression level patterns in hybrids. Complete silencing of one allele was also observed irrespective of ploidy level, but not transcriptome wide as previously speculated. Instead, it was found only in a restricted number of genes, particularly ones with functions related to mitochondria and ribosomes. This leads us to hypothesize that allelic silencing may be a way to overcome intergenomic gene expression interaction conflicts, and that homoeolog expression bias may be an important mechanism in the achievement of sustainable genomic interactions, mandatory to the success of allopolyploid systems, as in S. alburnoides.

Parent-of-Origin-Dependent Gene Expression in Male and Female Schistosome Parasites

Schistosomes are the causative agents of schistosomiasis, a neglected tropical disease affecting over 230 million people worldwide. Additionally to their major impact on human health, they are also models of choice in evolutionary biology. These parasitic flatworms are unique among the common hermaphroditic trematodes as they have separate sexes. This so-called “evolutionary scandal” displays a female heterogametic genetic sex-determination system (ZZ males and ZW females), as well as a pronounced adult sexual dimorphism. These phenotypic differences are determined by a shared set of genes in both sexes, potentially leading to intralocus sexual conflicts. To resolve these conflicts in sexually selected traits, molecular mechanisms such as sex-biased gene expression could occur, but parent-of-origin gene expression also provides an alternative. In this work we investigated the latter mechanism, that is, genes expressed preferentially from either the maternal or the paternal allele, in Schistosoma mansoni species. To this end, transcriptomes from male and female hybrid adults obtained by strain crosses were sequenced. Strain-specific single nucleotide polymorphism (SNP) markers allowed us to discriminate the parental origin, while reciprocal crosses helped to differentiate parental expression from strain-specific expression. We identified genes containing SNPs expressed in a parent-of-origin manner consistent with paternal and maternal imprints. Although the majority of the SNPs was identified in mitochondrial and Z-specific loci, the remaining SNPs found in male and female transcriptomes were situated in genes that have the potential to explain sexual differences in schistosome parasites. Furthermore, we identified and validated four new Z-specific scaffolds.

The extremely divergent maternally- and paternally-transmitted mitochondrial genomes are co-expressed in somatic tissues of two freshwater mussel species with doubly uniparental inheritance of mtDNA

Freshwater mussel species with doubly uniparental inheritance (DUI) of mtDNA are unique because they are naturally heteroplasmic for two extremely divergent mtDNAs with ~50% amino acid differences for protein-coding genes. The paternally-transmitted mtDNA (or M mtDNA) clearly functions in sperm in these species, but it is still unknown whether it is transcribed when present in male or female soma. In the present study, we used PCR and RT-PCR to detect the presence and expression of the M mtDNA in male and female somatic and gonadal tissues of the freshwater mussel species Venustaconcha ellipsiformis and Utterbackia peninsularis (Unionidae). This is the first study demonstrating that the M mtDNA is transcribed not only in male gonads, but also in male and female soma in freshwater mussels with DUI. Because of the potentially deleterious nature of heteroplasmy, we suggest the existence of different mechanisms in DUI species to deal with this possibly harmful situation, such as silencing mechanisms for the M mtDNA at the transcriptional, post-transcriptional and/or post-translational levels. These hypotheses will necessitate additional studies in distantly-related DUI species that could possess different mechanisms of action to deal with heteroplasmy.

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.

Testing the kinship theory of intragenomic conflict in honey bees (Apis mellifera)

Sexual reproduction brings genes from two parents (matrigenes and patrigenes) together into one individual. These genes, despite being unrelated, should show nearly perfect cooperation because each gains equally through the production of offspring. However, an individual's matrigenes and patrigenes can have different probabilities of being present in other relatives, so kin selection could act on them differently. Such intragenomic conflict could be implemented by partial or complete silencing (imprinting) of an allele by one of the parents. Evidence supporting this theory is seen in offspring-mother interactions, with patrigenes favoring acquisition of more of the mother's resources if some of the costs fall on half-siblings who do not share the patrigene. The kinship theory of intragenomic conflict is little tested in other contexts, but it predicts that matrigene-patrigene conflict may be rife in social insects. We tested the hypothesis that honey bee worker reproduction is promoted more by patrigenes than matrigenes by comparing across nine reciprocal crosses of two distinct genetic stocks. As predicted, hybrid workers show reproductive trait characteristics of their paternal stock, (indicating enhanced activity of the patrigenes on these traits), greater patrigenic than matrigenic expression, and significantly increased patrigenic-biased expression in reproductive workers. These results support both the general prediction that matrigene-patrigene conflict occurs in social insects and the specific prediction that honey bee worker reproduction is driven more by patrigenes. The success of these predictions suggests that intragenomic conflict may occur in many contexts where matrigenes and patrigenes have different relatednesses to affected kin.

Cytonuclear Coordination Is Not Immediate upon Allopolyploid Formation in Tragopogon miscellus (Asteraceae) Allopolyploids

Allopolyploids, formed by hybridization and chromosome doubling, face the immediate challenge of having duplicated nuclear genomes that interact with the haploid and maternally inherited cytoplasmic (plastid and mitochondrial) genomes. Most of our knowledge of the genomic consequences of allopolyploidy has focused on the fate of the duplicated nuclear genes without regard to their potential interactions with cytoplasmic genomes. As a step toward understanding the fates of nuclear-encoded subunits that are plastid-targeted, here we examine the retention and expression of the gene encoding the small subunit of Ribulose-1, 5-bisphosphate carboxylase/oxygenase (Rubisco; rbcS) in multiple populations of allotetraploid Tragopogon miscellus (Asteraceae). These polyploids formed recently (~80 years ago) and repeatedly from T. dubius and T. pratensis in the northwestern United States. Examination of 79 T. miscellus individuals from 10 natural populations, as well as 25 synthetic allotetraploids, including reciprocally formed plants, revealed a low percentage of naturally occurring individuals that show a bias in either gene (homeolog) loss (12%) or expression (16%), usually toward maintaining the maternal nuclear copy of rbcS. For individuals showing loss, seven retained the maternally derived rbcS homeolog only, while three had the paternally derived copy. All of the synthetic polyploid individuals examined (S0 and S1 generations) retained and expressed both parental homeologs. These results demonstrate that cytonuclear coordination does not happen immediately upon polyploid formation in Tragopogon miscellus.

Biased Allele Expression and Aggression in Hybrid Honeybees may be Influenced by Inappropriate Nuclear-Cytoplasmic Signaling

Hybrid effects are often exhibited asymmetrically between reciprocal families. One way this could happen is if silencing of one parent’s allele occurs in one lineage but not the other, which could affect the phenotypes of the hybrids asymmetrically by silencing that allele in only one of the hybrid families. We have previously tested for allele-specific expression biases in hybrids of European and Africanized honeybees and we found that there was an asymmetric overabundance of genes showing a maternal bias in the family with a European mother. Here, we further analyze allelic bias in these hybrids to ascertain whether they may underlie previously described asymmetries in metabolism and aggression in similar hybrid families and we speculate on what mechanisms may produce this biased allele usage. We find that there are over 500 genes that have some form of biased allele usage and over 200 of these are biased toward the maternal allele but only in the family with European maternity, mirroring the pattern observed for aggression and metabolic rate. This asymmetrically biased set is enriched for genes in loci associated with aggressive behavior and also for mitochondrial-localizing proteins. It contains many genes that play important roles in metabolic regulation. Moreover we find genes relating to the piwi-interacting RNA (piRNA) pathway, which is involved in chromatin modifications and epigenetic regulation and may help explain the mechanism underlying this asymmetric allele use. Based on these findings and previous work investigating aggression and metabolism in bees, we propose a novel hypothesis; that the asymmetric pattern of biased allele usage in these hybrids is a result of inappropriate use of piRNA-mediated nuclear-cytoplasmic signaling that is normally used to modulate aggression in honeybees. This is the first report of widespread asymmetric effects on allelic expression in hybrids and may represent a novel mechanism for gene regulation.

Whole genome resequencing of Botrytis cinerea isolates identifies high levels of standing diversity

How standing genetic variation within a pathogen contributes to diversity in host/pathogen interactions is poorly understood, partly because most studied pathogens are host-specific, clonally reproducing organisms which complicates genetic analysis. In contrast, Botrytis cinerea is a sexually reproducing, true haploid ascomycete that can infect a wide range of diverse plant hosts. While previous work had shown significant genomic variation between two isolates, we proceeded to assess the level and frequency of standing variation in a population of B. cinerea. To begin measuring standing genetic variation in B. cinerea, we re-sequenced the genomes of 13 different isolates and aligned them to the previously sequenced T4 reference genome. In addition one of these isolates was resequenced from four independently repeated cultures. A high level of genetic diversity was found within the 13 isolates. Within this variation, we could identify clusters of genes with major effect polymorphisms, i.e., polymorphisms that lead to a predicted functional knockout, that surrounded genes involved in controlling vegetative incompatibility. The genotype at these loci was able to partially predict the interaction of these isolates in vegetative fusion assays showing that these loci control vegetative incompatibility. This suggests that the vegetative incompatibility loci within B. cinerea are associated with regions of increased genetic diversity. The genome re-sequencing of four clones from the one isolate (Grape) that had been independently propagated over 10 years showed no detectable spontaneous mutation. This suggests that B. cinerea does not display an elevated spontaneous mutation rate. Future work will allow us to test if, and how, this diversity may be contributing to the pathogen's broad host range.

Nonadditive gene expression in polyploids

Allopolyploidy involves hybridization and duplication of divergent parental genomes and provides new avenues for gene expression. The expression levels of duplicated genes in polyploids can show deviation from parental additivity (the arithmetic average of the parental expression levels). Nonadditive expression has been widely observed in diverse polyploids and comprises at least three possible scenarios: (a) The total gene expression level in a polyploid is similar to that of one of its parents (expression-level dominance); (b) total gene expression is lower or higher than in both parents (transgressive expression); and (c) the relative contribution of the parental copies (homeologs) to the total gene expression is unequal (homeolog expression bias). Several factors may result in expression nonadditivity in polyploids, including maternal-paternal influence, gene dosage balance, cis- and/or trans-regulatory networks, and epigenetic regulation. As our understanding of nonadditive gene expression in polyploids remains limited, a new generation of investigators should explore additional phenomena (i.e., alternative splicing) and use other high-throughput "omics" technologies to measure the impact of nonadditive expression on phenotype, proteome, and metabolome.

A Search for Parent-of-Origin Effects on Honey Bee Gene Expression

Parent-specific gene expression (PSGE) is little known outside of mammals and plants. PSGE occurs when the expression level of a gene depends on whether an allele was inherited from the mother or the father. Kin selection theory predicts that there should be extensive PSGE in social insects because social insect parents can gain inclusive fitness benefits by silencing parental alleles in female offspring. We searched for evidence of PSGE in honey bees using transcriptomes from reciprocal crosses between European and Africanized strains. We found 46 transcripts with significant parent-of-origin effects on gene expression, many of which overexpressed the maternal allele. Interestingly, we also found a large proportion of genes showing a bias toward maternal alleles in only one of the reciprocal crosses. These results indicate that PSGE may occur in social insects. The nonreciprocal effects could be largely driven by hybrid incompatibility between these strains. Future work will help to determine if these are indeed parent-of-origin effects that can modulate inclusive fitness benefits.

Cytonuclear evolution of rubisco in four allopolyploid lineages

Allopolyploidization in plants entails the merger of two divergent nuclear genomes, typically with only one set (usually maternal) of parental plastidial and mitochondrial genomes and with an altered cytonuclear stoichiometry. Thus, we might expect cytonuclear coevolution to be an important dimension of allopolyploid evolution. Here, we investigate cytonuclear coordination for the key chloroplast protein rubisco (ribulose 1,5-bisphosphate carboxylase/oxygenase), which is composed of nuclear-encoded, small subunits (SSUs) and plastid-encoded, large subunits. By studying gene composition and diversity as well as gene expression in four model allopolyploid lineages, Arabidopsis, Arachis, Brassica, and Nicotiana, we demonstrate that paralogous nuclear-encoded rbcS genes within diploids are subject to homogenization via gene conversion and that such concerted evolution via gene conversion characterizes duplicated genes (homoeologs) at the polyploid level. Many gene conversions in the polyploids are intergenomic with respect to the diploid progenitor genomes, occur in functional domains of the homoeologous SSUs, and are directionally biased, such that the maternal amino acid states are favored. This consistent preferential maternal-to-paternal gene conversion is mirrored at the transcriptional level, with a uniform transcriptional bias of the maternal-like rbcS homoeologs. These data, repeated among multiple diverse angiosperm genera for an important photosynthetic enzyme, suggest that cytonuclear coevolution may be mediated by intergenomic gene conversion and altered transcription of duplicated, now homoeologous nuclear genes.

Non-conflict theories for the evolution of genomic imprinting

Theories focused on kinship and the genetic conflict it induces are widely considered to be the primary explanations for the evolution of genomic imprinting. However, there have appeared many competing ideas that do not involve kinship/conflict. These ideas are often overlooked because kinship/conflict is entrenched in the literature, especially outside evolutionary biology. Here we provide a critical overview of these non-conflict theories, providing an accessible perspective into this literature. We suggest that some of these alternative hypotheses may, in fact, provide tenable explanations of the evolution of imprinting for at least some loci.

Gene interactions in the evolution of genomic imprinting

Numerous evolutionary theories have been developed to explain the epigenetic phenomenon of genomic imprinting. Here, we explore a subset of theories wherein non-additive genetic interactions can favour imprinting. In the simplest genic interaction--the case of underdominance--imprinting can be favoured to hide effectively low-fitness heterozygous genotypes; however, as there is no asymmetry between maternally and paternally inherited alleles in this model, other means of enforcing monoallelic expression may be more plausible evolutionary outcomes than genomic imprinting. By contrast, more successful interaction models of imprinting rely on an asymmetry between the maternally and paternally inherited alleles at a locus that favours the silencing of one allele as a means of coordinating the expression of high-fitness allelic combinations. For example, with interactions between autosomal loci, imprinting functionally preserves high-fitness genotypes that were favoured by selection in the previous generation. In this scenario, once a focal locus becomes imprinted, selection at interacting loci favours a matching imprint. Uniparental transmission generates similar asymmetries for sex chromosomes and cytoplasmic factors interacting with autosomal loci, with selection favouring the expression of either maternal or paternally derived autosomal alleles depending on the pattern of transmission of the uniparentally inherited factor. In a final class of models, asymmetries arise when genes expressed in offspring interact with genes expressed in one of its parents. Under such a scenario, a locus evolves to have imprinted expression in offspring to coordinate the interaction with its parent's genome. We illustrate these models and explore key links and differences using a unified framework.

Analysis of genomic imprinting by quantitative allele-specific expression by Pyrosequencing(®)

Genomic imprinting is a parent-of-origin phenomenon whereby gene expression is restricted to the allele inherited from either the maternal or paternal parent. It has been described from flowering plants and eutherian mammals and may have evolved due to parental conflicts over resource allocation. In mammals, imprinted genes are responsible for ensuring correct rates of embryo development and for preventing parthenogenesis. The molecular basis of imprinting depends upon the presence of differential epigenetic marks on the alleles inherited from each parent, although in plants the exact mechanisms that control imprinting are still unclear in many cases. Recent studies have identified large numbers of candidate imprinted genes from Arabidopsis thaliana and other plants (see Chap. 7 by Köhler and colleagues elsewhere in this volume) providing the tools for more thorough investigation into how imprinted gene networks (IGNs) are regulated. Analysis of genomic imprinting in animals has revealed important information on how IGNs are regulated during development, which often involves intermediate levels of imprinting. In some instances, small but significant changes in the degree of parental bias in gene expression have been linked to developmental traits, livestock phenotypes, and human disease. As some of the imprinted genes recently reported from plants show differential rather than complete (binary) imprinting, there is a clear need for tools that can quantify the degree of allelic expression bias occurring at a transcribed locus. In this chapter, we describe the use of Quantification of Allele-Specific Expression by Pyrosequencing(®) (QUASEP) as a tool suitable for this challenge. We describe in detail the factors which ensure that a Pyrosequencing(®) assay will be suitable for giving robust QUASEP and the problems which may be encountered during the study of imprinted genes by Pyrosequencing(®), with particular reference to our work in A. thaliana and in cattle. We also discuss some considerations with respect to the statistical analysis of the resulting data. Finally, we provide a brief overview of the future possibility of adapting Pyrosequencing(®) for analyzing other aspects of imprinting including the analysis of methylated regions.

Cytoplasmic genetic variation and extensive cytonuclear interactions influence natural variation in the metabolome

Understanding genome to phenotype linkages has been greatly enabled by genomic sequencing. However, most genome analysis is typically confined to the nuclear genome. We conducted a metabolomic QTL analysis on a reciprocal RIL population structured to examine how variation in the organelle genomes affects phenotypic variation. This showed that the cytoplasmic variation had effects similar to, if not larger than, the largest individual nuclear locus. Inclusion of cytoplasmic variation into the genetic model greatly increased the explained phenotypic variation. Cytoplasmic genetic variation was a central hub in the epistatic network controlling the plant metabolome. This epistatic influence manifested such that the cytoplasmic background could alter or hide pairwise epistasis between nuclear loci. Thus, cytoplasmic genetic variation plays a central role in controlling natural variation in metabolomic networks. This suggests that cytoplasmic genomes must be included in any future analysis of natural variation. DOI: http://dx.doi.org/10.7554/eLife.00776.001.

Hierarchical nuclear and cytoplasmic genetic architectures for plant growth and defense within Arabidopsis

To understand how genetic architecture translates between phenotypic levels, we mapped the genetic architecture of growth and defense within the Arabidopsis thaliana Kas × Tsu recombinant inbred line population. We measured plant growth using traditional size measurements and size-corrected growth rates. This population contains genetic variation in both the nuclear and cytoplasmic genomes, allowing us to separate their contributions. The cytoplasmic genome regulated a significant variance in growth but not defense, which was due to cytonuclear epistasis. Furthermore, growth adhered to an infinitesimal model of genetic architecture, while defense metabolism was more of a moderate-effect model. We found a lack of concordance between quantitative trait loci (QTL) regulating defense and those regulating growth. Given the published evidence proving the link between glucosinolates and growth, this is likely a false negative result caused by the limited population size. This size limitation creates an inability to test the entire potential genetic landscape possible between these two parents. We uncovered a significant effect of glucosinolates on growth once we accounted for allelic differences in growth QTLs. Therefore, other growth QTLs can mask the effects of defense upon growth. Investigating direct links across phenotypic hierarchies is fraught with difficulty; we identify issues complicating this analysis.

The maternal cytoplasmic environment may be involved in the viability selection of gametes and zygotes

Segregation distortion is the phenomenon whereby the observed genotypic frequencies of a locus fall outside the expected Mendelian segregation ratio, and it is increasingly recognised as a potentially powerful evolutionary force. The main reason for segregation distortion is a difference in the viability of gametes and zygotes caused by viability loci in the segregating progeny. However, the maternal cytoplasm may also be involved in the viability selection of gametes and zygotes. The objectives of this study were to map the segregation distortion loci (SDL) in maize and to test the hypothesis that the viability of gametes and zygotes may also be associated with the maternal cytoplasmic environment. In the present study, a reciprocal mating design was conducted to generate an F2-segregating population. A linkage map was constructed with 126 microsatellite markers. A whole-genome scan was performed to detect the SDL in segregating populations with different maternal cytoplasm environments. Altogether, 14 SDL with strong LOD (logarithm (base 10) of odds) supports were identified in the specifically designed F2 populations. Interestingly, we found dramatic changes in the genotypic frequencies of the SDL in the two maternal cytoplasmic backgrounds, which indicated a change in the viability of gametes and zygotes in different cytoplasmic environments. Furthermore, in the JB cytoplasmic background, most of the detected SDL and complete distortion markers exhibited similar bias patterns favouring the Y53 alleles. These results suggested that selfish cytoplasmic elements may have an important role in shaping the patterns of segregation distortion in F2 populations through selective viability of gametes and zygotes.

Evolution of genomic imprinting as a coordinator of coadapted gene expression

Genomic imprinting is an epigenetic phenomenon in which the expression of a gene copy inherited from the mother differs from that of the copy inherited from the father. Many imprinted genes appear to be highly interconnected through interactions mediated by proteins, RNA, and DNA. These kinds of interactions often favor the evolution of genetic coadaptation, where beneficially interacting alleles evolve to become coinherited. Here I demonstrate theoretically that the presence of gene interactions that favor coadaptation can also favor the evolution of genomic imprinting. Selection favors genomic imprinting because it coordinates the coexpression of positively interacting alleles at different loci. Evolution is expected to proceed through a scenario where selection builds associations between beneficial combinations of alleles and, if one locus evolves to become imprinted, it leads to selection for its interacting partners to match its pattern of imprinting. This process should favor the evolution of physical linkage between interacting genes and therefore may help explain why imprinted genes tend to be found in clusters. The model suggests that, whereas some genes are expected to evolve their imprinting status because selection directly favors a specific pattern of parent-of-origin-dependent expression, other genes may evolve imprinting as a coevolutionary response to match the expression pattern of their interacting partners. As a result, some genes will show phenotypic effects consistent with the predictions of models for the evolution of genomic imprinting (e.g., conflict models), but other genes may not, having simply evolved imprinting to follow the lead of their interacting partners.

De novo assembly of the carrot mitochondrial genome using next generation sequencing of whole genomic DNA provides first evidence of DNA transfer into an angiosperm plastid genome

BACKGROUND

Sequence analysis of organelle genomes has revealed important aspects of plant cell evolution. The scope of this study was to develop an approach for de novo assembly of the carrot mitochondrial genome using next generation sequence data from total genomic DNA.

RESULTS

Sequencing data from a carrot 454 whole genome library were used to develop a de novo assembly of the mitochondrial genome. Development of a new bioinformatic tool allowed visualizing contig connections and elucidation of the de novo assembly. Southern hybridization demonstrated recombination across two large repeats. Genome annotation allowed identification of 44 protein coding genes, three rRNA and 17 tRNA. Identification of the plastid genome sequence allowed organelle genome comparison. Mitochondrial intergenic sequence analysis allowed detection of a fragment of DNA specific to the carrot plastid genome. PCR amplification and sequence analysis across different Apiaceae species revealed consistent conservation of this fragment in the mitochondrial genomes and an insertion in Daucus plastid genomes, giving evidence of a mitochondrial to plastid transfer of DNA. Sequence similarity with a retrotransposon element suggests a possibility that a transposon-like event transferred this sequence into the plastid genome.

CONCLUSIONS

This study confirmed that whole genome sequencing is a practical approach for de novo assembly of higher plant mitochondrial genomes. In addition, a new aspect of intercompartmental genome interaction was reported providing the first evidence for DNA transfer into an angiosperm plastid genome. The approach used here could be used more broadly to sequence and assemble mitochondrial genomes of diverse species. This information will allow us to better understand intercompartmental interactions and cell evolution.

Identification of imprinted genes subject to parent-of-origin specific expression in Arabidopsis thaliana seeds

BACKGROUND

Epigenetic regulation of gene dosage by genomic imprinting of some autosomal genes facilitates normal reproductive development in both mammals and flowering plants. While many imprinted genes have been identified and intensively studied in mammals, smaller numbers have been characterized in flowering plants, mostly in Arabidopsis thaliana. Identification of additional imprinted loci in flowering plants by genome-wide screening for parent-of-origin specific uniparental expression in seed tissues will facilitate our understanding of the origins and functions of imprinted genes in flowering plants.

RESULTS

cDNA-AFLP can detect allele-specific expression that is parent-of-origin dependent for expressed genes in which restriction site polymorphisms exist in the transcripts derived from each allele. Using a genome-wide cDNA-AFLP screen surveying allele-specific expression of 4500 transcript-derived fragments, we report the identification of 52 maternally expressed genes (MEGs) displaying parent-of-origin dependent expression patterns in Arabidopsis siliques containing F1 hybrid seeds (3, 4 and 5 days after pollination). We identified these MEGs by developing a bioinformatics tool (GenFrag) which can directly determine the identities of transcript-derived fragments from (i) their size and (ii) which selective nucleotides were added to the primers used to generate them. Hence, GenFrag facilitates increased throughput for genome-wide cDNA-AFLP fragment analyses. The 52 MEGs we identified were further filtered for high expression levels in the endosperm relative to the seed coat to identify the candidate genes most likely representing novel imprinted genes expressed in the endosperm of Arabidopsis thaliana. Expression in seed tissues of the three top-ranked candidate genes, ATCDC48, PDE120 and MS5-like, was confirmed by Laser-Capture Microdissection and qRT-PCR analysis. Maternal-specific expression of these genes in Arabidopsis thaliana F1 seeds was confirmed via allele-specific transcript analysis across a range of different accessions. Differentially methylated regions were identified adjacent to ATCDC48 and PDE120, which may represent candidate imprinting control regions. Finally, we demonstrate that expression levels of these three genes in vegetative tissues are MET1-dependent, while their uniparental maternal expression in the seed is not dependent on MET1.

CONCLUSIONS

Using a cDNA-AFLP transcriptome profiling approach, we have identified three genes, ATCDC48, PDE120 and MS5-like which represent novel maternally expressed imprinted genes in the Arabidopsis thaliana seed. The extent of overlap between our cDNA-AFLP screen for maternally expressed imprinted genes, and other screens for imprinted and endosperm-expressed genes is discussed.

Mating system and the critical migration rate for swamping selection

Crow et al. (1990) and Barton (1992) have examined the critical migration rate for swamping selection in the nuclear system. Here, I use the same methodology to examine the critical migration rate in the cytonuclear system for hermaphrodite plants with a mixed mating system. Two selection schemes for a nuclear gene (heterozygote disadvantage and directional selection) and the directional selection scheme for organelle genes are considered. Results show that under random mating, the previous results are applicable to plant species by appropriate re-parameterization of the migration rate for nuclear and paternal organelle genes. A simple complementary relationship exists between seed and pollen flow in contributing to the critical migration rate. Under the mixed mating system, the critical migration rate of seeds and pollen for nuclear and paternal organelle genes can be changed due to the effects of selection and the cytonuclear linkage disequilibrium generated by migration and inbreeding. A negative but not complementary relationship exists between seed and pollen flow in contributing to the critical migration rate, varying with the mating system. Partial selfing can also adjust the critical seed flow for the maternal organelle gene, with a small critical migration rate for species of a high selfing rate. Both concordance and discordance among cytonuclear genes can occur under certain conditions during the process of swamping selection. This theory predicts the presence of various contributions of seed versus pollen flow to genetic swamping for plants with diverse mating systems.

Matrisibs, patrisibs, and the evolution of imprinting on autosomes and sex chromosomes

The conflict theory of genomic imprinting argues that parent-of-origin effects on allelic expression evolve as a consequence of conflict between maternally and paternally derived genomes. I derive explicit population-genetic models of this theory when individuals in a cohort with an arbitrary and variable number of sires and dams interact. I show that the evolution of imprinting is governed by the reciprocal of the harmonic mean number of fathers but the reciprocal of the arithmetic mean number of mothers per cohort. Thus, a few monandrous females in a polyandrous population decrease the strength of the genetic conflict and the opportunity for conflict-driven paternal imprinting. In contrast, in populations in which few males control large harems, rare males with small harems do not have such a disproportionate effect on genetic conflicts and maternal imprinting. Additionally, I demonstrate that under the conflict theory, selection for imprinted expression on paternally derived X chromosomes is much weaker than it is on maternally derived X chromosomes or autosomes.

F(ST) in the cytonuclear system

Selection on nuclear (or organelle) sites inevitably affects the spatial distribution of a neutral organelle (or nuclear) allele via transient cytonuclear disequilibrium. Here I examine this effect in terms of F(st) for a neutral allele by bringing together cytonuclear genomes with contrasting modes of inheritance. The relationships between cytonuclear disequilibrium and increment in F(st) are explored and confirmed through Monte Carlo simulations. Results show that the transient increment in F(st) for a neutral allele is not only related to the vectors of seed and pollen dispersal but also to the mode of its inheritance. Such increments can be substantial under certain conditions. Seed dispersal is more effective than pollen dispersal in changing the transient increment. The cumulative effects from multiple selective nuclear sites can amplify the transient increment in F(st) for a neutral paternal or maternal organelle allele. Selection on selective organelle sites facilitates the transient increment in F(st) for a neutral nuclear allele. Partial selfing can significantly reinforce the transient increment in F(st). These theoretical insights highlight the roles of transient cytonuclear disequilibrium as a biological factor in evolving population differentiation and refine our practical interpretations of F(st) with cytonuclear markers.