Next generation synthetic memory via intercepting recombinase function.

Here we present a technology to facilitate synthetic memory in a living system via repurposing Transcriptional Programming (i.e., our decision-making technology) parts, to regulate (intercept) recombinase function post-translation. We show that interception synthetic memory can facilitate programmable loss-of-function via site-specific deletion, programmable gain-of-function by way of site-specific inversion, and synthetic memory operations with nested Boolean logical operations. We can expand interception synthetic memory capacity more than 5-fold for a single recombinase, with reconfiguration specificity for multiple sites in parallel. Interception synthetic memory is ~10-times faster than previous generations of recombinase-based memory. We posit that the faster recombination speed of our next-generation memory technology is due to the post-translational regulation of recombinase function. This iteration of synthetic memory is complementary to decision-making via Transcriptional Programming - thus can be used to develop intelligent synthetic biological systems for myriad applications.

Magnetic resonance imaging features of invasive breast cancer association with the tumour stromal ratio.

OBJECTIVE: To assess the association between breast cancer tumour stroma and magnetic resonance imaging (MRI) features. MATERIALS AND METHODS: A total of 84 patients with treatment-naïve invasive breast cancer were enrolled into this retrospective study. The tumour stroma ratio (TSR) was estimated from the amount of tumour stroma in the pathology specimen of the breast tumour. The MRI images of the patients were analysed based on Breast Imaging Reporting and Data Systems (ACR-BIRADS) for qualitative features which include T2- weighted, diffusion-weighted images (DWI) and dynamic contrast-enhanced (DCE) for kinetic features. The mean signal intensity (SI) of Short Tau Inversion Recovery (STIR), with the ratio of STIR of the lesion and pectoralis muscle (L/M ratio) and apparent diffusion coefficient (ADC) value, were measured for the quantitative features. Correlation tests were performed to assess the relationship between TSR and MRI features. RESULTS: There was a significant correlation between the margin of mass, enhancement pattern, and STIR signal intensity of breast cancer and TSR. There were 54.76% (n = 46) in the low stromal group and 45.24% (n = 38) in the high stromal group. A significant association were seen between the margin of the mass and TSR (p = 0.034) between the L/M ratio (p <0.001), and between STIR SI of the lesion and TSR (p<0.001). The median L/M ratio was significantly higher in the high TSR group as compared to the lower TSR group (p < 0.001). CONCLUSION: Breast cancer with high stroma had spiculated margins, lower STIR signal intensity, and a heterogeneous pattern of enhancement. Hence, in this preliminary study, certain MRI features showed a potential to predict TSR.

Long-molecule scars of backup DNA repair in BRCA1- and BRCA2-deficient cancers.

Homologous recombination (HR) deficiency is associated with DNA rearrangements and cytogenetic aberrations1. Paradoxically, the types of DNA rearrangements that are specifically associated with HR-deficient cancers only minimally affect chromosomal structure2. Here, to address this apparent contradiction, we combined genome-graph analysis of short-read whole-genome sequencing (WGS) profiles across thousands of tumours with deep linked-read WGS of 46 BRCA1- or BRCA2-mutant breast cancers. These data revealed a distinct class of HR-deficiency-enriched rearrangements called reciprocal pairs. Linked-read WGS showed that reciprocal pairs with identical rearrangement orientations gave rise to one of two distinct chromosomal outcomes, distinguishable only with long-molecule data. Whereas one (cis) outcome corresponded to the copying and pasting of a small segment to a distant site, a second (trans) outcome was a quasi-balanced translocation or multi-megabase inversion with substantial (10 kb) duplications at each junction. We propose an HR-independent replication-restart repair mechanism to explain the full spectrum of reciprocal pair outcomes. Linked-read WGS also identified single-strand annealing as a repair pathway that is specific to BRCA2 deficiency in human cancers. Integrating these features in a classifier improved discrimination between BRCA1- and BRCA2-deficient genomes. In conclusion, our data reveal classes of rearrangements that are specific to BRCA1 or BRCA2 deficiency as a source of cytogenetic aberrations in HR-deficient cells.

Single-cell analysis of multiple invertible promoters reveals differential inversion rates as a strong determinant of bacterial population heterogeneity.

Population heterogeneity can promote bacterial fitness in response to unpredictable environmental conditions. A major mechanism of phenotypic variability in the human gut symbiont Bacteroides spp. involves the inversion of promoters that drive the expression of capsular polysaccharides, which determine the architecture of the cell surface. High-throughput single-cell sequencing reveals substantial population heterogeneity generated through combinatorial promoter inversion regulated by a broadly conserved serine recombinase. Exploiting control over population diversification, we show that populations with different initial compositions converge to a similar composition over time. Combining our data with stochastic computational modeling, we demonstrate that the differential rates of promoter inversion are a major mechanism shaping population dynamics. More broadly, our approach could be used to interrogate single-cell combinatorial phase variable states of diverse microbes including bacterial pathogens.

Earthquake focal mechanisms with distributed acoustic sensing.

Earthquake focal mechanisms provide critical in-situ insights about the subsurface faulting geometry and stress state. For frequent small earthquakes (magnitude< 3.5), their focal mechanisms are routinely determined using first-arrival polarities picked on the vertical component of seismometers. Nevertheless, their quality is usually limited by the azimuthal coverage of the local seismic network. The emerging distributed acoustic sensing (DAS) technology, which can convert pre-existing telecommunication cables into arrays of strain/strain-rate meters, can potentially fill the azimuthal gap and enhance constraints on the nodal plane orientation through its long sensing range and dense spatial sampling. However, determining first-arrival polarities on DAS is challenging due to its single-component sensing and low signal-to-noise ratio for direct body waves. Here, we present a data-driven method that measures P-wave polarities on a DAS array based on cross-correlations between earthquake pairs. We validate the inferred polarities using the regional network catalog on two DAS arrays, deployed in California and each comprising ~ 5000 channels. We demonstrate that a joint focal mechanism inversion combining conventional and DAS polarity picks improves the accuracy and reduces the uncertainty in the focal plane orientation. Our results highlight the significant potential of integrating DAS with conventional networks for investigating high-resolution earthquake source mechanisms.

Factors contributing to mitogenome size variation and a recurrent intracellular DNA transfer in Melastoma.

BACKGROUND: Mitogenome sizes of seed plants vary substantially even among closely related species, which are often related to horizontal or intracellular DNA transfer (HDT or IDT) events. However, the mechanisms of this size variation have not been well characterized. RESULTS: Here we assembled and characterized the mitogenomes of three species of Melastoma, a tropical shrub genus experiencing rapid speciation. The mitogenomes of M. candidum (Mc), M. sanguineum (Ms) and M. dodecandrum (Md) were assembled to a circular mapping chromosome of 391,595 bp, 395,542 bp and 412,026 bp, respectively. While the mitogenomes of Mc and Ms showed good collinearity except for a large inversion of ~ 150 kb, there were many rearrangements in the mitogenomes between Md and either Mc or Ms. Most non-alignable sequences (> 80%) between Mc and Ms are from gain or loss of mitochondrial sequences. Whereas, between Md and either Mc or Ms, non-alignable sequences in Md are mainly chloroplast derived sequences (> 30%) and from putative horizontal DNA transfers (> 30%), and those in both Mc and Ms are from gain or loss of mitochondrial sequences (> 80%). We also identified a recurrent IDT event in another congeneric species, M. penicillatum, which has not been fixed as it is only found in one of the three examined populations. CONCLUSIONS: By characterizing mitochondrial genome sequences of Melastoma, our study not only helps understand mitogenome size evolution in closely related species, but also cautions different evolutionary histories of mitochondrial regions due to potential recurrent IDT events in some populations or species.

An algebraic model for inversion and deletion in bacterial genome rearrangement.

Inversions, also sometimes called reversals, are a major contributor to variation among bacterial genomes, with studies suggesting that those involving small numbers of regions are more likely than larger inversions. Deletions may arise in bacterial genomes through the same biological mechanism as inversions, and hence a model that incorporates both is desirable. However, while inversion distances between genomes have been well studied, there has yet to be a model which accounts for the combination of both deletions and inversions. To account for both of these operations, we introduce an algebraic model that utilises partial permutations. This leads to an algorithm for calculating the minimum distance to the most recent common ancestor of two bacterial genomes evolving by inversions (of adjacent regions) and deletions. The algebraic model makes the existing short inversion models more complete and realistic by including deletions, and also introduces new algebraic tools into evolutionary distance problems.

Chromosomal inversion polymorphisms shape human brain morphology.

The impact of chromosomal inversions on human brain morphology remains underexplored. We studied 35 common inversions classified from genotypes of 33,018 adults with European ancestry. The inversions at 2p22.3, 16p11.2, and 17q21.31 reach genome-wide significance, followed by 8p23.1 and 6p21.33, in their association with cortical and subcortical morphology. The 17q21.31, 8p23.1, and 16p11.2 regions comprise the LRRC37, OR7E, and NPIP duplicated gene families. We find the 17q21.31 MAPT inversion region, known for harboring neurological risk, to be the most salient locus among common variants for shaping and patterning the cortex. Overall, we observe the inverted orientations decreasing brain size, with the exception that the 2p22.3 inversion is associated with increased subcortical volume and the 8p23.1 inversion is associated with increased motor cortex. These significant inversions are in the genomic hotspots of neuropsychiatric loci. Our findings are generalizable to 3,472 children and demonstrate inversions as essential genetic variation to understand human brain phenotypes.

Recent Evolution of a Maternally Acting Sex-Determining Supergene in a Fly with Single-Sex Broods.

Sex determination is a key developmental process, yet it is remarkably variable across the tree of life. The dipteran family Sciaridae exhibits one of the most unusual sex determination systems in which mothers control offspring sex through selective elimination of paternal X chromosomes. Whereas in some members of the family females produce mixed-sex broods, others such as the dark-winged fungus gnat Bradysia coprophila are monogenic, with females producing single-sex broods. Female-producing females were previously found to be heterozygous for a large X-linked paracentric inversion (X'), which is maternally inherited and absent from male-producing females. Here, we assembled and characterized the X' sequence. As close sequence homology between the X and X' made identification of the inversion challenging, we developed a k-mer-based approach to bin genomic reads before assembly. We confirmed that the inversion spans most of the X' chromosome (∼55 Mb) and encodes ∼3,500 genes. Analysis of the divergence between the inversion and the homologous region of the X revealed that it originated very recently (<0.5 Ma). Surprisingly, we found that the X' is more complex than previously thought and is likely to have undergone multiple rearrangements that have produced regions of varying ages, resembling a supergene composed of evolutionary strata. We found functional degradation of ∼7.3% of genes within the region of recombination suppression, but no evidence of accumulation of repetitive elements. Our findings provide an indication that sex-linked inversions are driving turnover of the strange sex determination system in this family of flies.

Complex evolutionary processes maintain an ancient chromosomal inversion.

Genome re-arrangements such as chromosomal inversions are often involved in adaptation. As such, they experience natural selection, which can erode genetic variation. Thus, whether and how inversions can remain polymorphic for extended periods of time remains debated. Here we combine genomics, experiments, and evolutionary modeling to elucidate the processes maintaining an inversion polymorphism associated with the use of a challenging host plant (Redwood trees) in Timema stick insects. We show that the inversion is maintained by a combination of processes, finding roles for life-history trade-offs, heterozygote advantage, local adaptation to different hosts, and gene flow. We use models to show how such multi-layered regimes of balancing selection and gene flow provide resilience to help buffer populations against the loss of genetic variation, maintaining the potential for future evolution. We further show that the inversion polymorphism has persisted for millions of years and is not a result of recent introgression. We thus find that rather than being a nuisance, the complex interplay of evolutionary processes provides a mechanism for the long-term maintenance of genetic variation.

Recombination and sterility in inversion homo- and heterokaryotypes under a general counting model of chiasma interference.

It has long been known that the chiasmata are not independently distributed in most organisms, a phenomenon known as chiasma interference. In this paper, I suggest a model of chiasma interference that generalizes the Poisson model, the counting model, the Poisson-skip model, and the two-pathway counting model into a single framework, and use it to derive infinite series expressions for the sterility and recombination pattern probabilities in inversion homo- and heterokaryotypes, and a closed-form expression for the special case of the two-pathway counting model in homokaryotypes. I then use these expressions to perform maximum likelihood parameter estimations for recombination and tetrad data from various species. The results imply that the simpler counting models perform well compared to more complex ones, that interference works in a similar way in homo- and heterokaryotypes, and that the model fits well with data for the latter as well as the former. I also find evidence that the interference signal is broken by the centromere in some species, but not others, suggestions of negative interference in Aspergillus nidulans, and no consistent support for the theory that a second noninterfering chiasma pathway exists only in organisms that require double-strand break for synapsis. I suggest that the latter finding is at least partly due to issues involved in analyzing aggregate data from different experiments and individuals.

The effects of inversion polymorphisms on patterns of neutral genetic diversity.

The strong reduction in the frequency of recombination in heterozygotes for an inversion and a standard gene arrangement causes the arrangements to become partially isolated genetically, resulting in sequence divergence between them and changes in the levels of neutral variability at nucleotide sites within each arrangement class. Previous theoretical studies on the effects of inversions on neutral variability have assumed either that the population is panmictic or that it is divided into 2 populations subject to divergent selection. Here, the theory is extended to a model of an arbitrary number of demes connected by migration, using a finite island model with the inversion present at the same frequency in all demes. Recursion relations for mean pairwise coalescent times are used to obtain simple approximate expressions for diversity and divergence statistics for an inversion polymorphism at equilibrium under recombination and drift, and for the approach to equilibrium following the sweep of an inversion to a stable intermediate frequency. The effects of an inversion polymorphism on patterns of linkage disequilibrium are also examined. The reduction in effective recombination rate caused by population subdivision can have significant effects on these statistics. The theoretical results are discussed in relation to population genomic data on inversion polymorphisms, with an emphasis on Drosophila melanogaster. Methods are proposed for testing whether or not inversions are close to recombination-drift equilibrium, and for estimating the rate of recombinational exchange in heterozygotes for inversions; difficulties involved in estimating the ages of inversions are also discussed.

Long-read genome sequencing identifies cryptic structural variants in congenital aniridia cases.

BACKGROUND: Haploinsufficiency of the transcription factor PAX6 is the main cause of congenital aniridia, a genetic disorder characterized by iris and foveal hypoplasia. 11p13 microdeletions altering PAX6 or its downstream regulatory region (DRR) are present in about 25% of patients; however, only a few complex rearrangements have been described to date. Here, we performed nanopore-based whole-genome sequencing to assess the presence of cryptic structural variants (SVs) on the only two unsolved "PAX6-negative" cases from a cohort of 110 patients with congenital aniridia after unsuccessfully short-read sequencing approaches. RESULTS: Long-read sequencing (LRS) unveiled balanced chromosomal rearrangements affecting the PAX6 locus at 11p13 in these two patients and allowed nucleotide-level breakpoint analysis. First, we identified a cryptic 4.9 Mb de novo inversion disrupting intron 7 of PAX6, further verified by targeted polymerase chain reaction amplification and sequencing and FISH-based cytogenetic analysis. Furthermore, LRS was decisive in correctly mapping a t(6;11) balanced translocation cytogenetically detected in a second proband with congenital aniridia and considered non-causal 15 years ago. LRS resolved that the breakpoint on chromosome 11 was indeed located at 11p13, disrupting the DNase I hypersensitive site 2 enhancer within the DRR of PAX6, 161 Kb from the causal gene. Patient-derived RNA expression analysis demonstrated PAX6 haploinsufficiency, thus supporting that the 11p13 breakpoint led to a positional effect by cleaving crucial enhancers for PAX6 transactivation. LRS analysis was also critical for mapping the exact breakpoint on chromosome 6 to the highly repetitive centromeric region at 6p11.1. CONCLUSIONS: In both cases, the LRS-based identified SVs have been deemed the hidden pathogenic cause of congenital aniridia. Our study underscores the limitations of traditional short-read sequencing in uncovering pathogenic SVs affecting low-complexity regions of the genome and the value of LRS in providing insight into hidden sources of variation in rare genetic diseases.

Quantum oscillations in field-induced correlated insulators of a moiré superlattice.

We report an observation of quantum oscillations (QOs) in the correlated insulators with valley anisotropy of twisted double bilayer graphene (TDBG). The anomalous QOs are best captured in the magneto resistivity oscillations of the insulators at v = -2, with a period of 1/B and an oscillation amplitude as high as 150 kΩ. The QOs can survive up to ∼10 K, and above 12 K, the insulating behaviors are dominant. The QOs of the insulator are strongly D dependent: the carrier density extracted from the 1/B periodicity decreases almost linearly with D from -0.7 to -1.1 V/nm, suggesting a reduced Fermi surface; the effective mass from Lifshitz-Kosevich analysis depends nonlinearly on D, reaching a minimal value of 0.1 me at D = ∼ -1.0 V/nm. Similar observations of QOs are also found at v = 2, as well as in other devices without graphite gate. We interpret the D sensitive QOs of the correlated insulators in the picture of band inversion. By reconstructing an inverted band model with the measured effective mass and Fermi surface, the density of state at the gap, calculated from thermal broadened Landau levels, agrees qualitatively with the observed QOs in the insulators. While more theoretical understandings are needed in the future to fully account for the anomalous QOs in this moiré system, our study suggests that TDBG is an excellent platform to discover exotic phases where correlation and topology are at play.

The molecular mechanisms of recombinant chromosome 18 with parental pericentric inversions and a review of the literature.

Chromosomal rearrangements mostly result from non-allelic homologous recombination mediated by low-copy repeats (LCRs) or segmental duplications (SDs). Recent studies on recombinant chromosome 18 (rec (18)) have focused on diagnoses and clinical phenotypes. We diagnosed two cases of prenatal rec (18) and identified precise breakpoint intervals using karyotype and chromosomal microarray analyses. We analyzed the distribution characteristics of breakpoint repetitive elements to infer rearrangement mechanisms and reviewed relevant literature to identify genetic trends. Among the 12 families with 25 pregnancies analyzed, 68% rec (18), 24% spontaneous abortions, and 8% normal births were reported. In the 17 rec (18) cases, 65% presented maternal origin and 35% were paternal. Short-arm breakpoints at p11.31 were reported in 10 cases, whereas the long-arm breakpoints were located at q21.3 (6 cases) and q12 (4 cases). Breakpoints of pericentric inversions on chromosome 18 are concentrated in p11.31, q21.3, and q12 regions. Rearrangements at 18p11.31 are non-recurrent events. ALUs, LINE1s, and MIRs were enriched at the breakpoint regions (1.85 to 3.42-fold enrichment over the entire chromosome 18), while SDs and LCRs were absent. ALU subfamilies had sequence identities of 85.94% and 83.01% between two pair breakpoints. Small repetitive elements may mediate recombination-coupled DNA repair processes, facilitating rearrangements on chromosome 18. Maternal inversion carriers are more prone to abnormal recombination in prenatal families with rec (18). Recombinant chromosomes may present preferential segregation during gamete formation.

Genomic, functional and structural analyses elucidate evolutionary innovation within the sea anemone 8 toxin family.

BACKGROUND: The ShK toxin from Stichodactyla helianthus has established the therapeutic potential of sea anemone venom peptides, but many lineage-specific toxin families in Actiniarians remain uncharacterised. One such peptide family, sea anemone 8 (SA8), is present in all five sea anemone superfamilies. We explored the genomic arrangement and evolution of the SA8 gene family in Actinia tenebrosa and Telmatactis stephensoni, characterised the expression patterns of SA8 sequences, and examined the structure and function of SA8 from the venom of T. stephensoni. RESULTS: We identified ten SA8-family genes in two clusters and six SA8-family genes in five clusters for T. stephensoni and A. tenebrosa, respectively. Nine SA8 T. stephensoni genes were found in a single cluster, and an SA8 peptide encoded by an inverted SA8 gene from this cluster was recruited to venom. We show that SA8 genes in both species are expressed in a tissue-specific manner and the inverted SA8 gene has a unique tissue distribution. While the functional activity of the SA8 putative toxin encoded by the inverted gene was inconclusive, its tissue localisation is similar to toxins used for predator deterrence. We demonstrate that, although mature SA8 putative toxins have similar cysteine spacing to ShK, SA8 peptides are distinct from ShK peptides based on structure and disulfide connectivity. CONCLUSIONS: Our results provide the first demonstration that SA8 is a unique gene family in Actiniarians, evolving through a variety of structural changes including tandem and proximal gene duplication and an inversion event that together allowed SA8 to be recruited into the venom of T. stephensoni.

[Optical genome mapping analysis of a Chinese pedigree with a rare chromosome 17 paracentric inversion insertion].

OBJECTIVE: To carry out optical genome mapping (OGM) for a Chinese pedigree with a rare paracentric reverse insertion of chromosome 17. METHODS: A high-risk pregnant woman identified at the Prenatal Diagnosis Center of Hangzhou Women's Hospital in October 2021 and her family members were selected as the study subjects. Chromosome G banding analysis, fluorescence in situ hybridization (FISH), single nucleotide polymorphism array (SNP array) and OGM were applied to verify the balanced structural abnormality of chromosome 17 in the pedigree. RESULTS: Chromosomal karyotyping analysis and SNP array assay have identified a duplication of 17q23q25 in the fetus. Karyotyping analysis of the pregnant woman showed that the structure of chromosome 17 was abnormal, whilst SNP array has detected no abnormality. OGM revealed that the woman has carried a paracentric reverse insertion, which was confirmed by FISH. The karyotype of her husband was normal. CONCLUSION: The duplication of 17q23q25 in the fetus has derived from a paracentric reverse insertion of chromosome 17 in its mother. OGM has the advantage for delineating balanced chromosome structural abnormalities.

Deep learning models challenge the prevailing assumption that face-like effects for objects of expertise support domain-general mechanisms.

The question of whether task performance is best achieved by domain-specific, or domain-general processing mechanisms is fundemental for both artificial and biological systems. This question has generated a fierce debate in the study of expert object recognition. Because humans are experts in face recognition, face-like neural and cognitive effects for objects of expertise were considered support for domain-general mechanisms. However, effects of domain, experience and level of categorization, are confounded in human studies, which may lead to erroneous inferences. To overcome these limitations, we trained deep learning algorithms on different domains (objects, faces, birds) and levels of categorization (basic, sub-ordinate, individual), matched for amount of experience. Like humans, the models generated a larger inversion effect for faces than for objects. Importantly, a face-like inversion effect was found for individual-based categorization of non-faces (birds) but only in a network specialized for that domain. Thus, contrary to prevalent assumptions, face-like effects for objects of expertise do not support domain-general mechanisms but may originate from domain-specific mechanisms. More generally, we show how deep learning algorithms can be used to dissociate factors that are inherently confounded in the natural environment of biological organisms to test hypotheses about their isolated contributions to cognition and behaviour.

An Ancestral Balanced Inversion Polymorphism Confers Global Adaptation.

Since the pioneering work of Dobzhansky in the 1930s and 1940s, many chromosomal inversions have been identified, but how they contribute to adaptation remains poorly understood. In Drosophila melanogaster, the widespread inversion polymorphism In(3R)Payne underpins latitudinal clines in fitness traits on multiple continents. Here, we use single-individual whole-genome sequencing, transcriptomics, and published sequencing data to study the population genomics of this inversion on four continents: in its ancestral African range and in derived populations in Europe, North America, and Australia. Our results confirm that this inversion originated in sub-Saharan Africa and subsequently became cosmopolitan; we observe marked monophyletic divergence of inverted and noninverted karyotypes, with some substructure among inverted chromosomes between continents. Despite divergent evolution of this inversion since its out-of-Africa migration, derived non-African populations exhibit similar patterns of long-range linkage disequilibrium between the inversion breakpoints and major peaks of divergence in its center, consistent with balancing selection and suggesting that the inversion harbors alleles that are maintained by selection on several continents. Using RNA-sequencing, we identify overlap between inversion-linked single-nucleotide polymorphisms and loci that are differentially expressed between inverted and noninverted chromosomes. Expression levels are higher for inverted chromosomes at low temperature, suggesting loss of buffering or compensatory plasticity and consistent with higher inversion frequency in warm climates. Our results suggest that this ancestrally tropical balanced polymorphism spread around the world and became latitudinally assorted along similar but independent climatic gradients, always being frequent in subtropical/tropical areas but rare or absent in temperate climates.