Nuclear Integrants of Organellar DNA Contribute to Genome Structure and Evolution in Plants

The extent and characteristics of DNA transfer between plasmids and chromosomes

Plasmids are extrachromosomal genetic elements that reside in prokaryotes. The acquisition of plasmids encoding beneficial traits can facilitate short-term survival in harsh environmental conditions or long-term adaptation of new ecological niches. Due to their ability to transfer between cells, plasmids are considered agents of gene transfer. Nonetheless, the frequency of DNA transfer between plasmids and chromosomes remains understudied. Using a novel approach for detection of homologous loci between genome pairs, we uncover gene sharing with the chromosome in 1,974 (66%) plasmids residing in 1,016 (78%) taxonomically diverse isolates. The majority of homologous loci correspond to mobile elements, which may be duplicated in the host chromosomes in tens of copies. Neighboring shared genes often encode similar functional categories, indicating the transfer of multigene functional units. Rare transfer events of antibiotics resistance genes are observed mainly with mobile elements. The frequent erosion of sequence similarity in homologous regions indicates that the transferred DNA is often devoid of function. DNA transfer between plasmids and chromosomes thus generates genetic variation that is akin to workings of endosymbiotic gene transfer in eukaryotic evolution. Our findings imply that plasmid contribution to gene transfer most often corresponds to transfer of the plasmid entity rather than transfer of protein-coding genes between plasmids and chromosomes.

Residual Microglia Following Short-term PLX5622 Treatment in 5xFAD Mice Exhibit Diminished NLRP3 Inflammasome and mTOR Signaling, and Enhanced Autophagy

Chronic neuroinflammation represents a prominent hallmark of Alzheimer's disease (AD). While moderately activated microglia are pivotal in clearing amyloid beta (Aβ), hyperactivated microglia perpetuate neuroinflammation. Prior investigations have indicated that the elimination of ∼80% of microglia through a month-long inhibition of the colony-stimulating factor 1 receptor (CSF1R) during the advanced stage of neuroinflammation in 5xFamilial AD (5xFAD) mice mitigates synapse loss and neurodegeneration without impacting Aβ levels. Furthermore, prolonged CSF1R inhibition diminished the development of parenchymal plaques. Nonetheless, the immediate effects of short-term CSF1R inhibition during the early stages of neuroinflammation on residual microglial phenotype or metabolic fitness are unknown. Therefore, we investigated the effects of 10-day CSF1R inhibition in three-month-old female 5xFAD mice, a stage characterized by the onset of neuroinflammation and minimal Aβ plaques. We observed ∼65% microglia depletion in the hippocampus and cerebral cortex. The leftover microglia demonstrated a noninflammatory phenotype, with highly branched and ramified processes and reduced NOD-, LRR-, and pyrin domain-containing protein 3 (NLRP3) inflammasome complexes. Moreover, plaque-associated microglia were reduced in number with diminished Clec7a (dectin-1) expression. Additionally, both microglia and neurons displayed reduced mechanistic target of rapamycin (mTOR) signaling and autophagy. Biochemical assays validated the inhibition of NLRP3 inflammasome activation, decreased mTOR signaling, and enhanced autophagy. However, short-term CSF1R inhibition did not influence Aβ plaques, soluble Aβ-42 levels, or hippocampal neurogenesis. Thus, short-term CSF1R inhibition during the early stages of neuroinflammation in 5xFAD mice promotes the retention of homeostatic microglia with diminished inflammasome activation and mTOR signaling, alongside increased autophagy.

Unraveling the complex evolutionary features of the Cinnamomum camphora mitochondrial genome

KEY MESSAGE

We reported the mitochondrial genome of Cinnamomum camphora for the first time, revealing frequent rearrangement events in the non-coding regions of Magnoliids mitochondrial genomes. As one of the representative species in the Lauraceae family of Magnoliids, Cinnamomum camphora holds significant economic and ecological value. In this study, the mitochondrial genome (mitogenome) of C. camphora was complete assembled and annotated using PacBio HiFi sequencing. The C. camphora mitogenome is characterized by a branch structure, spans 900,894 bp, and contains 43 protein-coding genes (PCGs), 24 tRNAs, and 3 rRNAs. Most of these PCGs are under purifying selection, with only two (ccmFc and rps7) exhibiting signs of positive selection. The C. camphora mitogenome contains numerous repetitive sequences and intracellular gene transfers, with a total of 36 mitochondrial plastid DNAs, amounting to a combined length of 23,816 bp. Comparative analysis revealed that the non-coding regions of Magnoliids mitogenomes have undergone frequent rearrangements during evolution, but the coding sequences remain highly conserved (more than 98% similarity for protein-coding sequences). Furthermore, a maximum-likelihood phylogenetic tree was reconstructed based on 25 PCGs from 23 plant mitogenomes. The analysis supports the closest relationship between C. camphora and C. chekiangense, consistent with the APG IV classification system. This study elucidates the unique evolutionary features of the C. camphora mitogenome, which will provide valuable insights into the study of genetics and evolution of the family Lauraceae.

Mitochondrial genome complexity in Stemona sessilifolia: nanopore sequencing reveals chloroplast gene transfer and DNA rearrangements

Mitochondria are semi-autonomous organelles in eukaryotic cells with their own genome. Plant mitogenomes differ from animal mitogenomes in size, structure, and repetitive DNA sequences. Despite larger sizes, plant mitogenomes do not have significantly more genes. They exhibit diverse structures due to variations in size, repetitive DNA, recombination frequencies, low gene densities, and reduced nucleotide substitution rates. In this study, we analyzed the mitochondrial genome of Stemona sessilifolia using Nanopore and Illumina sequencing. De-novo assembly and annotation were conducted using Unicycler, Geseq, tRNAscan-SE and BLASTN, followed by codon usage, repeat sequence, RNA-editing, synteny, and phylogenetic analyses. S. sessilifolia's mitogenome consisted of one linear contig and six circular contigs totaling 724,751 bp. It had 39 protein-coding genes, 27 tRNA genes, and 3 rRNA genes. Transfer of chloroplast sequences accounted for 13.14% of the mitogenome. Various analyses provided insights into genetic characteristics, evolutionary dynamics, and phylogenetic placement. Further investigations can explore transferred genes' functions and RNA-editing's role in mitochondrial gene expression in S. sessilifolia.

Plastid DNA is a major source of nuclear genome complexity and of RNA genes in the orphan crop moringa

BACKGROUND

Unlike Transposable Elements (TEs) and gene/genome duplication, the role of the so-called nuclear plastid DNA sequences (NUPTs) in shaping the evolution of genome architecture and function remains poorly studied. We investigate here the functional and evolutionary fate of NUPTs in the orphan crop Moringa oleifera (moringa), featured by the highest fraction of plastid DNA found so far in any plant genome, focusing on (i) any potential biases in their distribution in relation to specific nuclear genomic features, (ii) their contribution to the emergence of new genes and gene regions, and (iii) their impact on the expression of target nuclear genes.

RESULTS

In agreement with their potential mutagenic effect, NUPTs are underrepresented among structural genes, although their overall transcription levels and broadness were only lower when involved exonic regions; the occurrence of plastid DNA generally did not result in a broader expression, except among those affected in introns by older NUPTs. In contrast, we found a strong enrichment of NUPTs among specific superfamilies of retrotransposons and several classes of RNA genes, including those participating in the protein biosynthetic machinery (i.e., rRNA and tRNA genes) and a specific class of regulatory RNAs. A significant fraction of NUPT RNA genes was found to be functionally expressed, thus potentially contributing to the nuclear pool.

CONCLUSIONS

Our results complete our view of the molecular factors driving the evolution of nuclear genome architecture and function, and support plastid DNA in moringa as a major source of (i) genome complexity and (ii) the nuclear pool of RNA genes.

Comprehensive analysis of the mitochondrial genome of Rehmannia glutinosa: insights into repeat-mediated recombinations and RNA editing-induced stop codon acquisition

Rehmannia glutinosa is an economically significant medicinal plant. Yet, the structure and sequence of its mitochondrial genome has not been published, which plays a crucial role in evolutionary analysis and regulating respiratory-related macromolecule synthesis. In this study, the R. glutinosa mitogenome was sequenced employing a combination of Illumina short reads and Nanopore long reads, with subsequent assembly using a hybrid strategy. We found that the predominant configuration of the R. glutinosa mitogenome comprises two circular chromosomes. The primary structure of the mitogenome encompasses two mitochondrial chromosomes corresponding to the two major configurations, Mac1-1 and Mac1-2. The R. glutinosa mitogenome encoded an angiosperm-typical set of 24 core genes, nine variable genes, three rRNA genes, and 15 tRNA genes. A phylogenetic analysis using the 16 shared protein-coding genes (PCG) yielded a tree consistent with the phylogeny of Lamiales species and two outgroup taxa. Mapping RNA-seq data to the coding sequences (CDS) of the PCGs revealed 507 C-to-U RNA editing sites across 31 PCGs of the R. glutinosa mitogenome. Furthermore, one start codon (nad4L) and two stop codons (rpl10 and atp6) were identified as products of RNA editing events in the R. glutinosa mitogenome.

Overnutrition is a risk factor for iron, but not for zinc or vitamin A deficiency in children and young people: a systematic review and meta-analysis

INTRODUCTION

Traditionally associated with undernutrition, increasing evidence suggests micronutrient deficiencies can coexist with overnutrition. Therefore, this work aimed to systematically review the associations between iron, zinc and vitamin A (VA) status and weight status (both underweight and overweight) in children and young people.

METHODS

Ovid Medline, Ovid Embase, Scopus and Cochrane databases were systematically searched for observational studies assessing micronutrient status (blood, serum or plasma levels of iron, zinc or VA biomarkers) and weight status (body mass index or other anthropometric measurement) in humans under 25 years of any ethnicity and gender. Risk of bias assessment was conducted using the American Dietetic Association Quality Criteria Checklist. Where possible, random effects restricted maximum likelihood meta-analyses were performed.

RESULTS

After screening, 83 observational studies involving 190 443 participants from 44 countries were identified, with many studies having reported on more than one micronutrient and/or weight status indicator. Iron was the most investigated micronutrient, with 46, 28 and 27 studies reporting data for iron, zinc and VA status, respectively. Synthesising 16 records of OR from seven eligible studies, overnutrition (overweight and obesity) increased odds of iron deficiency (ID) (OR (95% CI): 1.51 (1.20 to 1.82), p<0.0001, I2=40.7%). Odds appeared to be higher for children living with obesity (1.88 (1.33 to 2.43), p<0.0001, I2=20.6%) in comparison to those with overweight (1.31 (0.98 to 1.64), p<0.0001, I2=40.5%), although between group differences were not significant (p=0.08).

CONCLUSIONS

Overnutrition is associated with increased risk of ID, but not zinc or VA deficiencies, with an inverted U-shaped relationship observed between iron status and bodyweight. Our results highlight significant heterogeneity in the reporting of micronutrient biomarkers and how deficiencies were defined. Inflammation status was rarely adequately accounted for, and the burden of ID may well be under-recognised, particularly in children and young people living with overnutrition.

PROSPERO REGISTRATION NUMBER

CRD42020221523.

Evolutionary trajectory of organelle-derived nuclear DNAs in the Triticum/Aegilops complex species

Organelle-derived nuclear DNAs, nuclear plastid DNAs (NUPTs), and nuclear mitochondrial DNAs (NUMTs) have been identified in plants. Most, if not all, genes residing in NUPTs/NUMTs (NUPGs/NUMGs) are known to be inactivated and pseudogenized. However, the role of epigenetic control in silencing NUPGs/NUMGs and the dynamic evolution of NUPTs/NUMTs with respect to organismal phylogeny remain barely explored. Based on the available nuclear and organellar genomic resources of wheat (genus Triticum) and goat grass (genus Aegilops) within Triticum/Aegilops complex species, we investigated the evolutionary fates of NUPTs/NUMTs in terms of their epigenetic silencing and their dynamic occurrence rates in the nuclear diploid genomes and allopolyploid subgenomes. NUPTs and NUMTs possessed similar genomic atlas, including (i) predominantly located in intergenic regions and preferential integration to gene regulation regions and (ii) generating sequence variations in the nuclear genome. Unlike nuclear indigenous genes, the alien NUPGs/NUMGs were associated with repressive epigenetic signals, namely high levels of DNA methylation and low levels of active histone modifications. Phylogenomic analyses suggested that the species-specific and gradual accumulation of NUPTs/NUMTs accompanied the speciation processes. Moreover, based on further pan-genomic analyses, we found significant subgenomic asymmetry in the NUPT/NUMT occurrence, which accumulated during allopolyploid wheat evolution. Our findings provide insight into the dynamic evolutionary fates of organelle-derived nuclear DNA in plants.

Chromosome-Level Genome Assembly for the Angiosperm Silene conica

The angiosperm genus Silene has been the subject of extensive study in the field of ecology and evolution, but the availability of high-quality reference genome sequences has been limited for this group. Here, we report a chromosome-level assembly for the genome of Silene conica based on Pacific Bioscience HiFi, Hi-C, and Bionano technologies. The assembly produced 10 scaffolds (1 per chromosome) with a total length of 862 Mb and only ∼1% gap content. These results confirm previous observations that S. conica and its relatives have a reduced base chromosome number relative to the genus's ancestral state of 12. Silene conica has an exceptionally large mitochondrial genome (>11 Mb), predominantly consisting of sequence of unknown origins. Analysis of shared sequence content suggests that it is unlikely that transfer of nuclear DNA is the primary driver of this mitochondrial genome expansion. More generally, this assembly should provide a valuable resource for future genomic studies in Silene, including comparative analyses with related species that recently evolved sex chromosomes.

Chromosome-level genome assembly for the angiosperm Silene conica

The angiosperm genus Silene has been the subject of extensive study in the field of ecology and evolution, but the availability of high-quality reference genome sequences has been limited for this group. Here, we report a chromosome-level assembly for the genome of Silene conica based on PacBio HiFi, Hi-C and Bionano technologies. The assembly produced 10 scaffolds (one per chromosome) with a total length of 862 Mb and only ~1% gap content. These results confirm previous observations that S. conica and its relatives have a reduced base chromosome number relative to the genus's ancestral state of 12. Silene conica has an exceptionally large mitochondrial genome (>11 Mb), predominantly consisting of sequence of unknown origins. Analysis of shared sequence content suggests that it is unlikely that transfer of nuclear DNA is the primary driver of this mitochondrial genome expansion. More generally, this assembly should provide a valuable resource for future genomic studies in Silene, including comparative analyses with related species that recently evolved sex chromosomes.

Significance

Whole-genome sequences have been largely lacking for species in the genus Silene even though these flowering plants have been used for studying ecology, evolution, and genetics for over a century. Here, we address this gap by providing a high-quality nuclear genome assembly for S. conica, a species known to have greatly accelerated rates of sequence and structural divergence in its mitochondrial and plastid genomes. This resource will be valuable in understanding the coevolutionary interactions between nuclear and cytoplasmic genomes and in comparative analyses across this highly diverse genus.

Pangenome-based trajectories of intracellular gene transfers in Poaceae unveil high cumulation in Triticeae

Intracellular gene transfers (IGTs) between the nucleus and organelles, including plastids and mitochondria, constantly reshape the nuclear genome during evolution. Despite the substantial contribution of IGTs to genome variation, the dynamic trajectories of IGTs at the pangenomic level remain elusive. Here, we developed an approach, IGTminer, that maps the evolutionary trajectories of IGTs using collinearity and gene reannotation across multiple genome assemblies. We applied IGTminer to create a nuclear organellar gene (NOG) map across 67 genomes covering 15 Poaceae species, including important crops. The resulting NOGs were verified by experiments and sequencing data sets. Our analysis revealed that most NOGs were recently transferred and lineage specific and that Triticeae species tended to have more NOGs than other Poaceae species. Wheat (Triticum aestivum) had a higher retention rate of NOGs than maize (Zea mays) and rice (Oryza sativa), and the retained NOGs were likely involved in photosynthesis and translation pathways. Large numbers of NOG clusters were aggregated in hexaploid wheat during 2 rounds of polyploidization, contributing to the genetic diversity among modern wheat accessions. We implemented an interactive web server to facilitate the exploration of NOGs in Poaceae. In summary, this study provides resources and insights into the roles of IGTs in shaping interspecies and intraspecies genome variation and driving plant genome evolution.

Comparative analysis of organellar genomes between diploid and tetraploid Chrysanthemum indicum with its relatives

Chrysanthemum indicum, a species native to Eastern Asia is well known as one of the progenitor species of the cultivated Chrysanthemum which is grown for its ornamental and medicinal value. Previous genomic studies on Chrysanthemum have largely ignored the dynamics of plastid genome (plastome) and mitochondria genome (mitogenome) evolution when analyzing this plant lineage. In this study, we sequenced and assembled the plastomes and mitogenomes of diploid and tetraploid C. indicum as well as the morphologically divergent variety C. indicum var. aromaticum. We used published data from 27 species with both plastome and mitogenome complete sequences to explore differences in sequence evolution between the organellar genomes. The size and structure of organellar genome between diploid and tetraploid C. indicum were generally similar but the tetraploid C. indicum and C. indicum var. aromaticum were found to contain unique sequences in the mitogenomes which also contained previously undescribed open reading frames (ORFs). Across Chrysanthemum mitogenome structure varied greatly but sequences transferred from plastomes in to the mitogenomes were conserved. Finally, differences observed between mitogenome and plastome gene trees may be the result of the difference in the rate of sequence evolution between genes in these two genomes. In total the findings presented here greatly expand the resources for studying Chrysanthemum organellar genome evolution with possible applications to conservation, breeding, and gene banking in the future.

CRISPR-Cas9-based repeat depletion for high-throughput genotyping of complex plant genomes

High-throughput genotyping enables the large-scale analysis of genetic diversity in population genomics and genome-wide association studies that combine the genotypic and phenotypic characterization of large collections of accessions. Sequencing-based approaches for genotyping are progressively replacing traditional genotyping methods because of the lower ascertainment bias. However, genome-wide genotyping based on sequencing becomes expensive in species with large genomes and a high proportion of repetitive DNA. Here we describe the use of CRISPR-Cas9 technology to deplete repetitive elements in the 3.76-Gb genome of lentil (Lens culinaris), 84% consisting of repeats, thus concentrating the sequencing data on coding and regulatory regions (single-copy regions). We designed a custom set of 566,766 gRNAs targeting 2.9 Gbp of repeats and excluding repetitive regions overlapping annotated genes and putative regulatory elements based on ATAC-seq data. The novel depletion method removed ∼40% of reads mapping to repeats, increasing those mapping to single-copy regions by ∼2.6-fold. When analyzing 25 million fragments, this repeat-to-single-copy shift in the sequencing data increased the number of genotyped bases of ∼10-fold compared to nondepleted libraries. In the same condition, we were also able to identify ∼12-fold more genetic variants in the single-copy regions and increased the genotyping accuracy by rescuing thousands of heterozygous variants that otherwise would be missed because of low coverage. The method performed similarly regardless of the multiplexing level, type of library or genotypes, including different cultivars and a closely related species (L. orientalis). Our results showed that CRISPR-Cas9-driven repeat depletion focuses sequencing data on single-copy regions, thus improving high-density and genome-wide genotyping in large and repetitive genomes.

Polymeric structure of the Cannabis sativa L. mitochondrial genome identified with an assembly graph model

Cannabis sativa L. belongs to the family Cannabaceae in Rosales. It has been widely used as medicines, building materials, and textiles. Elucidating its genome is critical for molecular breeding and synthetic biology study. Many studies have shown that the mitochondrial genomes (mitogenomes) and even chloroplast genomes (plastomes) had complex polymeric structures. Using the Nanopore sequencing platform, we sequenced, assembled, and analyzed its mitogenome and plastome. The resulting unitig graph suggested that the mitogenome had a complex polymeric structure. However, a gap-free, circular sequence was further assembled from the unitig graph. In contrast, a circular sequence representing the plastome was obtained. The mitogenome major conformation was 415,837 bp long, and the plastome was 153,927 bp long. To test if the repeat sequences promote recombination, which corresponds to the branch points in the structure, we tested the sequences around repeats by long-read mapping. Among 208 pairs of predicted repeats, the mapping results supported the presence of cross-over around 25 pairs of repeats. Subsequent PCR amplification confirmed the presence of cross-over around 15 of the 25 repeats. By comparing the mitogenome and plastome sequences, we identified 19 mitochondria plastid DNAs, including seven complete genes (trnW-CCA, trnP-UGG, psbJ, trnN-GUU, trnD-GUC, trnH-GUG, trnM-CAU) and nine gene fragments. Furthermore, the selective pressure analysis results showed that five genes (atp1, ccmB, ccmC, cox1, nad7) had 19 positively selected sites. Lastly, we predicted 28 RNA editing sites. A total of 8 RNA editing sites located in the coding regions were successfully validated by PCR amplification and Sanger sequencing, of which four were synonymous, and four were nonsynonymous. In particular, the RNA editing events appeared to be tissue-specific in C. sativa mitogenome. In summary, we have confirmed the major confirmation of C. sativa mitogenome and characterized its structural features in detail. These results provide critical information for future variety breeding and resource development for C. sativa.

Variation in mitogenome structural conformation in wild and cultivated lineages of sorghum corresponds with domestication history and plastome evolution

BACKGROUND

Mitochondria are organelles within eukaryotic cells that are central to the metabolic processes of cellular respiration and ATP production. However, the evolution of mitochondrial genomes (mitogenomes) in plants is virtually unknown compared to animal mitogenomes or plant plastids, due to complex structural variation and long stretches of repetitive DNA making accurate genome assembly more challenging. Comparing the structural and sequence differences of organellar genomes within and between sorghum species is an essential step in understanding evolutionary processes such as organellar sequence transfer to the nuclear genome as well as improving agronomic traits in sorghum related to cellular metabolism.

RESULTS

Here, we assembled seven sorghum mitochondrial and plastid genomes and resolved reticulated mitogenome structures with multilinked relationships that could be grouped into three structural conformations that differ in the content of repeats and genes by contig. The grouping of these mitogenome structural types reflects the two domestication events for sorghum in east and west Africa.

CONCLUSIONS

We report seven mitogenomes of sorghum from different cultivars and wild sources. The assembly method used here will be helpful in resolving complex genomic structures in other plant species. Our findings give new insights into the structure of sorghum mitogenomes that provides an important foundation for future research into the improvement of sorghum traits related to cellular respiration, cytonuclear incompatibly, and disease resistance.

A genome resource for Acacia, Australia's largest plant genus

Acacia (Leguminosae, Caesalpinioideae, mimosoid clade) is the largest and most widespread genus of plants in the Australian flora, occupying and dominating a diverse range of environments, with an equally diverse range of forms. For a genus of its size and importance, Acacia currently has surprisingly few genomic resources. Acacia pycnantha, the golden wattle, is a woody shrub or tree occurring in south-eastern Australia and is the country's floral emblem. To assemble a genome for A. pycnantha, we generated long-read sequences using Oxford Nanopore Technology, 10x Genomics Chromium linked reads, and short-read Illumina sequences, and produced an assembly spanning 814 Mb, with a scaffold N50 of 2.8 Mb, and 98.3% of complete Embryophyta BUSCOs. Genome annotation predicted 47,624 protein-coding genes, with 62.3% of the genome predicted to comprise transposable elements. Evolutionary analyses indicated a shared genome duplication event in the Caesalpinioideae, and conflict in the relationships between Cercis (subfamily Cercidoideae) and subfamilies Caesalpinioideae and Papilionoideae (pea-flowered legumes). Comparative genomics identified a suite of expanded and contracted gene families in A. pycnantha, and these were annotated with both GO terms and KEGG functional categories. One expanded gene family of particular interest is involved in flowering time and may be associated with the characteristic synchronous flowering of Acacia. This genome assembly and annotation will be a valuable resource for all studies involving Acacia, including the evolution, conservation, breeding, invasiveness, and physiology of the genus, and for comparative studies of legumes.

Evaluation of Intracellular Gene Transfers from Plastome to Nuclear Genome across Progressively Improved Assemblies for Arabidopsis thaliana and Oryza sativa

DNA originating from organellar genomes are regularly discovered in nuclear sequences during genome assembly. Nevertheless, such insertions are sometimes omitted during the process of nuclear genome assembly because the inserted DNA is assigned to organellar genomes, leading to a systematic underestimation of their frequency. With the rapid development of high-throughput sequencing technology, more inserted fragments from organelle genomes can now be detected. Therefore, it is necessary to be aware of the insertion events from organellar genomes during nuclear genome assembly to properly attribute the impact and rate of such insertions in the evolution of nuclear genomes. Here, we investigated the impact of intracellular gene transfer (IGT) from the plastome to the nuclear genome using genome assemblies that were refined through time with technological improvements from two model species, Arabidopsis thaliana and Oryza sativa. We found that IGT from the plastome to the nuclear genome is a dynamic and ongoing process in both A. thaliana and O. sativa, and mostly occurred recently, as the majority of transferred sequences showed over 95% sequence similarity with plastome sequences of origin. Differences in the plastome-to-nuclear genome IGT between A. thaliana and O. sativa varied among the different assembly versions and were associated with the quality of the nuclear genome assembly. IGTs from the plastome to nuclear genome occurred more frequently in intergenic regions, which were often associated with transposable elements (TEs). This study provides new insights into intracellular genome evolution and nuclear genome assembly by characterizing and comparing IGT from the plastome into the nuclear genome for two model plant species.

Recovery of chloroplast genomes from medieval millet grains excavated from the Areni-1 cave in southern Armenia

Panicum miliaceum L. was domesticated in northern China at least 7000 years ago and was subsequentially adopted in many areas throughout Eurasia. One such locale is Areni-1 an archaeological cave site in Southern Armenia, where vast quantities archaeobotanical material were well preserved via desiccation. The rich botanical material found at Areni-1 includes P. miliaceum grains that were identified morphologically and¹⁴C dated to the medieval period (873 ± 36 CE and 1118 ± 35 CE). To investigate the demographic and evolutionary history of the Areni-1 millet, we used ancient DNA extraction, hybridization capture enrichment, and high throughput sequencing to assemble three chloroplast genomes from the medieval grains and then compared these sequences to 50 modern P. miliaceum chloroplast genomes. Overall, the chloroplast genomes contained a low amount of diversity with domesticated accessions separated by a maximum of 5 SNPs and little inference on demography could be made. However, in phylogenies the chloroplast genomes separated into two clades, similar to what has been reported for nuclear DNA from P. miliaceum. The chloroplast genomes of two wild (undomesticated) accessions of P. miliaceum contained a relatively large number of variants, 11 SNPs, not found in the domesticated accessions. These results demonstrate that P. miliaceum grains from archaeological sites can preserve DNA for at least 1000 years and serve as a genetic resource to study the domestication of this cereal crop.

Targeted base editing in the mitochondrial genome of Arabidopsis thaliana

SignificanceThe mitochondrial genomes of land plants encode genes for cellular energy production and agriculturally important traits, but modification of the genomes is still difficult. Targeted base editing is one of the best ways to modify genes and intergenic regions and thus understand their functions, without drastically changing genome structure. In this study, we succeeded in creating plantlets of the model plant Arabidopsis thaliana, in which all of the many copies of the mitochondrial genomes in each cell had a targeted C:G base pair converted to a T:A pair. Introduced mutations were stably inherited by the next generation. This method will help to unravel the mysteries of plant mitochondrial genomes and may also serve as a basis for increasing crop yields.

Born in the mitochondrion and raised in the nucleus: evolution of a novel tandem repeat family in Medicago polymorpha (Fabaceae)

Plant nuclear genomes harbor sequence elements derived from the organelles (mitochondrion and plastid) through intracellular gene transfer (IGT). Nuclear genomes also show a dramatic range of repeat content, suggesting that any sequence can be readily amplified. These two aspects of plant nuclear genomes are well recognized but have rarely been linked. Through investigation of 31 Medicago taxa we detected exceptionally high post-IGT amplification of mitochondrial (mt) DNA sequences containing rps10 in the nuclear genome of Medicago polymorpha and closely related species. The amplified sequences were characterized as tandem arrays of five distinct repeat motifs (2157, 1064, 987, 971, and 587 bp) that have diverged from the mt genome (mitogenome) in the M. polymorpha nuclear genome. The mt rps10-like arrays were identified in seven loci (six intergenic and one telomeric) of the nuclear chromosome assemblies and were the most abundant tandem repeat family, representing 1.6-3.0% of total genomic DNA, a value approximately three-fold greater than the entire mitogenome in M. polymorpha. Compared to a typical mt gene, the mt rps10-like sequence coverage level was 691.5-7198-fold higher in M. polymorpha and closely related species. In addition to the post-IGT amplification, our analysis identified the canonical telomeric repeat and the species-specific satellite arrays that are likely attributable to an ancestral chromosomal fusion in M. polymorpha. A possible relationship between chromosomal instability and the mt rps10-like tandem repeat family in the M. polymorpha clade is discussed.

A Chromosome-Level Genome Assembly of the European Beech (Fagus sylvatica) Reveals Anomalies for Organelle DNA Integration, Repeat Content and Distribution of SNPs

The European Beech is the dominant climax tree in most regions of Central Europe and valued for its ecological versatility and hardwood timber. Even though a draft genome has been published recently, higher resolution is required for studying aspects of genome architecture and recombination. Here, we present a chromosome-level assembly of the more than 300 year-old reference individual, Bhaga, from the Kellerwald-Edersee National Park (Germany). Its nuclear genome of 541 Mb was resolved into 12 chromosomes varying in length between 28 and 73 Mb. Multiple nuclear insertions of parts of the chloroplast genome were observed, with one region on chromosome 11 spanning more than 2 Mb which fragments up to 54,784 bp long and covering the whole chloroplast genome were inserted randomly. Unlike in Arabidopsis thaliana, ribosomal cistrons are present in Fagus sylvatica only in four major regions, in line with FISH studies. On most assembled chromosomes, telomeric repeats were found at both ends, while centromeric repeats were found to be scattered throughout the genome apart from their main occurrence per chromosome. The genome-wide distribution of SNPs was evaluated using a second individual from Jamy Nature Reserve (Poland). SNPs, repeat elements and duplicated genes were unevenly distributed in the genomes, with one major anomaly on chromosome 4. The genome presented here adds to the available highly resolved plant genomes and we hope it will serve as a valuable basis for future research on genome architecture and for understanding the past and future of European Beech populations in a changing climate.

A dated phylogeny shows Plio-Pleistocene climates spurred evolution of antibrowsing defences in the New Zealand flora

Some plant traits may be legacies of coevolution with extinct megafauna. One example is the convergent evolution of 'divaricate' cage architectures in many New Zealand lineages, interpreted as a response to recently extinct flightless avian browsers whose ancestors arrived during the Paleogene period. Although experiments have confirmed that divaricate habit deters extant browsers, its abundance on frosty, droughty sites appears consistent with an earlier interpretation as a response to cold, dry Plio-Pleistocene climates. We used 45 protein-coding sequences from plastid genomes to reconstruct the evolutionary history of the divaricate habit in extant New Zealand lineages. Our dated phylogeny of 215 species included 91% of New Zealand eudicot divaricate species. We show that 86% of extant divaricate plants diverged from non-divaricate sisters within the last 5 Ma, implicating Plio-Pleistocene climates in the proliferation of cage architectures in New Zealand. Our results, combined with other recent findings, are consistent with the synthetic hypothesis that the browser-deterrent effect of cage architectures was strongly selected only when Plio-Pleistocene climatic constraints prevented woody plants from growing quickly out of reach of browsers. This is consistent with the abundance of cage architectures in other regions where plant growth is restricted by aridity or short frost-free periods.

A Critical Review of the Current Global Ex Situ Conservation System for Plant Agrobiodiversity. II. Strengths and Weaknesses of the Current System and Recommendations for Its Improvement

In this paper, we review gene bank operations that have an influence on the global conservation system, with the intention to identify critical aspects that should be improved for optimum performance. We describe the role of active and base collections and the importance of linking germplasm conservation and use, also in view of new developments in genomics and phenomics that facilitate more effective and efficient conservation and use of plant agrobiodiversity. Strengths, limitations, and opportunities of the existing global ex situ conservation system are discussed, and measures are proposed to achieve a rational, more effective, and efficient global system for germplasm conservation and sustainable use. The proposed measures include filling genetic and geographic gaps in current ex situ collections; determining unique accessions at the global level for long-term conservation in virtual base collections; intensifying existing international collaborations among gene banks and forging collaborations with the botanic gardens community; increasing investment in conservation research and user-oriented supportive research; improved accession-level description of the genetic diversity of crop collections; improvements of the legal and policy framework; and oversight of the proposed network of global base collections.

Genomic landscape of mitochondrial DNA insertions in 23 bat genomes: characteristics, loci, phylogeny, and polymorphism

The transfer of mitochondrial DNA to the nuclear genome gives rise to the nuclear DNA sequences of mitochondrial origin (NUMTs), considered as a driving force in genome evolution. In this study, NUMTs in 23 bat genomes were investigated and compared systematically. The results showed that NUMTs existed in 22 genomes except for Noctilio leporinus, suggesting that mitochondrial fragment insertion in the nuclear genome was a common event in bat genomes. However, remarkable variations in NUMTs number, cumulative length, and proportion in the nuclear genome were discovered across bat species. Further orthologous NUMT loci analysis of the Phyllostomidae family indicated that the NUMTs insertion events in bat genomes were homoplasy-free. The NUMTs were mainly inserted into the intergenic regions, particularly, co-localized with repetitive sequences (especially transposable elements). However, several NUMTs were inserted into genes, some of which were in the exon region of functional genes. One NUMT in the genome of Pteropus alecto surprisingly matched with cDNA of ATP8B3 that provided evidence of NUMTs with coding function. Phylogenic analysis on NUMTs originating from COXI and COXII loci highlighted 2 clusters of Yinpterochiroptera and Yangochiroptera for Chiroptera. Seven NUMTs from Rhinolophus ferrumequinum were amplified, and the sequencing results confirmed the reliability of the NUMT analysis. One of them was polymorphic for the presence or absence of the NUMT insertion, and each genotype of NUMT loci showed a distinct regional distribution pattern. The information obtained in this study provides novel insights into the NUMT organization and features in bat genomes and establishes a basis for further studying of the evolution of bat species.

An Insight Into the Mechanism of Plant Organelle Genome Maintenance and Implications of Organelle Genome in Crop Improvement: An Update

Besides the nuclear genome, plants possess two small extra chromosomal genomes in mitochondria and chloroplast, respectively, which contribute a small fraction of the organelles’ proteome. Both mitochondrial and chloroplast DNA have originated endosymbiotically and most of their prokaryotic genes were either lost or transferred to the nuclear genome through endosymbiotic gene transfer during the course of evolution. Due to their immobile nature, plant nuclear and organellar genomes face continuous threat from diverse exogenous agents as well as some reactive by-products or intermediates released from various endogenous metabolic pathways. These factors eventually affect the overall plant growth and development and finally productivity. The detailed mechanism of DNA damage response and repair following accumulation of various forms of DNA lesions, including single and double-strand breaks (SSBs and DSBs) have been well documented for the nuclear genome and now it has been extended to the organelles also. Recently, it has been shown that both mitochondria and chloroplast possess a counterpart of most of the nuclear DNA damage repair pathways and share remarkable similarities with different damage repair proteins present in the nucleus. Among various repair pathways, homologous recombination (HR) is crucial for the repair as well as the evolution of organellar genomes. Along with the repair pathways, various other factors, such as the MSH1 and WHIRLY family proteins, WHY1, WHY2, and WHY3 are also known to be involved in maintaining low mutation rates and structural integrity of mitochondrial and chloroplast genome. SOG1, the central regulator in DNA damage response in plants, has also been found to mediate endoreduplication and cell-cycle progression through chloroplast to nucleus retrograde signaling in response to chloroplast genome instability. Various proteins associated with the maintenance of genome stability are targeted to both nuclear and organellar compartments, establishing communication between organelles as well as organelles and nucleus. Therefore, understanding the mechanism of DNA damage repair and inter compartmental crosstalk mechanism in various sub-cellular organelles following induction of DNA damage and identification of key components of such signaling cascades may eventually be translated into strategies for crop improvement under abiotic and genotoxic stress conditions. This review mainly highlights the current understanding as well as the importance of different aspects of organelle genome maintenance mechanisms in higher plants.

Organelle Genetics in Plants

Eleven published articles (4 reviews, 7 research papers) are collected in the Special Issue entitled “Organelle Genetics in Plants.” This selection of papers covers a wide range of topics related to chloroplasts and plant mitochondria research: (i) organellar gene expression (OGE) and, more specifically, chloroplast RNA editing in soybean, mitochondria RNA editing, and intron splicing in soybean during nodulation, as well as the study of the roles of transcriptional and posttranscriptional regulation of OGE in plant adaptation to environmental stress; (ii) analysis of the nuclear integrants of mitochondrial DNA (NUMTs) or plastid DNA (NUPTs); (iii) sequencing and characterization of mitochondrial and chloroplast genomes; (iv) recent advances in plastid genome engineering. Here we summarize the main findings of these works, which represent the latest research on the genetics, genomics, and biotechnology of chloroplasts and mitochondria.

Mitochondrial Fostering: The Mitochondrial Genome May Play a Role in Plant Orphan Gene Evolution

Plant mitochondrial genomes exhibit unique evolutionary patterns. They have a high rearrangement but low mutation rate, and a large size. Based on massive mitochondrial DNA transfers to the nucleus as well as the mitochondrial unique evolutionary traits, we propose a "Mitochondrial Fostering" theory where the organelle genome plays an integral role in the arrival and development of orphan genes (genes with no homologs in other lineages). Two approaches were used to test this theory: (1) bioinformatic analysis of nuclear mitochondrial DNA (Numts: mitochondrial originating DNA that migrated to the nucleus) at the genome level, and (2) bioinformatic analysis of particular orphan sequences present in both the mitochondrial genome and the nuclear genome of Arabidopsis thaliana. One study example is given about one orphan sequence that codes for two unique orphan genes: one in the mitochondrial genome and another one in the nuclear genome. DNA alignments show regions of this A. thaliana orphan sequence exist scattered throughout other land plant mitochondrial genomes. This is consistent with the high recombination rates of mitochondrial genomes in land plants. This may also enable the creation of novel coding sequences within the orphan loci, which can then be transferred to the nuclear genome and become exposed to new evolutionary pressures. Our study also reveals a high correlation between the amount of mitochondrial DNA transferred to the nuclear genome and the number of orphan genes in land plants. All the data suggests the mitochondrial genome may play a role in nuclear orphan gene evolution in land plants.

Whole-genome de novo assemblies reveal extensive structural variations and dynamic organelle-to-nucleus DNA transfers in African and Asian rice

Asian cultivated rice (Oryza sativa) and African cultivated rice (Oryza glaberrima) originated from the wild rice species Oryza rufipogon and Oryza barthii, respectively. The genomes of both cultivated species have undergone profound changes during domestication. Whole-genome de novo assemblies of O. barthii, O. glaberrima, O. rufipogon and Oryza nivara, produced using PacBio single-molecule real-time (SMRT) and next-generation sequencing (NGS) technologies, showed that Gypsy-like retrotransposons are the major contributors to genome size variation in African and Asian rice. Through the detection of genome-wide structural variations (SVs), we observed that besides 28 shared SV hot spots, another 67 hot spots existed in either the Asian or African rice genomes. Based on gene annotation information of the SVs, we established that organelle-to-nucleus DNA transfers resulted in numerous SVs that participated in the nuclear genome divergence of rice species and subspecies. We detected 52 giant nuclear integrants of organelle DNA (NORGs, defined as >10 kb) in six Oryza AA genomes. In addition, we developed an effective method to genotype giant NORGs, based on genome assembly, and first showed the dynamic change in the distribution of giant NORGs in rice natural population. Interestingly, 16 highly differentiated giant NORGs tended to accumulate in natural populations of Asian rice from higher latitude regions, grown at lower temperatures and light intensities. Our study provides new insight into the genome divergence of African and Asian rice, and establishes that organelle-to-nucleus DNA transfers, as potentially powerful contributors to environmental adaptation during rice evolution, play a major role in producing SVs in rice genomes.

Evolutionary analysis of the Moringa oleifera genome reveals a recent burst of plastid to nucleus gene duplications

It is necessary to identify suitable alternative crops to ensure the nutritional demands of a growing global population. The genome of Moringa oleifera, a fast-growing drought-tolerant orphan crop with highly valuable agronomical, nutritional and pharmaceutical properties, has recently been reported. We model here gene family evolution in Moringa as compared with ten other flowering plant species. Despite the reduced number of genes in the compact Moringa genome, 101 gene families, grouping 957 genes, were found as significantly expanded. Expanded families were highly enriched for chloroplastidic and photosynthetic functions. Indeed, almost half of the genes belonging to Moringa expanded families grouped with their Arabidopsis thaliana plastid encoded orthologs. Microsynteny analysis together with modeling the distribution of synonymous substitutions rates, supported most plastid duplicated genes originated recently through a burst of simultaneous insertions of large regions of plastid DNA into the nuclear genome. These, together with abundant short insertions of plastid DNA, contributed to the occurrence of massive amounts of plastid DNA in the Moringa nuclear genome, representing 4.71%, the largest reported so far. Our study provides key genetic resources for future breeding programs and highlights the potential of plastid DNA to impact the structure and function of nuclear genes and genomes.