Mitochondrial genome and transcriptome analysis of five alloplasmic male-sterile lines in Brassica juncea

Mitochondrial genome comparison and phylogenetic analysis of Dendrobium (Orchidaceae) based on whole mitogenomes

BACKGROUND

Mitochondrial genomes are essential for deciphering the unique evolutionary history of seed plants. However, the rules of their extreme variation in genomic size, multi-chromosomal structure, and foreign sequences remain unresolved in most plant lineages, which further hindered the application of mitogenomes in phylogenetic analyses.

RESULTS

Here, we took Dendrobium (Orchidaceae) which shows the great divergence of morphology and difficulty in species taxonomy as the study focus. We first de novo assembled two complete mitogenomes of Dendrobium wilsonii and Dendrobium henanense that were 763,005 bp and 807,551 bp long with multichromosomal structures. To understand the evolution of Dendrobium mitogenomes, we compared them with those of four other orchid species. The results showed great variations of repetitive and chloroplast-derived sequences in Dendrobium mitogenomes. Moreover, the intergenic content of Dendrobium mitogenomes has undergone expansion during evolution. We also newly sequenced mitogenomes of 26 Dendrobium species and reconstructed phylogenetic relationships of Dendrobium based on genomic mitochondrial and plastid data. The results indicated that the existence of chloroplast-derived sequences made the mitochondrial phylogeny display partial characteristics of the plastid phylogeny. Additionally, the mitochondrial phylogeny provided new insights into the phylogenetic relationships of Dendrobium species.

CONCLUSIONS

Our study revealed the evolution of Dendrobium mitogenomes and the potential of mitogenomes in deciphering phylogenetic relationships at low taxonomic levels.

Comparative Transcriptome Analysis Reveals a Potential Regulatory Network for Ogura Cytoplasmic Male Sterility in Cabbage (Brassica oleracea L.)

Ogura cytoplasmic male sterility (CMS) lines are widely used breeding materials in cruciferous crops and play important roles in heterosis utilization; however, the sterility mechanism remains unclear. To investigate the microspore development process and gene expression changes after the introduction of orf138 and Rfo, cytological observation and transcriptome analysis were performed using a maintainer line, an Ogura CMS line, and a restorer line. Semithin sections of microspores at different developmental stages showed that the degradation of tapetal cells began at the tetrad stage in the Ogura CMS line, while it occurred at the bicellular microspore stage to the tricellular microspore stage in the maintainer and restorer lines. Therefore, early degradation of tapetal cells may be the cause of pollen abortion. Transcriptome analysis results showed that a total of 1287 DEGs had consistent expression trends in the maintainer line and restorer line, but were significantly up- or down-regulated in the Ogura CMS line, indicating that they may be closely related to pollen abortion. Functional annotation showed that the 1287 core DEGs included a large number of genes related to pollen development, oxidative phosphorylation, carbohydrate, lipid, and protein metabolism. In addition, further verification elucidated that down-regulated expression of genes related to energy metabolism led to decreased ATP content and excessive ROS accumulation in the anthers of Ogura CMS. Based on these results, we propose a transcriptome-mediated induction and regulatory network for cabbage Ogura CMS. Our research provides new insights into the mechanism of pollen abortion and fertility restoration in Ogura CMS.

Complete mitochondrial genomes of three Mangifera species, their genomic structure and gene transfer from chloroplast genomes

Background

Among the Mangifera species, mango (Mangifera indica) is an important commercial fruit crop. However, very few studies have been conducted on the Mangifera mitochondrial genome. This study reports and compares the newly sequenced mitochondrial genomes of three Mangifera species.

Results

Mangifera mitochondrial genomes showed partial similarities in the overall size, genomic structure, and gene content. Specifically, the genomes are circular and contain about 63–69 predicted functional genes, including five ribosomal RNA (rRNA) genes and 24–27 transfer RNA (tRNA) genes. The GC contents of the Mangifera mitochondrial genomes are similar, ranging from 44.42–44.66%. Leucine (Leu) and serine (Ser) are the most frequently used, while tryptophan (Trp) and cysteine (Cys) are the least used amino acids among the protein-coding genes in Mangifera mitochondrial genomes. We also identified 7–10 large chloroplast genomic fragments in the mitochondrial genome, ranging from 1407 to 6142 bp. Additionally, four intact mitochondrial tRNAs genes (tRNA-Cys, tRNA-Trp, tRNA-Pro, and tRNA-Met) and intergenic spacer regions were identified. Phylogenetic analysis based on the common protein-coding genes of most branches provided a high support value.

Conclusions

We sequenced and compared the mitochondrial genomes of three Mangifera species. The results showed that the gene content and the codon usage pattern of Mangifera mitochondrial genomes is similar across various species. Gene transfer from the chloroplast genome to the mitochondrial genome were identified. This study provides valuable information for evolutionary and molecular studies of Mangifera and a basis for further studies on genomic breeding of mango.

Organelle genome assembly uncovers the dynamic genome reorganization and cytoplasmic male sterility associated genes in tomato

To identify cytoplasmic male sterility (CMS)-associated genes in tomato, we determined the genome sequences of mitochondria and chloroplasts in three CMS tomato lines derived from independent asymmetric cell fusions, their nuclear and cytoplasmic donors, and male fertile weedy cultivated tomato and wild relatives. The structures of the CMS mitochondrial genomes were highly divergent from those of the nuclear and cytoplasmic donors, and genes of the donors were mixed up in these genomes. On the other hand, the structures of CMS chloroplast genomes were moderately conserved across the donors, but CMS chloroplast genes were unexpectedly likely derived from the nuclear donors. Comparative analysis of the structures and contents of organelle genes and transcriptome analysis identified three genes that were uniquely present in the CMS lines, but not in the donor or fertile lines. RNA-sequencing analysis indicated that these three genes transcriptionally expressed in anther, and identified different RNA editing levels in one gene, orf265, that was partially similar to ATP synthase subunit 8, between fertile and sterile lines. The orf265 was a highly potential candidate for CMS-associated gene. This study suggests that organelle reorganization mechanisms after cell fusion events differ between mitochondria and chloroplasts, and provides insight into the development of new F1 hybrid breeding programs employing the CMS system in tomato.

A complete mitochondrial genome for fragrant Chinese rosewood (Dalbergia odorifera, Fabaceae) with comparative analyses of genome structure and intergenomic sequence transfers

BACKGROUND

Dalbergia odorifera is an economically and culturally important species in the Fabaceae because of the high-quality lumber and traditional Chinese medicines made from this plant, however, overexploitation has increased the scarcity of D. odorifera. Given the rarity and the multiple uses of this species, it is important to expand the genomic resources for utilizing in applications such as tracking illegal logging, determining effective population size of wild stands, delineating pedigrees in marker assisted breeding programs, and resolving gene networks in functional genomics studies. Even the nuclear and chloroplast genomes have been published for D. odorifera, the complete mitochondrial genome has not been assembled or assessed for sequence transfer to other genomic compartments until now. Such work is essential in understanding structural and functional genome evolution in a lineage (Fabaceae) with frequent intergenomic sequence transfers.

RESULTS

We integrated Illumina short-reads and PacBio CLR long-reads to assemble and annotate the complete mitochondrial genome of D. odorifera. The mitochondrial genome was organized as a single circular structure of 435 Kb in length containing 33 protein coding genes, 4 rRNA and 17 tRNA genes. Nearly 4.0% (17,386 bp) of the genome was annotated as repetitive DNA. From the sequence transfer analysis, it was found that 114 Kb of DNA originating from the mitochondrial genome has been transferred to the nuclear genome, with most of the transfer events having taken place relatively recently. The high frequency of sequence transfers from the mitochondria to the nuclear genome was similar to that of sequence transfer from the chloroplast to the nuclear genome.

CONCLUSION

For the first-time, the complete mitochondrial genome of D. odorifera was assembled in this study, which will provide a baseline resource in understanding genomic evolution in the highly specious Fabaceae. In particular, the assessment of intergenomic sequence transfer suggests that transfers have been common and recent indicating a possible role in environmental adaptation as has been found in other lineages. The high turnover rate of genomic colinearly and large differences in mitochondrial genome size found in the comparative analyses herein providing evidence for the rapid evolution of mitochondrial genome structure compared to chloroplasts in Faboideae. While phylogenetic analyses using functional genes indicate that mitochondrial genes are very slowly evolving compared to chloroplast genes.

Transcriptional Landscape and Splicing Efficiency in Arabidopsis Mitochondria

Plant mitochondrial transcription is initiated from multiple promoters without an apparent motif, which precludes their identification in other species based on sequence comparisons. Even though coding regions take up only a small fraction of plant mitochondrial genomes, deep RNAseq studies uncovered that these genomes are fully or nearly fully transcribed with significantly different RNA read depth across the genome. Transcriptomic analysis can be a powerful tool to understand the transcription process in diverse angiosperms, including the identification of potential promoters and co-transcribed genes or to study the efficiency of intron splicing. In this work, we analyzed the transcriptional landscape of the Arabidopsis mitochondrial genome (mtDNA) based on large-scale RNA sequencing data to evaluate the use of RNAseq to study those aspects of the transcription process. We found that about 98% of the Arabidopsis mtDNA is transcribed with highly different RNA read depth, which was elevated in known genes. The location of a sharp increase in RNA read depth upstream of genes matched the experimentally identified promoters. The continuously high RNA read depth across two adjacent genes agreed with the known co-transcribed units in Arabidopsis mitochondria. Most intron-containing genes showed a high splicing efficiency with no differences between cis and trans-spliced introns or between genes with distinct splicing mechanisms. Deep RNAseq analyses of diverse plant species will be valuable to recognize general and lineage-specific characteristics related to the mitochondrial transcription process.

Horizontal gene transfers dominate the functional mitochondrial gene space of a holoparasitic plant

Although horizontal gene transfer (HGT) is common in angiosperm mitochondrial DNAs (mtDNAs), few cases of functional foreign genes have been identified. The one outstanding candidate for large-scale functional HGT is the holoparasite Lophophytum mirabile, whose mtDNA has lost most native genes but contains intact foreign homologs acquired from legume host plants. To investigate the extent to which this situation results from functional replacement of native by foreign genes, functional mitochondrial gene transfer to the nucleus, and/or loss of mitochondrial biochemical function in the context of extreme parasitism, we examined the Lophophytum mitochondrial and nuclear transcriptomes by deep paired-end RNA sequencing. Most foreign mitochondrial genes in Lophophytum are highly transcribed, accurately spliced, and efficiently RNA edited. By contrast, we found no evidence for functional gene transfer to the nucleus or loss of mitochondrial functions in Lophophytum. Many functional replacements occurred via the physical replacement of native genes by foreign genes. Some of these events probably occurred as the final act of HGT itself. Lophophytum mtDNA has experienced an unprecedented level of functional replacement of native genes by foreign copies. This raises important questions concerning population-genetic and molecular regimes that underlie such a high level of foreign gene takeover.

Investigating the mitochondrial genomic landscape of Arabidopsis thaliana by long-read sequencing

Plant mitochondrial genomes have distinctive features compared to those of animals; namely, they are large and divergent, with sizes ranging from hundreds of thousands of to a few million bases. Recombination among repetitive regions is thought to produce similar structures that differ slightly, known as "multipartite structures," which contribute to different phenotypes. Although many reference plant mitochondrial genomes represent almost all the genes in mitochondria, the full spectrum of their structures remains largely unknown. The emergence of long-read sequencing technology is expected to yield this landscape; however, many studies aimed to assemble only one representative circular genome, because properly understanding multipartite structures using existing assemblers is not feasible. To elucidate multipartite structures, we leveraged the information in existing reference genomes and classified long reads according to their corresponding structures. We developed a method that exploits two classic algorithms, partial order alignment (POA) and the hidden Markov model (HMM) to construct a sensitive read classifier. This method enables us to represent a set of reads as a POA graph and analyze it using the HMM. We can then calculate the likelihood of a read occurring in a given cluster, resulting in an iterative clustering algorithm. For synthetic data, our proposed method reliably detected one variation site out of 9,000-bp synthetic long reads with a 15% sequencing-error rate and produced accurate clustering. It was also capable of clustering long reads from six very similar sequences containing only slight differences. For real data, we assembled putative multipartite structures of mitochondrial genomes of Arabidopsis thaliana from nine accessions sequenced using PacBio Sequel. The results indicated that there are recurrent and strain-specific structures in A. thaliana mitochondrial genomes.

The Mitochondrial Genome Assembly of Fennel (Foeniculum vulgare) Reveals Two Different atp6 Gene Sequences in Cytoplasmic Male Sterile Accessions

Cytoplasmic male sterility (CMS) has always aroused interest among researchers and breeders, being a valuable resource widely exploited not only to breed F1 hybrid varieties but also to investigate genes that control stamen and pollen development. With the aim of identifying candidate genes for CMS in fennel, we adopted an effective strategy relying on the comparison between mitochondrial genomes (mtDNA) of both fertile and sterile genotypes. mtDNA raw reads derived from a CMS genotype were assembled in a single molecule (296,483 bp), while a draft mtDNA assembly (166,124 nucleotides, 94 contigs) was performed using male fertile sample (MF) sequences. From their annotation and alignment, two atp6-like sequences were identified. atp6-, the putative mutant copy with a 300 bp truncation at the 5'-end, was found only in the mtDNA of CMS samples, while the wild type copy (atp6+) was detected only in the MF mtDNA. Further analyses (i.e., reads mapping and Sanger sequencing), revealed an atp6+ copy also in CMS samples, probably in the nuclear DNA. However, qPCRs performed on different tissues proved that, despite its availability, atp6+ is expressed only in MF samples, while apt6- mRNA was always detected in CMS individuals. In the light of these findings, the energy deficiency model could explain the pollen deficiency observed in male sterile flower. atp6- could represent a gene whose mRNA is translated into a not-fully functional protein leading to suboptimal ATP production that guarantees essential cellular processes but not a high energy demand process such as pollen development. Our study provides novel insights into the fennel mtDNA genome and its atp6 genes, and paves the way for further studies aimed at understanding their functional roles in the determination of male sterility.

Organelle genome composition and candidate gene identification for Nsa cytoplasmic male sterility in Brassica napus

Background

Nsa cytoplasmic male sterility (CMS) is a novel alloplasmic male sterility system derived from somatic hybridization between Brassica napus and Sinapis arvensis. Identification of the CMS-associated gene is a prerequisite for a better understanding of the origin and molecular mechanism of this CMS. With the development of genome sequencing technology, organelle genomes of Nsa CMS line and its maintainer line were sequenced by pyro-sequencing technology, and comparative analysis of the organelle genomes was carried out to characterize the organelle genome composition of Nsa CMS as well as to identify the candidate Nsa CMS-associated genes.

Results

Nsa CMS mitochondrial genome showed a higher collinearity with that of S. arvensis than B. napus, indicating that Nsa CMS mitochondrial genome was mainly derived from S. arvensis. However, mitochondrial genome recombination of parental lines was clearly detected. In contrast, the chloroplast genome of Nsa CMS was highly collinear with its B. napus parent, without any evidence of recombination of the two parental chloroplast genomes or integration from S. arvensis. There were 16 open reading frames (ORFs) specifically existed in Nsa CMS mitochondrial genome, which could not be identified in the maintainer line. Among them, three ORFs (orf224, orf309, orf346) possessing chimeric and transmembrane structure are most likely to be the candidate CMS genes. Sequences of all three candidate CMS genes in Nsa CMS line were found to be 100% identical with those from S. arvensis mitochondrial genome. Phylogenetic and homologous analysis showed that all the mitochondrial genes were highly conserved during evolution.

Conclusions

Nsa CMS contains a recombined mitochondrial genome of its two parental species with the majority form S. arvensis. Three candidate Nsa CMS genes were identified and proven to be derived from S. arvensis other than recombination of its two parental species. Further functional study of the candidate genes will help to identify the gene responsible for the CMS and the underlying molecular mechanism.

Organization Features of the Mitochondrial Genome of Sunflower (Helianthus annuus L.) with ANN2-Type Male-Sterile Cytoplasm

This study provides insights into the flexibility of the mitochondrial genome in sunflower (Helianthus annuus L.) as well as into the causes of ANN2-type cytoplasmic male sterility (CMS). De novo assembly of the mitochondrial genome of male-sterile HA89(ANN2) sunflower line was performed using high-throughput sequencing technologies. Analysis of CMS ANN2 mitochondrial DNA sequence revealed the following reorganization events: twelve rearrangements, seven insertions, and nine deletions. Comparisons of coding sequences from the male-sterile line with the male-fertile line identified a deletion of orf777 and seven new transcriptionally active open reading frames (ORFs): orf324, orf327, orf345, orf558, orf891, orf933, orf1197. Three of these ORFs represent chimeric genes involving atp6 (orf1197), cox2 (orf558), and nad6 (orf891). In addition, orf558, orf891, orf1197, as well as orf933, encode proteins containing membrane domain(s), making them the most likely candidate genes for CMS development in ANN2. Although the investigated CMS phenotype may be caused by simultaneous action of several candidate genes, we assume that orf1197 plays a major role in developing male sterility in ANN2. Comparative analysis of mitogenome organization in sunflower lines representing different CMS sources also allowed identification of reorganization hot spots in the mitochondrial genome of sunflower.