A genome‐scale analysis of mRNAs targeting to plant mitochondria: upstream AUGs in 5' untranslated regions reduce mitochondrial association

Intracellular sorting of mRNAs is an essential process for regulating gene expression and protein localization. Most mitochondrial proteins are nuclear‐encoded and imported into the mitochondria through post‐translational or co‐translational processes. In the latter case, mRNAs are found to be enriched in the vicinity of mitochondria. A genome‐scale analysis of mRNAs associated with mitochondria has been performed to determine plant cytosolic mRNAs targeted to the mitochondrial surface. Many messengers encoding mitochondrial proteins were found associated with mitochondria. These mRNAs correspond to particular functions and complexes, such as respiration or mitoribosomes, which indicates a coordinated control of mRNA localization within metabolic pathways. In addition, upstream AUGs in 5' untranslated regions (UTRs), which modulate the translation efficiency of downstream sequences, were found to negatively affect the association of mRNAs with mitochondria. A mutational approach coupled with in vivo mRNA visualization confirmed this observation. Moreover, this technique allowed the identification of 3'‐UTRs as another essential element for mRNA localization at the mitochondrial surface. Therefore, this work offers new insights into the mechanism, function and regulation of the association of cytosolic mRNAs with plant mitochondria.     

The first complete mitochondrial genome sequence of Nanorana parkeri and Nanorana ventripunctata (Amphibia: Anura: Dicroglossidae), with related phylogenetic analyses

Members of the Nanorana genus (family Dicroglossidae) are often referred to as excellent model species with which to study amphibian adaptations to extreme environments and also as excellent keystone taxa for providing insights into the evolution of the Dicroglossidae. However, a complete mitochondrial genome is currently only available for Nanorana pleskei. Thus, we analyzed the complete mitochondrial genomes of Nanorana parkeri and Nanorana ventripunctata to investigate their evolutionary relationships within Nanorana and their phylogenetic position in the family Dicroglossidae. Our results showed that the genomes of N. parkeri (17,837 bp) and N. ventripunctata (18,373 bp) encode 13 protein‐coding genes (PCGs), two ribosomal RNA genes, 23 transfer RNA (tRNA) genes, and a noncoding control region. Overall sequences and genome structure of the two species showed high degree of similarity with N. pleskei, although the motif structures and repeat sequences of the putative control region showed clear differences among these three Nanorana species. In addition, a tandem repeat of the tRNA‐Met gene was found located between the tRNA‐Gln and ND2 genes. On both the 5′ and 3′‐sides, the control region possessed distinct repeat regions; however, the CSB‐2 motif was not found in N. pleskei. Based on the nucleotide sequences of 13 PCGs, our phylogenetic analyses, using Bayesian inference and maximum‐likelihood methods, illustrate the taxonomic status of Nanorana with robust support showing that N. ventripunctata and N. pleskei are more closely related than they are to N. parkeri. In conclusion, our analyses provide a more robust and reliable perspective on the evolutionary history of Dicroglossidae than earlier analyses, which used only a single species (N. pleskei).     

Draft genome assembly and annotation of Glycyrrhiza uralensis,a medicinal legume

Chinese liquorice/licorice (Glycyrrhiza uralensis) is a leguminous plant species whose roots and rhizomes have been widely used as a herbal medicine and natural sweetener. Whole‐genome sequencing is essential for gene discovery studies and molecular breeding in liquorice. Here, we report a draft assembly of the approximately 379‐Mb whole‐genome sequence of strain 308‐19 of G. uralensis; this assembly contains 34 445 predicted protein‐coding genes. Comparative analyses suggested well‐conserved genomic components and collinearity of gene loci (synteny) between the genome of liquorice and those of other legumes such as Medicago and chickpea. We observed that three genes involved in isoflavonoid biosynthesis, namely, 2‐hydroxyisoflavanone synthase (CYP93C), 2,7,4′‐trihydroxyisoflavanone 4′‐O‐methyltransferase/isoflavone 4′‐O‐methyltransferase (HI4OMT) and isoflavone‐7‐O‐methyltransferase (7‐IOMT) formed a cluster on the scaffold of the liquorice genome and showed conserved microsynteny with Medicago and chickpea. Based on the liquorice genome annotation, we predicted genes in the P450 and UDP‐dependent glycosyltransferase (UGT) superfamilies, some of which are involved in triterpenoid saponin biosynthesis, and characterised their gene expression with the reference genome sequence. The genome sequencing and its annotations provide an essential resource for liquorice improvement through molecular breeding and the discovery of useful genes for engineering bioactive components through synthetic biology approaches.     

Gene family evolution in green plants with emphasis on the origination and evolution of Arabidopsis thaliana genes

Gene family size variation is an important mechanism that shapes the natural variation for adaptation in various species. Despite its importance, the pattern of gene family size variation in green plants is still not well understood. In particular, the evolutionary pattern of genes and gene families remains unknown in the model plant Arabidopsis thaliana in the context of green plants. In this study, eight representative genomes of green plants are sampled to study gene family evolution and characterize the origination of A. thaliana genes, respectively. Four important insights gained are that: (i) the rate of gene gains and losses is about 0.001359 per gene every million years, similar to the rate in yeast, Drosophila, and mammals; (ii) some gene families evolved rapidly with extreme expansions or contractions, and 2745 gene families present in all the eight species represent the ‘core’ proteome of green plants; (iii) 70% of A. thaliana genes could be traced back to 450 million years ago; and (iv) intriguingly, A. thaliana genes with early origination are under stronger purifying selection and more conserved. In summary, the present study provides genome‐wide insights into evolutionary history and mechanisms of genes and gene families in green plants and especially in A. thaliana.     

Gene targeting in polymerase theta-deficient Arabidopsis thaliana

Agrobacterium tumefaciens-mediated transformation has been for decades the preferred tool to generate transgenic plants. During this process, a T-DNA carrying transgenes is transferred from the bacterium to plant cells, where it randomly integrates into the genome via polymerase theta (Polθ)-mediated end joining (TMEJ). Targeting of the T-DNA to a specific genomic locus via homologous recombination (HR) is also possible, but such gene targeting (GT) events occur at low frequency and are almost invariably accompanied by random integration events. An additional complexity is that the product of recombination between T-DNA and target locus may not only map to the target locus (true GT), but also to random positions in the genome (ectopic GT). In this study, we have investigated how TMEJ functionality affects the biology of GT in plants, by using Arabidopsis thaliana mutated for the TEBICHI gene, which encodes for Polθ. Whereas in TMEJ-proficient plants we predominantly found GT events accompanied by random T-DNA integrations, GT events obtained in the teb mutant background lacked additional T-DNA copies, corroborating the essential role of Polθ in T-DNA integration. Polθ deficiency also prevented ectopic GT events, suggesting that the sequence of events leading up to this outcome requires TMEJ. Our findings provide insights that can be used for the development of strategies to obtain high-quality GT events in crop plants.     

Whole‐exome targeted sequencing of the uncharacterized pine genome

The large genome size of many species hinders the development and application of genomic tools to study them. For instance, loblolly pine (Pinus taeda L.), an ecologically and economically important conifer, has a large and yet uncharacterized genome of 21.7 Gbp. To characterize the pine genome, we performed exome capture and sequencing of 14 729 genes derived from an assembly of expressed sequence tags. Efficiency of sequence capture was evaluated and shown to be similar across samples with increasing levels of complexity, including haploid cDNA, haploid genomic DNA and diploid genomic DNA. However, this efficiency was severely reduced for probes that overlapped multiple exons, presumably because intron sequences hindered probe:exon hybridizations. Such regions could not be entirely avoided during probe design, because of the lack of a reference sequence. To improve the throughput and reduce the cost of sequence capture, a method to multiplex the analysis of up to eight samples was developed. Sequence data showed that multiplexed capture was reproducible among 24 haploid samples, and can be applied for high‐throughput analysis of targeted genes in large populations. Captured sequences were de novo assembled, resulting in 11 396 expanded and annotated gene models, significantly improving the knowledge about the pine gene space. Interspecific capture was also evaluated with over 98% of all probes designed from P. taeda that were efficient in sequence capture, were also suitable for analysis of the related species Pinus elliottii Engelm.     

Presence of nuclear copies of mitochondrial origin (NUMTs) in two related species of stingless bee genus Melipona (Hymenoptera: Meliponini)

Transferred copies of mitochondrial DNA (mtDNA) into the nuclear genome (numts) have been reported in several Hymenoptera species, even at a high density in the honey bee nuclear genome. The accidental amplification of numts in phylogenetic studies focused on mtDNA highlights the importance of a correct determination of numts and their related mtDNA sequences. We report here the presence of numts derived from a mitochondrial rDNA 16S gene in the genome of the stingless bee species Melipona colimana and M. fasciata (tribe Meliponini) from Western Mexico. PCR products were cloned in both species obtaining thirty paralogous numts. Numts were identified by the presence of insertions and deletions and the disruption of the 16S secondary structure. Further phylogenetic analyses including alternative mitochondrial cox1 and nuclear ITS1 genes have revealed the presence of another numt (cox1) in the nuclear genome of these two species, and place both as sister lineages within the subgenus Michmelia. This is one of the first studies reporting the presence of numts in Meliponini species, and supports previous studies suggesting frequent transfer of mtDNA to the nuclear genome in Hymenoptera.     

Nucleotide‐binding resistance gene signatures in sugar beet,insights from a new reference genome

Nucleotide‐binding (NB‐ARC), leucine‐rich‐repeat genes (NLRs) account for 60.8% of resistance (R) genes molecularly characterized from plants. NLRs exist as large gene families prone to tandem duplication and transposition, with high sequence diversity among crops and their wild relatives. This diversity can be a source of new disease resistance, but difficulty in distinguishing specific sequences from homologous gene family members hinders characterization of resistance for improving crop varieties. Current genome sequencing and assembly technologies, especially those using long‐read sequencing, are improving resolution of repeat‐rich genomic regions and clarifying locations of duplicated genes, such as NLRs. Using the conserved NB‐ARC domain as a model, 231 tentative NB‐ARC loci were identified in a highly contiguous genome assembly of sugar beet, revealing diverged and truncated NB‐ARC signatures as well as full‐length sequences. The NB‐ARC‐associated proteins contained NLR resistance gene domains, including TIR, CC and LRR, as well as other integrated domains. Phylogenetic relationships of partial and complete domains were determined, and patterns of physical clustering in the genome were evaluated. Comparison of sugar beet NB‐ARC domains to validated R‐genes from monocots and eudicots suggested extensive Beta vulgaris‐specific subfamily expansions. The NLR landscape in the rhizomania resistance conferring Rz region of Chromosome 3 was characterized, identifying 26 NLR‐like sequences spanning 20 MB. This work presents the first detailed view of NLR family composition in a member of the Caryophyllales, builds a foundation for additional disease resistance work in B. vulgaris, and demonstrates an additional nucleic‐acid‐based method for NLR prediction in non‐model plant species.     

Oligonucleotide‐directed mutagenesis for precision gene editing

Differences in gene sequences, many of which are single nucleotide polymorphisms, underlie some of the most important traits in plants. With humanity facing significant challenges to increase global agricultural productivity, there is an urgent need to accelerate the development of these traits in plants. oligonucleotide‐directed mutagenesis (ODM), one of the many tools of Cibus’ Rapid Trait Development System ( RTDS ^?) technology, offers a rapid, precise and non‐transgenic breeding alternative for trait improvement in agriculture to address this urgent need. This review explores the application of ODM as a precision genome editing technology, with emphasis on using oligonucleotides to make targeted edits in plasmid, episomal and chromosomal DNA of bacterial, fungal, mammalian and plant systems. The process of employing ODM by way of RTDS technology has been improved in many ways by utilizing a fluorescence conversion system wherein a blue fluorescent protein (BFP) can be changed to a green fluorescent protein (GFP) by editing a single nucleotide of the BFP gene (CAC→TAC; H66 to Y66). For example, dependent on oligonucleotide length, applying oligonucleotide‐mediated technology to target the BFP transgene in Arabidopsis thaliana protoplasts resulted in up to 0.05% precisely edited GFP loci. Here, the development of traits in commercially relevant plant varieties to improve crop performance by genome editing technologies such as ODM, and by extension RTDS , is reviewed.     

GABI‐DUPLO: a collection of double mutants to overcome genetic redundancy in Arabidopsis thaliana

Owing to duplication events in its progenitor, more than 90% of the genes in the Arabidopsis thaliana genome are members of multigene families. A set of 2108 gene families, each consisting of precisely two unlinked paralogous genes, was identified in the nuclear genome of A. thaliana on the basis of sequence similarity. A systematic method for the creation of double knock‐out lines for such gene pairs, designated as DUPLO lines, was established and 200 lines are now publicly available. Their initial phenotypic characterisation led to the identification of seven lines with defects that emerge only in the adult stage. A further six lines display seedling lethality and 23 lines were lethal before germination. Another 14 lines are known to show phenotypes under non‐standard conditions or at the molecular level. Knock‐out of gene pairs with very similar coding sequences or expression profiles is more likely to produce a mutant phenotype than inactivation of gene pairs with dissimilar profiles or sequences. High coding sequence similarity and highly similar expression profiles are only weakly correlated, implying that promoter and coding regions of these gene pairs display different degrees of diversification.     

Chromosome‐level reference genome of X12, a highly virulent race of the soybean cyst nematode Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is a major pest of soybean that is spreading across major soybean production regions worldwide. Increased SCN virulence has recently been observed in both the United States and China. However, no study has reported a genome assembly for H. glycines at the chromosome scale. Herein, the first chromosome‐level reference genome of X12, an unusual SCN race with high infection ability, is presented. Using whole‐genome shotgun (WGS) sequencing, Pacific Biosciences (PacBio) sequencing, Illumina paired‐end sequencing, 10X Genomics linked reads and high‐throughput chromatin conformation capture (Hi‐C) genome scaffolding techniques, a 141.01‐megabase (Mb) assembled genome was obtained with scaffold and contig N50 sizes of 16.27 Mb and 330.54 kilobases (kb), respectively. The assembly showed high integrity and quality, with over 90% of Illumina reads mapped to the genome. The assembly quality was evaluated using Core Eukaryotic Genes Mapping Approach and Benchmarking Universal Single‐Copy Orthologs. A total of 11,882 genes were predicted using de novo, homolog and RNAseq data generated from eggs, second‐stage juveniles (J2), third‐stage juveniles (J3) and fourth‐stage juveniles (J4) of X12, and 79.0% of homologous sequences were annotated in the genome. These high‐quality X12 genome data will provide valuable resources for research in a broad range of areas, including fundamental nematode biology, SCN–plant interactions and co‐evolution, and also contribute to the development of technology for overall SCN management.