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.

Coffea cytogenetics: from the first karyotypes to the meeting with genomics

Coffea karyotype organization and evolution has been uncovered by classical cytogenetics and cytogenomics. We revisit these discoveries and present new karyotype data. Coffea possesses ~ 124 species, including C. arabica and C. canephora responsible for commercial coffee production. We reviewed the Coffea cytogenetics, from the first chromosome counting, encompassing the karyotype characterization, chromosome DNA content, and mapping of chromosome portions and DNA sequences, until the integration with genomics. We also showed new data about Coffea karyotype. The 2n chromosome number evidenced the diploidy of almost all Coffea, and the C. arabica tetraploidy, as well as the polyploidy of other hybrids. Since then, other genomic similarities and divergences among the Coffea have been shown by karyotype morphology, nuclear and chromosomal C-value, AT and GC rich chromosome portions, and repetitive sequence and gene mapping. These cytogenomic data allowed us to know and understand the phylogenetic relations in Coffea, as well as their ploidy level and genomic origin, highlighting the relatively recent allopolyploidy. In addition to the euploidy, the role of the mobile elements in Coffea diversification is increasingly more evident, and the comparative analysis of their structure and distribution on the genome of different species is in the spotlight for future research. An integrative look at all these data is fundamental for a deeper understanding of Coffea karyotype evolution, including the key role of polyploidy in C. arabica origin. The 'Híbrido de Timor', a recent natural allotriploid, is also in the spotlight for its potential as a source of resistance genes and model for plant polyploidy research. Considering this, we also present some unprecedented results about the exciting evolutionary history of these polyploid Coffea.

Innate, translation-dependent silencing of an invasive transposon in Arabidopsis

Co‐evolution between hosts’ and parasites’ genomes shapes diverse pathways of acquired immunity based on silencing small (s)RNAs. In plants, sRNAs cause heterochromatinization, sequence degeneration, and, ultimately, loss of autonomy of most transposable elements (TEs). Recognition of newly invasive plant TEs, by contrast, involves an innate antiviral‐like silencing response. To investigate this response’s activation, we studied the single‐copy element EVADÉ (EVD), one of few representatives of the large Ty1/Copia family able to proliferate in Arabidopsis when epigenetically reactivated. In Ty1/Copia elements, a short subgenomic mRNA (shGAG) provides the necessary excess of structural GAG protein over the catalytic components encoded by the full‐length genomic flGAG‐POL. We show here that the predominant cytosolic distribution of shGAG strongly favors its translation over mostly nuclear flGAG‐POL. During this process, an unusually intense ribosomal stalling event coincides with mRNA breakage yielding unconventional 5’OH RNA fragments that evade RNA quality control. The starting point of sRNA production by RNA‐DEPENDENT‐RNA‐POLYMERASE‐6 (RDR6), exclusively on shGAG, occurs precisely at this breakage point. This hitherto‐unrecognized “translation‐dependent silencing” (TdS) is independent of codon usage or GC content and is not observed on TE remnants populating the Arabidopsis genome, consistent with their poor association, if any, with polysomes. We propose that TdS forms a primal defense against EVD de novo invasions that underlies its associated sRNA pattern.

The Dynamism of Transposon Methylation for Plant Development and Stress Adaptation

Plant development processes are regulated by epigenetic alterations that shape nuclear structure, gene expression, and phenotypic plasticity; these alterations can provide the plant with protection from environmental stresses. During plant growth and development, these processes play a significant role in regulating gene expression to remodel chromatin structure. These epigenetic alterations are mainly regulated by transposable elements (TEs) whose abundance in plant genomes results in their interaction with genomes. Thus, TEs are the main source of epigenetic changes and form a substantial part of the plant genome. Furthermore, TEs can be activated under stress conditions, and activated elements cause mutagenic effects and substantial genetic variability. This introduces novel gene functions and structural variation in the insertion sites and primarily contributes to epigenetic modifications. Altogether, these modifications indirectly or directly provide the ability to withstand environmental stresses. In recent years, many studies have shown that TE methylation plays a major role in the evolution of the plant genome through epigenetic process that regulate gene imprinting, thereby upholding genome stability. The induced genetic rearrangements and insertions of mobile genetic elements in regions of active euchromatin contribute to genome alteration, leading to genomic stress. These TE-mediated epigenetic modifications lead to phenotypic diversity, genetic variation, and environmental stress tolerance. Thus, TE methylation is essential for plant evolution and stress adaptation, and TEs hold a relevant military position in the plant genome. High-throughput techniques have greatly advanced the understanding of TE-mediated gene expression and its associations with genome methylation and suggest that controlled mobilization of TEs could be used for crop breeding. However, development application in this area has been limited, and an integrated view of TE function and subsequent processes is lacking. In this review, we explore the enormous diversity and likely functions of the TE repertoire in adaptive evolution and discuss some recent examples of how TEs impact gene expression in plant development and stress adaptation.

Phenotype study of multifoliolate leaf formation in Trifolium alexandrinum L

Background

The genus Trifolium is characterized by typical trifoliolate leaves. Alterations in leaf formats from trifoliolate to multifoliolate, i.e., individual plants bearing trifoliolate, quadrifoliolate, pentafoliolate or more leaflets, were previously reported among many species of the genus. The study is an attempt to develop pure pentafoliolate plants of T. alexandrinum and to understand its genetic control.

Methods

The experimental material consisted of two populations of T. alexandrinum with multifoliolate leaf expression, i.e.,interspecific hybrid progenies of T. alexandrinum with T. apertum, and T. alexandrinum genotype Penta-1. Penetrance of the multifoliolate trait was observed among multifoliolate and trifoliolate plant progenies. In vitro culture and regeneration of plantlets from the axillary buds from different plant sources was also attempted.

Results

The inheritance among a large number of plant progenies together with in vitro micro-propagation results did not establish a definite pattern. The multifoliolate leaf formation was of chimeric nature, i.e., more than one leaf format appearing on individual branches. Reversal to normal trifoliolate from multifoliolate was also quite common. Penetrance and expression of multifoliolate leaf formation was higher among the plants raised from multifoliolate plants. Multifoliolate and pure pentafoliolate plants were observed in the progenies of pure trifoliolate plants and vice-versa. There was an apparent increase in the pentafoliolate leaf formation frequency over the years due to targeted selection. A few progenies of the complete pentafoliolate plants in the first year were true breeding in the second year. Frequency of plantlets with multifoliolate leaf formation was also higher in in vitro axillary bud multiplication when the explant bud was excised from the multifoliolate leaf node.

Conclusion

Number of leaflets being a discrete variable, occurrence of multifoliolate leaves on individual branches, reversal of leaf formats on branches and developing true breeding pentafoliolates were the factors leading to a hypothesis beyond normal Mendelian inheritance. Transposable elements (TEs) involved in leaf development in combination with epigenetics were probably responsible for alterations in the expression of leaflet number. Putative TE’s movement owing to chromosomal rearrangements possibly resulted in homozygous pentafoliolate trait with evolutionary significance. The hypothesis provides a new insight into understanding the genetic control of this trait in T. alexandrinum and may also be useful in other Trifolium species where such observations are reported.

Evolutionary dynamics of transposable elements and satellite DNAs in polyploid Spartina species

Repeated sequences and polyploidy play a central role in plant genome dynamics. Here, we analyze the evolutionary dynamics of repeats in tetraploid and hexaploid Spartina species that diverged during the last 10 million years within the Chloridoideae, one of the poorest investigated grass lineages. From high-throughput genome sequencing, we annotated Spartina repeats and determined what sequence types account for the genome size variation among species. We examined whether differential genome size evolution correlated with ploidy levels and phylogenetic relationships. We also examined the tempo of repeat sequence dynamics associated with allopatric speciation over the last 3-6 million years between hexaploid species that diverged on the American and European Atlantic coasts and tetraploid species from North and South America. The tetraploid S. spartinae, whose phylogenetic placement has been debated, exhibits a similar repeat content as hexaploid species, suggesting common ancestry. Genome expansion or contraction resulting from repeat dynamics seems to be explained mostly by the contrasting divergence times between species, rather than by genome changes triggered by ploidy level change per se. One 370 bp satellite may be exhibiting 'meiotic drive' and driving chromosome evolution in S. alterniflora. Our results provide crucial insights for investigating the genetic and epigenetic consequences of such differential repeat dynamics on the ecology and distribution of the meso- and neopolyploid Spartina species.

Shaping Plant Adaptability, Genome Structure and Gene Expression through Transposable Element Epigenetic Control: Focus on Methylation

In plants, transposable elements (TEs) represent a large fraction of the genome, with potential to alter gene expression and produce genomic rearrangements. Epigenetic control of TEs is often used to stop unrestricted movement of TEs that would result in detrimental effects due to insertion in essential genes. The current review focuses on the effects of methylation on TEs and their genomic context, and how this type of epigenetic control affects plant adaptability when plants are faced with different stresses and changes. TEs mobilize in response to stress elicitors, including biotic and abiotic cues, but also developmental transitions and ‘genome shock’ events like polyploidization. These events transitionally lift TE repression, allowing TEs to move to new genomic locations. When TEs fall close to genes, silencing through methylation can spread to nearby genes, resulting in lower gene expression. The presence of TEs in gene promoter regions can also confer stress inducibility modulated through alternative methylation and demethylation of the TE. Bursts of transposition triggered by events of genomic shock can increase genome size and account for differences seen during polyploidization or species divergence. Finally, TEs have evolved several mechanisms to suppress their own repression, including the use of microRNAs to control genes that promote methylation. The interplay between silencing, transient TE activation, and purifying selection allows the genome to use TEs as a reservoir of potential beneficial modifications but also keeps TEs under control to stop uncontrolled detrimental transposition.

Unique features of the loblolly pine (Pinus taeda L.) megagenome revealed through sequence annotation

The largest genus in the conifer family Pinaceae is Pinus, with over 100 species. The size and complexity of their genomes (∼20-40 Gb, 2n = 24) have delayed the arrival of a well-annotated reference sequence. In this study, we present the annotation of the first whole-genome shotgun assembly of loblolly pine (Pinus taeda L.), which comprises 20.1 Gb of sequence. The MAKER-P annotation pipeline combined evidence-based alignments and ab initio predictions to generate 50,172 gene models, of which 15,653 are classified as high confidence. Clustering these gene models with 13 other plant species resulted in 20,646 gene families, of which 1554 are predicted to be unique to conifers. Among the conifer gene families, 159 are composed exclusively of loblolly pine members. The gene models for loblolly pine have the highest median and mean intron lengths of 24 fully sequenced plant genomes. Conifer genomes are full of repetitive DNA, with the most significant contributions from long-terminal-repeat retrotransposons. In depth analysis of the tandem and interspersed repetitive content yielded a combined estimate of 82%.

Insights into the loblolly pine genome: characterization of BAC and fosmid sequences

Despite their prevalence and importance, the genome sequences of loblolly pine, Norway spruce, and white spruce, three ecologically and economically important conifer species, are just becoming available to the research community. Following the completion of these large assemblies, annotation efforts will be undertaken to characterize the reference sequences. Accurate annotation of these ancient genomes would be aided by a comprehensive repeat library; however, few studies have generated enough sequence to fully evaluate and catalog their non-genic content. In this paper, two sets of loblolly pine genomic sequence, 103 previously assembled BACs and 90,954 newly sequenced and assembled fosmid scaffolds, were analyzed. Together, this sequence represents 280 Mbp (roughly 1% of the loblolly pine genome) and one of the most comprehensive studies of repetitive elements and genes in a gymnosperm species. A combination of homology and de novo methodologies were applied to identify both conserved and novel repeats. Similarity analysis estimated a repetitive content of 27% that included both full and partial elements. When combined with the de novo investigation, the estimate increased to almost 86%. Over 60% of the repetitive sequence consists of full or partial LTR (long terminal repeat) retrotransposons. Through de novo approaches, 6,270 novel, full-length transposable element families and 9,415 sub-families were identified. Among those 6,270 families, 82% were annotated as single-copy. Several of the novel, high-copy families are described here, with the largest, PtPiedmont, comprising 133 full-length copies. In addition to repeats, analysis of the coding region reported 23 full-length eukaryotic orthologous proteins (KOGS) and another 29 novel or orthologous genes. These discoveries, along with other genomic resources, will be used to annotate conifer genomes and address long-standing questions about gymnosperm evolution.

Identification, characterization and distribution of transposable elements in the flax (Linum usitatissimum L.) genome

BACKGROUND

Flax (Linum usitatissimum L.) is an important crop for the production of bioproducts derived from its seed and stem fiber. Transposable elements (TEs) are widespread in plant genomes and are a key component of their evolution. The availability of a genome assembly of flax (Linum usitatissimum) affords new opportunities to explore the diversity of TEs and their relationship to genes and gene expression.

RESULTS

Four de novo repeat identification algorithms (PILER, RepeatScout, LTR_finder and LTR_STRUC) were applied to the flax genome assembly. The resulting library of flax repeats was combined with the RepBase Viridiplantae division and used with RepeatMasker to identify TEs coverage in the genome. LTR retrotransposons were the most abundant TEs (17.2% genome coverage), followed by Long Interspersed Nuclear Element (LINE) retrotransposons (2.10%) and Mutator DNA transposons (1.99%). Comparison of putative flax TEs to flax transcript databases indicated that TEs are not highly expressed in flax. However, the presence of recent insertions, defined by 100% intra-element LTR similarity, provided evidence for recent TE activity. Spatial analysis showed TE-rich regions, gene-rich regions as well as regions with similar genes and TE density. Monte Carlo simulations for the 71 largest scaffolds (≥ 1 Mb each) did not show any regional differences in the frequency of TE overlap with gene coding sequences. However, differences between TE superfamilies were found in their proximity to genes. Genes within TE-rich regions also appeared to have lower transcript expression, based on EST abundance. When LTR elements were compared, Copia showed more diversity, recent insertions and conserved domains than the Gypsy, demonstrating their importance in genome evolution.

CONCLUSIONS

The calculated 23.06% TE coverage of the flax WGS assembly is at the low end of the range of TE coverages reported in other eudicots, although this estimate does not include TEs likely found in unassembled repetitive regions of the genome. Since enrichment for TEs in genomic regions was associated with reduced expression of neighbouring genes, and many members of the Copia LTR superfamily are inserted close to coding regions, we suggest Copia elements have a greater influence on recent flax genome evolution while Gypsy elements have become residual and highly mutated.

C. elegans germ cells show temperature and age-dependent expression of Cer1, a Gypsy/Ty3-related retrotransposon

Virus-like particles (VLPs) have not been observed in Caenorhabditis germ cells, although nematode genomes contain low numbers of retrotransposon and retroviral sequences. We used electron microscopy to search for VLPs in various wild strains of Caenorhabditis, and observed very rare candidate VLPs in some strains, including the standard laboratory strain of C. elegans, N2. We identified the N2 VLPs as capsids produced by Cer1, a retrotransposon in the Gypsy/Ty3 family of retroviruses/retrotransposons. Cer1 expression is age and temperature dependent, with abundant expression at 15°C and no detectable expression at 25°C, explaining how VLPs escaped detection in previous studies. Similar age and temperature-dependent expression of Cer1 retrotransposons was observed for several other wild strains, indicating that these properties are common, if not integral, features of this retroelement. Retrotransposons, in contrast to DNA transposons, have a cytoplasmic stage in replication, and those that infect non-dividing cells must pass their genomic material through nuclear pores. In most C. elegans germ cells, nuclear pores are largely covered by germline-specific organelles called P granules. Our results suggest that Cer1 capsids target meiotic germ cells exiting pachytene, when free nuclear pores are added to the nuclear envelope and existing P granules begin to be removed. In pachytene germ cells, Cer1 capsids concentrate away from nuclei on a subset of microtubules that are exceptionally resistant to microtubule inhibitors; the capsids can aggregate these stable microtubules in older adults, which exhibit a temperature-dependent decrease in egg viability. When germ cells exit pachytene, the stable microtubules disappear and capsids redistribute close to nuclei that have P granule-free nuclear pores. This redistribution is microtubule dependent, suggesting that capsids that are released from stable microtubules transfer onto new, dynamic microtubules to track toward nuclei. These studies introduce C. elegans as a model to study the interplay between retroelements and germ cell biology.

Retrotransposons and retroviruses pose enormous threats to animal and plants because of their ability to insert into host genes. Retroelements that replicate in germ cells can, if left unchecked, expand exponentially in the host genome. C. elegans has proven to be an exceptional model system for studying many facets of cell and molecular biology, and the genome contains both retrotransposon and retroviral sequences. However, no virus-like particles have been observed in C. elegans germ cells. We show here that Cer1, an endogenous Gypsy/Ty3 class retrotransposon, is expressed at very high levels in C. elegans germ cells, but escaped detection in previous studies because its expression is both temperature and age dependent. These studies reveal new aspects of microtubule regulation in C. elegans that the retroelement appears to exploit to navigate the germ cell cytoplasm, and demonstrate the power of C. elegans for studying host/pathogen interactions in germ cell biology.