RESUMEN

Polyploidy is a significant mechanism in eukaryotic evolution and is particularly prevalent in the plant kingdom. However, our knowledge about this phenomenon and its effects on evolution remains limited. A major obstacle to the study of polyploidy is the great difficulty in untangling the origins of allopolyploids. Due to the drastic genome changes and the erosion of allopolyploidy signals caused by the combined effects of hybridization and complex post-polyploid diploidization processes, resolving the origins of allopolyploids has long been a challenging task. Here we revisit this issue with the interesting case of subtribe Tussilagininae (Asteraceae: Senecioneae) and by developing HomeoSorter, a new pipeline for network inferences by phasing homeologs to parental subgenomes. The pipeline is based on the basic idea of a previous study but with major changes to address the scaling problem and implement some new functions. With simulated data, we demonstrate that HomeoSorter works efficiently on genome-scale data and has high accuracy in identifying polyploid patterns and assigning homeologs. Using HomeoSorter, the maximum pseudo-likelihood model of Phylonet, and genome-scale data, we further address the complex origin of Tussilagininae, a speciose group (ca. 45 genera and 710 species) characterized by having high base chromosome numbers (mainly x = 30, 40). In particular, the inferred patterns are strongly supported by the chromosomal evidence. Tussilagininae is revealed to comprise two large groups with successive allopolyploid origins: Tussilagininae s.s. (mainly x = 30) and the Gynoxyoid group (x = 40). Two allopolyploidy events first give rise to Tussilagininae s.s., with the first event occurring between the ancestor of subtribe Senecioninae (x = 10) and a lineage (highly probably with x = 10) related to the Brachyglottis alliance, and the resulting hybrid lineage crossing with the ancestor of Chersodoma (x = 10) and leading to Tussilagininae s.s. Then, after early diversification, the Central American group (mainly x = 30) of Tussilagininae s.s., is involved in a third allopolyploidy event with, again, the Chersodoma lineage and produces the Gynoxyoid group. Our study highlights the value of HomeoSorter and the homeolog-sorting approach in polyploid phylogenetics. With rich species diversity and clear evolutionary patterns, Tussilagininae s.s. and the Gynoxyoid group are also excellent models for future investigations of polyploidy.

RESUMEN

Camelina sativa (L.) Crantz, a member of the Brassicaceae, has potential as a biofuel feedstock which is attributable to the production of fatty acids in its seeds, its fast growth cycle, and low input requirements. While a genome assembly is available for camelina, it was generated from short sequence reads and is thus highly fragmented in nature. Using long read sequences, we generated a chromosome-scale, highly contiguous genome assembly (644,491,969 bp) for the spring biotype cultivar 'Suneson' with an N50 contig length of 12,031,512 bp and a scaffold N50 length of 32,184,682 bp. Annotation of protein-coding genes revealed 91,877 genes that encode 133,355 gene models. We identified a total of 4,467 genes that were significantly up-regulated under cold stress which were enriched in gene ontology terms associated with "response to cold" and "response to abiotic stress". Coexpression analyses revealed multiple coexpression modules that were enriched in genes differentially expressed following cold stress that had putative functions involved in stress adaptation, specifically within the plastid. With access to a highly contiguous genome assembly, comparative analyses with Arabidopsis thaliana revealed 23,625 A. thaliana genes syntenic with 45,453 Suneson genes. Of these, 24,960 Suneson genes were syntenic to 8,320 A. thaliana genes reflecting a 3 camelina homeolog to 1 Arabidopsis gene relationship and retention of all three homeologs. Some of the retained triplicated homeologs showed conserved gene expression patterns under control and cold-stressed conditions whereas other triplicated homeologs displayed diverged expression patterns revealing sub- and neo-functionalization of the homeologs at the transcription level. Access to the chromosome-scale assembly of Suneson will enable both basic and applied research efforts in the improvement of camelina as a sustainable biofuel feedstock.

RESUMEN

Peanut (Arachis) is a key oil and protein crop worldwide with large genome. The genomes of diploid and tetraploid peanuts have been sequenced, which were compared to decipher their genome structures, evolutionary, and life secrets. Genome sequencing efforts showed that different cultivars, although Bt homeologs being more privileged in gene retention and gene expression. This subgenome bias, extended to sequence variation and point mutation, might be related to the long terminal repeat (LTR) explosions after tetraploidization, especially in At subgenomes. Except that, whole-genome sequences revealed many important genes, for example, fatty acids and triacylglycerols pathway, NBS-LRR (nucleotide-binding site-leucine-rich repeats), and seed size decision genes, were enriched after recursive polyploidization. Each ancestral polyploidy, with old ones having occurred hundreds of thousand years ago, has thousands of duplicated genes in extant genomes, contributing to genetic novelty. Notably, although full genome sequences are available, the actual At subgenome ancestor has still been elusive, highlighted with new debate about peanut origin. Although being an orphan crop lagging behind other crops in genomic resources, the genome sequencing achievement has laid a solid foundation for advancing crop enhancement and system biology research of peanut.

Asunto(s)

RESUMEN

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.

Asunto(s)

RESUMEN

Research on the evolutionary fate of duplicated genes in recurrent polyploids is scarce due to the difficulties in disentangling the different homeologs and alleles of duplicated genes. This chapter describes the detailed procedures to identify different homeologs and alleles of duplicated genes, to analyze their molecular characteristics, and to reveal their functional divergence by gene editing with CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated system 9). Using the gene editing approach, we efficiently constructed multiple knockout mutant lines with single or simultaneously disrupted different homeologs or alleles in a recurrent polyploid fish, demonstrating its usability for targeting and mutating multiple divergent homeologs and alleles in recurrent duplicated genomes.

Asunto(s)

RESUMEN

The role of hybridization is significant in biological invasion, and thermotolerance is a trait critical to range expansions. The South American Sphagneticola trilobata is now widespread in South China, threatening the native S. calendulacea by competition and hybridization. Furthermore, upon formation, their F1 hybrid can quickly replace both parents. In this study, the three taxa were used as a model to investigate the consequences of hybridization on cold tolerance, particularly the effect of subgenome dominance in the hybrid. Upon chilling treatments, physiological responses and transcriptome profiles were compared across different temperature points to understand their differential responses to cold. While both parents showed divergent responses, the hybrid's responses showed an overall resemblance to S. calendulacea, but the contribution of homeolog expression bias to cold stress was not readily evident in the F1 hybrid possibly due to inherent bias that comes with the sampling location. Our findings provided insights into the role of gene expression in differential cold tolerance, and further contribute to predicting the invasive potential of other hybrids between S. trilobata and its congeners around the world.

RESUMEN

WRKY transcription factors play critical roles in plant growth and development or stress responses. Using up-to-date genomic data, a total of 64 and 257 WRKY genes have been identified in the diploid woodland strawberry, Fragaria vesca, and the more complex allo-octoploid commercial strawberry, Fragaria × ananassa cv. Camarosa, respectively. The completeness of the new genomes and annotations has enabled us to perform a more detailed evolutionary and functional study of the strawberry WRKY family members, particularly in the case of the cultivated hybrid, in which homoeologous and paralogous FaWRKY genes have been characterized. Analysis of the available expression profiles has revealed that many strawberry WRKY genes show preferential or tissue-specific expression. Furthermore, significant differential expression of several FaWRKY genes has been clearly detected in fruit receptacles and achenes during the ripening process and pathogen challenged, supporting a precise functional role of these strawberry genes in such processes. Further, an extensive analysis of predicted development, stress and hormone-responsive cis-acting elements in the strawberry WRKY family is shown. Our results provide a deeper and more comprehensive knowledge of the WRKY gene family in strawberry.

RESUMEN

BACKGROUND: Cultivated peanut (Arachis hypogaea, AABB genome), an allotetraploid from a cross between A. duranensis (AA genome) and A. ipaensis (BB genome), is an important oil and protein crop with released genome and RNA-seq sequence datasets. These datasets provide the molecular foundation for studying gene expression and evolutionary patterns. However, there are no reports on the proteomic data of A. hypogaea cv. Tifrunner, which limits understanding of its gene function and protein level evolution. RESULTS: This study sequenced the A. hypogaea cv. Tifrunner leaf and root proteome using the tandem mass tag technology. A total of 4803 abundant proteins were identified. The 364 differentially abundant proteins were estimated by comparing protein abundances between leaf and root proteomes. The differentially abundant proteins enriched the photosystem process. The number of biased abundant homeologs between the two sub-genomes A (87 homeologs in leaf and root) and B (69 and 68 homeologs in leaf and root, respectively) was not significantly different. However, homeologous proteins with biased abundances in different sub-genomes enriched different biological processes. In the leaf, homeologs biased to sub-genome A enriched biosynthetic and metabolic process, while homeologs biased to sub-genome B enriched iron ion homeostasis process. In the root, homeologs with biased abundance in sub-genome A enriched inorganic biosynthesis and metabolism process, while homeologs with biased abundance in sub-genome B enriched organic biosynthesis and metabolism process. Purifying selection mainly acted on paralogs and homeologs. The selective pressure values were negatively correlated with paralogous protein abundance. About 77.42% (24/31) homeologous and 80% (48/60) paralogous protein pairs had asymmetric abundance, and several protein pairs had conserved abundances in the leaf and root tissues. CONCLUSIONS: This study sequenced the proteome of A. hypogaea cv. Tifrunner using the leaf and root tissues. Differentially abundant proteins were identified, and revealed functions. Paralog abundance divergence and homeolog bias abundance was elucidated. These results indicate that divergent abundance caused retention of homologs in A. hypogaea cv. Tifrunner.

RESUMEN

Proton-pumping pyrophosphatases (H+-PPases) have been shown to enhance biomass and yield. However, to date, there has been little work towards identify genes encoding H+-PPases in bread wheat (Triticum aestivum) (TaVPs) and limited knowledge on how the expression of these genes varies across different growth stages and tissue types. In this study, the IWGSC database was used to identify two novel TaVP genes, TaVP4 and TaVP5, and elucidate the complete homeolog sequences of the three known TaVP genes, bringing the total number of bread wheat TaVPs from 9 to 15. Gene expression levels of each TaVP homeolog were assessed using quantitative real-time PCR (qRT-PCR) in four diverse wheat varieties in terms of phenotypic traits related to high vacuolar pyrophosphatase expression. Homeolog expression was analyzed across multiple tissue types and developmental stages. Expression levels of the TaVP homeologs were found to vary significantly between varieties, tissues and plant developmental stages. During early development (Z10 and Z13), expressions of TaVP1 and TaVP2 homeologs were higher in shoot tissue than root tissue, with both shoot and root expression increasing in later developmental stages (Z22). TaVP2-D was expressed in all varieties and tissue types and was the most highly expressed homeolog at all developmental stages. Expression of the TaVP3 homeologs was restricted to developing grain (Z75), while TaVP4 homeolog expression was higher at Z22 than earlier developmental stages. Variation in TaVP4B was detected among varieties at Z22 and Z75, with Buck Atlantico (high biomass) and Scout (elite Australian cultivar) having the highest levels of expression. These findings offer a comprehensive overview of the bread wheat H+-PPase family and identify variation in TaVP homeolog expression that will be of use to improve the growth, yield, and abiotic stress tolerance of bread wheat.

RESUMEN

Polyploidy and subsequent diploidization provide genomic opportunities for evolutionary innovations and adaptation. The researches on duplicated gene evolutionary fates in recurrent polyploids have seriously lagged behind that in paleopolyploids with diploidized genomes. Moreover, the antiviral mechanisms of Viperin remain largely unclear in fish. Here, we elaborate the distinct antiviral mechanisms of two viperin homeologs (Cgviperin-A and Cgviperin-B) in auto-allo-hexaploid gibel carp (Carassius gibelio). First, Cgviperin-A and Cgviperin-B showed differential and biased expression patterns in gibel carp adult tissues. Subsequently, using co-immunoprecipitation (Co-IP) screening analysis, both CgViperin-A and CgViperin-B were found to interact with crucian carp (C. auratus) herpesvirus (CaHV) open reading frame 46 right (ORF46R) protein, a negative herpesvirus regulator of host interferon (IFN) production, and to promote the proteasomal degradation of ORF46R via decreasing K63-linked ubiquitination. Additionally, CgViperin-B also mediated ORF46R degradation through autophagosome pathway, which was absent in CgViperin-A. Moreover, we found that the N-terminal α-helix domain was necessary for the localization of CgViperin-A and CgViperin-B at the endoplasmic reticulum (ER), and the C-terminal domain of CgViperin-A and CgViperin-B was indispensable for the interaction with degradation of ORF46R. Therefore, the current findings clarify the divergent antiviral mechanisms of the duplicated viperin homeologs in a recurrent polyploid fish, which will shed light on the evolution of teleost duplicated genes.

Asunto(s)

RESUMEN

BACKGROUND: Polyploid species often originate recurrently. While this is well known, there is little information on the extent to which distinct allotetraploid species formed from the same parent species differ in gene expression. The tetraploid yarrow species Achillea alpina and A. wilsoniana arose independently from allopolyploidization between diploid A. acuminata and A. asiatica. The genetics and geography of these origins are clear from previous studies, providing a solid basis for comparing gene expression patterns of sibling allopolyploid species that arose independently. RESULTS: We conducted comparative RNA-sequencing analyses on the two Achillea tetraploid species and their diploid progenitors to evaluate: 1) species-specific gene expression and coexpression across the four species; 2) patterns of inheritance of parental gene expression; 3) parental contributions to gene expression in the allotetraploid species, and homeolog expression bias. Diploid A. asiatica showed a higher contribution than diploid A. acuminata to the transcriptomes of both tetraploids and also greater homeolog bias in these transcriptomes, possibly reflecting a maternal effect. Comparing expressed genes in the two allotetraploids, we found expression of ca. 30% genes were species-specific in each, which were most enriched for GO terms pertaining to "defense response". Despite species-specific and differentially expressed genes between the two allotetraploids, they display similar transcriptome changes in comparison to their diploid progenitors. CONCLUSION: Two independently originated Achillea allotetraploid species exhibited difference in gene expression, some of which must be related to differential adaptation during their post-speciation evolution. On the other hand, they showed similar expression profiles when compared to their progenitors. This similarity might be expected when pairs of merged diploid genomes in tetraploids are similar, as is the case in these two particular allotetraploids.

Asunto(s)

RESUMEN

Chemokines, relatively small secreted proteins, are involved in cell migration and function in various biological events, including immunity, morphogenesis, and disease. Due to their nature, chemokines tend to be a target of hijacking of immunity by virus and therefore show an exceptionally high mutation rate. Xenopus laevis is considered an excellent model to investigate the effect of whole-genome duplication for gene family evolution. Because its allotetraploidization occurred around 17-18 million years ago, ancestral subgenomes L and S were well conserved. Based on the gene model of human and diploid frog Xenopus tropicalis, we identified 52 chemokine genes and 26 chemokine receptors in X. laevis. The retention rate of the gene in the X. laevis L and S subgenomes was 96% (45/47) and 68% (32/47), respectively. We conducted molecular phylogenetic analysis and found clear orthologies in all receptor genes but not in the ligand genes, suggesting rapid divergences of the ligand. dN/dS calculation demonstrated that dN/dS ratio greater than one was observed in the four ligand genes, cxcl8b.1.S, cxcl18.S, ccl21.S, and xcl1.L, but nothing in receptor genes. These results revealed that the whole-genome duplication promotes diversification of chemokine ligands in X. laevis while conserving the genes necessary for homeostasis, suggesting that selective pressure also supports a rapid divergence of the chemokines in amphibians.

RESUMEN

Polyploidization is pervasive in plants, but little is known about the niche divergence of wild allopolyploids (species that harbor polyploid genomes originating from different diploid species) relative to their diploid progenitor species and the gene expression patterns that may underlie such ecological divergence. We conducted a fine-scale empirical study on habitat and gene expression of an allopolyploid and its diploid progenitors. We quantified soil properties and light availability of habitats of an allotetraploid Cardamine flexuosa and its diploid progenitors Cardamine amara and Cardamine hirsuta in two seasons. We analyzed expression patterns of genes and homeologs (homeologous gene copies in allopolyploids) using RNA sequencing. We detected niche divergence between the allopolyploid and its diploid progenitors along water availability gradient at a fine scale: the diploids in opposite extremes and the allopolyploid in a broader range between diploids, with limited overlap with diploids at both ends. Most of the genes whose homeolog expression ratio changed among habitats in C. flexuosa varied spatially and temporally. These findings provide empirical evidence for niche divergence between an allopolyploid and its diploid progenitor species at a fine scale and suggest that divergent expression patterns of homeologs in an allopolyploid may underlie its persistence in diverse habitats.

Asunto(s)

RESUMEN

Contemporary speciation provides a unique opportunity to directly observe the traits and environmental responses of a new species. Cardamine insueta is an allotriploid species that appeared within the past 150 years in a Swiss village, Urnerboden. In contrast to its two progenitor species, Cardamine amara and Cardamine rivularis that live in wet and open habitats, respectively, C. insueta is found in-between their habitats with temporal water level fluctuation. This triploid species propagates clonally and serves as a triploid bridge to form higher ploidy species. Although niche separation is observed in field studies, the mechanisms underlying the environmental robustness of C. insueta are not clear. To characterize responses to a fluctuating environment, we performed a time-course analysis of homeolog gene expression in C. insueta in response to submergence treatment. For this purpose, the two parental (C. amara and C. rivularis) genome sequences were assembled with a reference-guided approach, and homeolog-specific gene expression was quantified using HomeoRoq software. We found that C. insueta and C. rivularis initiated vegetative propagation by forming ectopic meristems on leaves, while C. amara did not. We examined homeolog-specific gene expression of three species at nine time points during the treatment. The genome-wide expression ratio of homeolog pairs was 2:1 over the time-course, consistent with the ploidy number. By searching the genes with high coefficient of variation of expression over time-course transcriptome data, we found many known key transcriptional factors related to meristem development and formation upregulated in both C. rivularis and rivularis-homeolog of C. insueta, but not in C. amara. Moreover, some amara-homeologs of these genes were also upregulated in the triploid, suggesting trans-regulation. In turn, Gene Ontology analysis suggested that the expression pattern of submergence tolerant genes in the triploid was inherited from C. amara. These results suggest that the triploid C. insueta combined advantageous patterns of parental transcriptomes to contribute to its establishment in a new niche along a water-usage gradient.

RESUMEN

Empirical evidence is limited on whether allopolyploid species combine or merge parental adaptations to broaden habitats. The allopolyploid Arabidopsis kamchatica is a hybrid of the two diploid parents Arabidopsis halleri and Arabidopsis lyrata. A. halleri is a facultative heavy metal hyperaccumulator, and may be found in cadmium (Cd) and zinc (Zn) contaminated environments, as well as non-contaminated environments. A. lyrata is considered non-tolerant to these metals, but can be found in serpentine habitats. Therefore, the parents have adaptation to different environments. Here, we measured heavy metals in soils from native populations of A. kamchatica. We found that soil Zn concentration of nearly half of the sampled 40 sites was higher than the critical toxicity level. Many of the sites were near human construction, suggesting adaptation of A. kamchatica to artificially contaminated soils. Over half of the A. kamchatica populations had >1,000 µg g-1 Zn in leaf tissues. Using hydroponic treatments, most genotypes accumulated >3,000 µg g-1 Zn, with high variability among them, indicating substantial genetic variation in heavy metal accumulation. Genes involved in heavy metal hyperaccumulation showed an expression bias in the A. halleri-derived homeolog in widely distributed plant genotypes. We also found that two populations were found growing on serpentine soils. These data suggest that A. kamchatica can inhabit a range of both natural and artificial soil environments with high levels of ions that either of the parents specializes and that it can accumulate varying amount of heavy metals. Our field and experimental data provide a compelling example of combining genetic toolkits for soil adaptations to expand the habitat of an allopolyploid species.

RESUMEN

Polyploidy is an important evolutionary mechanism and is prevalent among land plants. Most polyploid species examined have multiple origins, which provide genetic diversity and may enhance the success of polyploids. In some polyploids, recurrent origins can result from reciprocal crosses between the same diploid progenitors. Although great progress has been made in understanding the genetic consequences of polyploidy, the genetic implications of reciprocal polyploidization remain poorly understood, especially in natural polyploids. Tragopogon (Asteraceae) has become an evolutionary model system for studies of recent and recurrent polyploidy. Allotetraploid T. miscellus has formed reciprocally in nature with resultant distinctive floral and inflorescence morphologies (i.e., short- vs. long-liguled forms). In this study, we performed comparative inflorescence transcriptome analyses of reciprocally formed T. miscellus and its diploid parents, T. dubius and T. pratensis. In both forms of T. miscellus, homeolog expression of ∼70% of the loci showed vertical transmission of the parental expression patterns (i.e., parental legacy), and ∼20% of the loci showed biased homeolog expression, which was unbalanced toward T. pratensis. However, 17.9% of orthologous pairs showed different homeolog expression patterns between the two forms of T. miscellus. No clear effect of cytonuclear interaction on biased expression of the maternal homeolog was found. In terms of the total expression level of the homeologs studied, 22.6% and 16.2% of the loci displayed non-additive expression in short- and long-liguled T. miscellus, respectively. Unbalanced expression level dominance toward T. pratensis was observed in both forms of T. miscellus. Significantly, genes annotated as being involved in pectin catabolic processes were highly expressed in long-liguled T. miscellus relative to the short-liguled form, and the majority of these differentially expressed genes were transgressively down-regulated in short-liguled T. miscellus. Given the known role of these genes in cell expansion, they may play a role in the differing floral and inflorescence morphologies of the two forms. In summary, the overall inflorescence transcriptome profiles are highly similar between reciprocal origins of T. miscellus. However, the dynamic homeolog-specific expression and non-additive expression patterns observed in T. miscellus emphasize the importance of reciprocal origins in promoting the genetic diversity of polyploids.

RESUMEN

Most organisms are diploid, meaning they only have two copies of each chromosome (one set inherited from each parent). Polyploid organisms have more than two paired (homologous) sets of chromosomes. Many plant species are polyploid. Polyploid species cope better with stresses thanks to the redundancy in the chromosome copy number and dispose in this way a greater flexibility in gene expression. Allopolyploid species are polyploids that contain an alternative set of chromosomes by the cross of two (or more) species. Gene variants unique for a preferential phenotype are most probable candidate markers controlling the observed phenotype. Organ or tissue-specific silencing or overexpression of one parental homeolog is quite common. It is very challenging to find those tissue-specific gene variants. High-throughput proteomics is a successful method to discover them. This chapter proposes two possible workflows depending on the available resources and the knowledge of the species. An example is given for an AAB hybrid and an ABB hybrid. Allele-specific gene responses are picked up in this workflow as gene loci displaying genotype-specific differential expression that often have single amino acid polymorphisms. If the resources are sufficient, a genotype-specific mRNAseq database is recommended where a link is made to the allele-specific transcription levels. If the resources are limited, allele-specific proteins can be detected by the detection of genotype-specific peptides and the identification against existing genomics libraries of the parents.

Asunto(s)

RESUMEN

Evolution of Brassica genome post-polyploidization reveals asymmetrical genome fractionation and copy number variation. Herein, we describe the impact of promoter divergence among SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1 (SOC1) homeologs on expression and function in Brassica spp. SOC1, a regulated floral pathway integrator, is conserved as 3 redundant homeologs in diploid Brassicas. Even with high sequence identity within coding regions (92.8-100%), the spatio-temporal expression patterns of 9 SOC1 homologs in B. juncea and B. nigra indicates regulatory divergence. While LF and MF2 SOC1 homeologs are upregulated during floral transition, MF1 is barely expressed. Also, MF2 homeolog levels do not decline post-flowering, unlike LF. To investigate the underlying source of divergence, we analyzed the sequence and phylogeny of all reported (22) and isolated (21) upstream regions of Brassica SOC1. Full length upstream regions (4712-19189 bp) reveal 5 ubiquitously conserved ancestral Blocks, harboring binding sites of 18 TFs (TFBSs) characterized in Arabidopsis thaliana. The orthologs of these TFBSs are differentially conserved among Brassica SOC1 homeologs, imparting expression divergence. No crucial TFBSs are exclusively lost from LF_SOC1 promoter, while MF1_SOC1 has lost NF-Y binding site crucial for SOC1 activation by CONSTANS. MF2_SOC1 homeologs have lost important TFBSs (SEP3, AP1 and SMZ), responsible for SOC1 repression post-flowering. BjuAALF_SOC1 promoter (proximal 2 kb) shows ubiquitous reporter expression in B. juncea cv. Varuna transgenics, while BjuAAMF1_SOC1 promoter shows absence of reporter expression, validating the impact of TFBS divergence. Conservation of the original primary protein sequence is discovered in B. rapa homeologs (46) of 18 TFs. Co-regulation pattern of these TFs appeared similar for B. rapa LF and MF2 SOC1 homeologs; MF1 shows significant variation. Strong regulatory association is recorded for AP1, AP2, SEP3, FLC and CONSTANS/NF-Y, highlighting their importance in homeolog-specific SOC1 regulation. Correlation of B. juncea AP1, AP2 and FLC expression with SOC1 homeologs also complies with the TFBS differences. We thus conclude that redundant SOC1 loci contribute differentially to cumulative expression of SOC1 due to divergent selection of ancestral TFBSs.

Asunto(s)

RESUMEN

Genome duplication with hybridization, or allopolyploidization, occurs in animals, fungi and plants, and is especially common in crop plants. There is an increasing interest in the study of allopolyploids because of advances in polyploid genome assembly; however, the high level of sequence similarity in duplicated gene copies (homeologs) poses many challenges. Here we compared standard RNA-seq expression quantification approaches used currently for diploid species against subgenome-classification approaches which maps reads to each subgenome separately. We examined mapping error using our previous and new RNA-seq data in which a subgenome is experimentally added (synthetic allotetraploid Arabidopsis kamchatica) or reduced (allohexaploid wheat Triticum aestivum versus extracted allotetraploid) as ground truth. The error rates in the two species were very similar. The standard approaches showed higher error rates (>10% using pseudo-alignment with Kallisto) while subgenome-classification approaches showed much lower error rates (<1% using EAGLE-RC, <2% using HomeoRoq). Although downstream analysis may partly mitigate mapping errors, the difference in methods was substantial in hexaploid wheat, where Kallisto appeared to have systematic differences relative to other methods. Only approximately half of the differentially expressed homeologs detected using Kallisto overlapped with those by any other method in wheat. In general, disagreement in low-expression genes was responsible for most of the discordance between methods, which is consistent with known biases in Kallisto. We also observed that there exist uncertainties in genome sequences and annotation which can affect each method differently. Overall, subgenome-classification approaches tend to perform better than standard approaches with EAGLE-RC having the highest precision.

Asunto(s)