RAD21 is a driver of chromosome 8 gain in Ewing sarcoma to mitigate replication stress.

Aneuploidy, defined as whole-chromosome gain or loss, causes cellular stress but, paradoxically, is a frequent occurrence in cancers. Here, we investigate why ∼50% of Ewing sarcomas, driven by the EWS-FLI1 fusion oncogene, harbor chromosome 8 gains. Expression of the EWS-FLI1 fusion in primary cells causes replication stress that can result in cellular senescence. Using an evolution approach, we show that trisomy 8 mitigates EWS-FLI1-induced replication stress through gain of a copy of RAD21. Low-level ectopic expression of RAD21 is sufficient to dampen replication stress and improve proliferation in EWS-FLI1-expressing cells. Conversely, deleting one copy in trisomy 8 cells largely neutralizes the fitness benefit of chromosome 8 gain and reduces tumorgenicity of a Ewing sarcoma cancer cell line in soft agar assays. We propose that RAD21 promotes tumorigenesis through single gene copy gain. Such genes may explain some recurrent aneuploidies in cancer.

Allopolyploid subgenome identification and implications for evolutionary analysis.

Whole-genome duplications (WGDs) are widespread genomic events in eukaryotes that are hypothesized to contribute to the evolutionary success of many lineages, including flowering plants, Saccharomyces yeast, and vertebrates. WGDs generally can be classified into autopolyploids (ploidy increase descended from one species) or allopolyploids (ploidy increase descended from multiple species). Assignment of allopolyploid progenitor species (called subgenomes in the polyploid) is important to understanding the biology and evolution of polyploids, including the asymmetric subgenome evolution following hybridization (biased fractionation). Here, I review the different methodologies used to identify the ancestors of allopolyploid subgenomes, discuss the advantages and disadvantages of these methods, and outline the implications of how these methods affect the subsequent evolutionary analysis of these genomes.

The cryptic complex rearrangements involving the DMD gene: etiologic clues about phenotypical differences revealed by optical genome mapping.

BACKGROUND: Deletion or duplication in the DMD gene is one of the most common causes of Duchenne and Becker muscular dystrophy (DMD/BMD). However, the pathogenicity of complex rearrangements involving DMD, especially segmental duplications with unknown breakpoints, is not well understood. This study aimed to evaluate the structure, pattern, and potential impact of rearrangements involving DMD duplication. METHODS: Two families with DMD segmental duplications exhibiting phenotypical differences were recruited. Optical genome mapping (OGM) was used to explore the cryptic pattern of the rearrangements. Breakpoints were validated using long-range polymerase chain reaction combined with next-generation sequencing and Sanger sequencing. RESULTS: A multi-copy duplication involving exons 64-79 of DMD was identified in Family A without obvious clinical symptoms. Family B exhibited typical DMD neuromuscular manifestations and presented a duplication involving exons 10-13 of DMD. The rearrangement in Family A involved complex in-cis tandem repeats shown by OGM but retained a complete copy (reading frame) of DMD inferred from breakpoint validation. A reversed insertion with a segmental repeat was identified in Family B by OGM, which was predicted to disrupt the normal structure and reading frame of DMD after confirming the breakpoints. CONCLUSIONS: Validating breakpoint and rearrangement pattern is crucial for the functional annotation and pathogenic classification of genomic structural variations. OGM provides valuable insights into etiological analysis of DMD/BMD and enhances our understanding for cryptic effects of complex rearrangements.

New Approaches for Inferring Phylogenies in the Presence of Paralogs.

The availability of whole genome sequences was expected to supply essentially unlimited data for phylogenetics. However, strict reliance on single-copy genes for this purpose has drastically limited the amount of data that can be used. Here, we review several approaches for increasing the amount of data used for phylogenetic inference, focusing on methods that allow for the inclusion of duplicated genes (paralogs). Recently developed methods that are robust to high levels of incomplete lineage sorting also appear to be robust to the inclusion of paralogs, suggesting a promising way to take full advantage of genomic data. We discuss the pitfalls of these approaches, as well as further avenues for research.

Genomic analyses of new genes and their phenotypic effects reveal rapid evolution of essential functions in Drosophila development.

It is a conventionally held dogma that the genetic basis underlying development is conserved in a long evolutionary time scale. Ample experiments based on mutational, biochemical, functional, and complementary knockdown/knockout approaches have revealed the unexpectedly important role of recently evolved new genes in the development of Drosophila. The recent progress in the genome-wide experimental testing of gene effects and improvements in the computational identification of new genes (< 40 million years ago, Mya) open the door to investigate the evolution of gene essentiality with a phylogenetically high resolution. These advancements also raised interesting issues in techniques and concepts related to phenotypic effect analyses of genes, particularly of those that recently originated. Here we reported our analyses of these issues, including reproducibility and efficiency of knockdown experiment and difference between RNAi libraries in the knockdown efficiency and testing of phenotypic effects. We further analyzed a large data from knockdowns of 11,354 genes (~75% of the Drosophila melanogaster total genes), including 702 new genes (~66% of the species total new genes that aged < 40 Mya), revealing a similarly high proportion (~32.2%) of essential genes that originated in various Sophophora subgenus lineages and distant ancestors beyond the Drosophila genus. The transcriptional compensation effect from CRISPR knockout were detected for highly similar duplicate copies. Knockout of a few young genes detected analogous essentiality in various functions in development. Taken together, our experimental and computational analyses provide valuable data for detection of phenotypic effects of genes in general and further strong evidence for the concept that new genes in Drosophila quickly evolved essential functions in viability during development.

The sunflower genome provides insights into oil metabolism, flowering and Asterid evolution.

The domesticated sunflower, Helianthus annuus L., is a global oil crop that has promise for climate change adaptation, because it can maintain stable yields across a wide variety of environmental conditions, including drought. Even greater resilience is achievable through the mining of resistance alleles from compatible wild sunflower relatives, including numerous extremophile species. Here we report a high-quality reference for the sunflower genome (3.6 gigabases), together with extensive transcriptomic data from vegetative and floral organs. The genome mostly consists of highly similar, related sequences and required single-molecule real-time sequencing technologies for successful assembly. Genome analyses enabled the reconstruction of the evolutionary history of the Asterids, further establishing the existence of a whole-genome triplication at the base of the Asterids II clade and a sunflower-specific whole-genome duplication around 29 million years ago. An integrative approach combining quantitative genetics, expression and diversity data permitted development of comprehensive gene networks for two major breeding traits, flowering time and oil metabolism, and revealed new candidate genes in these networks. We found that the genomic architecture of flowering time has been shaped by the most recent whole-genome duplication, which suggests that ancient paralogues can remain in the same regulatory networks for dozens of millions of years. This genome represents a cornerstone for future research programs aiming to exploit genetic diversity to improve biotic and abiotic stress resistance and oil production, while also considering agricultural constraints and human nutritional needs.

The clinical importance of tandem exon duplication-derived substitutions.

Most coding genes in the human genome are annotated with multiple alternative transcripts. However, clear evidence for the functional relevance of the protein isoforms produced by these alternative transcripts is often hard to find. Alternative isoforms generated from tandem exon duplication-derived substitutions are an exception. These splice events are rare, but have important functional consequences. Here, we have catalogued the 236 tandem exon duplication-derived substitutions annotated in the GENCODE human reference set. We find that more than 90% of the events have a last common ancestor in teleost fish, so are at least 425 million years old, and twenty-one can be traced back to the Bilateria clade. Alternative isoforms generated from tandem exon duplication-derived substitutions also have significantly more clinical impact than other alternative isoforms. Tandem exon duplication-derived substitutions have >25 times as many pathogenic and likely pathogenic mutations as other alternative events. Tandem exon duplication-derived substitutions appear to have vital functional roles in the cell and may have played a prominent part in metazoan evolution.

The fragility of a structurally diverse duplication block triggers recurrent genomic amplification.

The human genome contains hundreds of large, structurally diverse blocks that are insufficiently represented in the reference genome and are thus not amenable to genomic analyses. Structural diversity in the human population suggests that these blocks are unstable in the germline; however, whether or not these blocks are also unstable in the cancer genome remains elusive. Here we report that the 500 kb block called KRTAP_region_1 (KRTAP-1) on 17q12-21 recurrently demarcates the amplicon of the ERBB2 (HER2) oncogene in breast tumors. KRTAP-1 carries numerous tandemly-duplicated segments that exhibit diversity within the human population. We evaluated the fragility of the block by cytogenetically measuring the distances between the flanking regions and found that spontaneous distance outliers (i.e DNA breaks) appear more frequently at KRTAP-1 than at the representative common fragile site (CFS) FRA16D. Unlike CFSs, KRTAP-1 is not sensitive to aphidicolin. The exonuclease activity of DNA repair protein Mre11 protects KRTAP-1 from breaks, whereas CtIP does not. Breaks at KRTAP-1 lead to the palindromic duplication of the ERBB2 locus and trigger Breakage-Fusion-Bridge cycles. Our results indicate that an insufficiently investigated area of the human genome is fragile and could play a crucial role in cancer genome evolution.

Diversified regulation of circadian clock gene expression following whole genome duplication.

Across taxa, circadian control of physiology and behavior arises from cell-autonomous oscillations in gene expression, governed by a networks of so-called 'clock genes', collectively forming transcription-translation feedback loops. In modern vertebrates, these networks contain multiple copies of clock gene family members, which arose through whole genome duplication (WGD) events during evolutionary history. It remains unclear to what extent multiple copies of clock gene family members are functionally redundant or have allowed for functional diversification. We addressed this problem through an analysis of clock gene expression in the Atlantic salmon, a representative of the salmonids, a group which has undergone at least 4 rounds of WGD since the base of the vertebrate lineage, giving an unusually large complement of clock genes. By comparing expression patterns across multiple tissues, and during development, we present evidence for gene- and tissue-specific divergence in expression patterns, consistent with functional diversification of clock gene duplicates. In contrast to mammals, we found no evidence for coupling between cortisol and circadian gene expression, but cortisol mediated non-circadian regulated expression of a subset of clock genes in the salmon gill was evident. This regulation is linked to changes in gill function necessary for the transition from fresh- to sea-water in anadromous fish. Overall, this analysis emphasises the potential for a richly diversified clock gene network to serve a mixture of circadian and non-circadian functions in vertebrate groups with complex genomes.

The SNAP hypothesis: Chromosomal rearrangements could emerge from positive Selection during Niche Adaptation.

The relative linear order of most genes on bacterial chromosomes is not conserved over evolutionary timescales. One explanation is that selection is weak, allowing recombination to randomize gene order by genetic drift. However, most chromosomal rearrangements are deleterious to fitness. In contrast, we propose the hypothesis that rearrangements in gene order are more likely the result of selection during niche adaptation (SNAP). Partial chromosomal duplications occur very frequently by recombination between direct repeat sequences. Duplicated regions may contain tens to hundreds of genes and segregate quickly unless maintained by selection. Bacteria exposed to non-lethal selections (for example, a requirement to grow on a poor nutrient) can adapt by maintaining a duplication that includes a gene that improves relative fitness. Further improvements in fitness result from the loss or inactivation of non-selected genes within each copy of the duplication. When genes that are essential in single copy are lost from different copies of the duplication, segregation is prevented even if the original selection is lifted. Functional gene loss continues until a new genetic equilibrium is reached. The outcome is a rearranged gene order. Mathematical modelling shows that this process of positive selection to adapt to a new niche can rapidly drive rearrangements in gene order to fixation. Signature features (duplication formation and divergence) of the SNAP model were identified in natural isolates from multiple species showing that the initial two steps in the SNAP process can occur with a remarkably high frequency. Further bioinformatic and experimental analyses are required to test if and to which extend the SNAP process acts on bacterial genomes.

Genome-Wide Analysis and Expression Profiling of the Glutathione Peroxidase-like Enzyme Gene Family in Solanum tuberosum.

Glutathione peroxidase-like enzyme is an important enzymatic antioxidant in plants. It is involved in scavenging reactive oxygen species, which can effectively prevent oxidative damage and improve resistance. GPXL has been studied in many plants but has not been reported in potatoes, the world's fourth-largest food crop. This study identified eight StGPXL genes in potatoes for the first time through genome-wide bioinformatics analysis and further studied the expression patterns of these genes using qRT-PCR. The results showed that the expression of StGPXL1 was significantly upregulated under high-temperature stress, indicating its involvement in potato defense against high-temperature stress, while the expression levels of StGPXL4 and StGPXL5 were significantly downregulated. The expression of StGPXL1, StGPXL2, StGPXL3, and StGPXL6 was significantly upregulated under drought stress, indicating their involvement in potato defense against drought stress. After MeJA hormone treatment, the expression level of StGPXL6 was significantly upregulated, indicating its involvement in the chemical defense mechanism of potatoes. The expression of all StGPXL genes is inhibited under biotic stress, which indicates that GPXL is a multifunctional gene family, which may endow plants with resistance to various stresses. This study will help deepen the understanding of the function of the potato GPXL gene family, provide comprehensive information for the further analysis of the molecular function of the potato GPXL gene family as well as a theoretical basis for potato molecular breeding.

Protein kinase C family evolution in jawed vertebrates.

Protein kinase C (PKC) was one of the first kinases identified in human cells. It is now known to constitute a family of kinases that respond to diacylglycerol, phosphatidylserine and for some family members, Ca2+. They have a plethora of different functions, such as cell cycle regulation, immune response and memory formation. In mammals, 12 PKC family members have been described, usually divided into 4 different subfamilies. We present here a comprehensive evolutionary analysis of the PKC genes in jawed vertebrates with special focus on the impact of the two tetraploidizations (1R and 2R) before the radiation of jawed vertebrates and the teleost tetraploidization (3R), as illuminated by synteny and paralogon analysis including many neighboring gene families. We conclude that the vertebrate predecessor had five PKC genes, as tunicates and lancelets still do, and that the PKC family should therefore ideally be organized into five subfamilies. The 1R and 2R events led to a total of 12 genes distributed among these five subfamilies. All 12 genes are still present in some of the major lineages of jawed vertebrates, including mammals, whereas birds and cartilaginous fishes have lost one member. The 3R event added another nine genes in teleosts, bringing the total to 21 genes. The zebrafish, a common experimental model animal, has retained 19. We have found no independent gene duplications. Thus, the genome doublings completely account for the complexity of this gene family in jawed vertebrates and have thereby had a huge impact on their evolution.

Reconstructing the evolutionary history of the oxytocin and vasotocin receptor gene family: Insights on whole genome duplication scenarios.

Vertebrate genome evolution remains a hotly debated topic, specifically as regards the number and the timing of putative rounds of whole genome duplication events. In this study, I sought to shed light to this conundrum through assessing the evolutionary history of the oxytocin/vasotocin receptor family. I performed ancestral analyses of the genomic segments containing oxytocin and vasotocin receptors (OTR-VTRs) by mapping them back to the reconstructed ancestral vertebrate/chordate karyotypes reported in five independent studies (Nakatani et al., 2007; Putnam et al., 2008; Smith and Keinath, 2015; Smith et al., 2018; Simakov et al., 2020) and found that two alternative scenarios can account for their evolution: one consistent with one round of whole genome duplication in the common ancestor of lampreys and gnathostomes, followed by segmental duplications in both lineages, and another consistent with two rounds of whole genome duplication, with the first occurring in the gnathostome-lamprey ancestor and the second in the jawed vertebrate ancestor. Combining the data reported here with synteny and phylogeny data reported in our previous study (Theofanopoulou et al., 2021), I put forward that a single round of whole genome duplication scenario is more consistent with the synteny and evolution of chromosomes where OTR-VTRs are encountered, without excluding the possibility of a scenario including two rounds of whole genome duplication. Although the analysis of one gene family is not able to capture the full complexity of vertebrate genome evolution, this study can provide solid insight, since the gene family used here has been meticulously analyzed for its genes' orthologous and paralogous relationships across species using high quality genomes.

Altered chromatin conformation and transcriptional regulation in watermelon following genome doubling.

Polyploidy has played a crucial role in plant evolution, development and function. Synthetic autopolyploid represents an ideal system to investigate the effects of polyploidization on transcriptional regulation. In this study, we deciphered the impact of genome duplication at phenotypic and molecular levels in watermelon. Overall, 88% of the genes in tetraploid watermelon followed a >1:1 dosage effect, and accordingly, differentially expressed genes were largely upregulated. In addition, a great number of hypomethylated regions (1688) were identified in an isogenic tetraploid watermelon. These differentially methylated regions were localized in promoters and intergenic regions and near transcriptional start sites of the identified upregulated genes, which enhances the importance of methylation in gene regulation. These changes were reflected in sophisticated higher-order chromatin structures. The genome doubling caused switching of 108 A and 626 B compartments that harbored genes associated with growth, development and stress responses.

Whole-Genome Duplication and Yeast's Fruitful Way of Life.

Studies on the fate of Saccharomyces cerevisiae paralogous gene pairs that arose through a whole-genome duplication event have shown diversification of retained duplicated genes. Paralogous functional specialization often results in improved function and/or novel function that could contribute to the adaptation of the organism to a new lifestyle. Here, we analyze and discuss particular case studies of paralogous functional diversification that could have played a role in the acquisition of yeast fermentative metabolism.

Origin, conservation, and loss of alternative splicing events that diversify the proteome in Saccharomycotina budding yeasts.

Many eukaryotes use RNA processing, including alternative splicing, to express multiple gene products from the same gene. The budding yeast Saccharomyces cerevisiae has been successfully used to study the mechanism of splicing and the splicing machinery, but alternative splicing in yeast is relatively rare and has not been extensively studied. Alternative splicing of SKI7/HBS1 is widely conserved, but yeast and a few other eukaryotes have replaced this one alternatively spliced gene with a pair of duplicated, unspliced genes as part of a whole genome doubling (WGD). We show that other examples of alternative splicing known to have functional consequences are widely conserved within Saccharomycotina. A common mechanism by which alternative splicing has disappeared is by replacement of an alternatively spliced gene with duplicate unspliced genes. This loss of alternative splicing does not always take place soon after duplication, but can take place after sufficient time has elapsed for speciation. Saccharomycetaceae that diverged before WGD use alternative splicing more frequently than S. cerevisiae, suggesting that WGD is a major reason for infrequent alternative splicing in yeast. We anticipate that WGDs in other lineages may have had the same effect. Having observed that two functionally distinct splice-isoforms are often replaced by duplicated genes allowed us to reverse the reasoning. We thereby identify several splice isoforms that are likely to produce two functionally distinct proteins because we find them replaced by duplicated genes in related species. We also identify some alternative splicing events that are not conserved in closely related species and unlikely to produce functionally distinct proteins.

Performances of NIPT for copy number variations at different sequencing depths using the semiconductor sequencing platform.

OBJECTIVE: To evaluate the performance of noninvasive prenatal testing (NIPT) and NIPT-PLUS for the detection of genome-wide microdeletion and microduplication syndromes (MMSs) at different sequencing depths. The NIPT sequencing depth was 0.15X, and the data volume was 3 million reads; the NIPT-PLUS sequencing depth was 0.4X, and the data volume was 8 million reads. METHODS: A cohort of 50,679 pregnancies was recruited. A total of 42,969 patients opted for NIPT, and 7710 patients opted for NIPT-PLUS. All high-risk cases were advised to undergo invasive prenatal diagnosis and were followed up. RESULTS: A total of 373 cases had a high risk of a copy number variation (CNV) as predicted by NIPT and NIPT-PLUS: NIPT predicted 250 high-risk CNVs and NIPT-PLUS predicted 123. NIPT-PLUS increased the detection rate by 1.02% (0.58% vs 1.60%, p < 0.001). A total of 291 cases accepted noninvasive prenatal diagnosis, with 197 cases of NIPT and 94 cases of NIPT-PLUS. The PPV of CNV > 10 Mb for NIPT-PLUS was significantly higher than that for NIPT (p = 0.02). The total PPV of NIPT-PLUS was 12.56% higher than that of NIPT (43.61% vs 30.96%, p = 0.03). CONCLUSION: NIPT-PLUS had a better performance in detecting CNVs in terms of the total detection rate and total PPV. However, great care must be taken in presenting results and providing appropriate counseling to patients when deeper sequencing is performed in clinical practice.

Pathogenic convergence of CNVs in genes functionally associated to a severe neuromotor developmental delay syndrome.

BACKGROUND: Complex developmental encephalopathy syndromes might be the consequence of unknown genetic alterations that are likely to contribute to the full neurological phenotype as a consequence of pathogenic gene combinations. METHODS: To identify the additional genetic contribution to the neurological phenotype, we studied as a test case a boy, with a KCNQ2 exon-7 partial duplication, by single-nucleotide polymorphism (SNP) microarray to detect copy-number variations (CNVs). RESULTS: The proband presented a cerebral palsy like syndrome with a severe motor and developmental encephalopathy. The SNP array analysis detected in the proband several de novo CNVs, nine partial gene losses (LRRC55, PCDH9, NALCN, RYR3, ELAVL2, CDH13, ATP1A2, SLC17A5, ANO3), and two partial gene duplications (PCDH19, EFNA5). The biological functions of these genes are associated with ion channels such as calcium, chloride, sodium, and potassium with several membrane proteins implicated in neural cell-cell interactions, synaptic transmission, and axon guidance. Pathogenically, these functions can be associated to cerebral palsy, seizures, dystonia, epileptic crisis, and motor neuron dysfunction, all present in the patient. CONCLUSIONS: Severe motor and developmental encephalopathy syndromes of unknown origin can be the result of a phenotypic convergence by combination of several genetic alterations in genes whose physiological function contributes to the neurological pathogenic mechanism.

The impact of whole genome duplications on the human gene regulatory networks.

This work studies the effects of the two rounds of Whole Genome Duplication (WGD) at the origin of the vertebrate lineage on the architecture of the human gene regulatory networks. We integrate information on transcriptional regulation, miRNA regulation, and protein-protein interactions to comparatively analyse the role of WGD and Small Scale Duplications (SSD) in the structural properties of the resulting multilayer network. We show that complex network motifs, such as combinations of feed-forward loops and bifan arrays, deriving from WGD events are specifically enriched in the network. Pairs of WGD-derived proteins display a strong tendency to interact both with each other and with common partners and WGD-derived transcription factors play a prominent role in the retention of a strong regulatory redundancy. Combinatorial regulation and synergy between different regulatory layers are in general enhanced by duplication events, but the two types of duplications contribute in different ways. Overall, our findings suggest that the two WGD events played a substantial role in increasing the multi-layer complexity of the vertebrate regulatory network by enhancing its combinatorial organization, with potential consequences on its overall robustness and ability to perform high-level functions like signal integration and noise control. Lastly, we discuss in detail the RAR/RXR pathway as an illustrative example of the evolutionary impact of WGD duplications in human.

Formation of new chromatin domains determines pathogenicity of genomic duplications.

Chromosome conformation capture methods have identified subchromosomal structures of higher-order chromatin interactions called topologically associated domains (TADs) that are separated from each other by boundary regions. By subdividing the genome into discrete regulatory units, TADs restrict the contacts that enhancers establish with their target genes. However, the mechanisms that underlie partitioning of the genome into TADs remain poorly understood. Here we show by chromosome conformation capture (capture Hi-C and 4C-seq methods) that genomic duplications in patient cells and genetically modified mice can result in the formation of new chromatin domains (neo-TADs) and that this process determines their molecular pathology. Duplications of non-coding DNA within the mouse Sox9 TAD (intra-TAD) that cause female to male sex reversal in humans, showed increased contact of the duplicated regions within the TAD, but no change in the overall TAD structure. In contrast, overlapping duplications that extended over the next boundary into the neighbouring TAD (inter-TAD), resulted in the formation of a new chromatin domain (neo-TAD) that was isolated from the rest of the genome. As a consequence of this insulation, inter-TAD duplications had no phenotypic effect. However, incorporation of the next flanking gene, Kcnj2, in the neo-TAD resulted in ectopic contacts of Kcnj2 with the duplicated part of the Sox9 regulatory region, consecutive misexpression of Kcnj2, and a limb malformation phenotype. Our findings provide evidence that TADs are genomic regulatory units with a high degree of internal stability that can be sculptured by structural genomic variations. This process is important for the interpretation of copy number variations, as these variations are routinely detected in diagnostic tests for genetic disease and cancer. This finding also has relevance in an evolutionary setting because copy-number differences are thought to have a crucial role in the evolution of genome complexity.