Duplicate Gene Expression Research Articles

Translational Regulation of Duplicated Gene Expression Evolution in Allopolyploid Cotton.

Polyploidy, a prevalent event in plant evolution, drives phenotypic diversification and speciation. While transcriptional changes and regulation in polyploids have been extensively studied, the translational level impact remains largely unexplored. To address this gap, we conducted a comparative transcriptomic and translatomic analysis of cotton leaves from allopolyploid species G. hirsutum (AD1) and G. barbadense (AD2) relative to their model A-genome and D-genome diploid progenitors. Our data revealed that while allopolyploidization significantly affects the transcriptional landscape, its impact on translation was relatively modest, evidenced by a narrower expression range and fewer expression changes in ribosome-protected fragments than in mRNA levels. Allopolyploid-specific changes commonly identified in both AD1 and AD2 were observed in 7393 genes at either transcriptional or translational levels. Interestingly, the majority of translational changes exhibited concordant down-regulation in both ribosome-protected fragments and mRNA, particularly associated with terpenoid synthesis and metabolism (352 genes). Regarding translational efficiency (TE), at least one-fifth of cotton genes exhibit translational level regulation, with a general trend of more down-regulation (13.9-15.1%) than up-regulation (7.3-11.2%) of TE. The magnitude of translational regulation was slightly reduced in allopolyploids compared with diploids, and allopolyploidy tends to have a more profound impact on genes and functional associations with ultra-low TE. Moreover, we demonstrated a reduced extent of homeolog expression biases during translation compared with transcription. Our study provides insights into the regulatory consequences of allopolyploidy post-transcription, contributing to a comprehensive understanding of regulatory mechanisms of duplicated gene expression evolution.

Characteristics of duplicated gene expression and DNA methylation regulation in different tissues of allopolyploid Brassica napus

Plant polyploidization increases the complexity of epigenomes and transcriptional regulation, resulting in genome evolution and enhanced adaptability. However, few studies have been conducted on the relationship between gene expression and epigenetic modification in different plant tissues after allopolyploidization. In this study, we studied gene expression and DNA methylation modification patterns in four tissues (stems, leaves, flowers and siliques) of Brassica napusand its diploid progenitors. On this basis, the alternative splicing patterns and cis-trans regulation patterns of four tissues in B. napus and its diploid progenitors were also analyzed. It can be seen that the number of alternative splicing occurs in the B. napus is higher than that in the diploid progenitors, and the IR type increases the most during allopolyploidy. In addition, we studied the fate changes of duplicated genes after allopolyploidization in B. napus. We found that the fate of most duplicated genes is conserved, but the number of neofunctionalization and specialization is also large. The genetic fate of B. napus was classified according to five replication types (WGD, PD, DSD, TD, TRD). This study also analyzed generational transmission analysis of expression and DNA methylation patterns. Our study provides a reference for the fate differentiation of duplicated genes during allopolyploidization.

Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton.

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

Duplicate Gene Expression and Possible Mechanisms of Paralog Retention During Bacterial Genome Expansion.

Gene duplication contributes to the evolution of expression and the origin of new genes, but the relative importance of different patterns of duplicate gene expression and mechanisms of retention remains debated and particularly poorly understood in bacteria. Here, we investigated gene expression patterns for two lab strains of the cyanobacterium Acaryochloris marina with expanding genomes that contain about 10-fold more gene duplicates compared with most bacteria. Strikingly, we observed a generally stoichiometric pattern of greater combined duplicate transcript dosage with increased gene copy number, in contrast to the prevalence of expression reduction reported for many eukaryotes. We conclude that increased transcript dosage is likely an important mechanism of initial duplicate retention in these bacteria and may persist over long periods of evolutionary time. However, we also observed that paralog expression can diverge rapidly, including possible functional partitioning, for which different copies were respectively more highly expressed in at least one condition. Divergence may be promoted by the physical separation of most Acaryochloris duplicates on different genetic elements. In addition, expression pattern for ancestrally shared duplicates could differ between strains, emphasizing that duplicate expression fate need not be deterministic. We further observed evidence for context-dependent transcript dosage, where the aggregate expression of duplicates was either greater or lower than their single-copy homolog depending on physiological state. Finally, we illustrate how these different expression patterns of duplicated genes impact Acaryochloris biology for the innovation of a novel light-harvesting apparatus and for the regulation of recA paralogs in response to environmental change.

Genome-wide DNA methylation dynamics following recent polyploidy in the allotetraploid Tragopogon miscellus (Asteraceae).

Polyploidy is an important evolutionary force, yet epigenetic mechanisms, such as DNA methylation, that regulate genome-wide expression of duplicated genes remain largely unknown. Here, we use Tragopogon (Asteraceae) as a model system to discover patterns and temporal dynamics of DNA methylation in recently formed polyploids. The naturally occurring allotetraploid Tragopogon miscellus formed in the last 95-100 yr from parental diploids Tragopogon dubius and T. pratensis. We profiled the DNA methylomes of these three species using whole-genome bisulfite sequencing. Genome-wide methylation levels in T. miscellus were intermediate between its diploid parents. However, nonadditive CG and CHG methylation occurred in transposable elements (TEs), with variation among TE types. Most differentially methylated regions (DMRs) showed parental legacy, but some novel DMRs were detected in the polyploid. Differentially methylated genes (DMGs) were also identified and characterized. This study provides the first assessment of both overall and locus-specific patterns of DNA methylation in a recent natural allopolyploid and shows that novel methylation variants can be generated rapidly after polyploid formation. Together, these results demonstrate that mechanisms to regulate duplicate gene expression may arise soon after allopolyploid formation and that these mechanisms vary among genes.

Young duplicate genes show developmental stage- and cell type-specific expression and function in Caenorhabditis elegans

Gene duplication produces the material that fuels evolutionary innovation. The "out-of-testis" hypothesis suggests that sperm competition creates selective pressure encouraging the emergence of new genes in male germline, but the somatic expression and function of the newly evolved genes are not well understood. We systematically mapped the expression of young duplicate genes throughout development in Caenorhabditis elegans using both whole-organism and single-cell transcriptomic data. Based on the expression dynamics across developmental stages, young duplicate genes fall into three clusters that are preferentially expressed in early embryos, mid-stage embryos, and late-stage larvae. Early embryonic genes are involved in protein degradation and develop essentiality comparable to the genomic average. In mid-to-late embryos and L4-stage larvae, young genes are enriched in intestine, epidermal cells, coelomocytes, and amphid chemosensory neurons. Their molecular functions and inducible expression indicate potential roles in innate immune response and chemosensory perceptions, which may contribute to adaptation outside of the sperm.

Regulatory controls of duplicated gene expression during fiber development in allotetraploid cotton

Polyploidy complicates transcriptional regulation and increases phenotypic diversity in organisms. The dynamics of genetic regulation of gene expression between coresident subgenomes in polyploids remains to be understood. Here we document the genetic regulation of fiber development in allotetraploid cotton Gossypium hirsutum by sequencing 376 genomes and 2,215 time-series transcriptomes. We characterize 1,258 genes comprising 36 genetic modules that control staged fiber development and uncover genetic components governing their partitioned expression relative to subgenomic duplicated genes (homoeologs). Only about 30% of fiber quality-related homoeologs show phenotypically favorable allele aggregation in cultivars, highlighting the potential for subgenome additivity in fiber improvement. We envision a genome-enabled breeding strategy, with particular attention to 48 favorable alleles related to fiber phenotypes that have been subjected to purifying selection during domestication. Our work delineates the dynamics of gene regulation during fiber development and highlights the potential of subgenomic coordination underpinning phenotypes in polyploid plants.

Symmetry of duplicated gene expression in the Interferon Stimulated Pathways of the tetraploid common carp

Symmetry of duplicated gene expression in the Interferon Stimulated Pathways of the tetraploid common carp

RNA silencing proteins and small RNAs in oomycete plant pathogens and biocontrol agents.

Oomycetes cause several damaging diseases of plants and animals, and some species also act as biocontrol agents on insects, fungi, and other oomycetes. RNA silencing is increasingly being shown to play a role in the pathogenicity of Phytophthora species, either through trans-boundary movement of small RNAs (sRNAs) or through expression regulation of infection promoting effectors. To gain a wider understanding of RNA silencing in oomycete species with more diverse hosts, we mined genome assemblies for Dicer-like (DCL), Argonaute (AGO), and RNA dependent RNA polymerase (RDRP) proteins from Phytophthora plurivora, Ph. cactorum, Ph. colocasiae, Pythium oligandrum, Py. periplocum, and Lagenidium giganteum. Moreover, we sequenced small RNAs from the mycelium stage in each of these species. Each of the species possessed a single DCL protein, but they differed in the number and sequence of AGOs and RDRPs. SRNAs of 21nt, 25nt, and 26nt were prevalent in all oomycetes analyzed, but the relative abundance and 5' base preference of these classes differed markedly between genera. Most sRNAs mapped to transposons and other repeats, signifying that the major role for RNA silencing in oomycetes is to limit the expansion of these elements. We also found that sRNAs may act to regulate the expression of duplicated genes. Other sRNAs mapped to several gene families, and this number was higher in Pythium spp., suggesting a role of RNA silencing in regulating gene expression. Genes for most effector classes were the source of sRNAs of variable size, but some gene families showed a preference for specific classes of sRNAs, such as 25/26 nt sRNAs targeting RxLR effector genes in Phytophthora species. Novel miRNA-like RNAs (milRNAs) were discovered in all species, and two were predicted to target transcripts for RxLR effectors in Ph. plurivora and Ph. cactorum, indicating a putative role in regulating infection. Moreover, milRNAs from the biocontrol Pythium species had matches in the predicted transcriptome of Phytophthora infestans and Botrytis cinerea, and L. giganteum milRNAs matched candidate genes in the mosquito Aedes aegypti. This suggests that trans-boundary RNA silencing may have a role in the biocontrol action of these oomycetes.

Admixture of divergent genomes facilitates hybridization across species in the family Brassicaceae.

Summary Hybridization and polyploidization are pivotal to plant evolution. Genetic crosses between distantly related species are rare in nature due to reproductive barriers but how such hurdles can be overcome is largely unknown. Here we report the hybrid genome structure of xBrassicoraphanus, a synthetic allotetraploid of Brassica rapa and Raphanus sativus.We performed cytogenetic analysis and de novo genome assembly to examine chromosome behaviors and genome integrity in the hybrid. Transcriptome analysis was conducted to investigate expression of duplicated genes in conjunction with epigenome analysis to address whether genome admixture entails epigenetic reconfiguration.Allotetraploid xBrassicoraphanus retains both parental chromosomes without genome rearrangement. Meiotic synapsis formation and chromosome exchange are avoided between nonhomologous progenitor chromosomes. Reconfiguration of transcription network occurs, and less divergent cis‐elements of duplicated genes are associated with convergent expression. Genome‐wide DNA methylation asymmetry between progenitors is largely maintained but, notably, B. rapa‐originated transposable elements are transcriptionally silenced in xBrassicoraphanus through gain of DNA methylation.Our results demonstrate that hybrid genome stabilization and transcription compatibility necessitate epigenome landscape adjustment and rewiring of cis–trans interactions. Overall, this study suggests that a certain extent of genome divergence facilitates hybridization across species, which may explain the great diversification and expansion of angiosperms during evolution.

Developmental Evolution of Hypaxial Muscles: Insights From Cyclostomes and Chondrichthyans.

Jawed vertebrates possess two distinct groups of muscles in the trunk (epaxial and hypaxial muscles) primarily defined by the pattern of motor innervation from the spinal cord. Of these, the hypaxial group includes muscles with highly differentiated morphology and function, such as the muscles associated with paired limbs, shoulder girdles and tongue/infrahyoid (hypobranchial) muscles. Here we summarize the latest findings on the evolutionary mechanisms underlying the morphological variety of hypaxial musculature, with special reference to the molecular insights obtained from several living species that diverged early in vertebrate evolution. Lampreys, extant jawless vertebrates, lack many of derived traits characteristic of the gnathostomes, such as jaws, paired fins and epaxial/hypaxial distinction of the trunk skeletal musculatures. However, these animals possess the primitive form of the hypobranchial muscle. Of the gnathostomes, the elasmobranchs exhibit developmental mode of hypaxial muscles that is not identical to that of other gnathostomes in that the muscle primordia relocate as coherent cell aggregates. Comparison of expression of developmental genes, including Lbx genes, has delineated the temporal order of differentiation of various skeletal muscles, such as the hypobranchial, posterior pharyngeal and cucullaris (trapezius) muscles. We have proposed that the sequential addition of distal muscles, associated with expression of duplicated Lbx genes, promoted the elaboration of skeletal musculature. These analyses have revealed the framework of an evolutionary pathway that gave rise to the morphological complexity and diversity of vertebrate body patterns.

Gapless indica rice genome reveals synergistic contributions of active transposable elements and segmental duplications to rice genome evolution

The ultimate goal of genome assembly is a high-accuracy gapless genome. Here, we report a new assembly pipeline that is used to produce a gapless genome for the indica rice cultivar Minghui 63. The resulting 397.71-Mb final assembly is composed of 12 contigs with a contig N50 size of 31.93 Mb. Each chromosome is represented by a single contig and the genomic sequences of all chromosomes are gapless. Quality evaluation of this gapless genome assembly showed that gene regions in our assembly have the highest completeness compared with the other 15 reported high-quality rice genomes. Further comparison with the japonica rice genome revealed that the gapless indica genome assembly contains more transposable elements (TEs) and segmental duplications (SDs), the latter of which produce many duplicated genes that can affect agronomic traits through dose effect or sub-/neo-functionalization. The insertion of TEs can also affect the expression of duplicated genes, which may drive the evolution of these genes. Furthermore, we found the expansion of nucleotide-binding site with leucine-rich repeat disease-resistance genes and cis-zeatin-O-glucosyltransferase growth-related genes in SDs in the gapless indica genome assembly, suggesting that SDs contribute to the adaptive evolution of rice disease resistance and developmental processes. Collectively, our findings suggest that active TEs and SDs synergistically contribute to rice genome evolution.

Region-specific expression of young small-scale duplications in the human central nervous system

BackgroundThe duplication of genes is one of the main genetic mechanisms that led to the gain in complexity of biological tissue. Although the implication of duplicated gene expression in brain evolution was extensively studied through comparisons between organs, their role in the regional specialization of the adult human central nervous system has not yet been well described.ResultsOur work explored intra-organ expression properties of paralogs through multiple territories of the human central nervous system (CNS) using transcriptome data generated by the Genotype-Tissue Expression (GTEx) consortium. Interestingly, we found that paralogs were associated with region-specific expression in CNS, suggesting their involvement in the differentiation of these territories. Beside the influence of gene expression level on region-specificity, we observed the contribution of both duplication age and duplication type to the CNS region-specificity of paralogs. Indeed, we found that small scale duplicated genes (SSDs) and in particular ySSDs (SSDs younger than the 2 rounds of whole genome duplications) were more CNS region-specific than other paralogs. Next, by studying the two paralogs of ySSD pairs, we observed that when they were region-specific, they tend to be specific to the same region more often than for other paralogs, showing the high co-expression of ySSD pairs. The extension of this analysis to families of paralogs showed that the families with co-expressed gene members (i.e. homogeneous families) were enriched in ySSDs. Furthermore, these homogeneous families tended to be region-specific families, where the majority of their gene members were specifically expressed in the same region.ConclusionsOverall, our study suggests the involvement of ySSDs in the differentiation of human central nervous system territories. Therefore, we show the relevance of exploring region-specific expression of paralogs at the intra-organ level.

Bioinformatic analysis of chromatin organization and biased expression of duplicated genes between two poplars with a common whole-genome duplication

The nonrandom three-dimensional organization of chromatin plays an important role in the regulation of gene expression. However, it remains unclear whether this organization is conserved and whether it is involved in regulating gene expression during speciation after whole-genome duplication (WGD) in plants. In this study, high-resolution interaction maps were generated using high-throughput chromatin conformation capture (Hi-C) techniques for two poplar species, Populus euphratica and Populus alba var. pyramidalis, which diverged ~14 Mya after a common WGD. We examined the similarities and differences in the hierarchical chromatin organization between the two species, including A/B compartment regions and topologically associating domains (TADs), as well as in their DNA methylation and gene expression patterns. We found that chromatin status was strongly associated with epigenetic modifications and gene transcriptional activity, yet the conservation of hierarchical chromatin organization across the two species was low. The divergence of gene expression between WGD-derived paralogs was associated with the strength of chromatin interactions, and colocalized paralogs exhibited strong similarities in epigenetic modifications and expression levels. Thus, the spatial localization of duplicated genes is highly correlated with biased expression during the diploidization process. This study provides new insights into the evolution of chromatin organization and transcriptional regulation during the speciation process of poplars after WGD.

Copy number variations and young duplicate genes have high methylation levels in sticklebacks.

Gene duplication is an important driver of genomic diversity that can promote adaptive evolution. However, like most mutations, a newly duplicated gene is often deleterious and removed from the genome by drift or natural selection. The early molecular changes that occur soon after duplication therefore may influence the long-term survival of gene duplicates, but relatively little empirical data exist on the events near the onset of duplication before mutations have time to accumulate. In this study, we contrast gene expression and DNA methylation levels of duplicate genes in the threespine stickleback, Gasterosteus aculeatus, including recently emerged duplications that segregate as copy number variations (CNVs). We find that younger duplicate genes have higher levels of promoter methylation than older genes, and that gene CNVs have higher promoter methylation than non-CNVs. These results suggest preferential duplication of highly methylated genes or rapid methylation changes soon after duplication. We also find a negative association between methylation and expression, providing a putative role for methylation in suppressing transcription that compensates for increases in gene copy numbers and promoting paralog retention. We propose that methylation contributes to the longevity of young duplicate genes, extending the window of opportunity for functional divergence via mutation.

In situ dissecting the evolution of gene duplication with different histone modification patterns based on high-throughput data analysis in Arabidopsis thaliana.

BackgroundGenetic regulation is known to contribute to the divergent expression of duplicate genes; however, little is known about how epigenetic modifications regulate the expression of duplicate genes in plants.MethodsThe histone modification (HM) profile patterns of different modes of gene duplication, including the whole genome duplication, proximal duplication, tandem duplication and transposed duplication were characterized based on ChIP-chip or ChIP-seq datasets. In this study, 10 distinct HM marks including H2Bub, H3K4me1, H3K4me2, H3K4me3, H3K9ac, H3K9me2, H3K27me1, H3K27me3, H3K36me3 and H3K14ac were analyzed. Moreover, the features of gene duplication with different HM patterns were characterized based on 88 RNA-seq datasets of Arabidopsis thaliana.ResultsThis study showed that duplicate genes in Arabidopsis have a more similar HM pattern than single-copy genes in both their promoters and protein-coding regions. The evolution of HM marks is found to be coupled with coding sequence divergence and expression divergence after gene duplication. We found that functionally selective constraints may impose on epigenetic evolution after gene duplication. Furthermore, duplicate genes with distinct functions have more divergence in histone modification compared with the ones with the same function, while higher expression divergence is found with mutations of chromatin modifiers. This study shows the role of epigenetic marks in regulating gene expression and functional divergence after gene duplication in plants based on sequencing data.

Genome-wide identification, and phylogenetic and expression profiling analyses, of XTH gene families in Brassica rapa L. and Brassica oleracea L.

BackgroundXyloglucan endotransglucosylase/hydrolase genes (XTHs) are a multigene family and play key roles in regulating cell wall extensibility in plant growth and development. Brassica rapa and Brassica oleracea contain XTHs, but detailed identification and characterization of the XTH family in these species, and analysis of their tissue expression profiles, have not previously been carried out.ResultsIn this study, 53 and 38 XTH genes were identified in B. rapa and B. oleracea respectively, which contained some novel members not observed in previous studies. All XTHs of B. rapa, B. oleracea and Arabidopsis thaliana could be classified into three groups, Group I/II, III and the Early diverging group, based on phylogenetic relationships. Gene structures and motif patterns were similar within each group. All XTHs in this study contained two characteristic conserved domains (Glyco_hydro and XET_C). XTHs are located mainly in the cell wall but some are also located in the cytoplasm. Analyses of the mechanisms of gene family expansion revealed that whole-genome triplication (WGT) events and tandem duplication (TD) may have been the major mechanisms accounting for the expansion of the XTH gene family. Interestingly, TD genes all belonged to Group I/II, suggesting that TD was the main reason for the largest number of genes being in these groups. B. oleracea had lost more of the XTH genes, the conserved domain XET_C and the conserved active-site motif EXDXE compared with B. rapa, consistent with asymmetrical evolution between the two Brassica genomes. A majority of XTH genes exhibited different tissue-specific expression patterns based on RNA-seq data analyses. Moreover, there was differential expression of duplicated XTH genes in the two species, indicating that their functional differentiation occurred after B. rapa and B. oleracea diverged from a common ancestor.ConclusionsWe carried out the first systematic analysis of XTH gene families in B. rapa and B. oleracea. The results of this investigation can be used for reference in further studies on the functions of XTH genes and the evolution of this multigene family.

Global insights into duplicated gene expression and alternative splicing in polyploid Brassica napus under heat, cold, and drought stress.

Polyploidy has been a prevalent process during plant evolution and it has made a major impact on the structure and evolution of plant genomes. Many important crop plants are polyploid. There is considerable interest in expression patterns of duplicated genes in polyploids. Alternative splicing (AS) is a fundamental aspect of gene expression that produces multiple final transcript types from a single type of mRNAs. The effects of abiotic stress conditions on AS in polyploids has received little attention. We conducted a global transcriptome analysis of Brassica napus, an allotetraploid derived from B. rapa (AT ) and B. oleracea (CT ), by RNA-Seq of plants subjected to cold, heat, and drought stress treatments. Analyses of 27,360 pairs of duplicated genes revealed overall AT subgenome biases in gene expression and CT subgenome biases in the extent of alternative splicing under all three stress treatments. More genes increased in expression than decreased in response to the stresses. Negative correlations were found between expression levels and AS frequency for each type of AS. Cold stress produced the greatest changes in gene expression and AS. Cold-induced AS changes were more likely to be shared with those generated by drought than by heat stress. We used homeologs of FLC and CCA1 as case studies to show the dynamics of how duplicates in a polyploid respond to cold stress. Our results suggest that divergence in gene expression and AS patterns between duplicated genes may increase the flexibility of polyploids when responding to abiotic stressors.

Homoeologous gene expression and co-expression network analyses and evolutionary inference in allopolyploids.

Polyploidy is a widespread phenomenon throughout eukaryotes. Due to the coexistence of duplicated genomes, polyploids offer unique challenges for estimating gene expression levels, which is essential for understanding the massive and various forms of transcriptomic responses accompanying polyploidy. Although previous studies have explored the bioinformatics of polyploid transcriptomic profiling, the causes and consequences of inaccurate quantification of transcripts from duplicated gene copies have not been addressed. Using transcriptomic data from the cotton genus (Gossypium) as an example, we present an analytical workflow to evaluate a variety of bioinformatic method choices at different stages of RNA-seq analysis, from homoeolog expression quantification to downstream analysis used to infer key phenomena of polyploid expression evolution. In general, EAGLE-RC and GSNAP-PolyCat outperform other quantification pipelines tested, and their derived expression dataset best represents the expected homoeolog expression and co-expression divergence. The performance of co-expression network analysis was less affected by homoeolog quantification than by network construction methods, where weighted networks outperformed binary networks. By examining the extent and consequences of homoeolog read ambiguity, we illuminate the potential artifacts that may affect our understanding of duplicate gene expression, including an overestimation of homoeolog co-regulation and the incorrect inference of subgenome asymmetry in network topology. Taken together, our work points to a set of reasonable practices that we hope are broadly applicable to the evolutionary exploration of polyploids.

Aspontaneous complex structural variant in rcan-1 increases exploratory behavior and laboratory fitness of Caenorhabditis elegans.

Over long evolutionary timescales, major changes to the copy number, function, and genomic organization of genes occur, however, our understanding of the individual mutational events responsible for these changes is lacking. In this report, we study the genetic basis of adaptation of two strains of C. elegans to laboratory food sources using competition experiments on a panel of 89 recombinant inbred lines (RIL). Unexpectedly, we identified a single RIL with higher relative fitness than either of the parental strains. This strain also displayed a novel behavioral phenotype, resulting in higher propensity to explore bacterial lawns. Using bulk-segregant analysis and short-read resequencing of this RIL, we mapped the change in exploration behavior to a spontaneous, complex rearrangement of the rcan-1 gene that occurred during construction of the RIL panel. We resolved this rearrangement into five unique tandem inversion/duplications using Oxford Nanopore long-read sequencing. rcan-1 encodes an ortholog to human RCAN1/DSCR1 calcipressin gene, which has been implicated as a causal gene for Down syndrome. The genomic rearrangement in rcan-1 creates two complete and two truncated versions of the rcan-1 coding region, with a variety of modified 5' and 3' non-coding regions. While most copy-number variations (CNVs) are thought to act by increasing expression of duplicated genes, these changes to rcan-1 ultimately result in the reduction of its whole-body expression due to changes in the upstream regions. By backcrossing this rearrangement into a common genetic background to create a near isogenic line (NIL), we demonstrate that both the competitive advantage and exploration behavioral changes are linked to this complex genetic variant. This NIL strain does not phenocopy a strain containing an rcan-1 loss-of-function allele, which suggests that the residual expression of rcan-1 is necessary for its fitness effects. Our results demonstrate how colonization of new environments, such as those encountered in the laboratory, can create evolutionary pressure to modify gene function. This evolutionary mismatch can be resolved by an unexpectedly complex genetic change that simultaneously duplicates and diversifies a gene into two uniquely regulated genes. Our work shows how complex rearrangements can act to modify gene expression in ways besides increased gene dosage.