Modeling transposable elements dynamics during polyploidization in plants

In this work we study the proliferation of transposable elements (TEs) and the epigenetic response of plants during the process of polyploidization. Through a deterministic model, expanding on our previous work on TE proliferation under epigenetic regulation, we study the long-term TE distribution and TE stability in the subgenomes of both autopolyploids and allopolyploids. We also explore different small-interfering RNA (siRNA) action modes on the subgenomes, including a model where siRNAs are not directed to specific genomes and one where siRNAs are directed - i.e. more active - in subgenomes with higher TE loads. In the autopolyploid case, we find long-term stable equilbria that tend to equilibrate the number of active TEs between subgenomes. In the allopolyploid case, directed siRNA action is fundamental to avoid a "winner takes all" outcome of the competition between the TEs in the different subgenomes. We also show that decaying oscillations in the number of TEs occur naturally in all cases, perhaps explaining some of the observed features of 'genomic shock' after hybridization events, and that the balance in the dynamics of the different types of siRNA is determinant for the synchronization of these oscillations.

Genome Editing by CRISPR/Cas9 in Polyploids

CRISPR/Cas system has been widely used for genome editing in the past few years. Even though it has been performed in many polyploid species to date, its efficient accomplishment in these organisms is still a challenge. The presence of multiple homoeologous genes as targets for their editing requires more rigorous work and specific needs to assess successful genome editing. Here, we describe a general stepwise protocol to select target sites, design sgRNAs, indicate vector requirements, and screen CRISPR/Cas9-mediated genome editing in polyploid species.

Pangenome Analysis of Plant Transcripts and Coding Sequences

The pangenome of a species is the sum of the genomes of its individuals. As coding sequences often represent only a small fraction of each genome, analyzing the pangene set can be a cost-effective strategy for plants with large genomes or highly heterozygous species. Here, we describe a step-by-step protocol to analyze plant pangene sets with the software GET_HOMOLOGUES-EST . After a short introduction, where the main concepts are illustrated, the remaining sections cover the installation and typical operations required to analyze and annotate pantranscriptomes and gene sets of plants. The recipes include instructions on how to call core and accessory genes, how to compute a presence-absence pangenome matrix, and how to identify and analyze private genes, present only in some genotypes. Downstream phylogenetic analyses are also discussed.

Scripting Analyses of Genomes in Ensembl Plants

Ensembl Plants ( http://plants.ensembl.org ) offers genome-scale information for plants, with four releases per year. As of release 47 (April 2020) it features 79 species and includes genome sequence, gene models, and functional annotation. Comparative analyses help reconstruct the evolutionary history of gene families, genomes, and components of polyploid genomes. Some species have gene expression baseline reports or variation across genotypes. While the data can be accessed through the Ensembl genome browser, here we review specifically how our plant genomes can be interrogated programmatically and the data downloaded in bulk. These access routes are generally consistent across Ensembl for other non-plant species, including plant pathogens, pests, and pollinators.

Integrated transcriptome and miRNA sequencing approaches provide insights into salt tolerance in allotriploid Populus cathayana

Some salt-stress responsive DEGs, mainly involved in ion transmembrane transport, hormone regulation, antioxidant system, osmotic regulation, and some miRNA jointly regulated the salt response process in allotriploid Populus cathayana. The molecular mechanism of plant polyploid stress resistance has been a hot topic in biological research. In this study, Populus diploids and first division restitution (FDR) and second division restitution (SDR) triploids were selected as research materials. All materials were treated with 70 mM NaCl solutions for 30 days in the same pot environment. We observed the growth state of triploids and diploids and determined the ratio of potassium and sodium ions, peroxidase (POD) activity, proline content, and ABA and jasmonic acid (JA) hormone content in leaves in the same culture environment with the same concentration of NaCl solution treatment. In addition, RNA-seq technology was used to study the differential expression of mRNA and miRNA. The results showed that triploid Populus grew well and the K⁺ content and the K⁺/Na⁺ ratio in the salt treatment were significantly lower than those in the control. The contents of ABA, JA, POD, and proline were increased compared with contents in diploid under salt stress. The salt-stress responsive DEGs were mainly involved in ion transport, cell homeostasis, the MAPK signaling pathway, peroxisome, citric acid cycle, and other salt response and growth pathways. The transcription factors mainly included NAC, MYB, MYB_related and AP2/ERF. Moreover, the differentially expressed miRNAs involved 32 families, including 743 miRNAs related to predicted target genes, among which 22 miRNAs were significantly correlated with salt-stress response genes and related to the regulation of hormones, ion transport, reactive oxygen species (ROS) and other biological processes. Our results provided insights into the physiological and molecular aspects for further research into the response mechanisms of allotriploid Populus cathayana to salt stress. This study provided valuable information for the salt tolerance mechanism of allopolyploids.

Plastome-based phylogenomics elucidate relationships in rare Isoëtes species groups from the Neotropics

The genus Isoëtes is globally distributed. Within the Neotropics, Isoëtes occurs in various habitats and ecosystems, making it an interesting case study to address phylogenetic and biogeographic questions. We sequenced and assembled plastomes and ribosomal DNA (rDNA) sequences to reconstruct phylogenetic relationships in Isoëtes from tropical regions in the Neotropics. The ploidy level of nine taxa was established to address the potential source of phylogenetic incongruence in the genus. Node ages were estimated using MCMCTree. The ancestral range estimates were conducted in BioGeoBEARS. Plastome-based phylogenies were congruent throughout distinct matrices and partition schemes, exhibiting high support for almost all nodes. Whereas, we found incongruences between the rDNA and plastome datasets. Chromosome counts identified three diploids, five tetraploids and one likely hexaploid among Neotropical species. Plastome-based node age estimates showed that the radiation of the crown Isoëtes group occurred at 20 Ma, with the diversification of the tropical American (TAA) clade taking place in the Pleistocene at 1.7 Ma. Ancestral range estimates showed that the ancestor of the TAA clade may have evolved first in the dry diagonal area in South America before reaching more humid regions. In addition, the colonization of the Brazilian semiarid region occurred three times, while the occupation of the Cerrado and Amazon regions occurred twice and once, respectively. Our study showed a large unobserved diversity within the genus in warm-dry regions in the Neotropics. Plastomes provided sufficient genomic information to establish a robust phylogenetic framework to answer evolutionary questions in Isoëtes from the Neotropics.

ChIP-cloning analysis uncovers centromere-specific retrotransposons in Brassica nigra and reveals their rapid diversification in Brassica allotetraploids

Centromeres are indispensable functional units of chromosomes. The evolutionary mechanisms underlying the rapid evolution of centromeric repeats, especially those following polyploidy, remain unknown. In this study, we isolated centromeric sequences of Brassica nigra, a model diploid progenitor (B genome) of the allopolyploid species B. juncea (AB genome) and B. carinata (BC genome) by chromatin immunoprecipitation of nucleosomes containing the centromere-specific histone CENH3. Sequence analysis detected no centromeric satellite DNAs, and most B. nigra centromeric repeats were found to originate from Tyl/copia-class retrotransposons. In cytological analyses, six of the seven analyzed repeat clusters had no FISH signals in A or C genomes of the related diploid species B. rapa and B. oleracea. Notably, five repeat clusters had FISH signals in both A and B subgenomes in the tetraploid B. juncea. In the tetraploid B. carinata, only CL23 displayed three pairs of signals in terminal or interstitial regions of the C-derived chromosome, and no evidence of colonization of CLs onto C-subgenome centromeres was found in B. carinata. This observation suggests that centromeric repeats spread and proliferated between genomes after polyploidization. CL3 and CRB are likely ancient centromeric sequences arising prior to the divergence of diploid Brassica which have detected signals across the genus. And in allotetraploids B. juncea and B. carinata, the FISH signal intensity of CL3 and CRB differed among subgenomes. We discussed possible mechanisms for centromeric repeat divergence during Brassica speciation and polyploid evolution, thus providing insights into centromeric repeat establishment and targeting.

FitTetra 2.0 - improved genotype calling for tetraploids with multiple population and parental data support

BACKGROUND

Genetic studies in tetraploids are lagging behind in comparison with studies of diploids as the complex genetics of tetraploids require much more elaborated computational methodologies. Recent advancements in development of molecular techniques and computational tools facilitate new methods for automated, high-throughput genotype calling in tetraploid species. We report on the upgrade of the widely-used fitTetra software aiming to improve its accuracy, which to date is hampered by technical artefacts in the data.

RESULTS

Our upgrade of the fitTetra package is designed for a more accurate modelling of complex collections of samples. The package fits a mixture model where some parameters of the model are estimated separately for each sub-collection. When a full-sib family is analyzed, we use parental genotypes to predict the expected segregation in terms of allele dosages in the offspring. More accurate modelling and use of parental data increases the accuracy of dosage calling. We tested the package on data obtained with an Affymetrix Axiom 60 k array and compared its performance with the original version and the recently published ClusterCall tool, showing that at least 20% more SNPs could be called with our updated.

CONCLUSION

Our updated software package shows clearly improved performance in genotype calling accuracy. Estimation of mixing proportions of the underlying dosage distributions is separated for full-sib families (where mixture proportions can be estimated from the parental dosages and inheritance model) and unstructured populations (where they are based on the assumption of Hardy-Weinberg equilibrium). Additionally, as the distributions of signal ratios of the dosage classes can be assumed to be the same for all populations, including parental data for some subpopulations helps to improve fitting other populations as well. The R package fitTetra 2.0 is freely available under the GNU Public License as Additional file with this article.

pSBVB: A Versatile Simulation Tool To Evaluate Genomic Selection in Polyploid Species

Genomic Selection (GS) is the procedure whereby molecular information is used to predict complex phenotypes and it is standard in many animal and plant breeding schemes. However, only a small number of studies have been reported in horticultural crops, and in polyploid species in particular. In this paper, we have developed a versatile forward simulation tool, called polyploid Sequence Based Virtual Breeding (pSBVB), to evaluate GS strategies in polyploids; pSBVB is an efficient gene dropping software that can simulate any number of complex phenotypes, allowing a very flexible modeling of phenotypes suited to polyploids. As input, it takes genotype data from the founder population, which can vary from single nucleotide polymorphisms (SNP) chips up to sequence, a list of causal variants for every trait and their heritabilities, and the pedigree. Recombination rates between homeologous chromosomes can be specified, so that both allo- and autopolyploid species can be considered. The program outputs phenotype and genotype data for all individuals in the pedigree. Optionally, it can produce several genomic relationship matrices that consider exact or approximate genotype values. pSBVB can therefore be used to evaluate GS strategies in polyploid species (say varying SNP density, genetic architecture or population size, among other factors), or to optimize experimental designs for association studies. We illustrate pSBVB with SNP data from tetraploid potato and partial sequence data from octoploid strawberry, and we show that GS is a promising breeding strategy for polyploid species but that the actual advantage critically depends on the underlying genetic architecture. Source code, examples and a complete manual are freely available in GitHub https://github.com/lauzingaretti/pSBVB.

Meiotic abnormalities affect genetic constitution and pollen viability in dicots from Indian cold deserts

BACKGROUND

Meiotic abnormalities lead to morphological and genetic variations which caused not only to evolution but also intraspecific reproductive barriers. During present study of detailed meiotic course in dicotyledonous plants sampled from Indian cold deserts, various meiotic abnormalities have been detected. For this, the plant materials fixed in Carnoy's fixative and studied detailed meiotic course by standard squash method in 1% acetocarmine.

RESULTS

Meiotic abnormalities have been presently detected in 71 species which include multiple associations in diploids (Achillea millefolium L.), multivalents and univalents in polyploids (4 species), cytomixis (40 species), chromosome stickiness (20 species), nonsynchronous disjunction of bivalents (32 species), interbivalent connections (15 species), synaptic mutants (2 species), syncyte meiocytes (2 species), abnormal spindles (7 species), and fusion of pollen grains (1 species), laggards and chromatin bridges, hypo-, hyperploid PMCs, monads, dyads, triads, tetrads with micronuclei and polyads.

CONCLUSIONS

Consequently, variable sized apparently fertile pollen grains and considerable amount of sterile pollen grains are resulted as end products which lead to different genetic constitution (aneuploids and polyploids) and curtailed sexual reproductive success in these species.

Chromosome painting in meiosis reveals pairing of specific chromosomes in polyploid Solanum species

Analysis of chromosome pairing has been an important tool to assess the genetic similarity of homologous and homoeologous chromosomes in polyploids. However, it is technically challenging to monitor the pairing of specific chromosomes in polyploid species, especially for plant species with a large number of small chromosomes. We developed oligonucleotide-based painting probes for four different potato chromosomes. We demonstrate that these probes are robust enough to monitor a single chromosome throughout the prophase I of meiosis in polyploid Solanum species. Cultivated potato (Solanum tuberosum, 2n = 4x = 48) is an autotetraploid. We demonstrate that the four copies of each potato chromosome pair as a quadrivalent in 66-78% of the meiotic cells at the pachytene stage. Solanum demissum (2n = 6x = 72) is a hexaploid and has been controversial regarding its nature as an autopolyploid or allopolyploid. Interestingly, no hexavalent pairing was observed in meiosis. Instead, we observed three independent bivalents in 83-98% of the meiotic cells at late diakinesis and early metaphase I for the four chromosomes. These results suggest that S. demissum has evolved into a cytologically stable state with predominantly bivalent pairing in meiosis.

Tissue-specific expression of telomerase reverse transcriptase gene variants in Nicotiana tabacum

MAIN CONCLUSION

In tobacco, three sequence variants of the TERT gene have been described. We revealed unbalanced levels of TERT variant transcripts in vegetative tobacco tissues and enhanced TERT transcription and telomerase activity in reproductive tissues. Telomerase is a ribonucleoprotein complex responsible for the maintenance of telomeres, structures delimiting ends of linear eukaryotic chromosomes. In the Nicotiana tabacum (tobacco) allotetraploid plant, three sequence variants (paralogs) of the gene coding for the telomerase reverse transcriptase subunit (TERT) have been described, two of them derived from the maternal N. sylvestris genome (TERT_Cs, TERT_D) and one originated from the N. tomentosiformis paternal genome (TERT_Ct). In this work, we analyzed the transcription of TERT variants in correlation with telomerase activity in tobacco tissues. High and approximately comparable levels of TERT_Ct and TERT_Cs transcripts were detected in seedlings, roots, flower buds and leaves, while the transcript of the TERT_D variant was markedly underrepresented. Similarly, in N. sylvestris tissues, TERT_Cs transcript significantly predominated. A specific pattern of TERT transcripts was found in samples of tobacco pollen with the TERT_Cs variant clearly dominating particularly at the early stage of pollen development. Detailed analysis of TERT_C variants representation in functionally distinct fractions of pollen transcriptome revealed their prevalence in large ribonucleoprotein particles encompassing translationally silent mRNA; only a minority of TERT_Ct and TERT_Cs transcripts were localized in actively translated polysomes. Histones of the TERT_C chromatin were decorated predominantly with the euchromatin-specific epigenetic modification in both telomerase-positive and telomerase-negative tobacco tissues. We conclude that the existence and transcription pattern of tobacco TERT paralogs represents an interesting phenomenon and our results indicate its functional significance. Nicotiana species have again proved to be appropriate and useful model plants in telomere biology studies.

PERGOLA: fast and deterministic linkage mapping of polyploids

BACKGROUND

A large share of agriculturally and horticulturally important plant species are polyploid. Linkage maps are used to locate associations between genes and traits by breeders and geneticists. Linkage map creation for polyploid species is not supported by standard tools. We want to overcome this limitation and validate our results with simulation studies.

RESULTS

We developed PERGOLA, a deterministic and heuristic method that addresses this problem. We show that it creates correct linkage groups, marker orders and distances for simulated and real datasets. We compare it to existing tools and demonstrate that it overcomes limitations in ploidy and outperforms them in computational time and mapping accuracy. We represent linkage maps as dendrograms and show that this has advantages in the comparison of different maps.

CONCLUSIONS

PERGOLA can be used successfully to calculate linkage maps for diploid and polyploid species and outperforms existing tools.

A rise of ploidy level influences the rate of cytomixis in tobacco male meiosis

The effect of plant ploidy level on the rate of cytomixis in microsporogenesis has been analyzed with the help of a unique model, the collection of tobacco plants of different ploidies (2n = 2x = 24, 4x = 48, 6x = 72, and 8x = 96). As has been shown, the rate of cytomixis proportionally increases in 6x and 8x cytotypes, being rather similar in 2x and 4x plants. The rate of cytomixis is highly variable, differing even in the genetically identical plants grown under the same conditions. The cytological pattern of cytomixis in the microsporogenesis of control 4x plants has been compared with the corresponding patterns of 2x, 6x, and 8x plants. Involvement of cytomixis in production of unreduced gametes and stabilization of the newly formed hybrid and polyploidy genomes is discussed.

Advantages of continuous genotype values over genotype classes for GWAS in higher polyploids: a comparative study in hexaploid chrysanthemum

Background

Association studies are an essential part of modern plant breeding, but are limited for polyploid crops. The increased number of possible genotype classes complicates the differentiation between them. Available methods are limited with respect to the ploidy level or data producing technologies. While genotype classification is an established noise reduction step in diploids, it gains complexity with increasing ploidy levels. Eventually, the errors produced by misclassifications exceed the benefits of genotype classes. Alternatively, continuous genotype values can be used for association analysis in higher polyploids. We associated continuous genotypes to three different traits and compared the results to the output of the genotype caller SuperMASSA. Linear, Bayesian and partial least squares regression were applied, to determine if the use of continuous genotypes is limited to a specific method. A disease, a flowering and a growth trait with h² of 0.51, 0.78 and 0.91 were associated with a hexaploid chrysanthemum genotypes. The data set consisted of 55,825 probes and 228 samples.

Results

We were able to detect associating probes using continuous genotypes for multiple traits, using different regression methods. The identified probe sets were overlapping, but not identical between the methods. Baysian regression was the most restrictive method, resulting in ten probes for one trait and none for the others. Linear and partial least squares regression led to numerous associating probes. Association based on genotype classes resulted in similar values, but missed several significant probes. A simulation study was used to successfully validate the number of associating markers.

Conclusions

Association of various phenotypic traits with continuous genotypes is successful with both uni- and multivariate regression methods. Genotype calling does not improve the association and shows no advantages in this study. Instead, use of continuous genotypes simplifies the analysis, saves computational time and results more potential markers.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2926-5) contains supplementary material, which is available to authorized users.

Effects of genome structure variation, homeologous genes and repetitive DNA on polyploid crop research in the age of genomics

Compared to diploid species, allopolyploid crop species possess more complex genomes, higher productivity, and greater adaptability to changing environments. Next generation sequencing techniques have produced high-density genetic maps, whole genome sequences, transcriptomes and epigenomes for important polyploid crops. However, several problems interfere with the full application of next generation sequencing techniques to these crops. Firstly, different types of genomic variation affect sequence assembly and QTL mapping. Secondly, duplicated or homoeologous genes can diverge in function and then lead to emergence of many minor QTL, which increases difficulties in fine mapping, cloning and marker assisted selection. Thirdly, repetitive DNA sequences arising in polyploid crop genomes also impact sequence assembly, and are increasingly being shown to produce small RNAs to regulate gene expression and hence phenotypic traits. We propose that these three key features should be considered together when analyzing polyploid crop genomes. It is apparent that dissection of genomic structural variation, elucidation of the function and mechanism of interaction of homoeologous genes, and investigation of the de novo roles of repeat sequences in agronomic traits are necessary for genomics-based crop breeding in polyploids.