The role of organelle genomes in plant adaptation: time to get to work!

A Scalable and Robust Chloroplast Genotyping Solution: Development and Application of SNP and InDel Markers in the Maize Chloroplast Genome

Maize(Zea mays. L) is a globally important crop, and understanding its genetic diversity is crucial for plant breeding phylogenetic analyses and comparative genetics. While nuclear markers have been extensively used for mapping agriculturally important genes, they are limited in recognizing characteristics, such as cytoplasmic male sterility and reciprocal cross hybrids. In this study, we performed next-generation sequencing of 176samples, and the maize cultivars represented five distinct groups. A total of 89 single nucleotide polymorphisms (SNPs) and 11 insertion/deletion polymorphisms (InDels) were identified. To enable high-throughput detection, we successfully amplified and confirmed 49 SNP and InDel markers, which were defined as a Varietal Chloroplast Panel (VCP) using the Kompetitive Allele Specific PCR (KASP). The specific markers provided a valuable tool for identifying chloroplast groups. The verification experiment, focusing on the identification of reciprocal cross hybrids and cytoplasmic male sterility hybrids, demonstrated the significant advantages of VCP markers in maternal inheritance characterization. Furthermore, only a small subset of these markers is needed to provide useful information, showcasing the effectiveness of these markers in elucidating the artificial selection process of elite maize lines.

Comparative Analysis of Tylosema esculentum Mitochondrial DNA Revealed Two Distinct Genome Structures

Tylosema esculentum, commonly known as the marama bean, is an underutilized legume with nutritious seeds, holding potential to enhance food security in southern Africa due to its resilience to prolonged drought and heat. To promote the selection of this agronomically valuable germplasm, this study assembled and compared the mitogenomes of 84 marama individuals, identifying variations in genome structure, single-nucleotide polymorphisms (SNPs), insertions/deletions (indels), heteroplasmy, and horizontal transfer. Two distinct germplasms were identified, and a novel mitogenome structure consisting of three circular molecules and one long linear chromosome was discovered. The structural variation led to an increased copy number of specific genes, nad5, nad9, rrnS, rrn5, trnC, and trnfM. The two mitogenomes also exhibited differences at 230 loci, with only one notable nonsynonymous substitution in the matR gene. Heteroplasmy was concentrated at certain loci on chromosome LS1 (OK638188). Moreover, the marama mitogenome contained an over 9 kb insertion of cpDNA, originating from chloroplast genomes, but had accumulated mutations and lost gene functionality. The evolutionary and comparative genomics analysis indicated that mitogenome divergence in marama might not be solely constrained by geographical factors. Additionally, marama, as a member from the Cercidoideae subfamily, tends to possess a more complete set of mitochondrial genes than Faboideae legumes.

Comparative analysis of the organelle genomes of three Rhodiola species provide insights into their structural dynamics and sequence divergences

BACKGROUND

Plant organelle genomes are a valuable resource for evolutionary biology research, yet their genome architectures, evolutionary patterns and environmental adaptations are poorly understood in many lineages. Rhodiola species is a type of flora mainly distributed in highland habitats, with high medicinal value. Here, we assembled the organelle genomes of three Rhodiola species (R. wallichiana, R. crenulata and R. sacra) collected from the Qinghai-Tibet plateau (QTP), and compared their genome structure, gene content, structural rearrangements, sequence transfer and sequence evolution rates.

RESULTS

The results demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes in three Rhodiola species, with the former possessing more conserved genome structure but faster evolutionary rates of sequence, while the latter exhibiting structural diversity but slower rates of sequence evolution. Some lineage-specific features were observed in Rhodiola mitogenomes, including chromosome fission, gene loss and structural rearrangement. Repeat element analysis shows that the repeats occurring between the two chromosomes may mediate the formation of multichromosomal structure in the mitogenomes of Rhodiola, and this multichromosomal structure may have recently formed. The identification of homologous sequences between plastomes and mitogenomes reveals several unidirectional protein-coding gene transfer events from chloroplasts to mitochondria. Moreover, we found that their organelle genomes contained multiple fragments of nuclear transposable elements (TEs) and exhibited different preferences for TEs insertion type. Genome-wide scans of positive selection identified one gene matR from the mitogenome. Since the matR is crucial for plant growth and development, as well as for respiration and stress responses, our findings suggest that matR may participate in the adaptive response of Rhodiola species to environmental stress of QTP.

CONCLUSION

The study analyzed the organelle genomes of three Rhodiola species and demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes. Signals of positive selection were detected in the matR gene of Rhodiola mitogenomes, suggesting the potential role of this gene in Rhodiola adaptation to QTP. Together, the study is expected to enrich the genomic resources and provide valuable insights into the structural dynamics and sequence divergences of Rhodiola species.

Genomics Armed With Diversity Leads the Way in Brassica Improvement in a Changing Global Environment

Meeting the needs of a growing world population in the face of imminent climate change is a challenge; breeding of vegetable and oilseed Brassica crops is part of the race in meeting these demands. Available genetic diversity constituting the foundation of breeding is essential in plant improvement. Elite varieties, land races, and crop wild species are important resources of useful variation and are available from existing genepools or genebanks. Conservation of diversity in genepools, genebanks, and even the wild is crucial in preventing the loss of variation for future breeding efforts. In addition, the identification of suitable parental lines and alleles is critical in ensuring the development of resilient Brassica crops. During the past two decades, an increasing number of high-quality nuclear and organellar Brassica genomes have been assembled. Whole-genome re-sequencing and the development of pan-genomes are overcoming the limitations of the single reference genome and provide the basis for further exploration. Genomic and complementary omic tools such as microarrays, transcriptomics, epigenetics, and reverse genetics facilitate the study of crop evolution, breeding histories, and the discovery of loci associated with highly sought-after agronomic traits. Furthermore, in genomic selection, predicted breeding values based on phenotype and genome-wide marker scores allow the preselection of promising genotypes, enhancing genetic gains and substantially quickening the breeding cycle. It is clear that genomics, armed with diversity, is set to lead the way in Brassica improvement; however, a multidisciplinary plant breeding approach that includes phenotype = genotype × environment × management interaction will ultimately ensure the selection of resilient Brassica varieties ready for climate change.

Study of Androgenic Plant Families of Alloplasmic Introgression Lines (H. vulgare) -T. aestivum and the Use of Sister DH Lines in Breeding

One of the limitations in obtaining the genetic diversity of doubled haploid (DH) lines via anther culture is the development of families of regenerants, and each family represents a clone. This work examines the results of studying this phenomenon in anther culture of alloplasmic (H. vulgare)–T. aestivum and euplasmic lines with 1RS.1BL and 7DL-7Ai translocations and hybrids between them. Parameters of androgenesis such as the number of embryo-like structures, the total number of regenerants, and the number of green regenerants per 100 anthers varied depending on the genotype. In all genotypes from embryo-like structures, predominant development of families of plantlets rather than single plantlets was found. The source of family plantlets was polyembryos. About 75% of families consisted of regenerants at the same fertility level. On average, 37.74%4% of the R0 plants were fertile. The sister DH lines of three hybrid combinations were formed from seeds of R1 plants (2n = 42) with high fertility and in the presence of wheat–alien translocations. After four years of breeding trials, the sister DH lines of three families with fungal disease resistance increased yield, and some parameters of grain quality exceeding the controls were identified as promising for breeding.

Plastid-targeted forms of restriction endonucleases enhance the plastid genome rearrangement rate and trigger the reorganization of its genomic architecture

Plant cells have acquired chloroplasts (plastids) with a unique genome (ptDNA), which developed during the evolution of endosymbiosis. The gene content and genome structure of ptDNAs in land plants are considerably stable, although those of algal ptDNAs are highly varied. Plant cells seem, therefore, to be intolerant of any structural or organizational changes in the ptDNA. Genome rearrangement functions as a driver of genomic evolutionary divergence. Here, we aimed to create various types of rearrangements in the ptDNA of Arabidopsis genomes using plastid-targeted forms of restriction endonucleases (pREs). Arabidopsis plants expressing each of the three specific pREs, i.e., pTaqI, pHinP1I, and pMseI, were generated; they showed the leaf variegation phenotypes associated with impaired chloroplast development. We confirmed that these pREs caused double-stranded breaks (DSB) at their recognition sites in ptDNAs. Genome-wide analysis of ptDNAs revealed that the transgenic lines exhibited a large number of rearrangements such as inversions and deletions/duplications, which were dominantly repaired by microhomology-mediated recombination and microhomology-mediated end-joining, and less by non-homologous end-joining. Notably, pHinP1I, which recognized a small number of sites in ptDNA, induced drastic structural changes, including regional copy number variations throughout ptDNAs. In contrast, the transient expression of either pTaqI or pMseI, whose recognition site numbers were relatively larger, resulted in small-scale changes at the whole genome level. These results indicated that DSB frequencies and their distribution are major determinants in shaping ptDNAs.

Wide crossing diversify mitogenomes of rice

BACKGROUND

In most angiosperms, the inheritance of the mitochondria takes place in a typical maternal manner. However, very less information is available about if the existence of structural variations or not in mitochondrial genomes (mitogenomes) between maternal parents and their progenies.

RESULTS

In order to find the answer, a stable rice backcross inbred line (BIL) population was derived from the crosses of Oryza glaberrima/Oryza sativa//Oryza sativa. The current study presents a comparative analysis of the mitogenomes between maternal parents and five BILs. There were recorded universal structural variations such as reversal, translocation, fusion, and fission among the BILs. The repeat-mediated recombination and non-homologous end-joining contributed virtually equal to the rearrangement of mitogenomes. Similarly, the relative order, copy-number, expression level, and RNA-editing rate of mitochondrial genes were also extensively varied among BILs.

CONCLUSIONS

These novel findings unraveled an unusual mystery of the maternal inheritance and possible cause for heterogeneity of mitogenomes in rice population. The current piece of work will greatly develop our understanding of the plant nucleo-cytoplasmic interaction and their potential role in plant growth and developmental processes.

Cytonuclear Genetic Incompatibilities in Plant Speciation

Due to the endosymbiotic origin of organelles, a pattern of coevolution and coadaptation between organellar and nuclear genomes is required for proper cell function. In this review, we focus on the impact of cytonuclear interaction on the reproductive isolation of plant species. We give examples of cases where species exhibit barriers to reproduction which involve plastid-nuclear or mito-nuclear genetic incompatibilities, and describe the evolutionary processes at play. We also discuss potential mechanisms of hybrid fitness recovery such as paternal leakage. Finally, we point out the possible interplay between plant mating systems and cytonuclear coevolution, and its consequence on plant speciation.

Two reproductive traits show contrasting genetic architectures in Plantago lanceolata

In many species, temperature-sensitive phenotypic plasticity (i.e., an individual's phenotypic response to temperature) displays a positive correlation with latitude, a pattern presumed to reflect local adaptation. This geographical pattern raises two general questions: (a) Do a few large-effect genes contribute to latitudinal variation in a trait? (b) Is the thermal plasticity of different traits regulated pleiotropically? To address the questions, we crossed individuals of Plantago lanceolata derived from northern and southern European populations. Individuals naturally exhibited high and low thermal plasticity in floral reflectance and flowering time. We grew parents and offspring in controlled cool- and warm-temperature environments, mimicking what plants would encounter in nature. We obtained genetic markers via genotype-by-sequencing, produced the first recombination map for this ecologically important nonmodel species, and performed quantitative trait locus (QTL) mapping of thermal plasticity and single-environment values for both traits. We identified a large-effect QTL that largely explained the reflectance plasticity differences between northern and southern populations. We identified multiple smaller-effect QTLs affecting aspects of flowering time, one of which affected flowering time plasticity. The results indicate that the genetic architecture of thermal plasticity in flowering is more complex than for reflectance. One flowering time QTL showed strong cytonuclear interactions under cool temperatures. Reflectance and flowering plasticity QTLs did not colocalize, suggesting little pleiotropic genetic control and freedom for independent trait evolution. Such genetic information about the architecture of plasticity is environmentally important because it informs us about the potential for plasticity to offset negative effects of climate change.

Mitochondrial genetic effects on reproductive success: signatures of positive intrasexual, but negative intersexual pleiotropy

Theory predicts that maternal inheritance of mitochondria will facilitate the accumulation of mtDNA mutations that are male biased, or even sexually antagonistic, in effect. While there are many reported cases of mtDNA mutations conferring cytoplasmic male sterility in plants, historically it was assumed such mutations would not persist in the streamlined mitochondrial genomes of bilaterian metazoans. Intriguingly, recent cases of mitochondrial variants exerting male biases in effect have come to light in bilaterians. These cases aside, it remains unknown whether the mitochondrial genetic variation affecting phenotypic expression, and in particular reproductive performance, in bilaterians is routinely composed of sex-biased or sex-specific variation. If selection consistently favours mtDNA variants that augment female fitness, but at cost to males, this could shape patterns of pleiotropy and lead to negative intersexual correlations across mtDNA haplotypes. Here, we show that genetic variation across naturally occurring mitochondrial haplotypes affects components of reproductive success in both sexes, in the fruit fly Drosophila melanogaster We find that intrasexual correlations across mitochondrial haplotypes, for components of reproductive success, are generally positive, while intersexual correlations are negative. These results accord with theoretical predictions, suggesting that maternal inheritance has led to the fixation of numerous mutations of sexually antagonistic effect.

The plastid and mitochondrial genomes of Eucalyptus grandis

BACKGROUND

Land plant organellar genomes have significant impact on metabolism and adaptation, and as such, accurate assembly and annotation of plant organellar genomes is an important tool in understanding the evolutionary history and interactions between these genomes. Intracellular DNA transfer is ongoing between the nuclear and organellar genomes, and can lead to significant genomic variation between, and within, species that impacts downstream analysis of genomes and transcriptomes.

RESULTS

In order to facilitate further studies of cytonuclear interactions in Eucalyptus, we report an updated annotation of the E. grandis plastid genome, and the second sequenced and annotated mitochondrial genome of the Myrtales, that of E. grandis. The 478,813 bp mitochondrial genome shows the conserved protein coding regions and gene order rearrangements typical of land plants. There have been widespread insertions of organellar DNA into the E. grandis nuclear genome, which span 141 annotated nuclear genes. Further, we identify predicted editing sites to allow for the discrimination of RNA-sequencing reads between nuclear and organellar gene copies, finding that nuclear copies of organellar genes are not expressed in E. grandis.

CONCLUSIONS

The implications of organellar DNA transfer to the nucleus are often ignored, despite the insight they can give into the ongoing evolution of plant genomes, and the problems they can cause in many applications of genomics. Future comparisons of the transcription and regulation of organellar genes between Eucalyptus genotypes may provide insight to the cytonuclear interactions that impact economically important traits in this widely grown lignocellulosic crop species.

Novel Cytonuclear Combinations Modify Arabidopsis thaliana Seed Physiology and Vigor

Dormancy and germination vigor are complex traits of primary importance for adaptation and agriculture. Intraspecific variation in cytoplasmic genomes and cytonuclear interactions were previously reported to affect germination in Arabidopsis using novel cytonuclear combinations that disrupt co-adaptation between natural variants of nuclear and cytoplasmic genomes. However, specific aspects of dormancy and germination vigor were not thoroughly explored, nor the parental contributions to the genetic effects. Here, we specifically assessed dormancy, germination performance and longevity of seeds from Arabidopsis plants with natural and new genomic compositions. All three traits were modified by cytonuclear reshuffling. Both depth and release rate of dormancy could be modified by a changing of cytoplasm. Significant changes on dormancy and germination performance due to specific cytonuclear interacting combinations mainly occurred in opposite directions, consistent with the idea that a single physiological consequence of the new genetic combination affected both traits oppositely. However, this was not always the case. Interestingly, the ability of parental accessions to contribute to significant cytonuclear interactions modifying the germination phenotype was different depending on whether they provided the nuclear or cytoplasmic genetic compartment. The observed deleterious effects of novel cytonuclear combinations (in comparison with the nuclear parent) were consistent with a contribution of cytonuclear interactions to germination adaptive phenotypes. More surprisingly, we also observed favorable effects of novel cytonuclear combinations, suggesting suboptimal genetic combinations exist in natural populations for these traits. Reduced sensitivity to exogenous ABA and faster endogenous ABA decay during germination were observed in a novel cytonuclear combination that also exhibited enhanced longevity and better germination performance, compared to its natural nuclear parent. Taken together, our results strongly support that cytoplasmic genomes represent an additional resource of natural variation for breeding seed vigor traits.

Genetic diversity of submergence stress response in cytoplasms of the Triticum-Aegilops complex

Genetic diversity in cytoplasmic and nuclear genomes and their interaction affecting adaptive traits is an attractive research subject in plants. We addressed submergence stress response of wheat that has become increasingly important but remained largely uninvestigated. Our primary aim was to disclose cytoplasmic diversity using nucleus-cytoplasm (NC) hybrids possessing a series of heterologous cytoplasms in a common nuclear background. Effects of submergence on seedling emergence and growth from imbibed seeds were studied and compared with euplasmic lines. Marked phenotypic variabilities were observed among both lines, demonstrating divergent cytoplasmic and nuclear effects on submergence response. NC hybrids with cytoplasm of Aegilops mutica showed a less inhibition, indicative of their positive contribution to submergence tolerance, whereas cytoplasms of Aegilops umbellulata and related species caused a greater inhibition. Superoxide dismutase (SOD) activity showed a marked increase accompanied by retardation of seedling growth in a susceptible NC hybrid. The observation suggested that the elevated SOD activity was resulted from a high level of reactive oxygen species accumulated and remained in susceptible seedlings. Taken together, our results point to the usefulness of NC hybrids in further studies needed to clarify molecular mechanisms underlying the nucleus-cytoplasm interaction regulating submergence stress response in wheat.

Cytonuclear interactions affect adaptive traits of the annual plant Arabidopsis thaliana in the field

Although the contribution of cytonuclear interactions to plant fitness variation is relatively well documented at the interspecific level, the prevalence of cytonuclear interactions at the intraspecific level remains poorly investigated. In this study, we set up a field experiment to explore the range of effects that cytonuclear interactions have on fitness-related traits in Arabidopsis thaliana To do so, we created a unique series of 56 cytolines resulting from cytoplasmic substitutions among eight natural accessions reflecting within-species genetic diversity. An assessment of these cytolines and their parental lines scored for 28 adaptive whole-organism phenotypes showed that a large proportion of phenotypic traits (23 of 28) were affected by cytonuclear interactions. The effects of these interactions varied from slight but frequent across cytolines to strong in some specific parental pairs. Two parental pairs accounted for half of the significant pairwise interactions. In one parental pair, Ct-1/Sha, we observed symmetrical phenotypic responses between the two nuclear backgrounds when combined with specific cytoplasms, suggesting nuclear differentiation at loci involved in cytonuclear epistasis. In contrast, asymmetrical phenotypic responses were observed in another parental pair, Cvi-0/Sha. In the Cvi-0 nuclear background, fecundity and phenology-related traits were strongly affected by the Sha cytoplasm, leading to a modified reproductive strategy without penalizing total seed production. These results indicate that natural variation in cytoplasmic and nuclear genomes interact to shape integrative traits that contribute to adaptation, thereby suggesting that cytonuclear interactions can play a major role in the evolutionary dynamics ofA. thaliana.

Novel nuclear-cytoplasmic interaction in wheat (Triticum aestivum) induces vigorous plants

Interspecific hybridization can be considered an accelerator of evolution, otherwise a slow process, solely dependent on mutation and recombination. Upon interspecific hybridization, several novel interactions between nuclear and cytoplasmic genomes emerge which provide additional sources of diversity. The magnitude and essence of intergenomic interactions between nuclear and cytoplasmic genomes remain unknown due to the direction of many crosses. This study was conducted to address the role of nuclear-cytoplasmic interactions as a source of variation upon hybridization. Wheat (Triticum aestivum) alloplasmic lines carrying the cytoplasm of Aegilops mutica along with an integrated approach utilizing comparative quantitative trait locus (QTL) and epigenome analysis were used to dissect this interaction. The results indicate that cytoplasmic genomes can modify the magnitude of QTL controlling certain physiological traits such as dry matter weight. Furthermore, methylation profiling analysis detected eight polymorphic regions affected by the cytoplasm type. In general, these results indicate that novel nuclear-cytoplasmic interactions can potentially trigger an epigenetic modification cascade in nuclear genes which eventually change the genetic network controlling physiological traits. These modified genetic networks can serve as new sources of variation to accelerate the evolutionary process. Furthermore, this variation can synthetically be produced by breeders in their programs to develop epigenomic-segregating lines.

On the adaptive value of cytoplasmic genomes in plants

Is DNA variation maintained in organelle genomes selectively neutral? The answer to this question has important implications for many aspects of ecology and evolution. While traditionally the answer has been 'yes', recent studies in animals have shown that, on the contrary, mitochondrial DNA polymorphism is frequently adaptive. In plants, however, the neutrality assumption has not been strongly challenged. Here, we begin with a critical evaluation of arguments in favour of this long-held view. We then discuss the latest empirical evidence for the opposing prediction that sequence variation in plant cytoplasmic genomes is frequently adaptive. While outstanding research progress is being made towards understanding this fundamental topic, we highlight the need for studies that combine information ranging from field experiments to physiology to molecular evolutionary biology. Such an interdisciplinary approach provides a means for determining the frequency, drivers and evolutionary significance of adaptive organelle DNA variation.

Genetic signatures of natural selection in response to air pollution in red spruce (Picea rubens, Pinaceae)

One of the most important drivers of local adaptation for forest trees is climate. Coupled to these patterns, however, are human-induced disturbances through habitat modification and pollution. The confounded effects of climate and disturbance have rarely been investigated with regard to selective pressure on forest trees. Here, we have developed and used a population genetic approach to search for signals of selection within a set of 36 candidate genes chosen for their putative effects on adaptation to climate and human-induced air pollution within five populations of red spruce (Picea rubens Sarg.), distributed across its natural range and air pollution gradient in eastern North America. Specifically, we used FST outlier and environmental correlation analyses to highlight a set of seven single nucleotide polymorphisms (SNPs) that were overly correlated with climate and levels of sulphate pollution after correcting for the confounding effects of population history. Use of three age cohorts within each population allowed the effects of climate and pollution to be separated temporally, as climate-related SNPs (n = 7) showed the strongest signals in the oldest cohort, while pollution-related SNPs (n = 3) showed the strongest signals in the youngest cohorts. These results highlight the usefulness of population genetic scans for the identification of putatively nonneutral evolution within genomes of nonmodel forest tree species, but also highlight the need for the development and application of robust methodologies to deal with the inherent multivariate nature of the genetic and ecological data used in these types of analyses.

The evolutionary ecology of cytonuclear interactions in angiosperms

Interactions between cytoplasmic and nuclear genomes have significant evolutionary consequences. In angiosperms, the most common cytonuclear interaction is between mitochondrial genes that disrupt pollen production (cytoplasmic male sterility, CMS) and nuclear genes that restore it (nuclear male fertility restorers, Rf). The outcome of CMS/Rf interactions can depend on whether Rf alleles have negative pleiotropic effects on fitness. Although these fitness costs are often considered to be independent of the ecological context, we argue that the effects of Rf alleles on fitness should be context dependent. Thus, measuring the cost of restoration across a range of environments could help explain geographic and phylogenetic variation in the distribution of Rf alleles and the outcome of CMS/Rf interactions.

Genetic basis of local adaptation and flowering time variation in Arabidopsis lyrata

Understanding how genetic variation at individual loci contributes to adaptation of populations to different local environments is an important topic in modern evolutionary biology. To date, most evidence has pointed to conditionally neutral quantitative trait loci (QTL) showing fitness effects only in some environments, while there has been less evidence for single-locus fitness trade-offs. At QTL underlying local adaptation, alleles from the local population are expected to show a fitness advantage. Cytoplasmic genomes also can have a role in local adaptation, but the role of cytonuclear interactions in adaptive differentiation has remained largely unknown. We mapped genomic regions underlying adaptive differentiation in multiple fitness components and flowering time in diverged populations of a perennial plant Arabidopsis lyrata. Experimental hybrids for this purpose were grown in natural field conditions of the parental populations in Norway and North Carolina (NC), USA, and in the greenhouse. We found QTL where high fitness and early flowering were associated with local alleles, indicating a role of different selection pressures in phenotypic differentiation. At two QTL regions, a fitness component showing local adaptation between the parental populations also showed signs of putative fitness trade-offs. Beneficial dominance effects of conditionally neutral QTL for different fitness components resulted in hybrid vigour at the Norwegian site in the F(2) hybrids. We also found that cytoplasmic genomes contributed to local adaptation and hybrid vigour by interacting with nuclear QTL, but these interactions did not show evidence for cytonuclear coadaptation (high fitness of local alleles combined with the local cytoplasm).