Recurrent inversion polymorphisms in humans associate with genetic instability and genomic disorders.

Unlike copy number variants (CNVs), inversions remain an underexplored genetic variation class. By integrating multiple genomic technologies, we discover 729 inversions in 41 human genomes. Approximately 85% of inversions <2 kbp form by twin-priming during L1 retrotransposition; 80% of the larger inversions are balanced and affect twice as many nucleotides as CNVs. Balanced inversions show an excess of common variants, and 72% are flanked by segmental duplications (SDs) or retrotransposons. Since flanking repeats promote non-allelic homologous recombination, we developed complementary approaches to identify recurrent inversion formation. We describe 40 recurrent inversions encompassing 0.6% of the genome, showing inversion rates up to 2.7 × 10-4 per locus per generation. Recurrent inversions exhibit a sex-chromosomal bias and co-localize with genomic disorder critical regions. We propose that inversion recurrence results in an elevated number of heterozygous carriers and structural SD diversity, which increases mutability in the population and predisposes specific haplotypes to disease-causing CNVs.

Interfaces of Malignant and Immunologic Clonal Dynamics in Ovarian Cancer.

High-grade serous ovarian cancer (HGSC) exhibits extensive malignant clonal diversity with widespread but non-random patterns of disease dissemination. We investigated whether local immune microenvironment factors shape tumor progression properties at the interface of tumor-infiltrating lymphocytes (TILs) and cancer cells. Through multi-region study of 212 samples from 38 patients with whole-genome sequencing, immunohistochemistry, histologic image analysis, gene expression profiling, and T and B cell receptor sequencing, we identified three immunologic subtypes across samples and extensive within-patient diversity. Epithelial CD8+ TILs negatively associated with malignant diversity, reflecting immunological pruning of tumor clones inferred by neoantigen depletion, HLA I loss of heterozygosity, and spatial tracking between T cell and tumor clones. In addition, combinatorial prognostic effects of mutational processes and immune properties were observed, illuminating how specific genomic aberration types associate with immune response and impact survival. We conclude that within-patient spatial immune microenvironment variation shapes intraperitoneal malignant spread, provoking new evolutionary perspectives on HGSC clonal dispersion.

Why do sex chromosomes progressively lose recombination?

Progressive recombination loss is a common feature of sex chromosomes. Yet, the evolutionary drivers of this phenomenon remain a mystery. For decades, differences in trait optima between sexes (sexual antagonism) have been the favoured hypothesis, but convincing evidence is lacking. Recent years have seen a surge of alternative hypotheses to explain progressive extensions and maintenance of recombination suppression: neutral accumulation of sequence divergence, selection of nonrecombining fragments with fewer deleterious mutations than average, sheltering of recessive deleterious mutations by linkage to heterozygous alleles, early evolution of dosage compensation, and constraints on recombination restoration. Here, we explain these recent hypotheses and dissect their assumptions, mechanisms, and predictions. We also review empirical studies that have brought support to the various hypotheses.

The impact of inversions across 33,924 families with rare disease from a national genome sequencing project.

Detection of structural variants (SVs) is currently biased toward those that alter copy number. The relative contribution of inversions toward genetic disease is unclear. In this study, we analyzed genome sequencing data for 33,924 families with rare disease from the 100,000 Genomes Project. From a database hosting >500 million SVs, we focused on 351 genes where haploinsufficiency is a confirmed disease mechanism and identified 47 ultra-rare rearrangements that included an inversion (24 bp to 36.4 Mb, 20/47 de novo). Validation utilized a number of orthogonal approaches, including retrospective exome analysis. RNA-seq data supported the respective diagnoses for six participants. Phenotypic blending was apparent in four probands. Diagnostic odysseys were a common theme (>50 years for one individual), and targeted analysis for the specific gene had already been performed for 30% of these individuals but with no findings. We provide formal confirmation of a European founder origin for an intragenic MSH2 inversion. For two individuals with complex SVs involving the MECP2 mutational hotspot, ambiguous SV structures were resolved using long-read sequencing, influencing clinical interpretation. A de novo inversion of HOXD11-13 was uncovered in a family with Kantaputra-type mesomelic dysplasia. Lastly, a complex translocation disrupting APC and involving nine rearranged segments confirmed a clinical diagnosis for three family members and resolved a conundrum for a sibling with a single polyp. Overall, inversions play a small but notable role in rare disease, likely explaining the etiology in around 1/750 families across heterogeneous clinical cohorts.

Constraining Earth's nonlinear mantle viscosity using plate-boundary resolving global inversions.

Variable viscosity in Earth's mantle exerts a fundamental control on mantle convection and plate tectonics, yet rigorously constraining the underlying parameters has remained a challenge. Inverse methods have not been sufficiently robust to handle the severe viscosity gradients and nonlinearities (arising from dislocation creep and plastic failure) while simultaneously resolving the megathrust and bending slabs globally. Using global plate motions as constraints, we overcome these challenges by combining a scalable nonlinear Stokes solver that resolves the key tectonic features with an adjoint-based Bayesian approach. Assuming plate cooling, variations in the thickness of continental lithosphere, slabs, and broad scale lower mantle structure as well as a constant grain size through the bulk of the upper mantle, a good fit to global plate motions is found with a nonlinear upper mantle stress exponent of 2.43 [Formula: see text] 0.25 (mean [Formula: see text] SD). A relatively low yield stress of 151 [Formula: see text] 19 MPa is required for slabs to bend during subduction and transmit a slab pull that generates asymmetrical subduction. The recovered long-term strength of megathrusts (plate interfaces) varies between different subduction zones, with South America having a larger strength and Vanuatu and Central America having lower values with important implications for the stresses driving megathrust earthquakes.

Observation-derived 2010-2019 trends in methane emissions and intensities from US oil and gas fields tied to activity metrics.

The United States is the world's largest oil/gas methane emitter according to current national reports. Reducing these emissions is a top priority in the US government's climate action plan. Here, we use a 2010 to 2019 high-resolution inversion of surface and satellite observations of atmospheric methane to quantify emission trends for individual oil/gas production regions in North America and relate them to production and infrastructure. We estimate a mean US oil/gas methane emission of 14.8 (12.4 to 16.5) Tg a-1 for 2010 to 2019, 70% higher than reported by the US Environmental Protection Agency. While emissions in Canada and Mexico decreased over the period, US emissions increased from 2010 to 2014, decreased until 2017, and rose again afterward. Increases were driven by the largest production regions (Permian, Anadarko, Marcellus), while emissions in the smaller production regions generally decreased. Much of the year-to-year emission variability can be explained by oil/gas production rates, active well counts, and new wells drilled, with the 2014 to 2017 decrease driven by reduction in new wells and the 2017 to 2019 surge driven by upswing of production. We find a steady decrease in the oil/gas methane intensity (emission per unit methane gas production) for almost all major US production regions. The mean US methane intensity decreased from 3.7% in 2010 to 2.5% in 2019. If the methane intensity for the oil/gas supply chain continues to decrease at this pace, we may expect a 32% decrease in US oil/gas emissions by 2030 despite projected increases in production.

A novel application of data-consistent inversion to overcome spurious inference in genome-wide association studies.

The genome-wide association studies (GWAS) typically use linear or logistic regression models to identify associations between phenotypes (traits) and genotypes (genetic variants) of interest. However, the use of regression with the additive assumption has potential limitations. First, the normality assumption of residuals is the one that is rarely seen in practice, and deviation from normality increases the Type-I error rate. Second, building a model based on such an assumption ignores genetic structures, like, dominant, recessive, and protective-risk cases. Ignoring genetic variants may result in spurious conclusions about the associations between a variant and a trait. We propose an assumption-free model built upon data-consistent inversion (DCI), which is a recently developed measure-theoretic framework utilized for uncertainty quantification. This proposed DCI-derived model builds a nonparametric distribution on model inputs that propagates to the distribution of observed data without the required normality assumption of residuals in the regression model. This characteristic enables the proposed DCI-derived model to cover all genetic variants without emphasizing on additivity of the classic-GWAS model. Simulations and a replication GWAS with data from the COPDGene demonstrate the ability of this model to control the Type-I error rate at least as well as the classic-GWAS (additive linear model) approach while having similar or greater power to discover variants in different genetic modes of transmission.

Strong methane point sources contribute a disproportionate fraction of total emissions across multiple basins in the United States.

Understanding, prioritizing, and mitigating methane (CH4) emissions requires quantifying CH4 budgets from facility scales to regional scales with the ability to differentiate between source sectors. We deployed a tiered observing system for multiple basins in the United States (San Joaquin Valley, Uinta, Denver-Julesburg, Permian, Marcellus). We quantify strong point source emissions (>10 kg CH4 h-1) using airborne imaging spectrometers, attribute them to sectors, and assess their intermittency with multiple revisits. We compare these point source emissions to total basin CH4 fluxes derived from inversion of Sentinel-5p satellite CH4 observations. Across basins, point sources make up on average 40% of the regional flux. We sampled some basins several times across multiple months and years and find a distinct bimodal structure to emission timescales: the total point source budget is split nearly in half by short-lasting and long-lasting emission events. With the increasing airborne and satellite observing capabilities planned for the near future, tiered observing systems will more fully quantify and attribute CH4 emissions from facility to regional scales, which is needed to effectively and efficiently reduce methane emissions.

Neo-sex chromosome evolution shapes sex-dependent asymmetrical introgression barrier.

Sex chromosomes play a special role in the evolution of reproductive barriers between species. Here we describe conflicting roles of nascent sex chromosomes on patterns of introgression in an experimental hybrid swarm. Drosophila nasuta and Drosophila albomicans are recently diverged, fully fertile sister species that have different sex chromosome systems. The fusion between an autosome (Muller CD) with the ancestral X and Y gave rise to neo-sex chromosomes in D. albomicans, while Muller CD remains unfused in D. nasuta. We found that a large block containing overlapping inversions on the neo-sex chromosome stood out as the strongest barrier to introgression. Intriguingly, the neo-sex chromosome introgression barrier is asymmetrical and sex-dependent. Female hybrids showed significant D. albomicans­biased introgression on Muller CD (neo-X excess), while males showed heterosis with excessive (neo-X, D. nasuta Muller CD) genotypes. We used a population genetic model to dissect the interplay of sex chromosome drive, heterospecific pairing incompatibility between the neo-sex chromosomes and unfused Muller CD, neo-Y disadvantage, and neo-X advantage in generating the observed sex chromosome genotypes in females and males. We show that moderate neo-Y disadvantage and D. albomicans specific meiotic drive are required to observe female-specific D. albomicans­biased introgression in this system, together with pairing incompatibility and neo-X advantage. In conclusion, this hybrid swarm between a young species pair sheds light onto the multifaceted roles of neo-sex chromosomes in a sex-dependent asymmetrical introgression barrier at a species boundary.

Transient chirality inversion during racemization of a helical cobalt(III) complex.

SignificanceWe first observed a transient chirality inversion on a simple unimolecular platform during the racemization of a chiral helical complex [LCo3A6]3+, i.e., the helicity changed from P-rich (right-handed) to M-rich (left-handed), which then racemized to a P/M equimolar mixture in spite of the absence of a reagent that could induce the M helix. This transient chirality inversion was observed only in the forward reaction, whereas the reverse reaction showed a simple monotonic change with an induction time. Consequently, the M helicity appeared only in the forward reaction. These forward and reverse reactions constitute a hysteretic cycle. Compounds showing such unique time responses would be useful for developing time-programmable switchable materials that can control the physical/chemical properties in a time-dependent manner.

Controlling inversion disorder in a stoichiometric spinel magnet.

In the study of frustrated quantum magnets, it is essential to be able to control the nature and degree of site disorder during the growth process, as many measurement techniques are incapable of distinguishing between site disorder and frustration-induced spin disorder. Pyrochlore-structured spinel oxides can serve as model systems of geometrically frustrated three-dimensional quantum magnets; however, the nature of the magnetism in one well-studied spinel, ZnFe2O4, remains unclear. Here, we demonstrate simultaneous control of both stoichiometry and inversion disorder in the growth of ZnFe2O4 single crystals, directly yielding a revised understanding of both the collective spin behavior and lattice symmetry. Crystals grown in the stoichiometric limit with minimal site inversion disorder contravene all the previously suggested exotic spin phases in ZnFe2O4. Furthermore, the structure is confirmed on the [Formula: see text] space group with broken inversion symmetry that induces antiferroelectricity. The effective tuning of magnetic behavior by site disorder in the presence of robust antiferroelectricity makes ZnFe2O4 of special interest to multiferroic devices.

Unprecedented variation pattern of plastid genomes and the potential role in adaptive evolution in Poales.

BACKGROUND: The plastid is the photosynthetic organelle in plant cell, and the plastid genomes (plastomes) are generally conserved in evolution. As one of the most economically and ecologically important order of angiosperms, Poales was previously documented to exhibit great plastomic variation as an order of photoautotrophic plants. RESULTS: We acquired 93 plastomes, representing all the 16 families and 5 major clades of Poales to reveal the extent of their variation and evolutionary pattern. Extensive variation including the largest one in monocots with 225,293 bp in size, heterogeneous GC content, and a wide variety of gene duplication and loss were revealed. Moreover, rare occurrences of three inverted repeat (IR) copies in angiosperms and one IR loss were observed, accompanied by short IR (sIR) and small direct repeat (DR). Widespread structural heteroplasmy, diversified inversions, and unusual genomic rearrangements all appeared in Poales, occasionally within a single species. Extensive repeats in the plastomes were found to be positively correlated with the observed inversions and rearrangements. The variation all showed a "small-large-moderate" trend along the evolution of Poales, as well as for the sequence substitution rate. Finally, we found some positively selected genes, mainly in C4 lineages, while the closely related lineages of those experiencing gene loss tended to have undergone more relaxed purifying selection. CONCLUSIONS: The variation of plastomes in Poales may be related to its successful diversification into diverse habitats and multiple photosynthetic pathway transitions. Our order-scale analyses revealed unusual evolutionary scenarios for plastomes in the photoautotrophic order of Poales and provided new insights into the plastome evolution in angiosperms as a whole.

Unconventional Superconducting Diode Effects via Antisymmetry and Antisymmetry Breaking.

Symmetry breaking plays a pivotal role in unlocking intriguing properties and functionalities in material systems. For example, the breaking of spatial and temporal symmetries leads to a fascinating phenomenon: the superconducting diode effect. However, generating and precisely controlling the superconducting diode effect pose significant challenges. Here, we take a novel route with the deliberate manipulation of magnetic charge potentials to realize unconventional superconducting flux-quantum diode effects. We achieve this through suitably tailored nanoengineered arrays of nanobar magnets on top of a superconducting thin film. We demonstrate the vital roles of inversion antisymmetry and its breaking in evoking unconventional superconducting effects, namely a magnetically symmetric diode effect and an odd-parity magnetotransport effect. These effects are nonvolatilely controllable through in situ magnetization switching of the nanobar magnets. Our findings promote the use of antisymmetry (breaking) for initiating unconventional superconducting properties, paving the way for exciting prospects and innovative functionalities in superconducting electronics.

Anatomy of Hidden Dzyaloshinskii-Moriya Interactions and Topological Spin Textures in Centrosymmetric Crystals.

The Dzyaloshinskii-Moriya interaction (DMI) is understood to be forbidden by the symmetry of centrosymmetric systems, thus restricting the candidate types for investigating many correlated physical phenomena. Here, we report the hidden DMI existing in centrosymmetric magnets driven by the local inversion symmetry breaking of specific spin sublattices. The opposite DMI spatially localized on the inverse spin sublattice favors the separated spin spiral with opposite chirality. Furthermore, we elucidate that hidden DMI widely exists in many potential candidates, from the first-principles calculations on the mature crystal database. Interestingly, novel topological spin configurations, such as the anti-chirality-locked merons and antiferromagnetic-ferromagnetic meron chains, are stabilized as a consequence of hidden DMI. Our understanding enables the effective control of DMI by symmetry operations at the atomic level and enlarges the range of currently useful magnets for topological magnetism.

Electric Tuning of Vortex Ratchet Effect in NbSe2.

Inversion symmetry breaking has played an important role in recent discoveries of nonreciprocal charge transport. Niobium diselenide, for example, lacks an inversion center in the monolayer form and can host prominent nonreciprocal transport property. Here, however, we observe a nonreciprocal transport signal in the second-harmonic channel of bulk-like NbSe2, in which inversion symmetry of the lattice seems preserved. The second-harmonic signal occurs along different in-plane current orientations and appears not only in the vortex-liquid regime but also even in the superconducting fluctuation regime without an applied magnetic field. By adding a direct current (DC) bias, we quantify the symmetry breaking effect in the vortex-liquid regime. The DC bias also suggests that the rectification effect at the contacts may account for the seemingly nonreciprocal transport at zero magnetic field. Our results demonstrate that DC biasing is a useful knob for addressing nonreciprocal charge transport in a wide range of materials.

High-resolution sequencing of nine elite upland cotton cultivars uncovers genic variations and breeding improvement targets.

Structural variations (SVs) are critical factors affecting genome evolution and important traits. However, identification results and functional analyses of SVs in upland cotton are rare. Here, based on the genetic relationships, breeding history and cumulative planting area of upland cotton in China, nine predominant cultivars from the past 60 years (1950s-2010s) were selected for long read sequencing to uncover genic variations and breeding improvement targets for this crop. Based on the ZM24 reference genome, 0.88-1.47 × 104 SVs per cultivar were identified, and an SV set was constructed. SVs affected the expression of a large number of genes during fiber elongation, and a transposable element insertion resulted in the glandless phenotype in upland cotton. Six widespread inversions were identified based on nine draft genomes and high-throughput chromosome conformation capture data. Multiple haplotype blocks that were always associated with aggregated SVs were demonstrated to play a pivotal role in the agronomic traits of upland cotton and drove its adaptation to the northern planting region. Exotic introgression was the source of these haplotype blocks and increased the genetic diversity of upland cotton. Our results enrich the genome resources of upland cotton, and the identified SVs will promote genetic and breeding research in cotton.

A comparative plastome approach enhances the assessment of genetic variation in the Melilotus genus.

BACKGROUND: Melilotus, a member of the Fabaceae family, is a pivotal forage crop that is extensively cultivated in livestock regions globally due to its notable productivity and ability to withstand abiotic stress. However, the genetic attributes of the chloroplast genome and the evolutionary connections among different Melilotus species remain unresolved. RESULTS: In this study, we compiled the chloroplast genomes of 18 Melilotus species and performed a comprehensive comparative analysis. Through the examination of protein-coding genes, we successfully established a robust phylogenetic tree for these species. This conclusion is further supported by the phylogeny derived from single-nucleotide polymorphisms (SNPs) across the entire chloroplast genome. Notably, our findings revealed that M. infestus, M. siculus, M. sulcatus, and M. speciosus formed a distinct subgroup within the phylogenetic tree. Additionally, the chloroplast genomes of these four species exhibit two shared inversions. Moreover, inverted repeats were observed to have reemerged in six species within the IRLC. The distribution patterns of single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) within protein-coding genes indicated that ycf1 and ycf2 accumulated nonconservative alterations during evolutionary development. Furthermore, an examination of the evolutionary rate of protein-coding genes revealed that rps18, rps7, and rpl16 underwent positive selection specifically in Melilotus. CONCLUSIONS: We present a comparative analysis of the complete chloroplast genomes of Melilotus species. This study represents the most thorough and detailed exploration of the evolution and variability within the genus Melilotus to date. Our study provides valuable chloroplast genomic information for improving phylogenetic reconstructions and making biogeographic inferences about Melilotus and other Papilionoideae species.

Medial positioning of the hippocampus and hippocampal fissure volume in developmental topographical disorientation.

Developmental topographical disorientation (DTD) refers to the lifelong inability to orient by means of cognitive maps in familiar surroundings despite otherwise well-preserved general cognitive functions, and the absence of any acquired brain injury or neurological condition. While reduced functional connectivity between the hippocampus and other brain regions has been reported in DTD individuals, no structural differences in gray matter tissue for the whole brain neither for the hippocampus were detected. Considering that the human hippocampus is the main structure associated with cognitive map-based navigation, here, we investigated differences in morphological and morphometric hippocampal features between individuals affected by DTD (N = 20) and healthy controls (N = 238). Specifically, we focused on a developmental anomaly of the hippocampus that is characterized by the incomplete infolding of hippocampal subfields during fetal development, giving the hippocampus a more round or pyramidal shape, called incomplete hippocampal inversion (IHI). We rated IHI according to standard criteria and extracted hippocampal subfield volumes after FreeSurfer's automatic segmentation. We observed similar IHI prevalence in the group of individuals with DTD with respect to the control population. Neither differences in whole hippocampal nor major hippocampal subfield volumes have been observed between groups. However, when assessing the IHI independent criteria, we observed that the hippocampus in the DTD group is more medially positioned comparing to the control group. In addition, we observed bigger hippocampal fissure volume for the DTD comparing to the control group. Both of these findings were stronger for the right hippocampus comparing to the left. Our results provide new insights regarding the hippocampal morphology of individuals affected by DTD, highlighting the role of structural anomalies during early prenatal development in line with the developmental nature of the spatial disorientation deficit.

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.