Differences in REE patterns of calcite from extensional and shear veins of the Sestola Vidiciatico Tectonic Unit in the Northern Apennines suggest variations in fluid source during the seismic cycle in an ancient analogue of a shallow megathrust (T_max c. 100–150°C). In shear veins, a positive Eu anomaly suggests an exotic fluid source, probably hotter than the fault environment. Small-scale extensional veins were derived instead from a local fluid in equilibrium with the fault rocks. Mutually crosscutting relations between two extensional vein sets, parallel and perpendicular to the megathrust, suggest repeated shifting of the ₁ and ₃ stresses during the seismic cycle. This is consistent with: (1) a seismic phase, with brittle failure along the thrust, crystallization of shear veins from an exotic fluid, stress drop and stress rotation; (2) a post-seismic phase, with fault-normal compaction and formation of fault-normal extensional veins fed by local fluids; (3) a reloading phase, where shear stress and pore pressure are gradually restored and fault-parallel extensional veins form, until the thrust fails again. The combination of geochemical and structural analyses in veins from exhumed megathrust analogues represents a promising tool to better understand the interplay between stress state and fluids in modern subduction zones.

Supplementary material: Cathodoluminescence microphotographs, methodological details of the microstructural analysis, microphotographs of the location of analysed spots and a geochemical data table are available at https://doi.org/10.6084/m9.figshare.c.4842165

Thematic collection: This article is part of the Polygenetic mélanges collection available at: https://www.lyellcollection.org/cc/polygenetic-melanges

The authors present results of the first high-resolution deep seismic reflection survey in the Pieniny Klippen Belt (PKB) in Poland. This survey sheds new light on the matter of olistostromes and the mélange character of the PKB. The sedimentary mass-transport deposits represented by olistoliths and olistostromes manifest themselves by different petrophysical parameters of rocks (velocity, density and resistivity) and seismic attributes. Seismic attributes are very effective in the interpretation of the geology of complex mélanges. The authors used selected attributes: low-pass filter, energy, energy gradient, dip-steered median filter, Prewitt filter, Laplacian edge enhancing filter and square root of the energy gradient. These attributes emphasize changes of the seismic image inside mélange zones. The distinguished olistoliths are now inside imbricated thrust structures and they are tectonically rearranged. Polygenetic mélanges in the PKB originated as a result of sedimentary and tectonic processes. The PKB in the investigated area forms several north-vergent thrust sheets belonging to the Złatne and Hulina nappes. Both nappes contain large chaotic, non-reflective olistoliths as well as the smaller mainly high-reflective olistoliths. Olistoliths are arranged parallel to the flysch layering and thrusts. The results presented confirm the postulated two olistostrome belts within the PKB structure.

Thematic collection: This article is part of the Polygenetic mélanges collection available at: https://www.lyellcollection.org/cc/polygenetic-melanges

The Dukla Nappe in the Skrzydlna area exposes two types of mélange reflecting two different phases of basinal and tectonic evolution of the Outer West Carpathian orogen in its Polish sector. The Oligocene-age sedimentary mélange (olistostrome) is related to growth of the accretionary wedge, whereas the Miocene-age diapiric mélange postdates the orogenic thrusting. Textural and structural features of the very coarse-grained sedimentary mélange suggest non-cohesive debris flows and high-density turbidity currents as predominant emplacement mechanisms. Growth strata, associated with progressive unconformities, and facies contrast between the underlying fine-grained unit and the overlying olistostrome reflect a considerable uplift of the source area and rotation of the adjacent part of the basin floor. The olistostrome and the overlying turbidite succession form a retrogressive sequence interpreted as a submarine canyon infill grading to a small submarine fan. The diapiric mélange, injected into the Oligocene-age succession of the Dukla Nappe, contains the Early and Late Cretaceous-age blocks and matrix derived from the underlying Silesian Nappe. The features reflecting diapiric emplacement include matrix proportion increase and block content decrease towards the mélange margins, scaly fabric and shear zones. Both mélanges, interpreted in the past as chaotic bodies, upon detailed examination reveal genesis-related subtle internal organization.

The Lower and Middle Ordovician of the Yangtze Platform, China, is characterized by a sedimentary succession dominated by carbonate rocks. Three sections spanning the Nantsinkuan/Lunshan, Fenhsiang, Hunghuayuan, and Dawan/Zitai formations, corresponding to the Tremadocian–Dapingian in age, have been sampled for high-resolution ¹³C chemostratigraphy (542 samples in total). Our new ¹³C data reveal five tie-points with the potential for global correlation: (1) a positive ¹³C excursion in the lower Nantsinkuan Formation within the Tremadocian Rossodus manitouensis Zone; (2) an excursion with two peaks roughly within the late Tremadocian Paltodus ‘deltifer’ Zone; (3) a positive ¹³C shift in the lower Hunghuayuan Formation, within the early Floian Serratognathus diversus Zone; (4) a gradual positive ¹³C shift in the late Floian, ranging from the uppermost S. diversus Zone to the basal Oepikodus evae Zone; (5) a minor negative shift in the lower Dawan/Zitai Formation, within the early Dapingian Baltoniodus triangularis Zone. These excursions are herein used for correlation of the Yangtze Platform strata with successions from South China, North China, the Argentine Precordillera, North America and Baltica. From a palaeogeographical perspective the Gudongkou, Xiangshuidong and Daling sections represent depositional environments along an inner to outer ramp profile. ¹³C data from these sections show successively heavier (higher) ¹³C values with increasing depositional depth. This is interpreted as due to remineralization of organic carbon within the carbonate rocks.

Supplementary material: Carbon and oxygen isotope data are available at: https://doi.org/10.6084/m9.figshare.c.4767080

Laterites preserved on both sides of the Western Ghats Escarpment of Peninsular India have formed by long-term lateritic weathering essentially after India–Seychelles continental break-up following Deccan Traps emplacement (c. 63 myr ago). Supergene manganese ores of the Western Ghats were formed on Late Archean manganese protores. Among Mn oxides composing the ores, cryptomelane (K-rich Mn oxide) was characterized and dated by ⁴⁰Ar/³⁹Ar geochronology. Measured ages complement those previously obtained in other South Indian manganese ores from the hinterland plateau and further document three major weathering periods, c. 53–44, c. 39–22 and c. 14–10 Ma, the last being documented for the first time in India. These periods coincide with global palaeoclimatic proxies and date the lateritic weathering of three successive palaeolandscapes of the Western Ghats that evolved under slow denudation (c. 8 m Ma^–1) over the last 44 myr and were mostly incised during the Neogene (<22 Ma). This indicates that the Western Ghats are a relict of a South Indian plateau preserved at the headwaters of very long east-flowing river systems and above the Western Ghats Escarpment. Topography and denudation history of this landscape do not require Neogene tilt of the Peninsula as recently proposed.

Supplementary material: Full details of field and sample description, methods and analytical data including electron probe microanalyses of cryptomelane, and isotopic analyses and degassing spectra of irradiated cryptomelane grains are available at https://doi.org/10.6084/m9.figshare.c.4726661

In the Central Western Carpathians (CWC), most published paleomagnetic results from Permo-Mesozoic rocks document extensive remagnetizations and come from thin-skinned thrust units that have undergone multistage deformation. We present results from lower Triassic redbeds from the autochthonous cover overlying the basement that carry a primary magnetization. Petromagnetic results indicate that the dominant ferromagnetic carrier is hematite, while magnetic susceptibility and its anisotropy are controlled by both ferromagnetic and paramagnetic minerals. Magnetic fabrics document weak deformation related to Late Cretaceous shortening. The directions of the high unblocking temperature remanence components pass both reversal and fold tests, attesting to their primary nature. Paleomagnetic inclinations are flatter than expected from reference datasets, suggesting small latitudinal separation between the CWC and stable Europe. Paleomagnetic declinations are mostly clustered within individual mountain massifs, implying their tectonic coherence. They show only minor differences between the massifs, indicating a lack of significant vertical-axis tectonic rotations within the studied central parts of the CWC. The paleomagnetic declinations are therefore representative of the whole of the CWC in terms of regional paleogeographic interpretations, and imply moderate counterclockwise rotations (c. 26°) of the region with respect to stable Europe since the Early Triassic.

Mechanical heterogeneity of a sedimentary sequence exerts a primary control on the geometry of fault zones and the proportion of offset accommodated by folding. The Wildensbuch Fault Zone in the Swiss Molasse Basin, with a maximum throw of 40 m, intersects a Mesozoic section containing a thick (120 m) clay-dominated unit (Opalinus Clay) over- and underlain by more competent limestone units. Interpretation of a 3D seismic reflection survey indicates that the fault zone formed by upward propagation of an east–west-trending basement structure, through the Mesozoic section, in response to NE–SW Miocene extension. This configuration formed an array of left-stepping normal fault segments above and below the Opalinus Clay. In cross-section a broad monoclinal fold is observed in the Opalinus Clay. Folding, however, is not ubiquitous and occurs in the Opalinus Clay where fault segments above and below are oblique to one another; where they are parallel the fault passes through the Opalinus Clay with little folding. These observations demonstrate that, even in strongly heterogeneous sequences, here a four-fold difference in both Young's modulus and cohesion between layers, the occurrence of folding may depend on the local relationship between fault geometry and applied stress field rather than rheological properties alone.

The Northern Paraguai Belt, at the SE border of the Amazonian Craton, central Brazil, has been interpreted as a Brasiliano–Pan-African (c. 650–600 Ma) belt with a foreland basin, recording collisional polyphase tectonism and greenschist-facies metamorphism extending from the late Precambrian to the Cambrian–Ordovician. New structural investigations indicate that the older metasedimentary rocks of the Cuiabá Group represent a Tonian–Cryogenian basement assemblage deformed in two contemporaneous fault-bounded structural sub-domains of wrench-dominated (rake <10°) and contraction-dominated (rake ~30–40°) sinistral transpression, with tectonic vergence towards the SE. The younger late Cryogenian to early Cambrian sedimentary rocks lying to the NW of the Cuiabá Group are non-metamorphic and display only pervasive brittle transtension characterized by normal-oblique faults, fractures and forced drag folds with no consistent vergence pattern. Our analyses suggest that an unconformity separates the metasedimentary Cuiabá Group basement of the Northern Paraguai Belt from the unmetamorphosed sedimentary cover. It is proposed that the latter units were deposited during a post-glacial marine transgression (after c. 635 Ma) in a post-collisional basin. The Paraguai Belt basement and its post-collisional sedimentary cover share a number of significant geological similarities with sequences in the Bassarides Belt and Taoudéni Basin in the SW portion of the West African Craton.

The Tabaquito batholith (Frontal Cordillera, western Argentina), is mainly composed of shallowly emplaced granodiorite to minor monzogranite with abundant mafic microgranular enclaves. New sensitive high-resolution ion microprobe U–Pb zircon ages of c. 337 Ma (biotite granodiorite) and c. 284 Ma (mafic dyke) along with previously published geochronological data suggest that a long-lived magmatic system formed through at least two magmatic pulses at c. 337 and c. 322 Ma with later superimposition of Permian magmatism. The Tabaquito granitoids are metaluminous, calc-alkalic and magnesian with I-type affinity. Elevated Th/Nb, Y/Nb and La/Nb ratios along with negative Nb–Ta and positive Pb anomalies are consistent with a continental arc setting. Hf, Nd and Sr isotopic composition of the Tabaquito granitoids suggests that their source could result from mixing of an old felsic crustal component and a juvenile mafic to intermediate component. New geochronological and geochemical data together with published data reveal a continuous arc setting from the Carboniferous to the Permian in Argentina, and important magmatic compositional variations through time and space controlled by episodic fluctuations in the subduction angle of the oceanic plate. Reported and compiled data allow us to infer the continuity of the Carboniferous magmatic arc along the west margin of Gondwana.

Supplementary material: Detailed petrography, analytical methods and data, zircon cathodoluminescence images and supplementary figures are available at https://doi.org/10.6084/m9.figshare.c.4763993

Select agents (SA) pose unique challenges for licensing vaccines and therapies. In the case of toxin-mediated diseases, HHS assigns guidelines for SA use, oversees vaccine and therapy development, and approves animal models and approaches to identify mechanisms for toxin neutralization. In this commentary, we discuss next-generation vaccines and therapies against ricin toxin and botulinum toxin, which are regulated SA toxins that utilize structure-based approaches for countermeasures to guide rapid response to future biothreats.

Bovine tuberculosis (TB), caused by Mycobacterium bovis, remains an important zoonotic disease posing a serious threat to livestock and wildlife. The current TB tests relying on cell-mediated and humoral immune responses in cattle have performance limitations. To identify new serodiagnostic markers of bovine TB, we screened a panel of 101 recombinant proteins, including 10 polyepitope fusions, by a multiantigen print immunoassay (MAPIA) with well-characterized serum samples serially collected from cattle with experimental or naturally acquired M. bovis infection. A novel set of 12 seroreactive antigens was established. Evaluation of selected proteins in the dual-path platform (DPP) assay showed that the highest diagnostic accuracy (~95%) was achieved with a cocktail of five best-performing antigens, thus demonstrating the potential for development of an improved and more practical serodiagnostic test for bovine TB.

The kinetics, longevity, and breadth of antibodies to influenza virus neuraminidase (NA) in archival, sequential serum/plasma samples from influenza A virus (IAV) H5N1 infection survivors and from patients infected with the 2009 pandemic IAV (H1N1) virus were determined using an enzyme-linked lectin-based assay. The reverse-genetics-derived H4N1 viruses harboring a hemagglutinin (HA) segment from A/duck/Shan Tou/461/2000 (H4N9) and an NA segment derived from either IAV H5N1 clade 1, IAV H5N1 clade 2.3.4, the 2009 pandemic IAV (H1N1) (H1N1pdm), or A/Puerto Rico/8/1934 (H1N1) virus were used as the test antigens. These serum/plasma samples were also investigated by microneutralization (MN) and/or hemagglutination inhibition (HI) assays. Neuraminidase-inhibiting (NI) antibodies against N1 NA of both homologous and heterologous viruses were observed in H5N1 survivors and H1N1pdm patients. H5N1 survivors who were never exposed to H1N1pdm virus developed NI antibodies to H1N1pdm NA. Seroconversion of NI antibodies was observed in 65% of the H1N1pdm patients at day 7 after disease onset, but an increase in titer was not observed in serum samples obtained late in infection. On the other hand, an increase in seroconversion rate with the HI assay was observed in the follow-up series of sera obtained on days 7, 14, 28, and 90 after infection. The study also showed that NI antibodies are broadly reactive, while MN and HI antibodies are more strain specific.

RiVax is a promising recombinant ricin toxin A subunit (RTA) vaccine antigen that has been shown to be safe and immunogenic in humans and effective at protecting rhesus macaques against lethal-dose aerosolized toxin exposure. We previously used a panel of RTA-specific monoclonal antibodies (MAbs) to demonstrate, by competition enzyme-linked immunosorbent assay (ELISA), that RiVax elicits similar serum antibody profiles in humans and macaques. However, the MAb binding sites on RiVax have yet to be defined. In this study, we employed hydrogen exchange-mass spectrometry (HX-MS) to localize the epitopes on RiVax recognized by nine toxin-neutralizing MAbs and one nonneutralizing MAb. Based on strong protection from hydrogen exchange, the nine MAbs grouped into four spatially distinct epitope clusters (namely, clusters I to IV). Cluster I MAbs protected RiVax's α-helix B (residues 94 to 107), a protruding immunodominant secondary structure element known to be a target of potent toxin-neutralizing antibodies. Cluster II consisted of two subclusters located on the "back side" (relative to the active site pocket) of RiVax. One subcluster involved α-helix A (residues 14 to 24) and α-helices F-G (residues 184 to 207); the other encompassed β-strand d (residues 62 to 69) and parts of α-helices D-E (154 to 164) and the intervening loop. Cluster III involved α-helices C and G on the front side of RiVax, while cluster IV formed a sash from the front to back of RiVax, spanning strands b, c, and d (residues 35 to 59). Having a high-resolution B cell epitope map of RiVax will enable the development and optimization of competitive serum profiling assays to examine vaccine-induced antibody responses across species.

We previously produced a heavy-chain-only antibody (Ab) VH domain (V_HH)-displayed phage library from two alpacas that had been immunized with ricin toxoid and nontoxic mixtures of the enzymatic ricin toxin A subunit (RTA) and binding ricin toxin B subunit (RTB) (D. J. Vance, J. M. Tremblay, N. J. Mantis, and C. B. Shoemaker, J Biol Chem 288:36538–36547, 2013, https://doi.org/10.1074/jbc.M113.519207). Initial and subsequent screens of that library by direct enzyme-linked immunosorbent assay (ELISA) yielded more than two dozen unique RTA- and RTB-specific V_HHs, including 10 whose structures were subsequently solved in complex with RTA. To generate a more complete antigenic map of ricin toxin and to define the epitopes associated with toxin-neutralizing activity, we subjected the V_HH-displayed phage library to additional "pannings" on both receptor-bound ricin and antibody-captured ricin. We now report the full-length DNA sequences, binding affinities, and neutralizing activities of 68 unique V_HHs: 31 against RTA, 33 against RTB, and 4 against ricin holotoxin. Epitope positioning was achieved through cross-competition ELISAs performed with a panel of monoclonal antibodies (MAbs) and verified, in some instances, with hydrogen-deuterium exchange mass spectrometry. The 68 V_HHs grouped into more than 20 different competition bins. The RTA-specific V_HHs with strong toxin-neutralizing activities were confined to bins that overlapped two previously identified neutralizing hot spots, termed clusters I and II. The four RTB-specific V_HHs with potent toxin-neutralizing activity grouped within three adjacent bins situated at the RTA-RTB interface near cluster II. These results provide important insights into epitope interrelationships on the surface of ricin and delineate regions of vulnerability that can be exploited for the purpose of vaccine and therapeutic development.

The global burden of disease caused by extraintestinal pathogenic Escherichia coli (ExPEC) is increasing as the prevalence of multidrug-resistant strains rises. A multivalent ExPEC O-antigen bioconjugate vaccine could have a substantial impact in preventing bacteremia and urinary tract infections. Development of an ExPEC vaccine requires a readout to assess the functionality of antibodies. We developed an opsonophagocytic killing assay (OPA) for four ExPEC serotypes (serotypes O1A, O2, O6A, and O25B) based on methods established for pneumococcal conjugate vaccines. The performance of the assay was assessed with human serum by computing the precision, linearity, trueness, total error, working range, and specificity. Serotypes O1A and O6A met the acceptance criteria for precision (coefficient of variation for repeatability and intermediate precision, ≤50%), linearity (90% confidence interval of the slope of each strain, 0.80, 1.25), trueness (relative bias range, –30% to 30%), and total error (total error range, –65% to 183%) at five serum concentrations and serotypes O2 and O25B met the acceptance criteria at four concentrations (the lowest concentration for serotypes O2 and O25B did not meet the system suitability test of maximum killing of ≥85% of E. coli cells). All serotypes met the acceptance criteria for specificity (opsonization index value reductions of ≤20% for heterologous serum preadsorption and ≥70% for homologous serum preadsorption). The assay working range was defined on the basis of the lowest and highest concentrations at which the assay jointly fulfilled the target acceptance criteria for linearity, precision, and accuracy. An OPA suitable for multiple E. coli serotypes has been developed, qualified, and used to assess the immunogenicity of a 4-valent E. coli bioconjugate vaccine (ExPEC4V) administered to humans.

In the article "Neurologic syndromes related to anti-GAD65: Clinical and serologic response to treatment" by Muñoz-Lopetegi et al.,¹ published online March 2, 2020, the y-axis label for figure 5’s right graph should be "CSF anti-GAD65 concentration (IU/mL)." The editorial office regrets the error.

Objective

To describe the presentations, radiologic features, and outcomes of children with autoimmune encephalitis associated with myelin oligodendrocyte glycoprotein antibodies (MOG abs).

Objective

To determine whether live-attenuated yellow fever vaccine (YFV) was associated with MS relapse, we evaluated the clinical courses of 23 patients in the year before and the year after immunization at the university hospital of Geneva, Switzerland.

Objective

To assess possible adverse effects on breastfed infants of mothers receiving monoclonal antibodies (MAbs) during pregnancy and/or lactation.

Objective

To determine the characteristic clinical and spinal MRI phenotypes of sarcoidosis-associated myelopathy (SAM), we analyzed a large cohort of patients with this disorder.

Most British railway embankments were constructed between 120 and 180 years ago without the benefit of modern design and construction methods. This can result in undesirable load-deformation characteristics and consequent disruption to present-day railway operations, for which there is unprecedented demand. Annual rail passenger kilometres have approximately doubled in the last 20 years and freight has increased by 60% over the same period. Whereas elements such as rails or bridges can be refurbished or replaced to meet increasing demand, the same is not usually feasible for embankments. Development of techniques to assess embankment performance risks posed by operational capacity enhancements is therefore of increasing significance to railway geotechnical asset management. The two case studies presented in this paper demonstrate how geospatial analysis and data management techniques may be applied to this challenge at both strategic (regional or national) and tactical (site-specific) scales for embankments incorporating plastic clay fill. The case studies also demonstrate, in a world of ever more abundant data, the growing need for engineering geologists and geotechnical engineers to augment their traditional knowledge with comprehensive data management and geospatial analysis skills, these being essential for modern infrastructure asset management.

Thematic collection: This article is part of the ‘Ground-related risk to transportation infrastructure’ collection available at: https://www.lyellcollection.org/cc/Ground-related-risk-to-transportation-infrastructure

Transportation infrastructure owners manage an array of different asset types such as bridges, road pavements, earthworks and drainage. Currently, most organization management procedures are siloed by asset type; however, there are important interactions between these asset groups that need to be managed in a cross-asset way. Although these interactions are known, there is little or no quantification of these interactions. For the first time, this paper quantifies that 74% of Highways England's earthwork failures are a result of drainage-related problems, either the lack of drainage infrastructure or the poor performance of it. The analysis undertaken is an important first step not only in moving towards more connected asset management planning for earthworks and drainage, but to also provide guidance for other owners of earthwork infrastructure assets to improve their strategic asset management procedures.

Thematic collection: This article is part of the Ground-related risk to transportation infrastructure collection available at https://www.lyellcollection.org/cc/Ground-related-risk-to-transportation-infrastructure

South Africa is a mineral-rich country with a diverse geology and a long history of mining. The rich history of mining activities includes the extraction of coal from the Ecca Group Sediments of the Karoo Supergroup (250 Ma), gold and uranium from the Witwatersrand Supergroup (2900 Ma), as well as platinum, uranium, tin and lead from the layered Bushveld Igneous Complex (BIC) (2150 Ma). The extraction of gold, copper, tin, lead and rare earth minerals also took place in the Archean rocks of Swazium age (3500–3000 Ma). The historical mining records have either not been accurately recorded or have been lost over time. This has resulted in significant geohazard risk during infrastructure development, especially in and around historical mining towns, such as Johannesburg and Ermelo. These geohazard risks require careful appraisal and quantification prior to any infrastructure design or construction.

This case study aims to set out the development aspects of the Multi-Faceted Geophysical Modelling Systems approach, which was used by the South African National Roads Agency SOC Ltd (SANRAL) during an investigation of undermined ground for the historical coal-mining town of Ermelo in South Africa. The N11/N2 ring road was planned to go around Ermelo to ensure mobility between major routes, whilst still maintaining town access.

The systems approach used a combination of airborne geophysics (Versatile Time Domain Electromagnetic System (VTEM^TM) and magnetics), generally used in mining exploration, land-based and borehole geophysics, borehole water testing, and ground-truthing. The approach was continuous and iterative, building on the data at hand and reducing unnecessary investigations while eliminating the possibility of anomalies being missed, as in the case of conventional discrete drilling. The investigation ensured that 100% of the route was comprehensively investigated with a high confidence in the geological and geophysical data, and concomitant mitigation of infrastructure risk.

The Multi-Faceted Geophysical Modelling Systems approach was successfully used to identify a previously unknown 1 x 1 m mining stope cavity at 90 m depth and a 3 x 5 m access tunnel at 24 m depth in a timely and cost-effective manner. Seven reverse-circulation percussion boreholes confirmed the structural integrity of these underground cavities, as well as the structural geology along the centreline. Based on the great success achieved in identifying shallow anomalies, this Multi-Faceted Geophysical Modelling Systems approach is now being considered for field trails on the dolomitic formations and the Wild Coast greenfields road project where there are large historical slumps and many fault lines.

Thematic collection: This article is part of the Ground-related risk to transportation infrastructure collection available at https://www.lyellcollection.org/cc/Ground-related-risk-to-transportation-infrastructure

Strategic geotechnical asset management considers the whole of an organization's earthworks portfolio and is concerned with setting an overall earthworks asset management policy with long-term objectives related to asset performance, safety and condition, and identifying how those objectives can best be met, now and into the future. A risk-based approach is adopted that requires an understanding of the likelihood that any of the earthworks may fail, combined with a knowledge of the consequences should they fail. Procedures are required to identify those earthworks that are most vulnerable to failure under the influence of triggering events, such as extreme weather. The risks are managed through a mix of interventions to reduce the likelihood of failure and mitigations to reduce the impact of failure. Many of the challenges of implementing a strategic earthworks policy have, or are, being met by the main UK transportation infrastructure organizations.

Thematic collection: This article is part of the Ground-related risk to transportation infrastructure collection available at https://www.lyellcollection.org/cc/Ground-related-risk-to-transportation-infrastructure

This case study aims to identify the cubic capacity and geometry of the geological body of silt–organic sediments in the environment of a former gravel pit situated in a drainless depression of the alluvium of the Čierna voda River. It is located in the well-known geotourism locality of Slnečné Jazera in Senec, in the SW of Slovakia. To identify the body, electrical resistivity tomography was combined with the use of sonar. The research shows that this approach is appropriate for a number of activities that are subjects of engineering-geological investigations. The cubic capacity and geometry of specific aqueous engineering-geological environments must be determined in connection with the need for the management, control, quantification and extraction of selected sedimentary bodies. In addition, the management of sustainable development of reservoirs, sedimentation basins, industrial ponds, settling pits and natural pools for recreation (the subject of the case study) is important to control the limit amounts of sediments in such environments. The results of this study may be applied in analogous engineering-geological conditions. The drainless depression Slnečné Jazera contained a body of silt–organic sediments amounting to 23 000 m³ (41 Olympic-size pools of 25 m x 12.5 m x 1.8 m). The maximum thickness of the bottom sediments was about 6.3 m on the alluvium with an articulated morphology owing to the submerged digging of gravel. The proposed approach improved the control of extraction of the sedimentary body and optimized the engineering-geological conditions in this geotourism locality.

This paper uses a 3D model for stability assessment of Zhujiabaobao waste dump with ground cracks. The study data were gathered via reconnaissance, geomorphological analysis and laboratory experiment. A 3D finite extended element method model that can consider cracks was then used to calculate the factor of safety (FOS) of the waste dump via the strength reduction technique. The simulation shows the dump to have an FOS of 1.22 and both the position and depth of penetration of cracks in the waste dump have a crucial impact on the stability of the slope. Because the study area is located in a seismically active area, simulation and analysis of the dynamic response of the waste dump under different magnitudes of seismic waves (peak acceleration is 0.05, 0.15, 0.25 and 0.45g) were performed via an explicit dynamic model. The simulation shows that high steps in the slope are particularly responsive to earthquakes. The approach used here for analysing stability under static and dynamic loads is useful for hazard prevention and mitigation.

Groundwater recharge and runoff conditions are ascertained in the Suixiao coal-mining district using the hydrogen and oxygen isotopes and the trace elements in the unconsolidated pore aquifer of the Cenozoic group, the fissured sandstone aquifer of the Permian system, and the karst fissured limestone aquifer of the Carboniferous Taiyuan Formation and the Ordovician system, which are the main recharge aquifers during coal mining. The main water–rock interactions are pyrite oxidation, cation exchange and adsorption, and carbonate acidification, which are educed by principal component analysis of conventional ions. These results combined with geological conditions prove that hydraulic connection exists generally between the main recharge aquifers, and the groundwater circulation is controlled by faults in the sandstone and limestone aquifers. The water–rock interaction is very weak in the east of the district, which is proved to be a recharge area by Fisher discriminant analysis. This study provides the theoretical basis for the hydrochemistry exploration and the establishment of a water-inrush warning system in a concealed coalfield.

A logistic regression model and a bivariate statistical analysis were used in this paper to evaluate the groundwater recharge susceptibility. The approach is based on the assessment of the relationship involving groundwater recharge and parameters that influence this hydrological process. Surface parameters and aquifer-related parameters were evaluated as thematic map layers using ArcGIS. Then, a weighted-rating method was adopted to categorize each parameter's map. To assess the role of each parameter in the aquifer recharge, a logistic regression model and a bivariate statistical analysis were applied to the Guenniche phreatic aquifer (Tunisia). Models are explored to establish a map showing the aquifer recharge susceptibility. The code Modflow was used to simulate the consequence of the recharge. The recharge amount was introduced in the model and was tested to verify the recharge effect on the hydraulic head for the two models. The obtained results reveal that the recharge as mapped in the bivariate statistical model has a minor impact on the hydraulic head. Results of the logistic regression model are more significant as the hydraulic head is widely affected. This model provides good results in mapping the spatial distribution of the aquifer recharge susceptibility.

Routine identification of fungal pathogens from positive blood cultures by culture-based methods can be time-consuming, delaying treatment with appropriate antifungal agents. The GenMark Dx ePlex investigational use only blood culture identification fungal pathogen panel (BCID-FP) rapidly detects 15 fungal targets simultaneously in blood culture samples positive for fungi by Gram staining. We aimed to determine the performance of the BCID-FP in a multicenter clinical study. Blood culture samples collected at 10 United States sites and tested with BCID-FP at 4 sites were compared to the standard-of-care microbiological and biochemical techniques, fluorescence in situ hybridization using peptide nucleic acid probes (PNA-FISH) and matrix-assisted laser desorption ionization–time of flight mass spectrometry (MALDI-TOF MS). Discrepant results were analyzed by bi-directional PCR/sequencing of residual blood culture samples. A total of 866 clinical samples, 120 retrospectively and 21 prospectively collected, along with 725 contrived samples were evaluated. Sensitivity and specificity of detection of Candida species (C. albicans, C. auris, C. dubliniensis, C. famata, C. glabrata, C. guilliermondii, C. kefyr, C. krusei, C. lusitaniae, C. parapsilosis, and C. tropicalis) ranged from 97.1 to 100% and 99.8 to 100%, respectively. For the other organism targets, sensitivity and specificity were as follows: 100% each for Cryptococcus neoformans and C. gattii, 98.6% and 100% for Fusarium spp., and 96.2% and 99.9% for Rhodotorula spp., respectively. In 4 of the 141 clinical samples, the BCID-FP panel correctly identified an additional Candida species, undetected by standard-of-care methods. The BCID-FP panel offers a faster turnaround time for identification of fungal pathogens in positive blood cultures that may allow for earlier antifungal interventions and includes C. auris, a highly multidrug-resistant fungus.

Infections due to methicillin-resistant Staphylococcus aureus (MRSA) are present worldwide and represent a major public health concern. The capability of PCR followed by high-resolution melt (HRM) curve analysis for the detection of community-associated and livestock-associated MRSA strains and the identification of staphylococcal protein A (spa) locus was evaluated in 74 MRSA samples which were isolated from the environment, humans, and pigs on a single piggery. PCR-HRM curve analysis identified four spa types among MRSA samples and differentiated MRSA strains accordingly. A nonsubjective differentiation model was developed according to genetic confidence percentage values produced by tested samples, which did not require visual interpretation of HRM curve results. The test was carried out at different settings, and result data were reanalyzed and confirmed with DNA sequencing. PCR-HRM curve analysis proved to be a robust and reliable test for spa typing and can be used as a tool in epidemiological studies.

Q fever, caused by Coxiella burnetii, is a worldwide zoonotic disease that may cause severe forms in humans and requires a specific and prolonged antibiotic treatment. Although current serological and molecular detection tools allow a reliable diagnosis of the disease, culture of C. burnetii strains is mandatory to assess their susceptibility to antibiotics and sequence their genome in order to optimize patient management and epidemiological studies. However, cultivating this fastidious microorganism is difficult and restricted to reference centers, as it requires biosafety level 3 laboratories and relies on cell culture performed by experienced technicians. In addition, the culture yield is low, which results in a small number of isolates being available. In this work, we developed a novel high-content screening (HCS) isolation strategy based on optimized high-throughput cell culture and automated microscopic detection of infected cells with specifically designed algorithms targeting cytopathic effects. This method was more efficient than the shell vial assay, at the level of time dependency, when applied to both frozen specimens (7 isolates recovered by HCS only, sensitivity 91% versus 78% for shell vial) and fresh samples (1 additional isolate using HCS, sensitivity 7% versus 5% for shell vial), for which most strains were recovered more rapidly with the new technique. In addition, detecting positive cultures by an automated microscope reduced the need for expertise and saved 24% of technician working time. Application of HCS to antibiotic susceptibility testing of 12 strains demonstrated that it was as efficient as the standard procedure that combines shell vial culture and quantitative PCR.

Microscopy and rapid diagnostic tests (RDTs) are the main diagnostic tools for malaria but fail to detect low-density parasitemias that are important for maintaining malaria transmission. To complement existing diagnostic methods, an isothermal reverse transcription-recombinase polymerase amplification and lateral flow assay (RT-RPA) was developed. We compared the performance with that of ultrasensitive reverse transcription-quantitative PCR (uRT-qPCR) using nucleic acid extracts from blood samples (n = 114) obtained after standardized controlled human malaria infection (CHMI) with Plasmodium falciparum sporozoites. As a preliminary investigation, we also sampled asymptomatic individuals (n = 28) in an area of malaria endemicity (Lambaréné, Gabon) to validate RT-RPA and assess its performance with unprocessed blood samples (dbRT-RPA). In 114 samples analyzed from CHMI trials, the positive percent agreement to uRT-qPCR was 90% (95% confidence interval [CI], 80 to 96). The negative percent agreement was 100% (95% CI, 92 to 100). The lower limit of detection was 64 parasites/ml. In Gabon, RT-RPA was 100% accurate with asymptomatic volunteers (n = 28), while simplified dbRT-RPA showed 89% accuracy. In a subgroup analysis, RT-RPA detected 9/10 RT-qPCR-positive samples, while loop-mediated isothermal amplification (LAMP) detected 2/10. RT-RPA is a reliable diagnostic test for asymptomatic low-density infections. It is particularly useful in settings where uRT-qPCR is difficult to implement.

Human granulocytic anaplasmosis (HGA) is a tick-borne disease caused by the obligate intracellular Gram-negative bacterium Anaplasma phagocytophilum. The disease often presents with nonspecific symptoms with negative serology during the acute phase. Direct pathogen detection is the best approach for early confirmatory diagnosis. Over the years, PCR-based molecular detection methods have been developed, but optimal sensitivity is not achieved by conventional PCR while real-time PCR requires expensive and sophisticated instruments. To improve the sensitivity and also develop an assay that can be used in resource-limited areas, an isothermal DNA amplification assay based on recombinase polymerase amplification (RPA) was developed. To do this, we identified a 171-bp DNA sequence within multiple paralogous copies of msp2 within the genome of A. phagocytophilum. Our novel RPA assay targeting this sequence has an analytical limit of detection of one genome equivalent copy of A. phagocytophilum and can reliably detect 125 bacteria/ml in human blood. A high level of specificity was demonstrated by the absence of nonspecific amplification using genomic DNA from human or DNA from other closely-related pathogenic bacteria, such as Anaplasma platys, Ehrlichia chaffeensis, Orientia tsutsugamushi, and Rickettsia rickettsii, etc. When applied to patient DNA extracted from whole blood, this new RPA assay was able to detect 100% of previously diagnosed A. phagocytophilum cases. The sensitivity and rapidness of this assay represents a major improvement for early diagnosis of A. phagocytophilum in human patients and suggest a role for better surveillance in its reservoirs or vectors, especially in remote regions where resources are limited.

Limited treatment options contribute to high morbidity/mortality rates with carbapenem-resistant, Gram-negative bacterial infections. New approaches for carbapenemase-producing organism (CPO) detection may help inform clinician decision-making on patient treatment and infection control. BD Phoenix CPO detect (CPO detect) detects and classifies carbapenemases in Enterobacterales, Acinetobacter baumannii, and Pseudomonas aeruginosa during susceptibility testing. The clinical performance of CPO detect is reported here. Enterobacterales, Acinetobacter baumannii, and Pseudomonas aeruginosa isolates were evaluated across three sites using CPO detect and a composite reference method (RM); the latter was comprised of the modified carbapenem inactivation method and a MIC screen for ertapenem, imipenem, and meropenem. Multiplex PCR testing was also utilized for Ambler class determination. Positive and negative percentages of agreement (PPA and NPA, respectively) between CPO detect and the RM were determined. The PPA and NPA for Enterobacterales were 98.5% (confidence intervals, 96.6%, 99.4%) and 97.2% (95.8%, 98.2%), respectively. The A. baumannii PPA and NPA, respectively, were 97.1% (90.2%, 99.2%) and 97.1% (89.9%, 99.2%). The P. aeruginosa PPA and NPA, respectively, were 95.9% (88.6%, 98.6%) and 92.3% (86.7%, 95.6%). The PPA values for carbapenemase class designations for all organisms combined and Enterobacterales alone, respectively, were 95.3% (90.2%, 97.8%) and 94.6% (88.8%, 97.5%) for class A, 94.0% (88.7%, 96.6%) and 96.4% (90.0%, 98.8%) for class B, and 95.0% (90.1%, 97.6%) and 99.0% (94.4%, 99.8%) for class D carbapenemases. NPA values for all organisms and Enterobacterales alone ranged from 98.5% to 100%. CPO detect provided accurate detection and classification of CPOs for the majority of isolates of Enterobacterales, Acinetobacter baumannii, and Pseudomonas aeruginosa tested.

MALDI-TOF mass spectrometry (MS) identification of pathogenic filamentous fungi is often impaired by difficulties in harvesting hyphae embedded in the medium and long extraction protocols. The ID Fungi Plate (IDFP) is a novel culture method developed to address such difficulties and improve the identification of filamentous fungi by MALDI-TOF MS. We cultured 64 strains and 11 clinical samples on IDFP, Sabouraud agar-chloramphenicol (SAB), and ChromID Candida agar (CAN2). We then compared the three media for growth, ease of harvest, amount of material picked, and MALDI-TOF identification scores after either rapid direct transfer (DT) or a long ethanol-acetonitrile (EA) extraction protocol. Antifungal susceptibility testing and microscopic morphology after subculture on SAB and IDFP were also compared for ten molds. Growth rates and morphological aspects were similar for the three media. With IDFP, harvesting of fungal material for the extraction procedure was rapid and easy in 92.4% of cases, whereas it was tedious on SAB or CAN2 in 65.2% and 80.3% of cases, respectively. The proportion of scores above 1.7 (defined as acceptable identification) were comparable for both extraction protocols using IDFP (P = 0.256). Moreover, rates of acceptable identification after DT performed on IDFP (93.9%) were significantly higher than those obtained after EA extraction with SAB (69.7%) or CAN2 (71.2%) (P = <0.001 and P = 0.001, respectively). Morphological aspects and antifungal susceptibility testing were similar between IDFP and SAB. IDFP is a culture plate that facilitates and improves the identification of filamentous fungi, allowing accurate routine identification of molds with MALDI-TOF-MS using a rapid-extraction protocol.

Clostridioides difficile infection (CDI) is one of the most common health care-associated infections that can cause significant morbidity and mortality. CDI diagnosis involves laboratory testing in conjunction with clinical assessment. The objective of this study was to assess the performance of various C. difficile tests and to compare clinical characteristics, Xpert C. difficile/Epi (PCR) cycle threshold (C_T), and Singulex Clarity C. diff toxins A/B (Clarity) concentrations between groups with discordant test results. Unformed stool specimens from 200 hospitalized adults (100 PCR positive and 100 negative) were tested by cell cytotoxicity neutralization assay (CCNA), C. diff Quik Chek Complete (Quik Chek), Premier Toxins A and B, and Clarity. Clinical data, including CDI severity and CDI risk factors, were compared between discordant test results. Compared to CCNA, PCR had the highest sensitivity at 100% and Quik Chek had the highest specificity at 100%. Among clinical and laboratory data studied, prevalences of leukocytosis, prior antibiotic use, and hospitalizations were consistently higher across all subgroups in comparisons of toxin-positive to toxin-negative patients. Among PCR-positive samples, the median C_T was lower in toxin-positive samples than in toxin-negative samples; however, C_T ranges overlapped. Among Clarity-positive samples, the quantitative toxin concentration was significantly higher in toxin-positive samples than in toxin-negative samples as determined by CCNA and Quik Chek Toxin A and B. Laboratory tests for CDI vary in sensitivity and specificity. The quantitative toxin concentration may offer value in guiding CDI diagnosis and treatment. The presence of leukocytosis, prior antibiotic use, and previous hospitalizations may assist with CDI diagnosis, while other clinical parameters may not be consistently reliable.

The number of onsite clinical microbiology laboratories in hospitals is decreasing, likely related to the business model for laboratory consolidation and labor shortages, and this impacts a variety of clinical practices, including that of banking isolates for clinical or epidemiologic purposes. To determine the impact of these trends, infectious disease (ID) physicians were surveyed regarding their perceptions of offsite services. Clinical microbiology practices for retention of clinical isolates for future use were also determined. Surveys were sent to members of the Infectious Diseases Society of America’s (IDSA) Emerging Infections Network (EIN). The EIN is a sentinel network of ID physicians who care for adult and/or pediatric patients in North America and who are members of IDSA. The response rate was 763 (45%) of 1,680 potential respondents. Five hundred forty (81%) respondents reported interacting with the clinical microbiology laboratory. Eighty-six percent of respondents thought an onsite laboratory very important for timely diagnostic reporting and ongoing communication with the clinical microbiologist. Thirty-five percent practiced in institutions where the core microbiology laboratory has been moved offsite, and an additional 7% (n = 38) reported that movement of core laboratory functions offsite was being considered. The respondents reported that only 24% of laboratories banked all isolates, with the majority saving isolates for less than 30 days. Based on these results, the trend toward centralized core laboratories negatively impacts the practice of ID physicians, potentially delays effective implementation of prompt and targeted care for patients with serious infections, and similarly adversely impacts infection control epidemiologic investigations.

There are roughly 48,000 deaths caused by influenza annually and an estimated 200,000 people who have undiagnosed human immunodeficiency virus (HIV). These are examples of acute and chronic illnesses that can be identified by employing a CLIA-waived test. Pharmacies across the country have been incorporating CLIA-waived point-of-care tests (POCT) into disease screening and management programs offered in the pharmacy. The rationale behind these programs is discussed. Additionally, a summary of clinical data for some of these programs in the infectious disease arena is provided. Finally, we discuss the future potential for CLIA-waived POCT-based programs in community pharmacies.

Human herpesvirus 6 (HHV-6) is an important cause of meningitis and meningoencephalitis. As testing for HHV-6 in cerebrospinal fluid (CSF) is more readily available using the FilmArray Meningitis/Encephalitis panel (FA-ME; BioFire Diagnostics, Salt Lake City, UT), we aimed to determine the clinical significance of detecting HHV-6 in order to identify true infections and to ensure appropriate antiviral initiation. Chart review on 25 patients positive for HHV-6 by FA-ME was performed to determine clinical presentation, comorbidity, treatment, and outcome. The presence of chromosomally integrated HHV-6 (ciHHV-6) DNA was also investigated. Of 1,005 children tested by FA-ME, HHV-6 was detected in 25 (2.5%). Five patients were diagnosed with either HHV-6 meningitis or meningoencephalitis based on HHV-6 detection in CSF, clinical presentation, and radiographic findings. Detection of HHV-6 by FA-ME led to discontinuation of acyclovir within 12.0 h in all 12 patients empirically treated with acyclovir. Six of the 12 patients were started on ganciclovir therapy within 6.8 h; 4 of these were treated specifically for HHV-6 infection, whereas therapy was discontinued in the remaining 2 patients. CSF parameters were not generally predictive of HHV-6 positivity. The presence of ciHHV-6 was confirmed in 3 of 18 patients who could be tested. Five of the 25 patients included in the study were diagnosed with HHV-6 meningitis/meningoencephalitis. FA-ME results led to discontinuation of empirical antiviral treatment in 12 patients and appropriate initiation of ganciclovir in 4 patients. In our institution, detection of HHV-6 using FA-ME led to faster establishment of disease etiology and optimization of antimicrobial therapy.

On 31 December 2019, the World Health Organization was informed of a cluster of cases of pneumonia of unknown etiology in Wuhan, China. Subsequent investigations identified a novel coronavirus, now named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), from the affected patients. Highly sensitive and specific laboratory diagnostics are important for controlling the rapidly evolving SARS-CoV-2-associated coronavirus disease 2019 (COVID-19) epidemic. In this study, we developed and compared the performance of three novel real-time reverse transcription-PCR (RT-PCR) assays targeting the RNA-dependent RNA polymerase (RdRp)/helicase (Hel), spike (S), and nucleocapsid (N) genes of SARS-CoV-2 with that of the reported RdRp-P2 assay, which is used in >30 European laboratories. Among the three novel assays, the COVID-19-RdRp/Hel assay had the lowest limit of detection in vitro (1.8 50% tissue culture infective doses [TCID₅₀]/ml with genomic RNA and 11.2 RNA copies/reaction with in vitro RNA transcripts). Among 273 specimens from 15 patients with laboratory-confirmed COVID-19 in Hong Kong, 77 (28.2%) were positive by both the COVID-19-RdRp/Hel and RdRp-P2 assays. The COVID-19-RdRp/Hel assay was positive for an additional 42 RdRp-P2-negative specimens (119/273 [43.6%] versus 77/273 [28.2%]; P < 0.001), including 29/120 (24.2%) respiratory tract specimens and 13/153 (8.5%) non-respiratory tract specimens. The mean viral load of these specimens was 3.21 x 10⁴ RNA copies/ml (range, 2.21 x 10² to 4.71 x 10⁵ RNA copies/ml). The COVID-19-RdRp/Hel assay did not cross-react with other human-pathogenic coronaviruses and respiratory pathogens in cell culture and clinical specimens, whereas the RdRp-P2 assay cross-reacted with SARS-CoV in cell culture. The highly sensitive and specific COVID-19-RdRp/Hel assay may help to improve the laboratory diagnosis of COVID-19.

The IR Biotyper is a new automated typing system based on Fourier-transform infrared (FT-IR) spectroscopy that gives results within 4 h. We aimed (i) to use the IR Biotyper to retrospectively analyze an outbreak of extended-spectrum beta-lactamase-producing Klebsiella pneumoniae (ESBL-KP) in a neonatal intensive care unit and to compare results to BOX-PCR and whole-genome sequencing (WGS) results as the gold standard and (ii) to assess how the cutoff values used to define clusters affect the discriminatory power of the IR Biotyper. The sample consisted of 18 isolates from 14 patients. Specimens were analyzed in the IR Biotyper using the default analysis settings, and spectra were analyzed using OPUS 7.5 software. The software contains a feature that automatically proposes a cutoff value to define clusters; the cutoff value defines up to which distance the spectra are considered to be in the same cluster. Based on FT-IR, the outbreak represented 1 dominant clone, 1 secondary clone, and several unrelated clones. FT-IR results, using the cutoff value generated by the accompanying software after 4 replicates, were concordant with WGS for all but 1 isolate. BOX-PCR was underdiscriminatory compared to the other two methods. Using the cutoff value generated after 12 replicates, the results of FT-IR and WGS were completely concordant. The IR Biotyper can achieve the same typeability and discriminatory power as genome-based methods. However, to attain this high performance requires either previous, strain-dependent knowledge about the optimal technical parameters to be used or validation by a second method.

Serological testing for nasopharyngeal carcinoma (NPC) has recently been reinvigorated by the implementation of novel Epstein-Barr virus (EBV)-specific IgA and IgG antibodies from a proteome array. Although proteome arrays are well suited for comprehensive antigen selection, they are not applicable for large-scale studies. We adapted a 13-marker EBV antigen signature for NPC risk identified by proteome arrays to multiplex serology to establish an assay for large-scale studies. Taiwanese NPC cases (n = 175) and matched controls (n = 175) were used for assay validation. Spearman’s correlation was calculated, and the diagnostic value of all multiplex markers was assessed independently using the area under the receiver operating characteristic curve (AUC). Two refined signatures were identified using stepwise logistic regression and internally validated with 10-fold cross validation. Array and multiplex serology showed strong correlation for each individual EBV marker, as well as for a 13-marker combined model on continuous data. Two refined signatures with either four (LF2 and BGLF2 IgG, LF2 and BMRF1 IgA) or two (LF2 and BGLF2 IgG) antibodies on dichotomous data were identified as the most parsimonious set of serological markers able to distinguish NPC cases from controls with AUCs of 0.992 (95% confidence interval [CI], 0.983 to 1.000) and 0.984 (95% CI, 0.971 to 0.997), respectively. Neither differed significantly from the 13-marker model (AUC, 0.992; 95% CI, 0.982 to 1.000). All models were internally validated. Multiplex serology successfully validated the original EBV proteome microarray data. Two refined signatures of four and two antibodies were capable of detecting NPC with 99.2% and 98.4% accuracy.

In Arabidopsis, LTR retrotransposons are activated by mutations in the chromatin gene DECREASE in DNA METHYLATION 1 (DDM1), giving rise to 21- to 22-nt epigenetically activated siRNA (easiRNA) that depend on RNA DEPENDENT RNA POLYMERASE 6 (RDR6). We purified virus-like particles (VLPs) from ddm1 and ddm1rdr6 mutants in which genomic RNA is reverse transcribed into complementary DNA. High-throughput short-read and long-read sequencing of VLP DNA (VLP DNA-seq) revealed a comprehensive catalog of active LTR retrotransposons without the need for mapping transposition, as well as independent of genomic copy number. Linear replication intermediates of the functionally intact COPIA element EVADE revealed multiple central polypurine tracts (cPPTs), a feature shared with HIV in which cPPTs promote nuclear localization. For one member of the ATCOPIA52 subfamily (SISYPHUS), cPPT intermediates were not observed, but abundant circular DNA indicated transposon "suicide" by auto-integration within the VLP. easiRNA targeted EVADE genomic RNA, polysome association of GYPSY (ATHILA) subgenomic RNA, and transcription via histone H3 lysine-9 dimethylation. VLP DNA-seq provides a comprehensive landscape of LTR retrotransposons and their control at transcriptional, post-transcriptional, and reverse transcriptional levels.

Recent progress has been made in identifying genomic regions implicated in trait evolution on a microevolutionary scale in many species, but whether these are relevant over macroevolutionary time remains unclear. Here, we directly address this fundamental question using bird beak shape, a key evolutionary innovation linked to patterns of resource use, divergence, and speciation, as a model trait. We integrate class-wide geometric-morphometric analyses with evolutionary sequence analyses of 10,322 protein-coding genes as well as 229,001 genomic regions spanning 72 species. We identify 1434 protein-coding genes and 39,806 noncoding regions for which molecular rates were significantly related to rates of bill shape evolution. We show that homologs of the identified protein-coding genes as well as genes in close proximity to the identified noncoding regions are involved in craniofacial embryo development in mammals. They are associated with embryonic stem cell pathways, including BMP and Wnt signaling, both of which have repeatedly been implicated in the morphological development of avian beaks. This suggests that identifying genotype-phenotype association on a genome-wide scale over macroevolutionary time is feasible. Although the coding and noncoding gene sets are associated with similar pathways, the actual genes are highly distinct, with significantly reduced overlap between them and bill-related phenotype associations specific to noncoding loci. Evidence for signatures of recent diversifying selection on our identified noncoding loci in Darwin finch populations further suggests that regulatory rather than coding changes are major drivers of morphological diversification over macroevolutionary times.