RESUMO
Centriole biogenesis, as in most organelle assemblies, involves the sequential recruitment of sub-structural elements that will support its function. To uncover this process, we correlated the spatial location of 24 centriolar proteins with structural features using expansion microscopy. A time-series reconstruction of protein distributions throughout human procentriole assembly unveiled the molecular architecture of the centriole biogenesis steps. We found that the process initiates with the formation of a naked cartwheel devoid of microtubules. Next, the bloom phase progresses with microtubule blade assembly, concomitantly with radial separation and rapid cartwheel growth. In the subsequent elongation phase, the tubulin backbone grows linearly with the recruitment of the A-C linker, followed by proteins of the inner scaffold (IS). By following six structural modules, we modeled 4D assembly of the human centriole. Collectively, this work provides a framework to investigate the spatial and temporal assembly of large macromolecules.
Assuntos
Centríolos , Microtúbulos , Centríolos/metabolismo , Humanos , Microtúbulos/metabolismo , Tubulina (Proteína)/metabolismo , Proteínas de Ciclo Celular/metabolismoRESUMO
Duplication is a foundation of molecular evolution and a driver of genomic and complex diseases. Here, we develop a genome editing tool named Amplification Editing (AE) that enables programmable DNA duplication with precision at chromosomal scale. AE can duplicate human genomes ranging from 20 bp to 100 Mb, a size comparable to human chromosomes. AE exhibits activity across various cell types, encompassing diploid, haploid, and primary cells. AE exhibited up to 73.0% efficiency for 1 Mb and 3.4% for 100 Mb duplications, respectively. Whole-genome sequencing and deep sequencing of the junctions of edited sequences confirm the precision of duplication. AE can create chromosomal microduplications within disease-relevant regions in embryonic stem cells, indicating its potential for generating cellular and animal models. AE is a precise and efficient tool for chromosomal engineering and DNA duplication, broadening the landscape of precision genome editing from an individual genetic locus to the chromosomal scale.
Assuntos
Duplicação Gênica , Edição de Genes , Genoma Humano , Humanos , Edição de Genes/métodos , Sistemas CRISPR-Cas/genética , DNA/genética , Animais , Células-Tronco Embrionárias/metabolismo , Cromossomos Humanos/genéticaRESUMO
Duplicated genes are thought to follow one of three evolutionary trajectories that resolve their redundancy: neofunctionalization, subfunctionalization or pseudogenization. Differences in expression patterns have been documented for many duplicated gene pairs and interpreted as evidence of subfunctionalization and a loss of redundancy. However, little is known about the functional impact of such differences and about their molecular basis. Here, we investigate the genetic and molecular basis for the partial loss of redundancy between the two BLADE-ON-PETIOLE genes BOP1 and BOP2 in Red Shepherd's Purse (Capsella rubella) compared to Arabidopsis (Arabidopsis thaliana). While both genes remain almost fully redundant in A. thaliana, BOP1 in C. rubella can no longer ensure wild-type floral organ numbers and suppress bract formation, due to an altered expression pattern in the region of the cryptic bract primordium. We use two complementary approaches, transgenic rescue of A. thaliana atbop1 atbop2 double mutants and deletions in the endogenous AtBOP1 promoter, to demonstrate that several BOP1 promoter regions containing conserved non-coding sequences interact in a non-additive manner to control BOP1 expression in the bract primordium, and that changes in these interactions underlie the evolutionary divergence between C. rubella and A. thaliana BOP1 expression and activity. Similarly, altered interactions between cis-regulatory regions underlie the divergence in functional promoter architecture related to the control of floral organ abscission by BOP1. These findings highlight the complexity of promoter architecture in plants and suggest that changes in the interactions between cis-regulatory elements are key drivers for evolutionary divergence in gene expression and the loss of redundancy.
RESUMO
Following whole-genome duplication (WGD), duplicate gene pairs (homoeologs) can evolve varying degrees of expression divergence. However, the determinants influencing these relative expression level differences (RFPKM) between homoeologs remain elusive. Here, we analyzed the RFPKM between homoeologs in three angiosperms, Nymphaea colorata, Nelumbo nucifera, and Acorus tatarinowii, all having undergone a single WGD since the origin of angiosperms. Our results show significant positive correlations in RFPKM of homoeologs among tissues within the same species, and among orthologs across these three species, indicating convergent expression balance/bias between homoeologous gene copies following independent WGDs. We linked RFPKM between homoeologs to gene attributes associated with dosage balance constraints, such as protein-protein interactions, lethal-phenotype scores in Arabidopsis (Arabidopsis thaliana) orthologs, domain numbers, and expression breadth. Notably, homoeologs with lower RFPKM often had more interactions and higher lethal-phenotype scores, indicating selective pressures favoring balanced expression. Also, homoeologs with lower RFPKM were more likely to be retained after WGDs in angiosperms. Within Nelumbo, greater RFPKM between homoeologs correlated with increased cis- and trans-regulatory differentiation between species, highlighting the ongoing escalation of gene expression divergence. We further found that expression degeneration in one copy of homoeologs is inclined towards nonfunctionalization. Our research highlights the importance of balanced expression, shaped by dosage balance constraints, in the evolutionary retention of homoeologs in plants.
RESUMO
Hybridization blurs species boundaries and leads to intertwined lineages resulting in reticulate evolution. Polyploidy, the outcome of whole genome duplication (WGD), has more recently been implicated in promoting and facilitating hybridization between polyploid species, potentially leading to adaptive introgression. However, because polyploid lineages are usually ephemeral states in the evolutionary history of life it is unclear whether WGD-potentiated hybridization has any appreciable effect on their diploid counterparts. Here, we develop a model of cytotype dynamics within mixed-ploidy populations to demonstrate that polyploidy can in fact serve as a bridge for gene flow between diploid lineages, where introgression is fully or partially hampered by the species barrier. Polyploid bridges emerge in the presence of triploid organisms, which despite critically low levels of fitness, can still allow the transfer of alleles between diploid states of independently evolving mixed-ploidy species. Notably, while marked genetic divergence prevents polyploid-mediated interspecific gene flow, we show that increased recombination rates can offset these evolutionary constraints, allowing a more efficient sorting of alleles at higher-ploidy levels before introgression into diploid gene pools. Additionally, we derive an analytical approximation for the rate of gene flow at the tetraploid level necessary to supersede introgression between diploids with nonzero introgression rates, which is especially relevant for plant species complexes, where interspecific gene flow is ubiquitous. Altogether, our results illustrate the potential impact of polyploid bridges on the (re)distribution of genetic material across ecological communities during evolution, representing a potential force behind reticulation.
Assuntos
Fluxo Gênico , Hibridização Genética , Modelos Genéticos , Poliploidia , Evolução Molecular , Diploide , AlelosRESUMO
Homosporous lycophytes (Lycopodiaceae) are a deeply diverged lineage in the plant tree of life, having split from heterosporous lycophytes (Selaginella and Isoetes) ~400 Mya. Compared to the heterosporous lineage, Lycopodiaceae has markedly larger genome sizes and remains the last major plant clade for which no chromosome-level assembly has been available. Here, we present chromosomal genome assemblies for two homosporous lycophyte species, the allotetraploid Huperzia asiatica and the diploid Diphasiastrum complanatum. Remarkably, despite that the two species diverged ~350 Mya, around 30% of the genes are still in syntenic blocks. Furthermore, both genomes had undergone independent whole genome duplications, and the resulting intragenomic syntenies have likewise been preserved relatively well. Such slow genome evolution over deep time is in stark contrast to heterosporous lycophytes and is correlated with a decelerated rate of nucleotide substitution. Together, the genomes of H. asiatica and D. complanatum not only fill a crucial gap in the plant genomic landscape but also highlight a potentially meaningful genomic contrast between homosporous and heterosporous species.
Assuntos
Genoma de Planta , Genômica , Genoma de Planta/genética , Tamanho do Genoma , Filogenia , Evolução MolecularRESUMO
Oomycete protists share phenotypic similarities with fungi, including the ability to cause plant diseases, but branch in a distant region of the tree of life. It has been suggested that multiple horizontal gene transfers (HGTs) from fungi-to-oomycetes contributed to the evolution of plant-pathogenic traits. These HGTs are predicted to include secreted proteins that degrade plant cell walls, a barrier to pathogen invasion and a rich source of carbohydrates. Using a combination of phylogenomics and functional assays, we investigate the diversification of a horizontally transferred xyloglucanase gene family in the model oomycete species Phytophthora sojae. Our analyses detect 11 xyloglucanase paralogs retained in P. sojae. Using heterologous expression in yeast, we show consistent evidence that eight of these paralogs have xyloglucanase function, including variants with distinct protein characteristics, such as a long-disordered C-terminal extension that can increase xyloglucanase activity. The functional variants analyzed subtend a phylogenetic node close to the fungi-to-oomycete transfer, suggesting the horizontally transferred gene was a bona fide xyloglucanase. Expression of three xyloglucanase paralogs in Nicotiana benthamiana triggers high-reactive oxygen species (ROS) generation, while others inhibit ROS responses to bacterial immunogens, demonstrating that the paralogs differentially stimulate pattern-triggered immunity. Mass spectrometry of detectable enzymatic products demonstrates that some paralogs catalyze the production of variant breakdown profiles, suggesting that secretion of variant xyloglucanases increases efficiency of xyloglucan breakdown as well as diversifying the damage-associated molecular patterns released. We suggest that this pattern of neofunctionalization and the variant host responses represent an aspect of the Red Queen host-pathogen coevolutionary dynamic.
Assuntos
Transferência Genética Horizontal , Glicosídeo Hidrolases , Filogenia , Glicosídeo Hidrolases/metabolismo , Glicosídeo Hidrolases/genética , Phytophthora/patogenicidade , Phytophthora/genética , Doenças das Plantas/microbiologia , Doenças das Plantas/parasitologia , Evolução Molecular , Duplicação GênicaRESUMO
Genome duplications and ploidy transitions have occurred in nearly every major taxon of eukaryotes, but they are far more common in plants than in animals. Due to the conservation of the nuclear:cytoplasmic volume ratio increased DNA content results in larger cells. In plants, polyploid organisms are larger than diploids as cell number remains relatively constant. Conversely, vertebrate body size does not correlate with cell size and ploidy as vertebrates compensate for increased cell size to maintain tissue architecture and body size. This has historically been explained by a simple reduction in cell number that matches the increase in cell size maintaining body size as ploidy increases, but here we show that the compensatory mechanisms that maintain body size in triploid zebrafish are tissue-specific: A) erythrocytes respond in the classical pattern with a reduced number of larger erythrocytes in circulation, B) muscle, a tissue comprised of polynucleated muscle fibers, compensates by reducing the number of larger nuclei such that myofiber and myotome size in unaffected by ploidy, and C) vascular tissue compensates by thickening blood vessel walls, possibly at the expense of luminal diameter. Understanding the physiological implications of ploidy on tissue function requires a detailed description of the specific mechanisms of morphological compensation occurring in each tissue to understand how ploidy changes affect development and physiology.
Assuntos
Poliploidia , Peixe-Zebra , Animais , Peixe-Zebra/genética , Ploidias , Tamanho Celular , Tamanho CorporalRESUMO
Actinidia ('Mihoutao' in Chinese) includes species with complex ploidy, among which diploid Actinidia chinensis and hexaploid Actinidia deliciosa are economically and nutritionally important fruit crops. Actinidia deliciosa has been proposed to be an autohexaploid (2n = 174) with diploid A. chinensis (2n = 58) as the putative parent. A CCS-based assembly anchored to a high-resolution linkage map provided a chromosome-resolved genome for hexaploid A. deliciosa yielded a 3.91-Gb assembly of 174 pseudochromosomes comprising 29 homologous groups with 6 members each, which contain 39 854 genes with an average of 4.57 alleles per gene. Here we provide evidence that much of the hexaploid genome matches diploid A. chinensis; 95.5% of homologous gene pairs exhibited >90% similarity. However, intragenome and intergenome comparisons of synteny indicate chromosomal changes. Our data, therefore, indicate that if A. deliciosa is an autoploid, chromosomal rearrangement occurred following autohexaploidy. A highly diversified pattern of gene expression and a history of rapid population expansion after polyploidisation likely facilitated the adaptation and niche differentiation of A. deliciosa in nature. The allele-defined hexaploid genome of A. deliciosa provides new genomic resources to accelerate crop improvement and to understand polyploid genome evolution.
Assuntos
Actinidia , Actinidia/genética , Mapeamento Cromossômico , Genoma de Planta/genética , Ploidias , Cromossomos , Frutas/genéticaRESUMO
Restoring cytonuclear stoichiometry is necessary after whole-genome duplication (WGD) and interspecific/intergeneric hybridization in plants. We investigated this phenomenon in auto- and allopolyploids of the Festuca-Lolium complex providing insights into the mechanisms governing cytonuclear interactions in early polyploid and hybrid generations. Our study examined the main processes potentially involved in restoring the cytonuclear balance after WGD comparing diploids and new and well-established autopolyploids. We uncovered that both the number of chloroplasts and the number of chloroplast genome copies were significantly higher in the newly established autopolyploids and grew further in more established autopolyploids. The increase in the copy number of the chloroplast genome exceeded the rise in the number of chloroplasts and fully compensated for the doubling of the nuclear genome. In addition, changes in nuclear and organelle gene expression were insignificant. Allopolyploid Festuca × Lolium hybrids displayed potential structural conflicts in parental protein variants within the cytonuclear complexes. While biased maternal allele expression has been observed in numerous hybrids, our results suggest that its role in cytonuclear stabilization in the Festuca × Lolium hybrids is limited. This study provides insights into the restoration of the cytonuclear stoichiometry, yet it emphasizes the need for future research to explore post-transcriptional regulation and its impact on cytonuclear gene expression stoichiometry. Our findings may enhance the understanding of polyploid plant evolution, with broader implications for the study of cytonuclear interactions in diverse biological contexts.
Assuntos
Núcleo Celular , Festuca , Lolium , Poliploidia , Festuca/genética , Lolium/genética , Núcleo Celular/genética , Núcleo Celular/metabolismo , Genoma de Planta/genética , Genoma de Cloroplastos , Cloroplastos/genética , Cloroplastos/metabolismo , Hibridização Genética , Regulação da Expressão Gênica de PlantasRESUMO
Aneuploidy is common in eukaryotes, often leading to decreased fitness. However, evidence from fungi and human tumur cells suggests that specific aneuploidies can be beneficial under stressful conditions and facilitate adaptation. In a previous evolutionary experiment with yeast, populations evolving under heat stress became aneuploid, only to later revert to euploidy after beneficial mutations accumulated. It was therefore suggested that aneuploidy is a "stepping stone" on the path to adaptation. Here, we test this hypothesis. We use Bayesian inference to fit an evolutionary model with both aneuploidy and mutation to the experimental results. We then predict the genotype frequency dynamics during the experiment, demonstrating that most of the evolved euploid population likely did not descend from aneuploid cells, but rather from the euploid wild-type population. Our model shows how the beneficial mutation supply-the product of population size and beneficial mutation rate-determines the evolutionary dynamics: with low supply, much of the evolved population descends from aneuploid cells; but with high supply, beneficial mutations are generated fast enough to outcompete aneuploidy due to its inherent fitness cost. Our results suggest that despite its potential fitness benefits under stress, aneuploidy can be an evolutionary "diversion" rather than a "stepping stone": it can delay, rather than facilitate, the adaptation of the population, and cells that become aneuploid may leave less descendants compared to cells that remain diploid.
Assuntos
Aneuploidia , Fungos , Humanos , Teorema de Bayes , DiploideRESUMO
Gene duplication is a major force driving evolutionary innovation. A classic example is generating new animal toxins via duplication of physiological protein-encoding genes and recruitment into venom. While this process drives the innovation of many animal venoms, reverse recruitment of toxins into nonvenomous cells remains unresolved. Using comparative genomics, we find members of the Membrane Attack Complex and Perforin Family (MAC) have been recruited into venom-injecting cells (cnidocytes), in soft and stony corals and sea anemones, suggesting that the ancestral MAC was a cnidocyte expressed toxin. Further investigation into the model sea anemone Nematostella vectensis reveals that three members have undergone Nematostella-specific duplications leading to their reverse recruitment into endomesodermal cells. Furthermore, simultaneous knockdown of all three endomesodermally expressed MACs leads to mis-development, supporting that these paralogs have nonvenomous function. By resolving the evolutionary history and function of MACs in Nematostella, we provide the first proof for reverse recruitment from venom to organismal development.
Assuntos
Evolução Molecular , Perforina , Anêmonas-do-Mar , Animais , Anêmonas-do-Mar/genética , Perforina/metabolismo , Perforina/genética , Duplicação Gênica , Venenos de Cnidários/genética , Venenos de Cnidários/metabolismo , Filogenia , Família MultigênicaRESUMO
Selaginellaceae, originated in the Carboniferous and survived the Permian-Triassic mass extinction, is the largest family of lycophyte, which is sister to other tracheophytes. It stands out from tracheophytes by exhibiting extraordinary habitat diversity and lacking polyploidization. The organelle genome-based phylogenies confirmed the monophyly of Selaginella, with six or seven subgenera grouped into two superclades, but the phylogenetic positions of the enigmatic Selaginella sanguinolenta clade remained problematic. Here, we conducted a phylogenomic study on Selaginellaceae utilizing large-scale nuclear gene data from RNA-seq to elucidate the phylogeny and explore the causes of the phylogenetic incongruence of the S. sanguinolenta clade. Our phylogenetic analyses resolved three different positions of the S. sanguinolenta clade, which were supported by the sorted three nuclear gene sets, respectively. The results from the gene flow test, species network inference, and plastome-based phylogeny congruently suggested a probable hybrid origin of the S. sanguinolenta clade involving each common ancestor of the two superclades in Selaginellaceae. The hybrid hypothesis is corroborated by the evidence from rhizophore morphology and spore micromorphology. The chromosome observation and Ks distributions further suggested hybridization accompanied by polyploidization. Divergence time estimation based on independent datasets from nuclear gene sets and plastid genome data congruently inferred that allopolyploidization occurred in the Early Triassic. To our best knowledge, the allopolyploidization in the Mesozoic reported here represents the earliest record of tracheophytes. Our study revealed a unique triad of phylogenetic positions for a hybrid-originated group with comprehensive evidence and proposed a hypothesis for retaining both parental alleles through gene conversion.
Assuntos
Filogenia , Poliploidia , Selaginellaceae , Selaginellaceae/genética , Transcriptoma , Fluxo GênicoRESUMO
The molecular underpinnings and consequences of cycles of whole-genome duplication (WGD) and subsequent gene loss through subgenome fractionation remain largely elusive. Endogenous drivers, such as transposable elements (TEs), have been postulated to shape genome-wide dominance and biased fractionation, leading to a conserved least-fractionated (LF) subgenome and a degenerated most-fractionated (MF) subgenome. In contrast, the role of exogenous factors, such as those induced by environmental stresses, has been overlooked. In this study, a chromosome-scale assembly of the alpine buckler mustard (Biscutella laevigata; Brassicaceae) that underwent a WGD event about 11 million years ago is coupled with transcriptional responses to heat, cold, drought, and herbivory to assess how gene expression is associated with differential gene retention across the MF and LF subgenomes. Counteracting the impact of TEs in reducing the expression and retention of nearby genes across the MF subgenome, dosage balance is highlighted as a main endogenous promoter of the retention of duplicated gene products under purifying selection. Consistent with the "turn a hobby into a job" model, about one-third of environment-responsive duplicates exhibit novel expression patterns, with one copy typically remaining conditionally expressed, whereas the other copy has evolved constitutive expression, highlighting exogenous factors as a major driver of gene retention. Showing uneven patterns of fractionation, with regions remaining unbiased, but with others showing high bias and significant enrichment in environment-responsive genes, this mesopolyploid genome presents evolutionary signatures consistent with an interplay of endogenous and exogenous factors having driven gene content following WGD-fractionation cycles.
Assuntos
Genoma de Planta , Duplicação Gênica , Evolução Molecular , Elementos de DNA Transponíveis , Estresse Fisiológico , Brassicaceae/genética , Regulação da Expressão Gênica de PlantasRESUMO
Sex chromosomes display remarkable diversity and variability among vertebrates. Compared with research on the X/Y and Z/W chromosomes, which have long evolutionary histories in mammals and birds, studies on the sex chromosomes at early evolutionary stages are limited. Here, we precisely assembled the genomes of homozygous XX female and YY male Lanzhou catfish (Silurus lanzhouensis) derived from an artificial gynogenetic family and a self-fertilized family, respectively. Chromosome 24 (Chr24) was identified as the sex chromosome based on resequencing data. Comparative analysis of the X and Y chromosomes showed an approximate 320 kb Y-specific region with a Y-specific duplicate of anti-Mullerian hormone type-II receptor (amhr2y), which is consistent with findings in two other Silurus species but on different chromosomes (Chr24 of S. meridionalis and Chr5 of S. asotus). Deficiency of amhr2y resulted in male-to-female sex reversal, indicating that amhr2y plays a male-determining role in S. lanzhouensis. Phylogenetic analysis and comparative genomics revealed that the common sex-determining gene amhr2y was initially translocated to Chr24 of the Silurus ancestor along with the expansion of transposable elements. Chr24 was maintained as the sex chromosome in S. meridionalis and S. lanzhouensis, whereas a sex-determining region transition triggered sex chromosome turnover from Chr24 to Chr5 in S. asotus. Additionally, gene duplication, translocation, and degeneration were observed in the Y-specific regions of Silurus species. These findings present a clear case for the early evolutionary trajectory of sex chromosomes, including sex-determining gene origin, repeat sequence expansion, gene gathering and degeneration in sex-determining region, and sex chromosome turnover.
RESUMO
Many adhesion proteins, evolutionarily related through gene duplication, exhibit distinct and precise interaction preferences and affinities crucial for cell patterning. Yet, the evolutionary paths by which these proteins acquire new specificities and prevent cross-interactions within their family members remain unknown. To bridge this gap, this study focuses on Drosophila Down syndrome cell adhesion molecule-1 (Dscam1) proteins, which are cell adhesion proteins that have undergone extensive gene duplication. Dscam1 evolved under strong selective pressure to achieve strict homophilic recognition, essential for neuronal self-avoidance and patterning. Through a combination of phylogenetic analyses, ancestral sequence reconstruction, and cell aggregation assays, we studied the evolutionary trajectory of Dscam1 exon 4 across various insect lineages. We demonstrated that recent Dscam1 duplications in the mosquito lineage bind with strict homophilic specificities without any cross-interactions. We found that ancestral and intermediate Dscam1 isoforms maintained their homophilic binding capabilities, with some intermediate isoforms also engaging in promiscuous interactions with other paralogs. Our results highlight the robust selective pressure for homophilic specificity integral to the Dscam1 function within the process of neuronal self-avoidance. Importantly, our study suggests that the path to achieving such selective specificity does not introduce disruptive mutations that prevent self-binding but includes evolutionary intermediates that demonstrate promiscuous heterophilic interactions. Overall, these results offer insights into evolutionary strategies that underlie adhesion protein interaction specificities.
Assuntos
Moléculas de Adesão Celular , Proteínas de Drosophila , Evolução Molecular , Animais , Proteínas de Drosophila/genética , Proteínas de Drosophila/metabolismo , Moléculas de Adesão Celular/genética , Moléculas de Adesão Celular/metabolismo , Filogenia , Duplicação Gênica , Drosophila/genética , Culicidae/genéticaRESUMO
Biological systems can gain complexity over time. While some of these transitions are likely driven by natural selection, the extent to which they occur without providing an adaptive benefit is unknown. At the molecular level, one example is heteromeric complexes replacing homomeric ones following gene duplication. Here, we build a biophysical model and simulate the evolution of homodimers and heterodimers following gene duplication using distributions of mutational effects inferred from available protein structures. We keep the specific activity of each dimer identical, so their concentrations drift neutrally without new functions. We show that for more than 60% of tested dimer structures, the relative concentration of the heteromer increases over time due to mutational biases that favor the heterodimer. However, allowing mutational effects on synthesis rates and differences in the specific activity of homo- and heterodimers can limit or reverse the observed bias toward heterodimers. Our results show that the accumulation of more complex protein quaternary structures is likely under neutral evolution, and that natural selection would be needed to reverse this tendency.
Assuntos
Evolução Molecular , Duplicação Gênica , Mutação , Mapas de Interação de Proteínas , Seleção Genética , Mapas de Interação de Proteínas/genética , Multimerização Proteica , Modelos Genéticos , Proteínas/genética , Proteínas/metabolismo , Proteínas/química , Simulação por ComputadorRESUMO
The evolution of gene families is complex, involving gene-level evolutionary events such as gene duplication, horizontal gene transfer, and gene loss, and other processes such as incomplete lineage sorting (ILS). Because of this, topological differences often exist between gene trees and species trees. A number of models have been recently developed to explain these discrepancies, the most realistic of which attempts to consider both gene-level events and ILS. When unified in a single model, the interaction between ILS and gene-level events can cause polymorphism in gene copy number, which we refer to as copy number hemiplasy (CNH). In this paper, we extend the Wright-Fisher process to include duplications and losses over several species, and show that the probability of CNH for this process can be significant. We study how well two unified models-multilocus multispecies coalescent (MLMSC), which models CNH, and duplication, loss, and coalescence (DLCoal), which does not-approximate the Wright-Fisher process with duplication and loss. We then study the effect of CNH on gene family evolution by comparing MLMSC and DLCoal. We generate comparable gene trees under both models, showing significant differences in various summary statistics; most importantly, CNH reduces the number of gene copies greatly. If this is not taken into account, the traditional method of estimating duplication rates (by counting the number of gene copies) becomes inaccurate. The simulated gene trees are also used for species tree inference with the summary methods ASTRAL and ASTRAL-Pro, demonstrating that their accuracy, based on CNH-unaware simulations calibrated on real data, may have been overestimated.
Assuntos
Evolução Molecular , Dosagem de Genes , Modelos Genéticos , Duplicação Gênica , Família Multigênica , Filogenia , Classificação/métodos , Simulação por ComputadorRESUMO
Ovarian high-grade serous carcinoma (HGSC) originates in the fallopian tube, with secretory cells carrying a TP53 mutation, known as p53 signatures, identified as potential precursors. p53 signatures evolve into serous tubal intraepithelial carcinoma (STIC) lesions, which in turn progress into invasive HGSC, which readily spreads to the ovary and disseminates around the peritoneal cavity. We recently investigated the genomic landscape of early- and late-stage HGSC and found higher ploidy in late-stage (median 3.1) than early-stage (median 2.0) samples. Here, to explore whether the high ploidy and possible whole-genome duplication (WGD) observed in late-stage disease were determined early in the evolution of HGSC, we analysed archival formalin-fixed paraffin-embedded (FFPE) samples from five HGSC patients. p53 signatures and STIC lesions were laser-capture microdissected and sequenced using shallow whole-genome sequencing (sWGS), while invasive ovarian/fallopian tube and metastatic carcinoma samples underwent macrodissection and were profiled using both sWGS and targeted next-generation sequencing. Results showed highly similar patterns of global copy number change between STIC lesions and invasive carcinoma samples within each patient. Ploidy changes were evident in STIC lesions, but not p53 signatures, and there was a strong correlation between ploidy in STIC lesions and invasive ovarian/fallopian tube and metastatic samples in each patient. The reconstruction of sample phylogeny for each patient from relative copy number indicated that high ploidy, when present, occurred early in the evolution of HGSC, which was further validated by copy number signatures in ovarian and metastatic tumours. These findings suggest that aberrant ploidy, suggestive of WGD, arises early in HGSC and is detected in STIC lesions, implying that the trajectory of HGSC may be determined at the earliest stages of tumour development. © 2024 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.
Assuntos
Cistadenocarcinoma Seroso , Neoplasias das Tubas Uterinas , Neoplasias Ovarianas , Proteína Supressora de Tumor p53 , Humanos , Feminino , Neoplasias Ovarianas/genética , Neoplasias Ovarianas/patologia , Cistadenocarcinoma Seroso/genética , Cistadenocarcinoma Seroso/patologia , Neoplasias das Tubas Uterinas/genética , Neoplasias das Tubas Uterinas/patologia , Proteína Supressora de Tumor p53/genética , Carcinoma in Situ/genética , Carcinoma in Situ/patologia , Gradação de Tumores , Variações do Número de Cópias de DNA , Mutação , Genômica/métodos , Sequenciamento Completo do Genoma , Ploidias , Pessoa de Meia-Idade , Biomarcadores Tumorais/genética , Progressão da DoençaRESUMO
BACKGROUND: Tubulins are major components of the eukaryotic cytoskeletons that are crucial in many cellular processes. Ciliated protists comprise one of the oldest eukaryotic lineages possessing cilia over their cell surface and assembling many diverse microtubular structures. As such, ciliates are excellent model organisms to clarify the origin and evolution of tubulins in the early stages of eukaryote evolution. Nonetheless, the evolutionary history of the tubulin subfamilies within and among ciliate classes is unclear. RESULTS: We analyzed the evolutionary pattern of ciliate tubulin gene family based on genomes/transcriptomes of 60 species covering 10 ciliate classes. Results showed: (1) Six tubulin subfamilies (α_Tub, ß_Tub, γ_Tub, δ_Tub, ε_Tub, and ζ_Tub) originated from the last eukaryotic common ancestor (LECA) were observed within ciliates. Among them, α_Tub, ß_Tub, and γ_Tub were present in all ciliate species, while δ_Tub, ε_Tub, and ζ_Tub might be independently lost in some species. (2) The evolutionary history of the tubulin subfamilies varied. Evolutionary history of ciliate γ_Tub, δ_Tub, ε_Tub, and ζ_Tub showed a certain degree of consistency with the phylogeny of species after the divergence of ciliate classes, while the evolutionary history of ciliate α_Tub and ß_Tub varied among different classes. (3) Ciliate α- and ß-tubulin isoforms could be classified into an "ancestral group" present in LECA and a "divergent group" containing only ciliate sequences. Alveolata-specific expansion events probably occurred within the "ancestral group" of α_Tub and ß_Tub. The "divergent group" might be important for ciliate morphological differentiation and wide environmental adaptability. (4) Expansion events of the tubulin gene family appeared to be consistent with whole genome duplication (WGD) events in some degree. More Paramecium-specific tubulin expansions were detected than Tetrahymena-specific ones. Compared to other Paramecium species, the Paramecium aurelia complex underwent a more recent WGD which might have experienced more tubulin expansion events. CONCLUSIONS: Evolutionary history among different tubulin gene subfamilies seemed to vary within ciliated protists. And the complex evolutionary patterns of tubulins among different ciliate classes might drive functional diversification. Our investigation provided meaningful information for understanding the evolution of tubulin gene family in the early stages of eukaryote evolution.