MIPS PlantsDB: a database framework for comparative plant genome research

Adaptive evolution of the enigmatic Takakia now facing climate change in Tibet

The most extreme environments are the most vulnerable to transformation under a rapidly changing climate. These ecosystems harbor some of the most specialized species, which will likely suffer the highest extinction rates. We document the steepest temperature increase (2010-2021) on record at altitudes of above 4,000 m, triggering a decline of the relictual and highly adapted moss Takakia lepidozioides. Its de-novo-sequenced genome with 27,467 protein-coding genes includes distinct adaptations to abiotic stresses and comprises the largest number of fast-evolving genes under positive selection. The uplift of the study site in the last 65 million years has resulted in life-threatening UV-B radiation and drastically reduced temperatures, and we detected several of the molecular adaptations of Takakia to these environmental changes. Surprisingly, specific morphological features likely occurred earlier than 165 mya in much warmer environments. Following nearly 400 million years of evolution and resilience, this species is now facing extinction.

Comparative analysis of repeat content in plant genomes, large and small

The DNA Features pipeline is the analysis pipeline at EMBL-EBI that annotates repeat elements, including transposable elements. With Ensembl's goal to stay at the cutting edge of genome annotation, we proved that this pipeline needed an update. We then created a new analysis that allowed the Ensembl database to store the repeat classification from the PGSB repeat classification (Recat). This new dataset was then fetched using Perl scripts and used to prove that the pipeline modification induced a gain in sensitivity. Finally, we performed a comparative analysis of transposable element distribution in all plant species available, raising new questions about transposable elements in certain branches of the taxonomic tree.

A graph-based genome and pan-genome variation of the model plant Setaria

Setaria italica (foxtail millet), a founder crop of East Asian agriculture, is a model plant for C4 photosynthesis and developing approaches to adaptive breeding across multiple climates. Here we established the Setaria pan-genome by assembling 110 representative genomes from a worldwide collection. The pan-genome is composed of 73,528 gene families, of which 23.8%, 42.9%, 29.4% and 3.9% are core, soft core, dispensable and private genes, respectively; 202,884 nonredundant structural variants were also detected. The characterization of pan-genomic variants suggests their importance during foxtail millet domestication and improvement, as exemplified by the identification of the yield gene SiGW3, where a 366-bp presence/absence promoter variant accompanies gene expression variation. We developed a graph-based genome and performed large-scale genetic studies for 68 traits across 13 environments, identifying potential genes for millet improvement at different geographic sites. These can be used in marker-assisted breeding, genomic selection and genome editing to accelerate crop improvement under different climatic conditions.

Phylogenomics of Aralia sect. Aralia (Araliaceae): Signals of hybridization and insights into its species delimitations and intercontinental biogeography

Genome-scale data have significantly increased the number of informative characters for phylogenetic analyses and recent studies have also revealed widespread phylogenomic discordance in many plant lineages. Aralia sect. Aralia is a small plant lineage (14 spp.) of the ginseng family Araliaceae with a disjunct distribution between eastern Asia (11 spp.) and North America (3 spp.). We herein employ sequences of hundreds of nuclear loci and the complete plastomes using targeted sequence capture and genome skimming to reconstruct the phylogenetic and biogeographic history of this section. We detected substantial conflicts among nuclear genes, yet different analytical strategies generated largely congruent topologies from the nuclear data. Significant cytonuclear discordance was detected, especially concerning the positions of the three North American species. The phylogenomic results support two intercontinental disjunctions: (1) Aralia californica of western North America is sister to the eastern Asian clade consisting of A. cordata and A. continentalis in the nuclear tree, and (2) the eastern North American A. racemosa forms a clade with A. bicrenata from southwestern North America, and the North American A. racemosa - A. bicrenata clade is then sister to the eastern Asian clade consisting of A. glabra (Japan), A. fargesii (C China), and A. apioides and A. atropurpurea (the Hengduan Mountains). Aralia cordata is supported to be disjunctly distributed in Japan, Taiwan, the Ulleung island of Korea, and in Central, Southwest and South China, and Aralia continentalis is redefined with a narrower distribution in Northeast China, eastern Russia and peninsular Korea.

An improved assembly of the "Cascade" hop (Humulus lupulus) genome uncovers signatures of molecular evolution and refines time of divergence estimates for the Cannabaceae family

We present a chromosome-level assembly of the Cascade hop (Humulus lupulus L. var. lupulus) genome. The hop genome is large (2.8 Gb) and complex, and early attempts at assembly were fragmented. Recent advances have made assembly of the hop genome more tractable, transforming the extent of investigation that can occur. The chromosome-level assembly of Cascade was developed by scaffolding the previously reported Cascade assembly generated with PacBio long-read sequencing and polishing with Illumina short-read DNA sequencing. We developed gene models and repeat annotations and used a controlled bi-parental mapping population to identify significant sex-associated markers. We assessed molecular evolution in gene sequences, gene family expansion and contraction, and time of divergence from Cannabis sativa and other closely related plant species using Bayesian inference. We identified the putative sex chromosome in the female genome based on significant sex-associated markers from the bi-parental mapping population. While the estimate of repeat content (~64%) is similar to the estimate for the hemp genome, syntenic blocks in hop contain a greater percentage of LTRs. Hop is enriched for disease resistance-associated genes in syntenic gene blocks and expanded gene families. The Cascade chromosome-level assembly will inform cultivation strategies and serve to deepen our understanding of the hop genomic landscape, benefiting hop researchers and the Cannabaceae genomics community.

Annotation of Siberian Larch (Larix sibirica Ledeb.) Nuclear Genome—One of the Most Cold-Resistant Tree Species in the Only Deciduous GENUS in Pinaceae

The recent release of the nuclear, chloroplast and mitochondrial genome assemblies of Siberian larch (Larix sibirica Ledeb.), one of the most cold-resistant tree species in the only deciduous genus of Pinaceae, with seasonal senescence and a rot-resistant valuable timber widely used in construction, greatly contributed to the development of genomic resources for the larch genus. Here, we present an extensive repeatome analysis and the first annotation of the draft nuclear Siberian larch genome assembly. About 66% of the larch genome consists of highly repetitive elements (REs), with the likely wave of retrotransposons insertions into the larch genome estimated to occur 4–5 MYA. In total, 39,370 gene models were predicted, with 87% of them having homology to the Arabidopsis-annotated proteins and 78% having at least one GO term assignment. The current state of the genome annotations allows for the exploration of the gymnosperm and angiosperm species for relative gene abundance in different functional categories. Comparative analysis of functional gene categories across different angiosperm and gymnosperm species finds that the Siberian larch genome has an overabundance of genes associated with programmed cell death (PCD), autophagy, stress hormone biosynthesis and regulatory pathways; genes that may play important roles in seasonal senescence and stress response to extreme cold in larch. Despite being incomplete, the draft assemblies and annotations of the conifer genomes are at a point of development where they now represent a valuable source for further genomic, genetic and population studies.

The Larix kaempferi genome reveals new insights into wood properties

Here, through single-molecule real-time sequencing, we present a high-quality genome sequence of the Japanese larch (Larix kaempferi), a conifer species with great value for wood production and ecological afforestation. The assembled genome is 10.97 Gb in size, harboring 45,828 protein-coding genes. Of the genome, 66.8% consists of repeat sequences, of which long terminal repeat retrotransposons are dominant and make up 69.86%. We find that tandem duplications have been responsible for the expansion of genes involved in transcriptional regulation and stress responses, unveiling their crucial roles in adaptive evolution. Population transcriptome analysis reveals that lignin content in L. kaempferi is mainly determined by the process of monolignol polymerization. The expression values of six genes (LkCOMT7, LkCOMT8, LkLAC23, LkLAC102, LkPRX148, and LkPRX166) have significantly positive correlations with lignin content. These results indicated that the increased expression of these six genes might be responsible for the high lignin content of the larches' wood. Overall, this study provides new genome resources for investigating the evolution and biological function of conifer trees, and also offers new insights into wood properties of larches.

TransposonUltimate: software for transposon classification, annotation and detection

Most genomes harbor a large number of transposons, and they play an important role in evolution and gene regulation. They are also of interest to clinicians as they are involved in several diseases, including cancer and neurodegeneration. Although several methods for transposon identification are available, they are often highly specialised towards specific tasks or classes of transposons, and they lack common standards such as a unified taxonomy scheme and output file format. We present TransposonUltimate, a powerful bundle of three modules for transposon classification, annotation, and detection of transposition events. TransposonUltimate comes as a Conda package under the GPL-3.0 licence, is well documented and it is easy to install through https://github.com/DerKevinRiehl/TransposonUltimate. We benchmark the classification module on the large TransposonDB covering 891,051 sequences to demonstrate that it outperforms the currently best existing solutions. The annotation and detection modules combine sixteen existing softwares, and we illustrate its use by annotating Caenorhabditis elegans, Rhizophagus irregularis and Oryza sativa subs. japonica genomes. Finally, we use the detection module to discover 29 554 transposition events in the genomes of 20 wild type strains of C. elegans. Databases, assemblies, annotations and further findings can be downloaded from (https://doi.org/10.5281/zenodo.5518085).

Whole-genome sequence and methylome profiling of the almond [Prunus dulcis (Mill.) D.A. Webb] cultivar 'Nonpareil'

Almond [Prunus dulcis (Mill.) D.A. Webb] is an economically important, specialty nut crop grown almost exclusively in the United States. Breeding and improvement efforts worldwide have led to the development of key, productive cultivars, including 'Nonpareil,' which is the most widely grown almond cultivar. Thus far, genomic resources for this species have been limited, and a whole-genome assembly for 'Nonpareil' is not currently available despite its economic importance and use in almond breeding worldwide. We generated a 571X coverage genome sequence using Illumina, PacBio, and optical mapping technologies. Gene prediction revealed 49,321 putative genes using MinION Oxford nanopore and Illumina RNA sequencing, and genome annotation found that 68% of predicted models are associated with at least one biological function. Furthermore, epigenetic signatures of almond, namely DNA cytosine methylation, have been implicated in a variety of phenotypes including self-compatibility, bud dormancy, and development of noninfectious bud failure. In addition to the genome sequence and annotation, this report also provides the complete methylome of several almond tissues, including leaf, flower, endocarp, mesocarp, exocarp, and seed coat. Comparisons between methylation profiles in these tissues revealed differences in genome-wide weighted % methylation and chromosome-level methylation enrichment.

Exploration of glutathione reductase for abiotic stress response in bread wheat (Triticum aestivum L.)

A total of seven glutathione reductase (GR) genes were identified in Triticum aestivum, which were used for comparative structural characterization, phylogenetic analysis and expression profiling with the GR genes of other cereal plants. The modulated gene expression and enzyme activity revealed the role of GRs in abiotic stress response in T. aestivum. Glutathione reductase (GR) is an enzymatic antioxidant that converts oxidized glutathione (GSSG) into reduced glutathione (GSH) through the ascorbate-glutathione cycle. In this study, a total of seven GR genes forming two homeologous groups were identified in the allohexaploid genome of bread wheat (Triticum aestivum). Besides, we identified three GR genes in each Aegilops tauschii, Brachypodium distachyon, Triticum urartu and Sorghum bicolor, which were used for comparative characterization. Phylogenetic analysis revealed the clustering of GR proteins into two groups; class I and class II, which were predicted to be localized in cytoplasm and chloroplast, respectively. The exon-intron and conserved motif patterns were almost conserved in each group, in which a maximum of 10 and 17 exons were present in chloroplastic and cytoplasmic GRs, respectively. The protein structure analysis confirmed the occurrence of conserved pyridine nucleotide disulfide oxidoreductase (Pyr_redox) and pyridine nucleotide disulfide oxidoreductase dimerization (Pyr_redox_dim) domains in each GR. The active site of GR proteins consisted of two conserved cysteine residues separated by four amino acid residues. Promoter analysis revealed the occurrence of growth and stress-related cis-active elements. Tissue-specific expression profiling suggested the involvement of GRs in both vegetative and reproductive tissue development in various plants. The differential expression of TaGR genes and enhanced GR enzyme activity suggested their roles under drought, heat, salt and arsenic stress. Interaction of GRs with other proteins and chemical compounds of the ascorbate-glutathione cycle revealed their coordinated functioning. The current study will provide a foundation for the validation of the precise role of each GR gene in future studies.

Finding and Characterizing Repeats in Plant Genomes

Plant genomes contain a particularly high proportion of repeated structures of various types. This chapter proposes a guided tour of the available software that can help biologists to scan automatically for these repeats in sequence data or check hypothetical models intended to characterize their structures. Since transposable elements (TEs) are a major source of repeats in plants, many methods have been used or developed for this broad class of sequences. They are representative of the range of tools available for other classes of repeats and we have provided two sections on this topic (for the analysis of genomes or directly of sequenced reads), as well as a selection of the main existing software. It may be hard to keep up with the profusion of proposals in this dynamic field and the rest of the chapter is devoted to the foundations of an efficient search for repeats and more complex patterns. We first introduce the key concepts of the art of indexing and mapping or querying sequences. We end the chapter with the more prospective issue of building models of repeat families. We present the Machine Learning approach first, seeking to build predictors automatically for some families of ET, from a set of sequences known to belong to this family. A second approach, the linguistic (or syntactic) approach, allows biologists to describe themselves and check the validity of models of their favorite repeat family.

Ensembl Genomes 2022: an expanding genome resource for non-vertebrates

Ensembl Genomes (https://www.ensemblgenomes.org) provides access to non-vertebrate genomes and analysis complementing vertebrate resources developed by the Ensembl project (https://www.ensembl.org). The two resources collectively present genome annotation through a consistent set of interfaces spanning the tree of life presenting genome sequence, annotation, variation, transcriptomic data and comparative analysis. Here, we present our largest increase in plant, metazoan and fungal genomes since the project's inception creating one of the world's most comprehensive genomic resources and describe our efforts to reduce genome redundancy in our Bacteria portal. We detail our new efforts in gene annotation, our emerging support for pangenome analysis, our efforts to accelerate data dissemination through the Ensembl Rapid Release resource and our new AlphaFold visualization. Finally, we present details of our future plans including updates on our integration with Ensembl, and how we plan to improve our support for the microbial research community. Software and data are made available without restriction via our website, online tools platform and programmatic interfaces (available under an Apache 2.0 license). Data updates are synchronised with Ensembl's release cycle.

CicerSpTEdb: A web-based database for high-resolution genome-wide identification of transposable elements in Cicer species

Recently, Cicer species have experienced increased research interest due to their economic importance, especially in genetics, genomics, and crop improvement. The Cicer arietinum, Cicer reticulatum, and Cicer echinospermum genomes have been sequenced and provide valuable resources for trait improvement. Since the publication of the chickpea draft genome, progress has been made in genome assembly, functional annotation, and identification of polymorphic markers. However, work is still needed to identify transposable elements (TEs) and make them available for researchers. In this paper, we present CicerSpTEdb, a comprehensive TE database for Cicer species that aims to improve our understanding of the organization and structural variations of the chickpea genome. Using structure and homology-based methods, 3942 C. echinospermum, 3579 C. reticulatum, and 2240 C. arietinum TEs were identified. Comparisons between Cicer species indicate that C. echinospermum has the highest number of LTR-RT and hAT TEs. C. reticulatum has more Mutator, PIF Harbinger, Tc1 Mariner, and CACTA TEs, while C. arietinum has the highest number of Helitron. CicerSpTEdb enables users to search and visualize TEs by location and download their results. The database will provide a powerful resource that can assist in developing TE target markers for molecular breeding and answer related biological questions. Database URL: http://cicersptedb.easyomics.org/index.php.

Evolutionary history and pan-genome dynamics of strawberry (Fragaria spp.)

Strawberry (Fragaria spp.) has emerged as a model system for various fundamental and applied research in recent years. In total, the genomes of five different species have been sequenced over the past 10 y. Here, we report chromosome-scale reference genomes for five strawberry species, including three newly sequenced species' genomes, and genome resequencing data for 128 additional accessions to estimate the genetic diversity, structure, and demographic history of key Fragaria species. Our analyses obtained fully resolved and strongly supported phylogenies and divergence times for most diploid strawberry species. These analyses also uncovered a new diploid species (Fragaria emeiensis Jia J. Lei). Finally, we constructed a pan-genome for Fragaria and examined the evolutionary dynamics of gene families. Notably, we identified multiple independent single base mutations of the MYB10 gene associated with white pigmented fruit shared by different strawberry species. These reference genomes and datasets, combined with our phylogenetic estimates, should serve as a powerful comparative genomic platform and resource for future studies in strawberry.

K-mer counting and curated libraries drive efficient annotation of repeats in plant genomes

The annotation of repetitive sequences within plant genomes can help in the interpretation of observed phenotypes. Moreover, repeat masking is required for tasks such as whole-genome alignment, promoter analysis, or pangenome exploration. Although homology-based annotation methods are computationally expensive, k-mer strategies for masking are orders of magnitude faster. Here, we benchmarked a two-step approach, where repeats were first called by k-mer counting and then annotated by comparison to curated libraries. This hybrid protocol was tested on 20 plant genomes from Ensembl, with the k-mer-based Repeat Detector (Red) and two repeat libraries (REdat, last updated in 2013, and nrTEplants, curated for this work). Custom libraries produced by RepeatModeler were also tested. We obtained repeated genome fractions that matched those reported in the literature but with shorter repeated elements than those produced directly by sequence homology. Inspection of the masked regions that overlapped genes revealed no preference for specific protein domains. Most Red-masked sequences could be successfully classified by sequence similarity, with the complete protocol taking less than 2 h on a desktop Linux box. A guide to curating your own repeat libraries and the scripts for masking and annotating plant genomes can be obtained at https://github.com/Ensembl/plant-scripts.

Two gap-free reference genomes and a global view of the centromere architecture in rice

Rice (Oryza sativa), a major staple throughout the world and a model system for plant genomics and breeding, was the first crop genome sequenced almost two decades ago. However, reference genomes for all higher organisms to date contain gaps and missing sequences. Here, we report the assembly and analysis of gap-free reference genome sequences for two elite O. sativa xian/indica rice varieties, Zhenshan 97 and Minghui 63, which are being used as a model system for studying heterosis and yield. Gap-free reference genomes provide the opportunity for a global view of the structure and function of centromeres. We show that all rice centromeric regions share conserved centromere-specific satellite motifs with different copy numbers and structures. In addition, the similarity of CentO repeats in the same chromosome is higher than across chromosomes, supporting a model of local expansion and homogenization. Both genomes have over 395 non-TE genes located in centromere regions, of which ∼41% are actively transcribed. Two large structural variants at the end of chromosome 11 affect the copy number of resistance genes between the two genomes. The availability of the two gap-free genomes lays a solid foundation for further understanding genome structure and function in plants and breeding climate-resilient varieties.

Patterns of stability and change in the maize genome: a case study of small RNA transcriptomes in two recombinant inbred lines and their progenitors

Small RNAs (sRNAs) are epigenetic regulators of eukaryotic genes and transposable elements (TEs). Diverse sRNA expression patterns exist within a species, but how this diversity arises is not well understood. To provide a window into the dynamics of maize sRNA patterning, sRNA and mRNA transcriptomes were examined in two related Zea mays recombinant inbred lines (RILs) and their inbred parents. Analysis of these RILs revealed that most clusters of sRNA expression retained the parental sRNA expression level. However, expression states that differ from the parental allele were also observed, predominantly reflecting decreases in sRNA expression. When RIL sRNA expression differed from the parental allele, the new state was frequently similar between the two RILs, and similar to the expression state found at the allele in the other parent. Novel sRNA expression patterns, distinct from those of either parent, were rare. Additionally, examination of sRNA expression over TEs revealed one TE family, Gyma, which showed consistent enrichment for RIL sRNA expression differences compared to those found in parental alleles. These findings provide insights into how sRNA silencing might evolve over generations and suggest that further investigation into the molecular nature of sRNA trans regulators is warranted.

Concerted genomic and epigenomic changes accompany stabilization of Arabidopsis allopolyploids

During evolution successful allopolyploids must overcome 'genome shock' between hybridizing species but the underlying process remains elusive. Here, we report concerted genomic and epigenomic changes in resynthesized and natural Arabidopsis suecica (TTAA) allotetraploids derived from Arabidopsis thaliana (TT) and Arabidopsis arenosa (AA). A. suecica shows conserved gene synteny and content with more gene family gain and loss in the A and T subgenomes than respective progenitors, although A. arenosa-derived subgenome has more structural variation and transposon distributions than A. thaliana-derived subgenome. These balanced genomic variations are accompanied by pervasive convergent and concerted changes in DNA methylation and gene expression among allotetraploids. The A subgenome is hypomethylated rapidly from F1 to resynthesized allotetraploids and convergently to the T-subgenome level in natural A. suecica, despite many other methylated loci being inherited from F1 to all allotetraploids. These changes in DNA methylation, including small RNAs, in allotetraploids may affect gene expression and phenotypic variation, including flowering, silencing of self-incompatibility and upregulation of meiosis- and mitosis-related genes. In conclusion, concerted genomic and epigenomic changes may improve stability and adaptation during polyploid evolution.

Plant genomics in Africa: present and prospects

Plants are the world's most consumed goods. They are of high economic value and bring many health benefits. In most countries in Africa, the supply and quality of food will rise to meet the growing population's increasing demand. Genomics and other biotechnology tools offer the opportunity to improve subsistence crops and medicinal herbs in the continent. Significant advances have been made in plant genomics, which have enhanced our knowledge of the molecular processes underlying both plant quality and yield. The sequencing of complex genomes of African plant species, facilitated by the continuously evolving next-generation sequencing technologies and advanced bioinformatics approaches, has provided new opportunities for crop improvement. This review summarizes the achievements of genome sequencing projects of endemic African plants in the last two decades. We also present perspectives and challenges for future plant genomic studies that will accelerate important plant breeding programs for African communities. These challenges include a lack of basic facilities, a lack of sequencing and bioinformatics facilities, and a lack of skills to design genomics studies. However, it is imperative to state that African countries have become key players in the plant genome revolution and genome derived-biotechnology. Therefore, African governments should invest in public plant genomics research and applications, establish bioinformatics platforms and training programs, and stimulate university and industry partnerships to fully deploy plant genomics, particularly in the fields of agriculture and medicine.

Introgressed DNA within a Zea mays near isogenic line displays lower levels of 24nt sRNA expression than the homologous region from the recurrent parent

Near isogenic lines (NILs) are a classical genetic tool used to dissect the actions of an allele when placed in a uniform genetic background. Although the goal of NIL creation is to examine the effects of a single allele in isolation, DNA linked to the allele is invariably retained and can confound any allele specific effects. In addition to genetic variation, highly polymorphic species like Zea mays will contain introgressed polymorphisms encompassing transposable elements (TEs) and the cis acting small RNA (sRNA) that represses them. Through transcriptomics, we described the sRNA and TE transcriptional expression differences between a W22-derived introgression and its homologous B73 region. As anticipated, many sRNA expression differences were found. Unexpectedly, however, 24nt sRNA expression over the introgressed region was low overall compared to both the homologous B73 region and the rest of the genome. Across the introgression, low sRNA expression was accompanied by increased TE transcription. Possible explanations for the observed trends in sRNA and TE expression across the introgression are discussed. These findings support the notion that any introgressed allele is in an epigenetic environment distinct from that found at the allele from the recurrent parent. Additionally, these results suggest that further study of sRNA expression levels during the introgression process is warranted.

A Practical Guide on Computational Tools and Databases for Transposable Elements in Plants

In the age of big data, obtaining precise information about the research topic of interesting is extremely important. Keeping this in mind, this chapter focuses on providing a practical knowledge guide about computational tools and databases of transposable elements (TE) in plants. For that, we organize and present this text in three sections: (1) a discussion about tools and databases on this theme; (2) hands-on of how to use a few of them; (3) an exploratory data analysis on public TE data. Finally, we are going deep to present the main challenges and possible solutions to improve resources and tools.

Cis-regulatory evolution spotlights species differences in the adaptive potential of gene expression plasticity

Phenotypic plasticity is the variation in phenotype that a single genotype can produce in different environments and, as such, is an important component of individual fitness. However, whether the effect of new mutations, and hence evolution, depends on the direction of plasticity remains controversial. Here, we identify the cis-acting modifications that have reshaped gene expression in response to dehydration stress in three Arabidopsis species. Our study shows that the direction of effects of most cis-regulatory variants differentiating the response between A. thaliana and the sister species A. lyrata and A. halleri depends on the direction of pre-existing plasticity in gene expression. A comparison of the rate of cis-acting variant accumulation in each lineage indicates that the selective forces driving adaptive evolution in gene expression favors regulatory changes that magnify the stress response in A. lyrata. The evolutionary constraints measured on the amino-acid sequence of these genes support this interpretation. In contrast, regulatory changes that mitigate the plastic response to stress evolved more frequently in A. halleri. Our results demonstrate that pre-existing plasticity may be a stepping stone for adaptation, but its selective remodeling differs between lineages.

Sequence of the supernumerary B chromosome of maize provides insight into its drive mechanism and evolution

B chromosomes are enigmatic elements in thousands of plant and animal genomes that persist in populations despite being nonessential. They circumvent the laws of Mendelian inheritance but the molecular mechanisms underlying this behavior remain unknown. Here we present the sequence, annotation, and analysis of the maize B chromosome providing insight into its drive mechanism. The sequence assembly reveals detailed locations of the elements involved with the cis and trans functions of its drive mechanism, consisting of nondisjunction at the second pollen mitosis and preferential fertilization of the egg by the B-containing sperm. We identified 758 protein-coding genes in 125.9 Mb of B chromosome sequence, of which at least 88 are expressed. Our results demonstrate that transposable elements in the B chromosome are shared with the standard A chromosome set but multiple lines of evidence fail to detect a syntenic genic region in the A chromosomes, suggesting a distant origin. The current gene content is a result of continuous transfer from the A chromosomal complement over an extended evolutionary time with subsequent degradation but with selection for maintenance of this nonvital chromosome.

A chromosome-scale genome assembly of European hazel (Corylus avellana L.) reveals targets for crop improvement

The European hazelnut (Corylus avellana L.) is a tree crop of economic importance worldwide, but especially for northern Turkey, where the majority of production takes place. Hazelnut production is currently challenged by environmental stresses, such as a recent outbreak of severe powdery mildew disease; furthermore, allergy to hazelnuts is an increasing health concern in some regions. In order to provide a foundation for using the available hazelnut genetic resources for crop improvement, we produced a fully assembled genome sequence and annotation for a hazelnut species, from C. avellana cv. 'Tombul', one of the most important Turkish varieties. A hybrid sequencing strategy, combining short reads, long reads and proximity ligation methods, enabled us to resolve heterozygous regions and produce a high-quality 370-Mb assembly that agrees closely with cytogenetic studies and genetic maps of the 11 C. avellana chromosomes, and covers 97.8% of the estimated genome size. The genome includes 27 270 high-confidence protein-coding genes, over 20 000 of which were functionally annotated based on homology with known plant proteins. We focused particularly on gene families encoding hazelnut allergens, and the Mildew resistance Locus O (MLO) proteins that are an important susceptibility factor for powdery mildew. The complete assembly enabled us to differentiate between members of these families and to identify homologues that may be important in mildew disease and hazelnut allergy. These findings provide examples of how the genome can be used to guide research and to develop effective strategies for crop improvement in C. avellana.

A draft phased assembly of the diploid Cascade hop (Humulus lupulus) genome

Hop (Humulus lupulus L. var Lupulus) is a diploid, dioecious plant with a history of cultivation spanning more than one thousand years. Hop cones are valued for their use in brewing and contain compounds of therapeutic interest including xanthohumol. Efforts to determine how biochemical pathways responsible for desirable traits are regulated have been challenged by the large (2.8 Gb), repetitive, and heterozygous genome of hop. We present a draft haplotype-phased assembly of the Cascade cultivar genome. Our draft assembly and annotation of the Cascade genome is the most extensive representation of the hop genome to date. PacBio long-read sequences from hop were assembled with FALCON and partially phased with FALCON-Unzip. Comparative analysis of haplotype sequences provides insight into selective pressures that have driven evolution in hop. We discovered genes with greater sequence divergence enriched for stress-response, growth, and flowering functions in the draft phased assembly. With improved resolution of long terminal retrotransposons (LTRs) due to long-read sequencing, we found that hop is over 70% repetitive. We identified a homolog of cannabidiolic acid synthase (CBDAS) that is expressed in multiple tissues. The approaches we developed to analyze the draft phased assembly serve to deepen our understanding of the genomic landscape of hop and may have broader applicability to the study of other large, complex genomes.

The genetic basis of cytoplasmic male sterility and fertility restoration in wheat

Hybrid wheat varieties give higher yields than conventional lines but are difficult to produce due to a lack of effective control of male fertility in breeding lines. One promising system involves the Rf1 and Rf3 genes that restore fertility of wheat plants carrying Triticum timopheevii-type cytoplasmic male sterility (T-CMS). Here, by genetic mapping and comparative sequence analyses, we identify Rf1 and Rf3 candidates that can restore normal pollen production in transgenic wheat plants carrying T-CMS. We show that Rf1 and Rf3 bind to the mitochondrial orf279 transcript and induce cleavage, preventing expression of the CMS trait. The identification of restorer genes in wheat is an important step towards the development of hybrid wheat varieties based on a CMS-Rf system. The characterisation of their mode of action brings insights into the molecular basis of CMS and fertility restoration in plants.

Epigenetic Variability Among Saffron Crocus (Crocus sativus L.) Accessions Characterized by Different Phenotypes

This work represents the first epigenomic study carried out on saffron crocus. Five accessions of saffron, showing differences in tepal pigmentation, yield of saffron and flowering time, were analyzed at the epigenetic level by applying a methylation-sensitive restriction enzyme-sequencing (MRE-seq) approach. Five accession-specific hypomethylomes plus a reference hypomethylome, generated by combining the sequence data from the single accessions, were obtained. Assembled sequences were annotated against existing online databases. In the absence of the Crocus genome, the rice genome was mainly used as the reference as it is the best annotated genome among monocot plants. Comparison of the hypomethylomes revealed many differentially methylated regions, confirming the high epigenetic variability present among saffron accessions, including sequences encoding for proteins that could be good candidates to explain the accessions' alternative phenotypes. In particular, transcription factors involved in flowering process (MADS-box and TFL) and for the production of pigments (MYB) were detected. Finally, by comparing the generated sequences of the different accessions, a high number of SNPs, likely having arisen as a consequence of the prolonged vegetative propagation, were detected, demonstrating surprisingly high genetic variability. Gene ontology (GO) was performed to map and visualize sequence polymorphisms located within the GOs and to compare their distributions among different accessions. As well as suggesting the possible existence of alternative phenotypes with a genetic basis, a clear difference in polymorphic GO is present among accessions based on their geographic origin, supporting a possible signature of selection in the Indian accession with respect to the Spanish ones.

Molecular Characterization Revealed the Role of Thaumatin-Like Proteins of Bread Wheat in Stress Response

Thaumatin-like proteins (TLPs) are related to pathogenesis-related-5 (PR-5) family and involved in stress response. Herein, a total of 93 TLP genes were identified in the genome of Triticum aestivum. Further, we identified 26, 27, 39, and 37 TLP genes in the Brachypodium distachyon, Oryza sativa, Sorghum bicolor, and Zea mays genomes for comparative characterization, respectively. They could be grouped into small and long TLPs with conserved thaumatin signature motif. Tightly clustered genes exhibited conserved gene and protein structure. The physicochemical analyses suggested significant differences between small and long TLPs. Evolutionary analyses suggested the role of duplication events and purifying selection in the expansion of the TLP gene family. Expression analyses revealed the possible roles of TLPs in plant development and abiotic and fungal stress response. Recombinant expression of TaTLP2-B in Saccharomyces cerevisiae provided significant tolerance against cold, heat, osmotic, and salt stresses. The results depicted the importance of TLPs in cereal crops that would be highly useful in future crop improvement programs.

TE-greedy-nester: structure-based detection of LTR retrotransposons and their nesting

Motivation

Transposable elements (TEs) in eukaryotes often get inserted into one another, forming sequences that become a complex mixture of full-length elements and their fragments. The reconstruction of full-length elements and the order in which they have been inserted is important for genome and transposon evolution studies. However, the accumulation of mutations and genome rearrangements over evolutionary time makes this process error-prone and decreases the efficiency of software aiming to recover all nested full-length TEs.

Results

We created software that uses a greedy recursive algorithm to mine increasingly fragmented copies of full-length LTR retrotransposons in assembled genomes and other sequence data. The software called TE-greedy-nester considers not only sequence similarity but also the structure of elements. This new tool was tested on a set of natural and synthetic sequences and its accuracy was compared to similar software. We found TE-greedy-nester to be superior in a number of parameters, namely computation time and full-length TE recovery in highly nested regions.

Availability and implementation

http://gitlab.fi.muni.cz/lexa/nested.

Supplementary information

Supplementary data are available at Bioinformatics online.

NGOME-Lite: Proteome-wide prediction of spontaneous protein deamidation highlights differences between taxa

Asparagines in proteins deamidate spontaneously, which changes the chemical structure of a protein and often affects its function. Current prediction algorithms for asparagine deamidation require a structure as an input or are too slow to be applied at a proteomic scale. We present NGOME-Lite, a new version of our sequence-based predictor for spontaneous asparagine deamidation that is faster by over two orders of magnitude at a similar degree of accuracy. The algorithm takes into account intrinsic sequence propensities and slowing down of deamidation by local structure. NGOME-Lite can run in a proteomic analysis mode that provides the half-time of the intact form of each protein, predicted by taking into account sequence propensities and structural protection or sequence propensities only, and a structure protection factor. The detailed analysis mode also provides graphical output for all Asn residues in the query sequence. We applied NGOME-Lite to over 257,000 sequences in 38 proteomes and found that different taxa differ in their predicted deamidation dynamics. Spontaneous protein deamidation is faster in Eukarya than in Bacteria because of a higher degree of structural protection in the latter. Predicted protein deamidation half-lifes correlate with protein turnover in human, mouse, rat, C. elegans and budding yeast but not in two plants and two bacteria. NGOME-Lite is implemented in a docker container available at https://ngome.proteinphysiologylab.org.

Genome biology of the paleotetraploid perennial biomass crop Miscanthus

Miscanthus is a perennial wild grass that is of global importance for paper production, roofing, horticultural plantings, and an emerging highly productive temperate biomass crop. We report a chromosome-scale assembly of the paleotetraploid M. sinensis genome, providing a resource for Miscanthus that links its chromosomes to the related diploid Sorghum and complex polyploid sugarcanes. The asymmetric distribution of transposons across the two homoeologous subgenomes proves Miscanthus paleo-allotetraploidy and identifies several balanced reciprocal homoeologous exchanges. Analysis of M. sinensis and M. sacchariflorus populations demonstrates extensive interspecific admixture and hybridization, and documents the origin of the highly productive triploid bioenergy crop M. × giganteus. Transcriptional profiling of leaves, stem, and rhizomes over growing seasons provides insight into rhizome development and nutrient recycling, processes critical for sustainable biomass accumulation in a perennial temperate grass. The Miscanthus genome expands the power of comparative genomics to understand traits of importance to Andropogoneae grasses.

The perennial grass Miscanthus is a promising biomass crop. Here, via genomics and transcriptomics, the authors reveal its allotetraploid origin, characterize gene expression associated with rhizome development and nutrient recycling, and describe the hybrid origin of the triploid M. x giganteus.

Siberian larch (Larix sibirica Ledeb.) mitochondrial genome assembled using both short and long nucleotide sequence reads is currently the largest known mitogenome

Background

Plant mitochondrial genomes (mitogenomes) can be structurally complex while their size can vary from ~ 222 Kbp in Brassica napus to 11.3 Mbp in Silene conica. To date, in comparison with the number of plant species, only a few plant mitogenomes have been sequenced and released, particularly for conifers (the Pinaceae family). Conifers cover an ancient group of land plants that includes about 600 species, and which are of great ecological and economical value. Among them, Siberian larch (Larix sibirica Ledeb.) represents one of the keystone species in Siberian boreal forests. Yet, despite its importance for evolutionary and population studies, the mitogenome of Siberian larch has not yet been assembled and studied.

Results

Two sources of DNA sequences were used to search for mitochondrial DNA (mtDNA) sequences: mtDNA enriched samples and nucleotide reads generated in the de novo whole genome sequencing project, respectively. The assembly of the Siberian larch mitogenome contained nine contigs, with the shortest and the largest contigs being 24,767 bp and 4,008,762 bp, respectively. The total size of the genome was estimated at 11.7 Mbp. In total, 40 protein-coding, 34 tRNA, and 3 rRNA genes and numerous repetitive elements (REs) were annotated in this mitogenome. In total, 864 C-to-U RNA editing sites were found for 38 out of 40 protein-coding genes. The immense size of this genome, currently the largest reported, can be partly explained by variable numbers of mobile genetic elements, and introns, but unlikely by plasmid-related sequences. We found few plasmid-like insertions representing only 0.11% of the entire Siberian larch mitogenome.

Conclusions

Our study showed that the size of the Siberian larch mitogenome is much larger than in other so far studied Gymnosperms, and in the same range as for the annual flowering plant Silene conica (11.3 Mbp). Similar to other species, the Siberian larch mitogenome contains relatively few genes, and despite its huge size, the repeated and low complexity regions cover only 14.46% of the mitogenome sequence.

Genome wide annotation and characterization of young, intact long terminal repeat retrotransposons (In-LTR-RTs) of seven legume species

Availability of genome sequence of different legume species has provided an opportunity to characterize the abundance, distribution, and divergence of canonical intact long terminal retrotransposons (In-LTR-RT) superfamilies. Among seven legume species, Arachis ipaensis (Aip) showed the highest number of full-length canonical In-LTR-RTs (3325), followed by Glycine max (Gma, 2328), Vigna angularis (Van, 1625), Arachis durensis (Adu, 1348), Lotus japonicus (Lja, 1294), Medicago truncatula (Mtr, 788), and Circer arietinum (Car, 124). Divergence time analysis demonstrated that the amplification timeframe of LTR-RTs dramatically varied in different families. The average insertion time of Copia element varied from 0.51 (Van) to 1.37 million years ago (Mya) (Adu, and Aip), whereas that of Gypsy was between 0.22 (Mtr) and 1.82 Mya (Adu). Bayesian phylogenetic tree analysis suggested that the 1397 and 1917 reverse transcriptase (RT) domains of Copia and Gypsy families of the seven legume species were clustered into 7 and 14 major groups, respectively. The highest proportion (approximately 94.79-100%) of transposable element (TE)-associated genes assigned to pathways was mapped to metabolism-related pathways in all species. The results enabled the structural understanding of full-length In-LTR-RTs and will be valuable resource for the further study of the impact of TEs on gene structure and expression in legume species.

Oxford Nanopore sequencing: new opportunities for plant genomics?

DNA sequencing was dominated by Sanger's chain termination method until the mid-2000s, when it was progressively supplanted by new sequencing technologies that can generate much larger quantities of data in a shorter time. At the forefront of these developments, long-read sequencing technologies (third-generation sequencing) can produce reads that are several kilobases in length. This greatly improves the accuracy of genome assemblies by spanning the highly repetitive segments that cause difficulty for second-generation short-read technologies. Third-generation sequencing is especially appealing for plant genomes, which can be extremely large with long stretches of highly repetitive DNA. Until recently, the low basecalling accuracy of third-generation technologies meant that accurate genome assembly required expensive, high-coverage sequencing followed by computational analysis to correct for errors. However, today's long-read technologies are more accurate and less expensive, making them the method of choice for the assembly of complex genomes. Oxford Nanopore Technologies (ONT), a third-generation platform for the sequencing of native DNA strands, is particularly suitable for the generation of high-quality assemblies of highly repetitive plant genomes. Here we discuss the benefits of ONT, especially for the plant science community, and describe the issues that remain to be addressed when using ONT for plant genome sequencing.

Genomic re-assessment of the transposable element landscape of the potato genome

We provide a comprehensive and reliable potato TE landscape, based on a wide variety of identification tools and integrative approaches, producing clear and ready-to-use outputs for the scientific community. Transposable elements (TEs) are DNA sequences with the ability to autoreplicate and move throughout the host genome. TEs are major drivers in stress response and genome evolution. Given their significance, the development of clear and efficient TE annotation pipelines has become essential for many species. The latest de novo TE discovery tools, along with available TEs from Repbase and sRNA-seq data, allowed us to perform a reliable potato TEs detection, classification and annotation through an open-source and freely available pipeline ( https://github.com/DiegoZavallo/TE_Discovery ). Using a variety of tools, approaches and rules, we were able to provide a clearly annotated of characterized TEs landscape. Additionally, we described the distribution of the different types of TEs across the genome, where LTRs and MITEs present a clear clustering pattern in pericentromeric and subtelomeric/telomeric regions respectively. Finally, we analyzed the insertion age and distribution of LTR retrotransposon families which display a distinct pattern between the two major superfamilies. While older Gypsy elements concentrated around heterochromatic regions, younger Copia elements located predominantly on euchromatic regions. Overall, we delivered not only a reliable, ready-to-use potato TE annotation files, but also all the necessary steps to perform de novo detection for other species.

Molecular characterization of ascorbate peroxidase (APX) and APX-related (APX-R) genes in Triticum aestivum L

Ascorbate peroxidases (APXs) are heme-dependent H₂O₂ scavenging enzymes involved in myriad biological processes. Herein, a total of 21 TaAPX and six TaAPX-R genes were identified from the A, B and D sub-genomes of Triticum aestivum L. The occurrence of three paralogous gene pairs with unequal evolutionary rate suggested functional divergence. The phylogenetic analysis formed four distinct clades having conserved gene and protein architecture, and sub-cellular localization. The tertiary structure analysis revealed the presence of helices and coils and residues involved in ligand binding. Transcriptional profiling of each TaAPX and TaAPX-R gene suggested their specific role during development and stress response. Modulated transcript expression and APX enzyme activity during various stress conditions indicated their role in stress response. Interaction analyses suggested their association with other genes, miRNAs and various legends. The present study reported numerous features of these genes, and may provide a platform for their detailed functional characterization in future studies.

Mobile genomics: tools and techniques for tackling transposons

Next-generation sequencing approaches have fundamentally changed the types of questions that can be asked about gene function and regulation. With the goal of approaching truly genome-wide quantifications of all the interaction partners and downstream effects of particular genes, these quantitative assays have allowed for an unprecedented level of detail in exploring biological interactions. However, many challenges remain in our ability to accurately describe and quantify the interactions that take place in those hard to reach and extremely repetitive regions of our genome comprised mostly of transposable elements (TEs). Tools dedicated to TE-derived sequences have lagged behind, making the inclusion of these sequences in genome-wide analyses difficult. Recent improvements, both computational and experimental, allow for the better inclusion of TE sequences in genomic assays and a renewed appreciation for the importance of TE biology. This review will discuss the recent improvements that have been made in the computational analysis of TE-derived sequences as well as the areas where such analysis still proves difficult. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Miniature inverted-repeat transposable elements (MITEs), derived insertional polymorphism as a tool of marker systems for molecular plant breeding

Plant molecular breeding is expected to give significant gains in cultivar development through development and utilization of suitable molecular marker systems for genetic diversity analysis, rapid DNA fingerprinting, identification of true hybrids, trait mapping and marker-assisted selection. Transposable elements (TEs) are the most abundant component in a genome and being used as genetic markers in the plant molecular breeding. Here, we review on the high copious transposable element belonging to class-II DNA TEs called "miniature inverted-repeat transposable elements" (MITEs). MITEs are ubiquitous, short and non-autonomous DNA transposable elements which have a tendency to insert into genes and genic regions have paved a way for the development of functional DNA marker systems in plant genomes. This review summarises the characteristics of MITEs, principles and methodologies for development of MITEs based DNA markers, bioinformatics tools and resources for plant MITE discovery and their utilization in crop improvement.

Whole-genome sequencing and analysis of the Chinese herbal plant Gelsemium elegans

BACKGROUND

Gelsemium elegans (G. elegans) (2n = 2x = 16) is genus of flowering plants belonging to the Gelsemicaeae family.

METHOD

Here, a high-quality genome assembly using the Oxford Nanopore Technologies (ONT) platform and high-throughput chromosome conformation capture techniques (Hi-C) were used.

RESULTS

A total of 56.11 Gb of raw GridION X5 platform ONT reads (6.23 Gb per cell) were generated. After filtering, 53.45 Gb of clean reads were obtained, giving 160 × coverage depth. The de novo genome assemblies 335.13 Mb, close to the 338 Mb estimated by k-mer analysis, was generated with contig N50 of 10.23 Mb. The vast majority (99.2%) of the G. elegans assembled sequence was anchored onto 8 pseudo-chromosomes. The genome completeness was then evaluated and 1338 of the 1440 conserved genes (92.9%) could be found in the assembly. Genome annotation revealed that 43.16% of the G. elegans genome is composed of repetitive elements and 23.9% is composed of long terminal repeat elements. We predicted 26,768 protein-coding genes, of which 84.56% were functionally annotated.

CONCLUSION

The genomic sequences of G. elegans could be a valuable source for comparative genomic analysis in the Gelsemicaeae family and will be useful for understanding the phylogenetic relationships of the indole alkaloid metabolism.

Chromosomal-level assembly of Juglans sigillata genome using Nanopore, BioNano, and Hi-C analysis

BACKGROUND

Juglans sigillata, or iron walnut, belonging to the order Juglandales, is an economically important tree species in Asia, especially in the Yunnan province of China. However, little research has been conducted on J. sigillata at the molecular level, which hinders understanding of its evolution, speciation, and synthesis of secondary metabolites, as well as its wide adaptability to its plateau environment. To address these issues, a high-quality reference genome of J. sigillata would be useful.

FINDINGS

To construct a high-quality reference genome for J. sigillata, we first generated 38.0 Gb short reads and 66.31 Gb long reads using Illumina and Nanopore sequencing platforms, respectively. The sequencing data were assembled into a 536.50-Mb genome assembly with a contig N50 length of 4.31 Mb. Additionally, we applied BioNano technology to identify contacts among contigs, which were then used to assemble contigs into scaffolds, resulting in a genome assembly with scaffold N50 length of 16.43 Mb and contig N50 length of 4.34 Mb. To obtain a chromosome-level genome assembly, we constructed 1 Hi-C library and sequenced 79.97 Gb raw reads using the Illumina HiSeq platform. We anchored ∼93% of the scaffold sequences into 16 chromosomes and evaluated the quality of our assembly using the high contact frequency heat map. Repetitive elements account for 50.06% of the genome, and 30,387 protein-coding genes were predicted from the genome, of which 99.8% have been functionally annotated. The genome-wide phylogenetic tree indicated an estimated divergence time between J. sigillata and Juglans regia of 49 million years ago on the basis of single-copy orthologous genes.

CONCLUSIONS

We provide the first chromosome-level genome for J. sigillata. It will lay a valuable foundation for future research on the genetic improvement of J. sigillata.

Thaumatin-like protein kinases: Molecular characterization and transcriptional profiling in five cereal crops

Thaumatin-like protein kinases (TLPKs) are defense related proteins having antimicrobial property. Herein, we identified two TLPKs in the genome of Brachypodium distachyon and Oryza sativa, four in Hordeum vulgare and Sorghum bicolor, and 16 in Triticum aestivum. All the TLPKs were located at only one chromosome in each plant except T. aestivum, where they were located on chromosome 2 and chromosome 3. Paralogous analysis suggested the occurrence of one duplication event (DE) in each B. distachyon and O. sativa, two in H. vulgare while four DEs in T. aestivum genome during the evolution of TLPKs. The majority of TLPKs were intron less, while a few contains one or two introns. The introns were found in each 0, 1 and 2 phase. Protein structure analysis suggested the occurrence of a thaumatin and a kinase domain with a transmembrane (TM) helix in each TLPK. Further, a thaumatin family signature motif "GX[GF]XCXT[GA]DCX(1,2)GX(2,3)C", a "REDDD" motif and 16 cysteine residues were found conserved in the majority of TLPKs. Expression analysis indicated variable expression of TLPKs in various tissues of different cereal crops. They were high expressing in reproductive tissues in B. distachyon, while in leaves in T. aestivum. Modulated expression of TaTLPKs in the presence of fungal pathogen, and heat, drought and salt stress in T. aestivum suggested their roles in stress response. Co-expression analysis showed interaction of TLPKs with various development and stress related genes. The results indicated diverse roles of TLPKs, which can be utilized for the development of eco-friendly pest resistant crops in future.

Thaumatin-like protein kinases: Molecular characterization and transcriptional profiling in five cereal crops

Thaumatin-like protein kinases (TLPKs) are defense related proteins having antimicrobial property. Herein, we identified two TLPKs in the genome of Brachypodium distachyon and Oryza sativa, four in Hordeum vulgare and Sorghum bicolor, and 16 in Triticum aestivum. All the TLPKs were located at only one chromosome in each plant except T. aestivum, where they were located on chromosome 2 and chromosome 3. Paralogous analysis suggested the occurrence of one duplication event (DE) in each B. distachyon and O. sativa, two in H. vulgare while four DEs in T. aestivum genome during the evolution of TLPKs. The majority of TLPKs were intron less, while a few contains one or two introns. The introns were found in each 0, 1 and 2 phase. Protein structure analysis suggested the occurrence of a thaumatin and a kinase domain with a transmembrane (TM) helix in each TLPK. Further, a thaumatin family signature motif "GX[GF]XCXT[GA]DCX(1,2)GX(2,3)C", a "REDDD" motif and 16 cysteine residues were found conserved in the majority of TLPKs. Expression analysis indicated variable expression of TLPKs in various tissues of different cereal crops. They were high expressing in reproductive tissues in B. distachyon, while in leaves in T. aestivum. Modulated expression of TaTLPKs in the presence of fungal pathogen, and heat, drought and salt stress in T. aestivum suggested their roles in stress response. Co-expression analysis showed interaction of TLPKs with various development and stress related genes. The results indicated diverse roles of TLPKs, which can be utilized for the development of eco-friendly pest resistant crops in future.

In silico determination of transposon-derived miRNAs and targets in Aegilops species

Transposable elements (TEs) are found almost in all living organism, shaping organisms' genomes. miRNAs are noncoding RNA types which are especially important in gene expression regulations. Many previously determined plant miRNAs are identical/homologous to transposons (TE-MIR). The aim of this study was computational characterization of novel TE-related miRNAs and their targets in Aegilops genome by using stringent criteria. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed by BLAST2GO. Seventeen novel TE-related miRNAs in Aegilops genome were identified for the first time. GO analyses indicated that 40 targets played different roles in biological processes, cellular components and molecular functions. Moreover, these genes were involved in 10 metabolic pathways such as purine metabolism, nitrogen metabolism, oxidative phosphorylation, etc. as a result of KEGG analyses. Identification of miRNAs and their targets are significant to understand miRNA-TEs relationships and even how TEs affect plant growth and development. Obtaining results of this study are expected to provide possible new insight into Aegilops and its related species, wheat, with respect to miRNAs evolution and domestication.Communicated by Ramaswamy H. Sarma.

Structural Characterization of ABCB1, the Gene Underlying the d2 Dwarf Phenotype in Pearl Millet, Cenchrus Americanus (L.) Morrone

Pearl millet is an important food crop in arid and semi-arid regions of South Asia and sub-Saharan Africa and is grown in Australia and the United States as a summer fodder crop. The d2 dwarf germplasm has been widely used in the last half-century to develop high-performing pearl millet hybrids. We previously mapped the d2 phenotype to a 1.6 cM region in linkage group (LG) 4 and identified the ABCB1 gene as a candidate underlying the trait. Here, we report the sequence, structure and expression of ABCB1 in tall (D2D2) and d2 dwarf (d2d2) germplasm. The ABCB1 allele in d2 dwarfs differs from that in tall inbreds by the presence of two different high copy transposable elements, one in the coding region and the second located 664 bp upstream of the ATG start codon. These transposons were present in all d2 dwarfs tested that were reported to be of independent origin and absent in the analyzed wild-type tall germplasm. We also compared the expression profile of this gene in different organs of multiple tall and d2 dwarf inbreds, including the near-isogenic inbreds at the d2 locus, Tift 23B (D2D2) and Tift 23DB (d2d2). Heterologous transformation of the tall (Ca_ABCB1) and the d2 dwarf (Ca_abcb1) pearl millet alleles in the Arabidopsis double mutant abcb1abcb19 showed that the pearl millet D2 but not the d2 allele complements the Arabidopsis abcb1 mutation. Our studies also show the importance of the COOH-terminal 22 amino acids of the ABCB1 protein in either protein function or stability.

Interspecies association mapping links reduced CG to TG substitution rates to the loss of gene-body methylation

Comparative genomics can unravel the genetic basis of species differences; however, successful reports on quantitative traits are still scarce. Here we present genome assemblies of 31 so-far unassembled Brassicaceae plant species and combine them with 16 previously published assemblies to establish the Brassicaceae Diversity Panel. Using a new interspecies association strategy for quantitative traits, we found a so-far unknown association between the unexpectedly high variation in CG to TG substitution rates in genes and the absence of CHROMOMETHYLASE3 (CMT3) orthologues. Low substitution rates were associated with the loss of CMT3, while species with conserved CMT3 orthologues showed high substitution rates. Species without CMT3 also lacked gene-body methylation (gbM), suggesting an evolutionary trade-off between the unknown function of gbM and low substitution rates in Brassicaceae, possibly due to low mutability of non-methylated cytosines.

Molecular characterization of bacterial leaf streak resistance in hard winter wheat

Bacterial leaf streak (BLS) caused by Xanthomonas campestris pv. translucens is one of the major bacterial diseases threatening wheat production in the United States Northern Great Plains (NGP) region. It is a sporadic but widespread wheat disease that can cause significant loss in grain yield and quality. Identification and characterization of genomic regions in wheat that confer resistance to BLS will help track resistance genes/QTLs in future wheat breeding. In this study, we evaluated a hard winter wheat association mapping panel (HWWAMP) containing 299 hard winter wheat lines from the US hard winter wheat growing region for their reactions to BLS. We observed a range of BLS responses among the lines, importantly, we identified ten genotypes that showed a resistant reaction both in greenhouse and field evaluation. -Genome-wide association analysis with 15,990 SNPs was conducted using an exponentially compressed mixed linear model. Five genomic regions (p < 0.001) that regulate the resistance to BLS were identified on chromosomes 1AL, 1BS, 3AL, 4AL, and 7AS. The QTLs Q.bls.sdsu-1AL, Q.bls.sdsu-1BS, Q.bls.sdsu-3AL, Q.bls.sdsu-4AL, and Q.bls.sdsu-7AS explain a total of 42% of the variation. In silico analysis of sequences in the candidate regions on chromosomes 1AL, 1BS, 3AL, 4AL, and 7AS identified 10, 25, 22, eight, and nine genes, respectively with known plant defense-related functions. Comparative analysis with rice showed two syntenic regions in rice that harbor genes for bacterial leaf streak resistance. The ten BLS resistant genotypes and SNP markers linked to the QTLs identified in our study could facilitate breeding for BLS resistance in winter wheat.

Wavelet-Based Genomic Signal Processing for Centromere Identification and Hypothesis Generation

Various ‘omics data types have been generated for Populus trichocarpa, each providing a layer of information which can be represented as a density signal across a chromosome. We make use of genome sequence data, variants data across a population as well as methylation data across 10 different tissues, combined with wavelet-based signal processing to perform a comprehensive analysis of the signature of the centromere in these different data signals, and successfully identify putative centromeric regions in P. trichocarpa from these signals. Furthermore, using SNP (single nucleotide polymorphism) correlations across a natural population of P. trichocarpa, we find evidence for the co-evolution of the centromeric histone CENH3 with the sequence of the newly identified centromeric regions, and identify a new CENH3 candidate in P. trichocarpa.

A high-quality apple genome assembly reveals the association of a retrotransposon and red fruit colour

A complete and accurate genome sequence provides a fundamental tool for functional genomics and DNA-informed breeding. Here, we assemble a high-quality genome (contig N50 of 6.99 Mb) of the apple anther-derived homozygous line HFTH1, including 22 telomere sequences, using a combination of PacBio single-molecule real-time (SMRT) sequencing, chromosome conformation capture (Hi-C) sequencing, and optical mapping. In comparison to the Golden Delicious reference genome, we identify 18,047 deletions, 12,101 insertions and 14 large inversions. We reveal that these extensive genomic variations are largely attributable to activity of transposable elements. Interestingly, we find that a long terminal repeat (LTR) retrotransposon insertion upstream of MdMYB1, a core transcriptional activator of anthocyanin biosynthesis, is associated with red-skinned phenotype. This finding provides insights into the molecular mechanisms underlying red fruit coloration, and highlights the utility of this high-quality genome assembly in deciphering agriculturally important trait in apple.

Existing apple genome assemblies all derive from Golden Delicious. Here, the authors combine different sequencing technologies to assemble a high quality genome of an anther-derived homozygous genotype HFTH1 and find the association of a retrotransposon and red fruit colour.

Gossypium barbadense and Gossypium hirsutum genomes provide insights into the origin and evolution of allotetraploid cotton

Allotetraploid cotton is an economically important natural-fiber-producing crop worldwide. After polyploidization, Gossypium hirsutum L. evolved to produce a higher fiber yield and to better survive harsh environments than Gossypium barbadense, which produces superior-quality fibers. The global genetic and molecular bases for these interspecies divergences were unknown. Here we report high-quality de novo-assembled genomes for these two cultivated allotetraploid species with pronounced improvement in repetitive-DNA-enriched centromeric regions. Whole-genome comparative analyses revealed that species-specific alterations in gene expression, structural variations and expanded gene families were responsible for speciation and the evolutionary history of these species. These findings help to elucidate the evolution of cotton genomes and their domestication history. The information generated not only should enable breeders to improve fiber quality and resilience to ever-changing environmental conditions but also can be translated to other crops for better understanding of their domestication history and use in improvement.

Gossypium barbadense and Gossypium hirsutum genomes provide insights into the origin and evolution of allotetraploid cotton

Allotetraploid cotton is an economically important natural-fiber-producing crop worldwide. After polyploidization, Gossypium hirsutum L. evolved to produce a higher fiber yield and to better survive harsh environments than Gossypium barbadense, which produces superior-quality fibers. The global genetic and molecular bases for these interspecies divergences were unknown. Here we report high-quality de novo-assembled genomes for these two cultivated allotetraploid species with pronounced improvement in repetitive-DNA-enriched centromeric regions. Whole-genome comparative analyses revealed that species-specific alterations in gene expression, structural variations and expanded gene families were responsible for speciation and the evolutionary history of these species. These findings help to elucidate the evolution of cotton genomes and their domestication history. The information generated not only should enable breeders to improve fiber quality and resilience to ever-changing environmental conditions but also can be translated to other crops for better understanding of their domestication history and use in improvement.