Mitochondrial genome sequences from wild and cultivated barley (Hordeum vulgare)

Advancing the Conservation and Utilization of Barley Genetic Resources: Insights into Germplasm Management and Breeding for Sustainable Agriculture

Barley is a very important crop particularly in marginal dry areas, where it often serves as the most viable option for farmers. Additionally, barley carries great significance in the Western world, serving not only as a fundamental crop for animal feed and malting but also as a nutritious food source. The broad adaptability of barley and its ability to withstand various biotic and abiotic stresses often make this species the sole cereal that can be cultivated in arid regions. The collection and utilization of barley genetic resources are crucial for identifying valuable traits to enhance productivity and mitigate the adverse effects of climate change. This review aims to provide an overview of the management and exploitation of barley genetic resources. Furthermore, the review explores the relationship between gene banks and participatory breeding, offering insights into the diversity and utilization of barley genetic resources through some examples such as the initiatives undertaken by ICARDA. Finally, this contribution highlights the importance of these resources for boosting barley productivity, addressing climate change impacts, and meeting the growing food demands in a rapidly changing agriculture. The understanding and utilizing the rich genetic diversity of barley can contribute to sustainable agriculture and ensure the success of this vital crop for future generations globally.

Complete mitochondrial genome of Agrostis stolonifera: insights into structure, Codon usage, repeats, and RNA editing

BACKGROUND

Plants possess mitochondrial genomes that are large and complex compared to animals. Despite their size, plant mitochondrial genomes do not contain significantly more genes than their animal counterparts. Studies into the sequence and structure of plant mitochondrial genomes heavily imply that the main mechanism driving replication of plant mtDNA, and offer valuable insights into plant evolution, energy production, and environmental adaptation.

RESULTS

This study presents the first comprehensive analysis of Agrostis stolonifera's mitochondrial genome, characterized by a branched structure comprising three contiguous chromosomes, totaling 560,800 bp with a GC content of 44.07%. Annotations reveal 33 unique protein-coding genes (PCGs), 19 tRNA genes, and 3 rRNA genes. The predominant codons for alanine and glutamine are GCU and CAA, respectively, while cysteine and phenylalanine exhibit weaker codon usage biases. The mitogenome contains 73, 34, and 23 simple sequence repeats (SSRs) on chromosomes 1, 2, and 3, respectively. Chromosome 1 exhibits the most frequent A-repeat monomeric SSR, whereas chromosome 2 displays the most common U-repeat monomeric SSR. DNA transformation analysis identifies 48 homologous fragments between the mitogenome and chloroplast genome, representing 3.41% of the mitogenome's total length. The PREP suite detects 460 C-U RNA editing events across 33 mitochondrial PCGs, with the highest count in the ccmFn gene and the lowest in the rps7 gene. Phylogenetic analysis confirms A. stolonifera's placement within the Pooideae subfamily, showing a close relationship to Lolium perenne, consistent with the APG IV classification system. Numerous homologous co-linear blocks are observed in A. stolonifera's mitogenomes and those of related species, while certain regions lack homology.

CONCLUSIONS

The unique features and complexities of the A. stolonifera mitochondrial genome, along with its similarities and differences to related species, provide valuable insights into plant evolution, energy production, and environmental adaptation. The findings from this study significantly contribute to the growing body of knowledge on plant mitochondrial genomes and their role in plant biology.

Complete chloroplast genomes of three wild perennial Hordeum species from Central Asia: genome structure, mutation hotspot, phylogenetic relationships, and comparative analysis

Hordeum L. is widely distributed in mountain or plateau of subtropical and warm temperate regions around the world. Three wild perennial Hordeum species, including H. bogdanii, H. brevisubulatum, and H. violaceum, have been used as forage and for grassland ecological restoration in high-altitude areas in recent years. To date, the degree of interspecies sequence variation in the three Hordeum species within existing gene pools is still not well-defined. Herein, we sequenced and assembled chloroplast (cp) genomes of the three species. The results revealed that the cp genome of H. bogdanii showed certain sequence variations compared with the cp genomes of the other two species (H. brevisubulatum and H. violaceum), and the latter two were characterized by a higher relative affinity. Parity rule 2 plot (PR2) analysis illuminated that most genes of all ten Hordeum species were concentrated in nucleotide T and G. Numerous single nucleotide polymorphism (SNP) and insertion/deletion (In/Del) events were detected in the three Hordeum species. A series of hotspots regions (tRNA-GGU ~ tRNA-GCA, tRNA-UGU ~ ndhJ, psbE ~ rps18, ndhF ~ tRNA-UAG, etc.) were identified by mVISTA procedures, and the five highly polymorphic genes (tRNA-UGC, tRNA-UAA, tRNA-UUU, tRNA-UAC, and ndhA) were proved by the nucleotide diversity (Pi). Although the distribution and existence of cp simple sequence repeats (cpSSRs) were predicted in the three Hordeum cp genomes, no rearrangement was found between them. A similar phenomenon has been found in the cp genome of the other seven Hordeum species, which has been published so far. In addition, evolutionary relationships were reappraised based on the currently reported cp genome of Hordeum L. This study offers a framework for gaining a better understanding of the evolutionary history of Hordeum species through the re-examination of their cp genomes, and by identifying highly polymorphic genes and hotspot regions that could provide important insights into the genetic diversity and differentiation of these species.

Deciphering the enigma of RNA editing in the ATP1_alpha subunit of ATP synthase in Triticum aestivum

There is evidence that RNA editing is related to plant cellular stress as well as electron transport organelles, such as mitochondria. The mitochondrial atp1 gene encodes the alpha-subunit of Atp synthase. Control as well as two periods of drought stress treatments were analyzed in the cDNAs generated from the mitochondrial atp1 gene of two cultivars of Triticum aestivum [Giza 168 (G168) and Gemmiza 10 (GM10)]. Following RNA-seq data assembly, atp1 cDNAs from the control (acc. no. OQ129415), 2-hour (acc. no. OQ129416), and 12-hour (acc. no. OQ129417) time points of the T. aestivum cultivar G168 were obtained. Control (acc. no. OQ129419), 2-hour (acc. no. OQ129420), and 12-hour (acc. no. OQ129421) samples all included reconstructed atp1 transcripts from Gemmiza 10. Atp1 transcripts were assembled using the wheat atp1 gene (acc. no. NC_036024). RNA-seq raw data was utilized to identify 11 RNA editing sites in atp1 in the tolerant cultivar Giza168 and 6 in the sensitive cultivar Gemmiza10. The significant difference in RNA editing observed between control and drought stress conditions in sites led to synonymous amino acids. This led to no change in tertiary structure between tolerant and sensitive cultivars. But the change was focused between produced protein and its correspondence sequence on DNA.

The mitochondrial genome of the diploid oat Avena longiglumis

BACKGROUND

Avena longiglumis Durieu (2n = 2x = 14) is a wild relative of cultivated oat (Avena sativa, 2n = 6x = 42) with good agronomic and nutritional traits. The plant mitochondrial genome has a complex organization and carries genetic traits of value in exploiting genetic resources, not least male sterility alleles used to generate F₁ hybrid seeds. Therefore, we aim to complement the chromosomal-level nuclear and chloroplast genome assemblies of A. longiglumis with the complete assembly of the mitochondrial genome (mitogenome) based on Illumina and ONT long reads, comparing its structure with Poaceae species.

RESULTS

The complete mitochondrial genome of A. longiglumis can be represented by one master circular genome being 548,445 bp long with a GC content of 44.05%. It can be represented by linear or circular DNA molecules (isoforms or contigs), with multiple alternative configurations mediated by long (4,100-31,235 bp) and medium (144-792 bp) size repeats. Thirty-five unique protein-coding genes, three unique rRNA genes, and 11 unique tRNA genes are identified. The mitogenome is rich in duplications (up to 233 kb long) and multiple tandem or simple sequence repeats, together accounting for more than 42.5% of the total length. We identify homologous sequences between the mitochondrial, plastid and nuclear genomes, including the exchange of eight plastid-derived tRNA genes, and nuclear-derived retroelement fragments. At least 85% of the mitogenome is duplicated in the A. longiglumis nuclear genome. We identify 269 RNA editing sites in mitochondrial protein-coding genes including stop codons truncating ccmFC transcripts.

CONCLUSIONS

Comparative analysis with Poaceae species reveals the dynamic and ongoing evolutionary changes in mitochondrial genome structure and gene content. The complete mitochondrial genome of A. longiglumis completes the last link of the oat reference genome and lays the foundation for oat breeding and exploiting the biodiversity in the genus.

Sequencing and analysis of the complete mitochondrial genomes of Toona sinensis and Toona ciliata reveal evolutionary features of Toona

BACKGROUND

Toona is a critical genus in the Meliaceae, and the plants of this group are an asset for both restorative and restorative purposes, the most flexible of which are Toona sinensis and Toona ciliata. To concentrate on the advancement of mitochondrial(Mt) genome variety in T.sinensis and T.ciliata, the Mt genomes of the two species were sequenced in high throughput independently, after de novo assembly and annotation to construct a Mt genome map for comparison in genome structure. Find their repetitive sequences and analyze them in comparison with the chloroplast genome, along with Maximum-likelihood(ML) phylogenetic analysis with 16 other relatives.

RESULTS

(1) T. sinensis and T.ciliata are both circular structures with lengths of 683482 bp and 68300 bp, respectively. They share a high degree of similarity in encoding genes and have AT preferences. All of them have the largest Phe concentration and are the most frequently used codons. (2) Both of their Mt genome are highly preserved in terms of structural and functional genes, while the main variability is reflected in the length of tRNA, the number of genes, and the value of RSCU. (3) T. siniensis and T. ciliata were detected to have 94 and 87 SSRs, respectively, of which mononucleotides accounted for the absolute proportion. Besides, the vast majority of their SSRs were found to be poly-A or poly-T. (4)10 and 11 migrating fragments were identified in the comparison with the chloroplast genome, respectively. (5) In the ML evolutionary tree, T.sinensis and T.ciliata clustered individually into a small branch with 100% support, reflecting two species of Toona are very similarly related to each other.

CONCLUSIONS

This research provides a basis for the exploitation of T.sinensis and T.ciliata in terms of medicinal, edible, and timber resources to avoid confusion; at the same time, it can explore the evolutionary relationship between the Toona and related species, which does not only have an important practical value, but also provides a theoretical basis for future hybrid breeding of forest trees, molecular markers, and evolutionary aspects of plants, which has great scientific significance.

Effect of salinity on ccmfn gene RNA editing of mitochondria in wild barley and uncommon types of RNA editing

The primary function of mitochondria is cellular respiration and energy production. Cytochrome C complex is an essential complex that transports electrons in the respiratory chain between complex III and complex IV. One of this complex's main subunits is CcmFN, which is believed to be crucial for holocytochrome assembly. In wild-type plant Hordeum vulgare subsp. spontaneum, four ccmfn cDNAs are subjected to high salt stress (500 mM salinity), 0 h (or control) (GenBank accession no. ON764850), after 2 h (GenBank accession no. ON7648515), after 12 h (GenBank accession no. ON764852), and after 24 h (GenBank accession no. ON764853) and mtDNA of ccmfn gene (GenBank accession no. ON764854). Using raw data from RNA-seq, 47 sites with nucleotide and amino acid modifications were detected. There were ten different RNA editing types, with most of them are C to U. Unusual editing types in plants have also been found, such as A to C, C to A, A to G, A to U, T to A, T to C, C to G, G to C, and T to G. High levels of editing were observed in control as well as treatments of salinity stress. Amino acid changes were found in 43 sites; nearly all showed hydrophilic to hydrophilic alterations. Only C749 showed regulation under salinity stress.

Characterization and Phylogenetic Analyses of the Complete Mitochondrial Genome of Sugarcane (Saccharum spp. Hybrids) Line A1

Modern sugarcane cultivars are highly polyploid with complex nuclear genomic genetic background, while their mitochondrion (mt) genomes are much simpler, smaller and more manageable and could provide useful phylogenetic information. In this study, the mt genome of a modern commercial cultivar A1 was sequenced via Illumina Hiseq XTen and PacBio Sequel platform. The assembled and annotated mitochondrial genomes of A1 were composed of two circular DNA molecules, one large and one small, which were named Chromosome 1 and Chromosome 2. The two distinct circular chromosomes of mitogenome construct is consisted with other sugarcane cultivars i.e., Saccharum officinarum Khon Kaen 3 and Saccharum spp. hybrids ROC22 and FN15. The Chromosome 1 of A1 mitogenome is 300,822 bp in length with the GC content of 43.94%, and 7.14% of Chromosome 1 sequences (21,468 nucleotides) are protein coding genes (PCGs) while 92.86% (279,354 nucleotides) are intergenic region. The length of Chromosome 2 is 144,744 bp with the GC content of 43.57%, and 8.20% of Chromosome 2 sequences (11,865 nucleotides) are PCGs while 91.80% (132,879 nucleotides) are intergenic region. A total of 43 genes are located on Chromosome 1, which contains 22 PCGs (six nad genes, four rps genes, four atp genes, three ccm genes, three cox genes, one mat gene and one mtt gene) and 21 non-coding genes including 15 tRNAs and 6 rRNAs. Chromosome 2 includes 18 genes in total, which contains 13 PCGs (four nad genes, three rps genes, two atp genes, one ccm gene, one cob gene, one cox gene and one rpl gene) and five non-coding genes (tRNA genes). Analysis of codon usage of 35 PCGs showed that codon ending in A/U was preferred. Investigation of gene composition indicated that the types and copy numbers of CDS genes, tRNAs and rRNAs of A1 and FN15 were identical. The cox1 gene has two copies and the trnP gene has one copy in A1, FN15 and ROC22 three lines, while there is only one copy of cox1 and two copies of trnP in S. officinarum Khon Kaen 3. In addition, S. officinarum Khon Kaen 3 have no nad1 gene and rps7 gene. 100 sequence repeats, 38 SSRs and 444 RNA editing sites in A1 mt genome were detected. Moreover, the maximum likelihood phylogenetic analysis found that A1 were more closely related to S. spp. hybrid (ROC22 and FN15) and S. officinarum (Khon Kaen 3). Herein, the complete mt genome of A1 will provide essential DNA molecular information for further phylogenetic and evolutionary analysis for Saccharum and Poaceae.

Chromosome-scale assembly of barley cv. 'Haruna Nijo' as a resource for barley genetics

Cultivated barley (Hordeum vulgare ssp. vulgare) is used for food, animal feed, and alcoholic beverages and is widely grown in temperate regions. Both barley and its wild progenitor (H. vulgare ssp. spontaneum) have large 5.1-Gb genomes. High-quality chromosome-scale assemblies for several representative barley genotypes, both wild and domesticated, have been constructed recently to populate the nascent barley pan-genome infrastructure. Here, we release a chromosome-scale assembly of the Japanese elite malting barley cultivar 'Haruna Nijo' using a similar methodology as in the barley pan-genome project. The 4.28-Gb assembly had a scaffold N50 size of 18.9 Mb. The assembly showed high collinearity with the barley reference genome 'Morex' cultivar, with some inversions. The pseudomolecule assembly was characterized using transcript evidence of gene projection derived from the reference genome and de novo gene annotation achieved using published full-length cDNA sequences and RNA-Seq data for 'Haruna Nijo'. We found good concordance between our whole-genome assembly and the publicly available BAC clone sequence of 'Haruna Nijo'. Interesting phenotypes have since been identified in Haruna Nijo; its genome sequence assembly will facilitate the identification of the underlying genes.

The first report of RNA U to C or G editing in the mitochondrial NADH dehydrogenase subunit 5 (Nad5) transcript of wild barley

BACKGROUND

Nad dehydrogenase complex in mtDNA has a significant role in cellular respiration. One of the largest subunits in the complex is subunit 5 (Nad5).

METHODS AND RESULTS

Four cDNAs of the Hordeum vulgare subsp. spontaneum nad5 gene have been characterized and subjected to four phases of 0.5 M salinity, at 0 h (control, accession no. MT235236), after 2 h (acc. no. MT235237), after 12 h (acc. no. MT235238) and after 24 h (acc. no. MT235239). Utilizing raw data from RNA-seq, ten RNA editing sites were reported. Seven sites have common editing from C to U in positions (C1490, C1859, C1895, C1900, C1901, C1916, C1918). A rare editing event U to C was detected in two positions (U1650 and U1652) and a novel editing event U to G was for the first time in positions nad5-U231. The highest editing level was shown in 2 and 12 h after salinity exposure. After 24 h, these edits were disrupted, possibly due to the launch of the programed cell death mechanism. However, the RNA editing in positions U1650, U1652 and U231 was fixed at all exposure times.

CONCLUSIONS

Although study clarified the role of salinity stress in nad5 RNA editing sites, the main achievements are first report of U to G RNA editing in plants at position U231 and first report of U to C editing in the nad5 gene at U1650 and U1652.

Deciphering the organelle genomes and transcriptomes of a common ornamental plant Ligustrum quihoui reveals multiple fragments of transposable elements in the mitogenome

Ligustrum quihoui (L. quihoui) is an important hedge material for landscaping and also possesses medicinal value. To generate genomic resources for better understanding the evolutionary history of this important plant, the organelle genomes of L. quihoui are de novo assembled and functionally annotated. Compared with other Oleaceae species, the 163,069 bp chloroplast genome of L. quihoui exhibits a typical quadripartite structure with highly conserved gene content and gene order, while the 848,451 bp mitochondrial genome of L. quihoui exhibits highly divergent genome size and gene content. Codon usage analyses show that genes related with photosynthesis and mitochondrial respiratory chain show inconsistent codon biases. A total of 48,760 bp transposable elements (TEs) fragments and 41,887 bp chloroplast-like sequences are found in the L. quihoui mitochondrial genome. A striking discrepancy of RNA editing between the two organelle genomes is found in L. quihoui, in which 146 mitochondrial editing sites coexist with only 43 such sites in chloroplast. Based on DNA and RNA-Seq data, we propose that GTG may act as the start codon of mitochondrial rpl16 in Oleaceae species. Phylogenetic analysis based on chloroplast genome shows that L. quihoui and L. japonicum form a sister clade within the genus Ligustrum.

History and future perspectives of barley genomics

Barley (Hordeum vulgare), one of the most widely cultivated cereal crops, possesses a large genome of 5.1 Gbp. Through various international collaborations, the genome has recently been sequenced and assembled at the chromosome-scale by exploiting available genetic and genomic resources. Many wild and cultivated barley accessions have been collected and preserved around the world. These accessions are crucial to obtain diverse natural and induced barley variants. The barley bioresource project aims to investigate the diversity of this crop based on purified seed and DNA samples of a large number of collected accessions. The long-term goal of this project is to analyse the genome sequences of major barley accessions worldwide. In view of technical limitations, a strategy has been employed to establish the exome structure of a selected number of accessions and to perform high-quality chromosome-scale assembly of the genomes of several major representative accessions. For the future project, an efficient annotation pipeline is essential for establishing the function of genomes and genes as well as for using this information for sequence-based digital barley breeding. In this article, the author reviews the existing barley resources along with their applications and discuss possible future directions of research in barley genomics.

Salinity effects on nad3 gene RNA editing of wild barley mitochondria

Nad complex plays a very important role during cellular respiration. nad3 (nad dehydrogenase subunit 3) is one of the biggest subunits in this complex. Four cDNAs of nad3 gene were characterized in Hordeum vulgare subsp. spontaneum at exposed to four periods of 500 mM salinity, 0 h or control (accession no. MN066165), after 2 h (accession no. MN066166), after 12 h (accession no. MN066167) and after 24 h (accession no. MN066168) using RNA-seq raw data. Seventeen RNA editing sites were found in positions (or nucleotide nos. C5, C39, C44, C61, C62, C79, C80, C147, C185, C190, C191, C208, C209, C275, C317, C344, C349) within the nad3 coding region. These alterations represent differential editing at four exposure times. The maximum editing rate was revealed 2 and 12 h after salinity exposure. However, these edits were disrupted after 24 h probably due to the initiation of program cell death machinery. We found that RNA editing not only improved protein function but also may improve codon bias by altering the nucleotide without any change in amino acid. Characterization of pentatricopeptide repeat-containing protein At4g13650 (PPRSp1) in wild barley helped us to understand the behavior of editing sites C190 and C191 under salinity. Position - 6 in cis-element upstream editing sites of C155, C190 and C191 may be vital to the editing process in these sites by PPRSp1 protein. The differential editing of this gene under salinity led to a relationship between RNA editing and cellular respiration regulation.

The study of organelle DNA variability in alloplasmic barley lines in the NGS era

Alloplasmic lines are a suitable model for studying molecular coevolution and interrelations between genetic systems of plant cells. Whole chloroplast (cp) and mitochondrial (mt) genome sequences were obtained by the MiSeq System (Illumina). Organelle DNA samples were prepared from a set of 12 alloplasmic barley lines with different cytoplasms of Hordeum vulgare ssp. spontaneum and H. vulgare ssp. vulgare, as well as from their paternal varieties. A bioinformatic approach for analysis of NGS data obtained on an organellar DNA mix has been developed and verified. A comparative study of Hordeum organelle genomes’ variability and disposition of polymorphic loci was conducted. Eight types of chloroplast DNA and 5 types of mitochondrial DNA were distinguished for the barley sample set examined. These results were compared with the previous data of a restriction fragment length polymorphism (RFLP) study of organelle DNAs for the same material. Formerly established data about a field evaluation of alloplasmic barley lines were revised in the light of information about organelle genomes gained after NGS. Totally 17 polymorphic loci were found at exons of chloroplast genomes. Seven of the SNPs were located in the genes of the Ndh complex. The nonsynonymous changes of nucleotides were detected in the matK, rpoC1, ndhK, ndhG and infA genes. Some of the SNPs detected are very similar in codon position and in the type of amino acid substitution to the places where RNA editing can occur. Thus, these results outline new perspectives for the future study of nuclear-cytoplasmic interactions in alloplasmic lines.

The draft genome of a wild barley genotype reveals its enrichment in genes related to biotic and abiotic stresses compared to cultivated barley

Wild barley (Hordeum spontaneum) is the progenitor of cultivated barley (Hordeum vulgare) and provides a rich source of genetic variations for barley improvement. Currently, the genome sequences of wild barley and its differences with cultivated barley remain unclear. In this study, we report a high-quality draft assembly of wild barley accession (AWCS276; henceforth named as WB1), which consists of 4.28 Gb genome and 36 395 high-confidence protein-coding genes. BUSCO analysis revealed that the assembly included full lengths of 95.3% of the 956 single-copy plant genes, illustrating that the gene-containing regions have been well assembled. By comparing with the genome of the cultivated genotype Morex, it is inferred that the WB1 genome contains more genes involved in resistance and tolerance to biotic and abiotic stresses. The presence of the numerous WB1-specific genes indicates that, in addition to enhance allele diversity for genes already existing in the cultigen, exploiting the wild barley taxon in breeding should also allow the incorporation of novel genes. Furthermore, high levels of genetic variation in the pericentromeric regions were detected in chromosomes 3H and 5H between the wild and cultivated genotypes, which may be the results of domestication. This H. spontaneum draft genome assembly will help to accelerate wild barley research and be an invaluable resource for barley improvement and comparative genomics research.

Methods and Tools for Plant Organelle Genome Sequencing, Assembly, and Downstream Analysis

Organelles play an important role in a eukaryotic cell. Among them, the two organelles, chloroplast and mitochondria, are responsible for the critical function of photosynthesis and aerobic respiration. Organellar genomes are also very important for plant systematic studies. Here we have described the methods for isolation of the mitochondrial and plastid DNA and its subsequent sequencing with the help of NGS technology. We have also discussed in detail the various tools available for assembly, annotation, and visualization of the organelle genome sequence.

The Role of Transposable Elements in Pongamia Unigenes and Protein Diversity

Pongamia pinnata (also called Millettia pinnata), a non-edible oil yielding tree, is well known for its multipurpose benefits and acts as a potential source for medicine and biodiesel preparation. Due to increase in demand for cultivation, understanding of genetic diversity is an important parameter for further breeding and cultivation programme. Transposable elements (TEs) are a major component of plant genome but still, their evolutionary significance in Pongamia remains unexplored. In view to understand the role of TEs in genome diversity, Pongamia unigenes were screened for the presence of TE cassettes. Our analysis showed the presence of all categories of TE cassettes in unigenes with major contribution of long terminal repeat-retrotransposons towards unigene diversity. Interestingly, the insertion of some TEs was also observed in both organellar genomes. The study of insertion of TEs in coding sequence is of great interest as they may be responsible for protein diversity thereby influencing the phenotype. The present investigation confirms the exaptation phenomenon in pyruvate decarboxylase (PDC) gene where the entire exon sequence was derived from Ty3-gypsy like retrotransposon. The study of PDC protein revealed the translation of gypsy element into protein. Furthermore, the phylogenetic study confirmed the diversity in PDC gene due to insertion of the gypsy element, where the PDC genes with and without gypsy insertion were clustered separately.

Thermal plasticity of the circadian clock is under nuclear and cytoplasmic control in wild barley

Temperature compensation, expressed as the ability to maintain clock characteristics (mainly period) in face of temperature changes, that is, robustness, is considered a key feature of circadian clock systems. In this study, we explore the genetic basis for lack of robustness, that is, plasticity, of circadian clock as reflected by photosynthesis rhythmicity. The clock rhythmicity of a new wild barley reciprocal doubled haploid population was analysed with a high temporal resolution of pulsed amplitude modulation of chlorophyll fluorescence under optimal (22°C) and high (32°C) temperature. This comparison between two environments pointed to the prevalence of clock acceleration under heat. Genotyping by sequencing of doubled haploid lines indicated a rich recombination landscape with minor fixation (less than 8%) for one of the parental alleles. Quantitative genetic analysis included genotype by environment interactions and binary-threshold models. Variation in the circadian rhythm plasticity phenotypes, expressed as change (delta) of period and amplitude under two temperatures, was associated with maternal organelle genome (the plasmotype), as well as with several nuclear loci. This first reported rhythmicity driven by nuclear loci and plasmotype with few identified variants, paves the way for studying impact of cytonuclear variations on clock robustness and on plant adaptation to changing environments.

Pacbio Sequencing Reveals Identical Organelle Genomes between American Cranberry (Vaccinium macrocarpon Ait.) and a Wild Relative

Breeding efforts in the American cranberry (Vaccinium macrocarpon Ait.), a North American perennial fruit crop of great importance, have been hampered by the limited genetic and phenotypic variability observed among cultivars and experimental materials. Most of the cultivars commercially used by cranberry growers today were derived from a few wild accessions bred in the 1950s. In different crops, wild germplasm has been used as an important genetic resource to incorporate novel traits and increase the phenotypic diversity of breeding materials. Vaccinium microcarpum (Turcz. ex Rupr.) Schmalh. and V. oxycoccos L., two closely related species, may be cross-compatible with the American cranberry, and could be useful to improve fruit quality such as phytochemical content. Furthermore, given their northern distribution, they could also help develop cold hardy cultivars. Although these species have previously been analyzed in diversity studies, genomic characterization and comparative studies are still lacking. In this study, we sequenced and assembled the organelle genomes of the cultivated American cranberry and its wild relative, V. microcarpum. PacBio sequencing technology allowed us to assemble both mitochondrial and plastid genomes at very high coverage and in a single circular scaffold. A comparative analysis revealed that the mitochondrial genome sequences were identical between both species and that the plastids presented only two synonymous single nucleotide polymorphisms (SNPs). Moreover, the Illumina resequencing of additional accessions of V. microcarpum and V. oxycoccos revealed high genetic variation in both species. Based on these results, we provided a hypothesis involving the extension and dynamics of the last glaciation period in North America, and how this could have shaped the distribution and dispersal of V. microcarpum. Finally, we provided important data regarding the polyploid origin of V. oxycoccos.

High intraspecific diversity of Restorer-of-fertility-like genes in barley

Nuclear restorer of fertility (Rf) genes suppress the effects of mitochondrial genes causing cytoplasmic male sterility (CMS), a condition in which plants fail to produce viable pollen. Rf genes, many of which encode RNA-binding pentatricopeptide repeat (PPR) proteins, are applied in hybrid breeding to overcome CMS used to block self-pollination of the seed parent. Here, we characterise the repertoire of restorer-of-fertility-like (RFL) PPR genes in barley (Hordeum vulgare). We found 26 RFL genes in the reference genome ('Morex') and an additional 51 putative orthogroups (POGs) in a re-sequencing data set from 262 barley genotypes and landraces. Whereas the sequences of some POGs are highly conserved across hundreds of barley accessions, the sequences of others are much more variable. High sequence variation strongly correlates with genomic location - the most variable genes are found in a cluster on chromosome 1H. A much higher likelihood of diversifying selection was found for genes within this cluster than for genes present as singlets. This work includes a comprehensive analysis of the patterns of intraspecific variation of RFL genes. The RFL sequences characterised in this study will be useful for the development of new markers for fertility restoration loci.

Next generation sequencing reveals packaging of host RNAs by brome mosaic virus

Although RNA viruses evolved the mechanisms of specific encapsidation, miss-packaging of cellular RNAs has been reported in such RNA virus systems as flock house virus or cucumber necrosis virus. To find out if brome mosaic virus (BMV), a tripartite RNA virus, can package cellular RNAs, BMV was propagated in barley and in Nicotiana benthamiana hosts, purified by cesium chloride (CsCl) gradient ultracentrifugation followed by nuclease treatment to remove any contaminating cellular (host) RNAs. The extracted virion RNA was then sequenced by using next-generation sequencing (NGS RNA-Seq) with the Illumina protocol. Bioinformatic analysis revealed the content of host RNAs ranging from 0.07% for BMV extracted from barley to 0.10% for the virus extracted from N. benthamiana. The viruses from two sources appeared to co-encapsidate different patterns of host-RNAs, including ribosomal RNAs (rRNAs), messenger RNAs (mRNAs) but also mitochondrial and plastid RNAs and, interestingly, transposable elements, both transposons and retrotransposons. Our data reveal that BMV virions can carry host RNAs, having a potential to mediate horizontal gene transfer (HGT) in plants.

Complete mitochondrial genome sequences of Brassica rapa (Chinese cabbage and mizuna), and intraspecific differentiation of cytoplasm in B. rapa and Brassica juncea

The complete sequence of the mitochondrial genome was determined for two cultivars of Brassica rapa. After determining the sequence of a Chinese cabbage variety, 'Oushou hakusai', the sequence of a mizuna variety, 'Chusei shiroguki sensuji kyomizuna', was mapped against the sequence of Chinese cabbage. The precise sequences where the two varieties demonstrated variation were ascertained by direct sequencing. It was found that the mitochondrial genomes of the two varieties are identical over 219,775 bp, with a single nucleotide polymorphism (SNP) between the genomes. Because B. rapa is the maternal species of an amphidiploid crop species, Brassica juncea, the distribution of the SNP was observed both in B. rapa and B. juncea. While the mizuna type SNP was restricted mainly to cultivars of mizuna (japonica group) in B. rapa, the mizuna type was widely distributed in B. juncea. The finding that the two Brassica species have these SNP types in common suggests that the nucleotide substitution occurred in wild B. rapa before both mitotypes were domesticated. It was further inferred that the interspecific hybridization between B. rapa and B. nigra took place twice and resulted in the two mitotypes of cultivated B. juncea.