The genomic landscape of 2,023 colorectal cancers.

Colorectal carcinoma (CRC) is a common cause of mortality1, but a comprehensive description of its genomic landscape is lacking2-9. Here we perform whole-genome sequencing of 2,023 CRC samples from participants in the UK 100,000 Genomes Project, thereby providing a highly detailed somatic mutational landscape of this cancer. Integrated analyses identify more than 250 putative CRC driver genes, many not previously implicated in CRC or other cancers, including several recurrent changes outside the coding genome. We extend the molecular pathways involved in CRC development, define four new common subgroups of microsatellite-stable CRC based on genomic features and show that these groups have independent prognostic associations. We also characterize several rare molecular CRC subgroups, some with potential clinical relevance, including cancers with both microsatellite and chromosomal instability. We demonstrate a spectrum of mutational profiles across the colorectum, which reflect aetiological differences. These include the role of Escherichia colipks+ colibactin in rectal cancers10 and the importance of the SBS93 signature11-13, which suggests that diet or smoking is a risk factor. Immune-escape driver mutations14 are near-ubiquitous in hypermutant tumours and occur in about half of microsatellite-stable CRCs, often in the form of HLA copy number changes. Many driver mutations are actionable, including those associated with rare subgroups (for example, BRCA1 and IDH1), highlighting the role of whole-genome sequencing in optimizing patient care.

The co-evolution of the genome and epigenome in colorectal cancer.

Colorectal malignancies are a leading cause of cancer-related death1 and have undergone extensive genomic study2,3. However, DNA mutations alone do not fully explain malignant transformation4-7. Here we investigate the co-evolution of the genome and epigenome of colorectal tumours at single-clone resolution using spatial multi-omic profiling of individual glands. We collected 1,370 samples from 30 primary cancers and 8 concomitant adenomas and generated 1,207 chromatin accessibility profiles, 527 whole genomes and 297 whole transcriptomes. We found positive selection for DNA mutations in chromatin modifier genes and recurrent somatic chromatin accessibility alterations, including in regulatory regions of cancer driver genes that were otherwise devoid of genetic mutations. Genome-wide alterations in accessibility for transcription factor binding involved CTCF, downregulation of interferon and increased accessibility for SOX and HOX transcription factor families, suggesting the involvement of developmental genes during tumourigenesis. Somatic chromatin accessibility alterations were heritable and distinguished adenomas from cancers. Mutational signature analysis showed that the epigenome in turn influences the accumulation of DNA mutations. This study provides a map of genetic and epigenetic tumour heterogeneity, with fundamental implications for understanding colorectal cancer biology.

Phenotypic plasticity and genetic control in colorectal cancer evolution.

Genetic and epigenetic variation, together with transcriptional plasticity, contribute to intratumour heterogeneity1. The interplay of these biological processes and their respective contributions to tumour evolution remain unknown. Here we show that intratumour genetic ancestry only infrequently affects gene expression traits and subclonal evolution in colorectal cancer (CRC). Using spatially resolved paired whole-genome and transcriptome sequencing, we find that the majority of intratumour variation in gene expression is not strongly heritable but rather 'plastic'. Somatic expression quantitative trait loci analysis identified a number of putative genetic controls of expression by cis-acting coding and non-coding mutations, the majority of which were clonal within a tumour, alongside frequent structural alterations. Consistently, computational inference on the spatial patterning of tumour phylogenies finds that a considerable proportion of CRCs did not show evidence of subclonal selection, with only a subset of putative genetic drivers associated with subclone expansions. Spatial intermixing of clones is common, with some tumours growing exponentially and others only at the periphery. Together, our data suggest that most genetic intratumour variation in CRC has no major phenotypic consequence and that transcriptional plasticity is, instead, widespread within a tumour.

Assessing Outcomes of Patients Subject to Intensive Care to Facilitate Organ Donation: A Spanish Multicenter Prospective Study.

Intensive Care to facilitate Organ Donation (ICOD) consists of the initiation or continuation of intensive care measures in patients with a devastating brain injury (DBI) in whom curative treatment is deemed futile and death by neurological criteria (DNC) is foreseen, to incorporate organ donation into their end-of-life plans. In this study we evaluate the outcomes of patients subject to ICOD and identify radiological and clinical factors associated with progression to DNC. In this first prospective multicenter study we tested by multivariate regression the association of clinical and radiological severity features with progression to DNC. Of the 194 patients, 144 (74.2%) patients fulfilled DNC after a median of 25 h (95% IQR: 17-44) from ICOD onset. Two patients (1%) shifted from ICOD to curative treatment, both were alive at discharge. Factors associated with progression to DNC included: age below 70 years, clinical score consistent with severe brain injury, instability, intracranial hemorrhage, midline shift ≥5 mm and certain types of brain herniation. Overall 151 (77.8%) patients progressed to organ donation. Based on these results, we conclude that ICOD is a beneficial and efficient practice that can contribute to the pool of deceased donors.

CRISPRi-induced transcriptional regulation of IAH1 gene and its influence on volatile compounds profile in Kluyveromyces marxianus DU3.

Mezcal is a traditional Mexican distilled beverage, known for its marked organoleptic profile, which is influenced by several factors, such as the fermentation process, where a wide variety of microorganisms are present. Kluyveromyces marxianus is one of the main yeasts isolated from mezcal fermentations and has been associated with ester synthesis, contributing to the flavors and aromas of the beverage. In this study, we employed CRISPR interference (CRISPRi) technology, using dCas9 fused to the Mxi1 repressor factor domain, to down-regulate the expression of the IAH1 gene, encoding for an isoamyl acetate-hydrolyzing esterase, in K. marxianus strain DU3. The constructed CRISPRi plasmid successfully targeted the IAH1 gene, allowing for specific gene expression modulation. Through gene expression analysis, we assessed the impact of IAH1 down-regulation on the metabolic profile of volatile compounds. We also measured the expression of other genes involved in volatile compound biosynthesis, including ATF1, EAT1, ADH1, and ZWF1 by RT-qPCR. Results demonstrated successful down-regulation of IAH1 expression in K. marxianus strain DU3 using the CRISPRi system. The modulation of IAH1 gene expression resulted in alterations in the production of volatile compounds, specifically ethyl acetate, which are important contributors to the beverage's aroma. Changes in the expression levels of other genes involved in ester biosynthesis, suggesting that the knockdown of IAH1 may generate intracellular alterations in the balance of these metabolites, triggering a regulatory response. The application of CRISPRi technology in K. marxianus opens the possibility of targeted modulation of gene expression, metabolic engineering strategies, and synthetic biology in this yeast strain.

Current research on viral proteins that interact with fibrillarin.

The nucleolus is a multifunctional nuclear domain primarily dedicated to ribosome biogenesis. Certain viruses developed strategies to manipulate host nucleolar proteins to facilitate their replication by modulating ribosomal RNA (rRNA) processing. This association interferes with nucleolar functions resulting in overactivation or arrest of ribosome biogenesis, induction or inhibition of apoptosis, and affecting stress response. The nucleolar protein fibrillarin (FBL) is an important target of some plant and animal viruses. FBL is an essential and highly conserved S-adenosyl methionine (SAM) dependent methyltransferase, capable of rRNA degradation by its intrinsically disordered region (IDR), the glycine/arginine-rich (GAR) domain. It forms a ribonucleoprotein complex that directs 2'-O-methylations in more than 100 sites of pre-rRNAs. It is involved in multiple cellular processes, including initiation of transcription, oncogenesis, and apoptosis, among others. The interaction with animal viruses, including human viruses, triggered its redistribution to the nucleoplasm and cytoplasm, interfering with its role in pre-rRNA processing. Viral-encoded proteins with IDRs as nucleocapsids, matrix, Tat protein, and even a viral snoRNA, can associate with FBL, forcing the nucleolar protein to undergo atypical functions. Here we review the molecular mechanisms employed by animal and human viruses to usurp FBL functions and the effect on cellular processes, particularly in ribosome biogenesis.

LIPRNAseq: a method to discover lipid interacting RNAs by sequencing.

BACKGROUND: Current biological research extensively describes the interactions of molecules such as RNA with other nucleic acids or proteins. However, the relatively recent discovery of nuclear phospholipids playing biologically relevant processes outside membranes, as well as, RNA-lipid interactions shows the need for new methods to explore the identity of these RNAs. METHODS AND RESULTS: In this study, we describe the method for LIPID-RNA isolation followed by sequencing and analysis of the RNA that has the ability to interact with the selected lipids. Here we utilized specific phospholipid coated beads for selective RNA binding. We tested RNA from organisms belonging to different realms (human, plant, and yeast), and tested their ability to bind a specific lipid. CONCLUSIONS: The results show several RNAs differentially enriched in the pull-down of phosphatidyl Inositol 4,5 bisphosphate coated beads. This method is helpful to screen lipid-binding RNA, which may have relevant biological functions. The method can be used with different lipids and comparison of pull-downs and can narrow the selection of RNAs that interact with a particular lipid for further studies.

Binaural Beat Effects on Attention: A Study Based on the Oddball Paradigm.

The impact of binaural beats (BBs) on human cognition and behavior remains and various methods have been used to measure their effect, including neurophysiological, psychometric, and human performance evaluations. The few approaches where the level of neural synchronicity and connectivity were measured by neuroimaging techniques have only been undertaken in spontaneous mode. The present research proposes an approach based on the oddball paradigm to study BB effect by estimating the level of attention induced by BBs. Evoked activity of 25 young adults between 19 and 24 years old with no hearing impairments nor clinical neurological history were analyzed. The experiment was conducted in two different sessions of 24.5 min. The first part consisted of 20-min BB stimulation in either theta (BBθ) or beta (BBß). After the BB stimulation, an oddball paradigm was applied in each BB condition to assess the attentional effect induced by BBs. Attention enhancement is expected for BBß with respect to BBθ. Target event related potentials (ERPs) were mainly analyzed in the time and time-frequency domains. The frequency analysis was based on continuous wavelet transform (CWT), event-related spectral perturbation (ERSP), and inter-trial phase coherence (ITPC). The study revealed that the P300 component was not significantly different between conditions (BBθ vs. BBß). However, the target grand average ERP in BBθ condition was mainly composed of 8 Hz-frequency components, appearing before 400 ms post-stimulus, and mainly on the centro-parietal regions. In contrast, the target grand average ERP in BBß condition was mainly composed of frequency components below 6 Hz, mainly appearing at 400 ms post-stimulus on the parieto-occipital regions. Furthermore, ERPs in the BBθ condition were more phase locked than the BBß condition.

Parameter sensitivities in tilted-pulse-front based terahertz setups and their implications for high-energy terahertz source design and optimization.

Despite the popularity and ubiquity of the tilted-pulse-front technique for single-cycle terahertz (THz) pulse generation, there is a deficit of experimental studies comprehensively mapping out the dependence of the performance on key setup parameters. The most critical parameters include the pulse-front tilt, the effective length of the pump pulse propagation within the crystal as well as effective length over which the THz beam interacts with the pump before it spatially walks off. Therefore, we investigate the impact of these parameters on the conversion efficiency and the shape of the THz beam via systematically scanning the 5D parameter space spanned by pump fluence, pulse-front-tilt, crystal-position (2D), and the pump size experimentally. We verify predictions so far only made by theory regarding the optimum interaction lengths and map out the impact of cascading on the THz radiation generation process. Furthermore, distortions imposed on the spatial THz beam profile for larger than optimum interaction lengths are observed. Finally, we identify the most sensitive parameters and, based on our findings, propose a robust optimization strategy for tilted-pulse-front THz setups. These findings are relevant for all THz strong-field applications in high demand of robust high-energy table-top single-cycle THz sources such as THz plasmonics, high-harmonic generation in solids as well as novel particle accelerators and beam manipulators.

Fibrillarin evolution through the Tree of Life: Comparative genomics and microsynteny network analyses provide new insights into the evolutionary history of Fibrillarin.

Fibrillarin (FIB), a methyltransferase essential for life in the vast majority of eukaryotes, is involved in methylation of rRNA required for proper ribosome assembly, as well as methylation of histone H2A of promoter regions of rRNA genes. RNA viral progression that affects both plants and animals requires FIB proteins. Despite the importance and high conservation of fibrillarins, there little is known about the evolutionary dynamics of this small gene family. We applied a phylogenomic microsynteny-network approach to elucidate the evolutionary history of FIB proteins across the Tree of Life. We identified 1063 non-redundant FIB sequences across 1049 completely sequenced genomes from Viruses, Bacteria, Archaea, and Eukarya. FIB is a highly conserved single-copy gene through Archaea and Eukarya lineages, except for plants, which have a gene family expansion due to paleopolyploidy and tandem duplications. We found a high conservation of the FIB genomic context during plant evolution. Surprisingly, FIB in mammals duplicated after the Eutheria split (e.g., ruminants, felines, primates) from therian mammals (e.g., marsupials) to form two main groups of sequences, the FIB and FIB-like groups. The FIB-like group transposed to another genomic context and remained syntenic in all the eutherian mammals. This transposition correlates with differences in the expression patterns of FIB-like proteins and with elevated Ks values potentially due to reduced evolutionary constraints of the duplicated copy. Our results point to a unique evolutionary event in mammals, between FIB and FIB-like genes, that led to non-redundant roles of the vital processes in which this protein is involved.

Plant viral proteins and fibrillarin: the link to complete the infective cycle.

The interaction between viruses with the nucleolus is already a well-defined field of study in plant virology. This interaction is not restricted to those viruses that replicate in the nucleus, in fact, RNA viruses that replicate exclusively in the cytoplasm express proteins that localize in the nucleolus. Some positive single stranded RNA viruses from animals and plants have been reported to interact with the main nucleolar protein, Fibrillarin. Among nucleolar proteins, Fibrillarin is an essential protein that has been conserved in sequence and function throughout evolution. Fibrillarin is a methyltransferase protein with more than 100 methylation sites in the pre-ribosomal RNA, involved in multiple cellular processes, including initiation of transcription, oncogenesis, and apoptosis, among others. Recently, it was found that AtFib2 shows a ribonuclease activity. In plant viruses, Fibrillarin is involved in long-distance movement and cell-to-cell movement, being two highly different processes. The mechanism that Fibrillarin performs is still unknown. However, and despite belonging to very different viral families, the majority comply with the following. (1) They are positive single stranded RNA viruses; (2) encode different types of viral proteins that partially localize in the nucleolus; (3) interacts with Fibrillarin exporting it to the cytoplasm; (4) the viral protein-Fibrillarin interaction forms an RNP complex with the viral RNA and; (5) Fibrillarin depletion affects the infective cycle of the virus. Here we review the relationship of those plant viruses with Fibrillarin interaction, with special focus on the molecular processes of the virus to sequester Fibrillarin to complete its infective cycle.

Phase Separation of Intrinsically Disordered Nucleolar Proteins Relate to Localization and Function.

The process of phase separation allows for the establishment and formation of subcompartmentalized structures, thus enabling cells to perform simultaneous processes with precise organization and low energy requirements. Chemical modifications of proteins, RNA, and lipids alter the molecular environment facilitating enzymatic reactions at higher concentrations in particular regions of the cell. In this review, we discuss the nucleolus as an example of the establishment, dynamics, and maintenance of a membraneless organelle with a high level of organization.

Phylogenomic and Microsynteny Analysis Provides Evidence of Genome Arrangements of High-Affinity Nitrate Transporter Gene Families of Plants.

Nitrate transporter 2 (NRT2) and NRT3 or nitrate-assimilation-related 2 (NAR2) proteins families form a two-component, high-affinity nitrate transport system, which is essential for the acquisition of nitrate from soils with low N availability. An extensive phylogenomic analysis across land plants for these families has not been performed. In this study, we performed a microsynteny and orthology analysis on the NRT2 and NRT3 genes families across 132 plants (Sensu lato) to decipher their evolutionary history. We identified significant differences in the number of sequences per taxonomic group and different genomic contexts within the NRT2 family that might have contributed to N acquisition by the plants. We hypothesized that the greater losses of NRT2 sequences correlate with specialized ecological adaptations, such as aquatic, epiphytic, and carnivory lifestyles. We also detected expansion on the NRT2 family in specific lineages that could be a source of key innovations for colonizing contrasting niches in N availability. Microsyntenic analysis on NRT3 family showed a deep conservation on land plants, suggesting a high evolutionary constraint to preserve their function. Our study provides novel information that could be used as guide for functional characterization of these gene families across plant lineages.

Prediction of Benefit from Checkpoint Inhibitors in Mismatch Repair Deficient Metastatic Colorectal Cancer: Role of Tumor Infiltrating Lymphocytes.

BACKGROUND: Immunotherapy with immune checkpoint inhibitors (ICIs) is highly effective in microsatellite instability-high (MSI-H) metastatic colorectal cancer (mCRC); however, specific predictive biomarkers are lacking. PATIENTS AND METHODS: Data and samples from 85 patients with MSI-H mCRC treated with ICIs were gathered. Tumor infiltrating lymphocytes (TILs) and tumor mutational burden (TMB) were analyzed in an exploratory cohort of "super" responders and "clearly" refractory patients; TILs were then evaluated in the whole cohort of patients. Primary objectives were the correlation between the number of TILs and TMB and their role as biomarkers of ICI efficacy. Main endpoints included response rate (RR), progression-free survival (PFS), and overall survival (OS). RESULTS: In the exploratory cohort, an increasing number of TILs correlated to higher TMB (Pearson's test, p = .0429). In the whole cohort, median number of TILs was 3.6 in responders compared with 1.8 in nonresponders (Mann-Whitney test, p = .0448). RR was 70.6% in patients with high number of TILs (TILs-H) compared with 42.9% in patients with low number of TILs (odds ratio = 3.20, p = .0291). Survival outcomes differed significantly in favor of TILs-H (PFS: hazard ratio [HR] = 0.42, p = .0278; OS: HR = 0.41, p = .0463). CONCLUSION: A significant correlation between higher TMB and increased number of TILs was shown. A significantly higher activity and better PFS and OS with ICI in MSI-H mCRC were reported in cases with high number of TILs, thus supporting further studies of TIL count as predictive biomarker of ICI efficacy. IMPLICATIONS FOR PRACTICE: Microsatellite instability is the result of mismatch repair protein deficiency, caused by germline mutations or somatic modifications in mismatch repair genes. In metastatic colorectal cancer (mCRC), immunotherapy (with immune checkpoint inhibitors [ICIs]) demonstrated remarkable clinical benefit in microsatellite instability-high (MSI-H) patients. ICI primary resistance has been observed in approximately 25% of patients with MSI-H mCRC, underlining the need for predictive biomarkers. In this study, tumor mutational burden (TMB) and tumor infiltrating lymphocyte (TIL) analyses were performed in an exploratory cohort of patients with MSI-H mCRC treated with ICIs, demonstrating a significant correlation between higher TMB and increased number of TILs. Results also demonstrated a significant correlation between high number of TILs and clinical responses and survival benefit in a large data set of patients with MSI-H mCRC treated with ICI. TMB and TILs could represent predictive biomarkers of ICI efficacy in MSI-H mCRC and should be incorporated in future trials testing checkpoint inhibitors in colorectal cancer.

The SoNAP gene from sugarcane (Saccharum officinarum) encodes a senescence-associated NAC transcription factor involved in response to osmotic and salt stress.

Climate change has caused serious problems related to the productivity of agricultural crops directly affecting human well-being. Plants have evolved to produce molecular mechanisms in response to environmental stresses, such as transcription factors (TFs), to cope with abiotic stress. The NAC proteins constitute a plant-specific family of TFs involved in plant development processes and tolerance to biotic and abiotic stress. Sugarcane is a perennial grass that accumulates a large amount of sucrose and is a crucial agro-industry crop in tropical regions. Our previous transcriptome analyses on sugarcane that were exposed to drought conditions revealed significant increases in the expression of several NAC TFs through all of the time-point stress conditions. In this work, we characterize all previously detected sugarcane NAC genes, utilizing phylogenetics and expression analyses. Furthermore, we characterized a sugarcane NAC gene orthologous to the senescence-associated genes AtNAP and OsNAP via transient expression in tobacco calluses, from Arabidopsis and rice respectively, thus we named it the SoNAP gene. Transient localization assays on onion epidermal cells confirmed the nuclear localization of the SoNAP. Expression analysis showed that the SoNAP gene was induced by high salinity, drought, and abscisic acid treatments. The overexpression of the SoNAP gene in tobacco calluses caused a senescence associated phenotype. Overall, our results indicated that the SoNAP gene from sugarcane is transcriptionally induced under abiotic stress conditions and conserved the predicted senescence-associated functions when it was overexpressed in a heterologous plant model. This work provides key insights about the senescence mechanisms related to abiotic stress and it provides a benchmark for future work on the improvement of this economically important crop.

Allele balance bias identifies systematic genotyping errors and false disease associations.

In recent years, next-generation sequencing (NGS) has become a cornerstone of clinical genetics and diagnostics. Many clinical applications require high precision, especially if rare events such as somatic mutations in cancer or genetic variants causing rare diseases need to be identified. Although random sequencing errors can be modeled statistically and deep sequencing minimizes their impact, systematic errors remain a problem even at high depth of coverage. Understanding their source is crucial to increase precision of clinical NGS applications. In this work, we studied the relation between recurrent biases in allele balance (AB), systematic errors, and false positive variant calls across a large cohort of human samples analyzed by whole exome sequencing (WES). We have modeled the AB distribution for biallelic genotypes in 987 WES samples in order to identify positions recurrently deviating significantly from the expectation, a phenomenon we termed allele balance bias (ABB). Furthermore, we have developed a genotype callability score based on ABB for all positions of the human exome, which detects false positive variant calls that passed state-of-the-art filters. Finally, we demonstrate the use of ABB for detection of false associations proposed by rare variant association studies. Availability: https://github.com/Francesc-Muyas/ABB.

Single-cell genetic analysis of clonal dynamics in colorectal adenomas indicates CDX2 gain as a predictor of recurrence.

Colorectal adenomas are common precancerous lesions with the potential for malignant transformation to colorectal adenocarcinoma. Endoscopic polypectomy provides an opportunity for cancer prevention; however, recurrence rates are high. We collected formalin-fixed paraffin-embedded tissue of 15 primary adenomas with recurrence, 15 adenomas without recurrence, and 14 matched pair samples (primary adenoma and the corresponding recurrent adenoma). The samples were analysed by array-comparative genomic hybridisation (aCGH) and single-cell multiplex interphase fluorescence in situ hybridisation (miFISH) to understand clonal evolution, to examine the dynamics of copy number alterations (CNAs) and to identify molecular markers for recurrence prediction. The miFISH probe panel consisted of 14 colorectal carcinogenesis-relevant genes (COX2, PIK3CA, APC, CLIC1, EGFR, MYC, CCND1, CDX2, CDH1, TP53, HER2, SMAD7, SMAD4 and ZNF217), and a centromere probe (CEP10). The aCGH analysis confirmed the genetic landscape typical for colorectal tumorigenesis, that is, CNAs of chromosomes 7, 13q, 18 and 20q. Focal aberrations (≤10 Mbp) were mapped to chromosome bands 6p22.1-p21.33 (33.3%), 7q22.1 (31.4%) and 16q21 (29.4%). MiFISH detected gains of EGFR (23.6%), CDX2 (21.8%) and ZNF217 (18.2%). Most adenomas exhibited a major clone population which was accompanied by multiple smaller clone populations. Gains of CDX2 were exclusively seen in primary adenomas with recurrence (25%) compared to primary adenomas without recurrence (0%). Generation of phylogenetic trees for matched pair samples revealed four distinct patterns of clonal dynamics. In conclusion, adenoma development and recurrence are complex genetic processes driven by multiple CNAs whose evaluations by miFISH, with emphasis on CDX2, might serve as a predictor of recurrence.

The effect of plant stress on phosphoinositides.

Phosphoinositides are very versatile molecules with a plethora of functions such as cytokinesis, chemotaxis, cell survival, and cell death. Their functions depend on the proteins with which they interact. Thus, when interacting with phospholipases, phosphatases, or kinases, they can be precursors of second messengers in different signalling pathways. They could be second messengers themselves and interact directly with other proteins to modulate their functions trough changing its localization and activity or enhancing its synthesis rate. Because they are more abundant in animal cells and their importance in diseases such as cancer has taken priority, the study of the phosphoinositides in plants has not evolved to the same extent. Nevertheless, several studies have shown the significance of these lipids in plant cells viability and environmental response. This review focuses on phosphoinositides response to abiotic and biotic stress, showing their implication in plant survival during different stages of development. SIGNIFICANCE OF THE STUDY: This review is focused on plant PIPs functions in stress, highlighting in the main differences between plant and mammal PIPs and the novel interactions that could be extrapolated to animal models to contribute in a better understanding of these pivotal molecules.

Chromosome-level assembly of Arabidopsis thaliana Ler reveals the extent of translocation and inversion polymorphisms.

Resequencing or reference-based assemblies reveal large parts of the small-scale sequence variation. However, they typically fail to separate such local variation into colinear and rearranged variation, because they usually do not recover the complement of large-scale rearrangements, including transpositions and inversions. Besides the availability of hundreds of genomes of diverse Arabidopsis thaliana accessions, there is so far only one full-length assembled genome: the reference sequence. We have assembled 117 Mb of the A. thaliana Landsberg erecta (Ler) genome into five chromosome-equivalent sequences using a combination of short Illumina reads, long PacBio reads, and linkage information. Whole-genome comparison against the reference sequence revealed 564 transpositions and 47 inversions comprising ∼3.6 Mb, in addition to 4.1 Mb of nonreference sequence, mostly originating from duplications. Although rearranged regions are not different in local divergence from colinear regions, they are drastically depleted for meiotic recombination in heterozygotes. Using a 1.2-Mb inversion as an example, we show that such rearrangement-mediated reduction of meiotic recombination can lead to genetically isolated haplotypes in the worldwide population of A. thaliana Moreover, we found 105 single-copy genes, which were only present in the reference sequence or the Ler assembly, and 334 single-copy orthologs, which showed an additional copy in only one of the genomes. To our knowledge, this work gives first insights into the degree and type of variation, which will be revealed once complete assemblies will replace resequencing or other reference-dependent methods.

New insights into the phylogeny of the TMBIM superfamily across the tree of life: Comparative genomics and synteny networks reveal independent evolution of the BI and LFG families in plants.

The Transmembrane BAX Inhibitor Motif containing (TMBIM) superfamily, divided into BAX Inhibitor (BI) and Lifeguard (LFG) families, comprises a group of cytoprotective cell death regulators conserved in prokaryotes and eukaryotes. However, no research has focused on the evolution of this superfamily in plants. We identified 685 TMBIM proteins in 171 organisms from Archaea, Bacteria, and Eukarya, and provided a phylogenetic overview of the whole TMBIM superfamily. Then, we used orthology and synteny network analyses to further investigate the evolution and expansion of the BI and LFG families in 48 plants from diverse taxa. Plant BI family forms a single monophyletic group; however, monocot BI sequences transposed to another genomic context during evolution. Plant LFG family, which expanded trough whole genome and tandem duplications, is subdivided in LFG I, LFG IIA, and LFG IIB major phylogenetic groups, and retains synteny in angiosperms. Moreover, two orthologous groups (OGs) are shared between bryophytes and seed plants. Other several lineage-specific OGs are present in plants. This work clarifies the phylogenetic classification of the TMBIM superfamily across the three domains of life. Furthermore, it sheds new light on the evolution of the BI and LFG families in plants providing a benchmark for future research.