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

The presence of large palaeo-landslides can hinder the construction of railway lines if they cause an alteration of the natural balance of the slope, significantly increasing the cost of the project. During the construction works in a section of the Madrid–Asturias high-speed railway line (Spain), a large-scale hillside instability affected 460 m of the section. The ground movement began 10 months after the start of the excavations and remained active throughout the observation period (2008–2010). Data provided by fieldwork, boreholes and instrumental monitoring have allowed the investigation of the geological units involved, together with the geometry and the kinematics of the mass movement. This landslide involves a Paleozoic basement with an estimated volume of 4 400 000 m³. It shows low displacement rates (<45 mm a^–1) and was accommodated on a single surface of rupture, which reaches more than 60 m depth. The kinematics is favoured by structural and lithological conditions, there being a strong rheological contrast in the lower part of the mobilized ground. This movement is linked to a large palaeo-landslide that was partially reactivated when excavations undermined the lower slopes of the hillside. Technical solutions for the stabilization have already caused an additional cost of 17% over the initial budget.

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.

This study addresses the major geo-environmental hazards caused by shallow coal mining in China's western eco-environment frangible area. These hazards are related to the high overburden pressure, surface subsidence, soil and water losses, and land desertification, with consequent vegetation and wildlife losses. To mitigate these hazards, three alternative backfill mining methods are proposed, for three typical shallow coal mining conditions, using aeolian sand-based backfilling materials, which are readily available in this area. The main influencing factor is the backfill material compaction ratio. Its effect on aquiclude deformation and water-conducting fracture evolution are assessed by numerical and physical simulation methods. The potential application of the proposed backfill coal mining alternatives is evaluated and discussed in detail. The results obtained are considered to be valuable for developing a strategy for the coordinated exploitation of coal resources and environmental protection in China's western frangible eco-environment area.

The need to characterize and distinguish geographically adjacent aggregate quarry sources prompted the SEM-EDS analysis of pyrite (FeS₂) within fill material taken from eight different quarry sources. This experiment was undertaken to investigate the possibility of geochemically separating these quarry sources based on the major element concentration of their pyrite. The results show that median values for Fe and S vary by up to 7.6 and 8.55 wt% respectively. By implementing statistical methods, including k-means clustering and principal component analysis, it is possible to geochemically discriminate three of the eight sources.

The Clay-with-flints Formation outcrops on the high chalk plateaux and interfluves of the chalk downs in southern England. Both current and historical definitions of the Clay-with-flints are detailed and important distinctions are identified with other deposits that appear identical but are formed in different ways. Historically pedological or geomorphological studies have been carried out on the deposit. Engineering studies are only carried out where the deposit is crossed by infrastructure. The physical and chemical processes acting on the deposit and the resulting effects on the physical properties are discussed.

In this paper, a decision tree is presented, constructed on the basis of hydrogeological characteristics (water table depth, freshwater thickness, surface area required and distance between wells), to choose the optimal groundwater extraction method in the case of a coastal unconfined aquifer. A comparison is made of the groundwater extraction methods in a freshwater aquifer of limited thickness occurring in coastal dunes in the eastern region of the Province of Buenos Aires (Argentina). The negative effects brought about by the wrong use of the groundwater extraction methods are analysed, because, as a result of excessive extraction, such methods lead to the dramatic decrease of the freshwater reserves. The decision tree is a useful tool to assist decision-makers as it suggests the most suitable groundwater extraction method options (vertical wells or wellpoints), as well as identifying areas that are unsuitable for sustainable groundwater extraction.

Groundwater is best described as an open access good: as a common pool resource others cannot be excluded from using and as a finite resource, its consumption leaves less for others to enjoy. In addition to its obvious benefits to humans and animals, it contributes in sustaining the hydro-environment and life on the planet. Scarcity causes groundwater to become an economic good and have an economic value. By and large, economic valuations respond to changes in the quantity and quality of groundwater. They are based on willingness to pay for maintaining its benefits or willingness to accept compensation for giving them up. There have been some 50 published groundwater valuation studies worldwide in the last 30 years, mostly focusing on groundwater quality and contamination. There is sparsity of valuations in Africa and Asia. The results suggest median economic values of 60–160 and mean values of 70–480 US dollars per household per year. Generally, values are higher in the USA followed by Europe, and lowest in Asia. People's income is a major factor affecting values, but this does not mean that they value groundwater less. Economic valuations and cost–benefit analysis are useful in highlighting people's priorities and choice of options. However, economics alone should not dictate actions for protecting groundwater resources from degradation and depletion.

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.

In 2015, the M_w 7.8 Gorkha earthquake struck Nepal, triggering thousands of landslides across the central and eastern Himalayas. These landslides had many adverse effects, including causing widespread damage to low-grade transport routes (e.g. tracks, footpaths) in rural regions that depend on tourism for survival. Langtang Valley is a glacial–periglacial landscape located 60 km north of Kathmandu. It is one of the most popular trekking regions in Nepal and has been severely affected by Gorkha earthquake-triggered and monsoon-triggered landsliding. Here, qualitative and quantitative observations from fieldwork and remote sensing are used to describe the materials and geomorphology of the landslides across Langtang Valley, and to quantify the extent to which coseismic and monsoon-triggered landslides have affected Langtang's trekking infrastructure. The dominant bedrock materials involved within Langtang landslides are found to be a range of gneisses and intruded leucogranites. In total, 64 landslides are found to have intersected trekking paths across Langtang, with coseismic and monsoon-triggered landslides having an impact on c. 3 km and 0.8 km of path respectively. It is observed that the practice of reconstructing paths through unstable landslide deposits is leaving the trekking infrastructure across Langtang increasingly vulnerable to future failure.

Leprosy is caused by Mycobacterium leprae, and it remains underdiagnosed in Burkina Faso. We investigated the use of fluorescent in situ hybridization (FISH) for detecting M. leprae in 27 skin samples (skin biopsy samples, slit skin samples, and skin lesion swabs) collected from 21 patients from Burkina Faso and three from Côte d’Ivoire who were suspected of having cutaneous leprosy. In all seven Ziehl-Neelsen-positive skin samples (four skin biopsy samples and three skin swabs collected from the same patient), FISH specifically identified M. leprae, including one FISH-positive skin biopsy sample that remained negative after testing with PCR targeting the rpoB gene and with the GenoType LepraeDR assay. Twenty other skin samples and three negative controls all remained negative for Ziehl-Neelsen staining, FISH, and rpoB PCR. These data indicate the usefulness of a microscopic examination of skin samples after FISH for first-line diagnosis of cutaneous leprosy. Accordingly, FISH represents a potentially useful point-of-care test for the diagnosis of cutaneous leprosy.

Accurate detection of influenza A virus (IAV) is crucial for patient management, infection control, and epidemiological surveillance. The World Health Organization and the Centers for Disease Control and Prevention have recommended using the M gene as the diagnostic gene target for reverse-transcription-PCR (RT-PCR). However, M gene RT-PCR has reduced sensitivity for recent IAV due to novel gene mutations. Here, we sought to identify novel diagnostic targets for the molecular detection of IAV using long-read third-generation sequencing. Direct nanopore sequencing from 18 nasopharyngeal specimens and one saliva specimen showed that the 5' and 3' ends of the PB2 gene and the entire NS gene were highly abundant. Primers selected for PB2 and NS genes were well matched with seasonal or avian IAV gene sequences. Our novel PB2 and NS gene real-time RT-PCR assays showed limits of detection similar to or lower than that of M gene RT-PCR and achieved 100% sensitivity and specificity in the detection of A(H1N1), A(H3N2), and A(H7N9) in nasopharyngeal and saliva specimens. For 10 patients with IAV detected by M gene RT-PCR conversion in sequentially collected specimens, NS and/or PB2 gene RT-PCR was positive in 2 (20%) of the initial specimens that were missed by M gene RT-PCR. In conclusion, we have shown that PB2 or NS gene RT-PCRs are suitable alternatives to the recommended M gene RT-PCR for diagnosis of IAV. Long-read nanopore sequencing facilitates the identification of novel diagnostic targets.

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.

Quantitative bacterial culture of bronchoalveolar lavage fluids (BALF) is labor-intensive, and the delay involved in performing culture, definitive identification, and susceptibility testing often results in prolonged use of broad-spectrum antibiotics. The Unyvero lower respiratory tract (LRT) panel (Curetis, Holzgerlingen, Germany) allows the multiplexed rapid detection and identification of 20 potential etiologic agents of pneumonia within 5 h of collection. In addition, the assay includes detection of gene sequences that confer antimicrobial resistance. We retrospectively compared the performance of the molecular panel to routine quantitative bacterial culture methods on remnant BALF. Upon testing 175 BALF, we were able to analyze positive agreement of 181 targets from 129 samples, and 46 samples were negative. The positive percent agreement (PPA) among the microbial targets was 96.5%, and the negative percent agreement (NPA) was 99.6%. The targets with a PPA of <100% were Staphylococcus aureus (34/37 [91.9%]), Streptococcus pneumoniae (10/11 [90.9%]), and Enterobacter cloacae complex (2/4 [50%]). For the analyzable resistance targets, concordance with phenotypic susceptibility testing was 79% (14/18). This study found the Unyvero LRT panel largely concordant with culture results; however, no outcome or clinical impact studies were performed.

Klebsiella species are problematic pathogens in neonatal units and may cause outbreaks, for which the sources of transmission may be challenging to elucidate. We describe the use of whole-genome sequencing (WGS) to investigate environmental sources of transmission during an outbreak of extended-spectrum-β-lactamase (ESBL)-producing Klebsiella michiganensis colonizing neonates. Ceftriaxone-resistant Klebsiella spp. isolated from neonates (or their mothers) and the hospital environment were included. Short-read sequencing (Illumina) and long-read sequencing (MinION; Oxford Nanopore Technologies) were used to confirm species taxonomy, to identify antimicrobial resistance genes, and to determine phylogenetic relationships using single-nucleotide polymorphism profiling. A total of 21 organisms (10 patient-derived isolates and 11 environmental isolates) were sequenced. Standard laboratory methods identified the outbreak strain as an ESBL-producing Klebsiella oxytoca, but taxonomic assignment from WGS data suggested closer identity to Klebsiella michiganensis. Strains isolated from multiple detergent-dispensing bottles were either identical or closely related by single-nucleotide polymorphism comparison. Detergent bottles contaminated by K. michiganensis had been used for washing milk expression equipment. No new cases were identified once the detergent bottles were removed. Environmental reservoirs may be an important source in outbreaks of multidrug-resistant organisms. WGS, in conjunction with traditional epidemiological investigation, can be instrumental in revealing routes of transmission and guiding infection control responses.

Screening for Chlamydia trachomatis and Neisseria gonorrhoeae at the pharyngeal, urogenital, and anorectal sites is recommended for men who have sex with men (MSM). Combining the three individual-site samples into a single pooled sample could result in significant cost savings, provided there is no significant sensitivity reduction. The aim of this study was to examine the sensitivity of pooled samples for detecting chlamydia and gonorrhea in asymptomatic MSM using a nucleic acid amplification test. Asymptomatic MSM who tested positive for chlamydia or gonorrhoea were invited to participate. Paired samples were obtained from participants prior to administration of treatment. To form the pooled sample, the anorectal swab was agitated in the urine specimen transport tube and then discarded. The pharyngeal swab and 2 ml of urine sample were then added to the tube. The difference in sensitivity between testing of pooled samples and individual-site testing was calculated against an expanded gold standard, where an individual is considered positive if either pooled-sample or individual-site testing returns a positive result. All samples were tested using the Aptima Combo 2 assay. A total of 162 MSM were enrolled in the study. Sensitivities of pooled-sample testing were 86% (94/109; 95% confidence interval [CI], 79 to 92%]) for chlamydia and 91% (73/80; 95% CI, 83 to 96%) for gonorrhea. The sensitivity reduction was significant for chlamydia (P = 0.02) but not for gonorrhea (P = 0.34). Pooling caused 22 infections (15 chlamydia and 7 gonorrhoea) to be missed, and the majority were single-site infections (19/22). Pooling urogenital and extragenital samples from asymptomatic MSM reduced the sensitivity of detection by approximately 10% for chlamydia but not for gonorrhea.

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.

Lyme borreliosis is a tick-borne disease caused by the Borrelia burgdorferi sensu lato complex. Bio-Rad Laboratories has developed a fully automated multiplex bead-based assay for the detection of IgM and IgG antibodies to B. burgdorferi. The BioPlex 2200 Lyme Total assay exhibits an improved rate of seropositivity in patients with early Lyme infection. Asymptomatic subjects from endemic and nonendemic origins demonstrated a seroreactivity rate of approximately 4% that was similar to other commercial assays evaluated in this study. Coupled to this result was the observation that the Lyme Total assay retained a high first-tier specificity of 96% while demonstrating a relatively high sensitivity of 91% among a well-characterized CDC Premarketing Lyme serum panel. The Lyme Total assay also performs well under a modified two-tier algorithm (sensitivity, 84.4 to 88.9%; specificity, 98.4 to 99.5%). Furthermore, the new assay is able to readily detect early Lyme infection in patient samples from outside North America.

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.

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.

Prosthetic joint infections are difficult to diagnose and treat due to biofilm formation by the causative pathogens. Pathogen identification relies on microbial culture that requires days to weeks, and in the case of chronic biofilm infections, lacks sensitivity. Diagnosis of infection is often delayed past the point of effective treatment such that only the removal of the implant is curative. Early diagnosis of an infection based on antibody detection might lead to less invasive, early interventions. Our study examined antibody-based assays against the Staphylococcus aureus biofilm-upregulated antigens SAOCOL0486 (a lipoprotein), glucosaminidase (a domain of SACOL1062), and SACOL0688 (the manganese transporter MntC) for detection of chronic S. aureus infection. We evaluated these antigens by enzyme-linked immunosorbent assay (ELISA) using sera from naive rabbits and rabbits with S. aureus-mediated osteomyelitis, and then we validated a proof of concept for the lateral flow assay (LFA). The SACOL0688 LFA demonstrated 100% specificity and 100% sensitivity. We demonstrated the clinical diagnostic utility of the SACOL0688 antigen using synovial fluid (SF) from humans with orthopedic implant infections. Elevated antibody levels to SACOL0688 in clinical SF specimens correlated with 91% sensitivity and 100% specificity for the diagnosis of S. aureus infection by ELISA. We found measuring antibodies levels to SACOL0688 in SF using ELISA or LFA provides a tool for the sensitive and specific diagnosis of S. aureus prosthetic joint infection. Development of the LFA diagnostic modality is a desirable, cost-effective option, potentially providing rapid readout in minutes for chronic biofilm infections.

Pyrazinamide (PZA) is considered the pivot drug in all tuberculosis treatment regimens due to its particular action on the persistent forms of Mycobacterium tuberculosis. However, no drug susceptibility test (DST) is considered sufficiently reliable for routine application. Although molecular tests are endorsed, their application is limited to known PZA resistance associated mutations. Microbiological DSTs for PZA have been restricted by technical limitations, especially the necessity for an acidic pH. Here, for the first time, MODS culture at neutral pH was evaluated using high PZA concentrations (400 and 800 μg/ml) to determine PZA susceptibility directly from sputum samples. Sputum samples were cultured with PZA for up to 21 days at 37°C. Plate reading was performed at two time points: R1 (mean, 10 days) and R2 (mean, 13 days) for each PZA concentration. A consensus reference test, composed of MGIT-PZA, pncA sequencing, and the classic Wayne test, was used. A total of 182 samples were evaluated. The sensitivity and specificity for 400 μg/ml ranged from 76.9 to 89.7 and from 93.0 to 97.9%, respectively, and for 800 μg/ml ranged from 71.8 to 82.1 and from 95.8 to 98.6%, respectively. Compared to MGIT-PZA, our test showed a similar turnaround time (medians of 10 and 12 days for PZA-sensitive and -resistant isolates, respectively). In conclusion, MODS-PZA is presented as a fast, simple, and low-cost DST that could complement the MODS assay to evaluate resistance to the principal first-line antituberculosis drugs. Further optimization of test conditions would be useful in order to increase its performance.

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.

This minireview focuses on the microbiologic evaluation of patients with asymptomatic bacteriuria, as well as indications for antibiotic treatment. Asymptomatic bacteriuria is defined as two consecutive voided specimens (preferably within 2 weeks) with the same bacterial species, isolated in quantitative counts of ≥10⁵ CFU/ml in women, including pregnant women; a single voided urine specimen with one bacterial species isolated in a quantitative count ≥10⁵ CFU/ml in men; and a single catheterized urine specimen with one or more bacterial species isolated in a quantitative count of ≥10⁵ CFU/ml in either women or men (or ≥10² CFU/ml of a single bacterial species from a single catheterized urine specimen). Any urine specimen with ≥10⁴ CFU/ml group B Streptococcus is significant for asymptomatic bacteriuria in a pregnant woman. Asymptomatic bacteriuria occurs, irrespective of pyuria, in the absence of signs or symptoms of a urinary tract infection. The two groups with the best evidence of adverse outcomes in the setting of untreated asymptomatic bacteriuria include pregnant women and patients who undergo urologic procedures with risk of mucosal injury. Screening and treatment of asymptomatic bacteriuria is not recommended in the following patient populations: pediatric patients, healthy nonpregnant women, older patients in the inpatient or outpatient setting, diabetic patients, patients with an indwelling urethral catheter, patients with impaired voiding following spinal cord injury, patients undergoing nonurologic surgeries, and nonrenal solid-organ transplant recipients. Renal transplant recipients beyond 1 month posttransplant should not undergo screening and treatment for asymptomatic bacteriuria. There is insufficient evidence to recommend for or against screening of renal transplant recipients within 1 month, patients with high-risk neutropenia, or patients with indwelling catheters at the time of catheter removal. Unwarranted antibiotics place patients at increased risk of adverse effects (including Clostridioides difficile diarrhea) and contribute to antibiotic resistance. Methods to reduce unnecessary screening for and treatment of asymptomatic bacteriuria aid in antibiotic stewardship.

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 new decade of the 21^st century (2020) started with the emergence of a novel coronavirus known as SARS-CoV-2 that caused an epidemic of coronavirus disease (COVID-19) in Wuhan, China. It is the third highly pathogenic and transmissible coronavirus after severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) emerged in humans. The source of origin, transmission to humans, and mechanisms associated with the pathogenicity of SARS-CoV-2 are not yet clear, however, its resemblance to SARS-CoV and several other bat coronaviruses was recently confirmed through genome sequencing-related studies. The development of therapeutic strategies is necessary in order to prevent further epidemics and cure infections. In this review, we summarize current information about the emergence, origin, diversity, and epidemiology of three pathogenic coronaviruses with a specific focus on the current outbreak in Wuhan, China. Furthermore, we discuss the clinical features and potential therapeutic options that may be effective against SARS-CoV-2.

The QIAstat-Dx Respiratory Panel (QIAstat-Dx RP) is a multiplex in vitro diagnostic test for the qualitative detection of 20 pathogens directly from nasopharyngeal swab (NPS) specimens. The assay is performed using a simple sample-to-answer platform with results available in approximately 69 min. The pathogens identified are adenovirus, coronavirus 229E, coronavirus HKU1, coronavirus NL63, coronavirus OC43, human metapneumovirus A and B, influenza A, influenza A H1, influenza A H3, influenza A H1N1/2009, influenza B, parainfluenza virus 1, parainfluenza virus 2, parainfluenza virus 3, parainfluenza virus 4, rhinovirus/enterovirus, respiratory syncytial virus A and B, Bordetella pertussis, Chlamydophila pneumoniae, and Mycoplasma pneumoniae. This multicenter evaluation provides data obtained from 1,994 prospectively collected and 310 retrospectively collected (archived) NPS specimens with performance compared to that of the BioFire FilmArray Respiratory Panel, version 1.7. The overall percent agreement between QIAstat-Dx RP and the comparator testing was 99.5%. In the prospective cohort, the QIAstat-Dx RP demonstrated a positive percent agreement of 94.0% or greater for the detection of all but four analytes: coronaviruses 229E, NL63, and OC43 and rhinovirus/enterovirus. The test also demonstrated a negative percent agreement of ≥97.9% for all analytes. The QIAstat-Dx RP is a robust and accurate assay for rapid, comprehensive testing for respiratory pathogens.

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.

Large-scale metagenomic and metatranscriptomic data analyses are often restricted by their gene-centric approach, limiting the ability to understand organismal and community biology. De novo assembly of large and mosaic eukaryotic genomes from complex meta-omics data remains a challenging task, especially in comparison with more straightforward bacterial and archaeal systems. Here, we use a transcriptome reconstruction method based on clustering co-abundant genes across a series of metagenomic samples. We investigated the co-abundance patterns of ~37 million eukaryotic unigenes across 365 metagenomic samples collected during the Tara Oceans expeditions to assess the diversity and functional profiles of marine plankton. We identified ~12,000 co-abundant gene groups (CAGs), encompassing ~7 million unigenes, including 924 metagenomics-based transcriptomes (MGTs, CAGs larger than 500 unigenes). We demonstrated the biological validity of the MGT collection by comparing individual MGTs with available references. We identified several key eukaryotic organisms involved in dimethylsulfoniopropionate (DMSP) biosynthesis and catabolism in different oceanic provinces, thus demonstrating the potential of the MGT collection to provide functional insights on eukaryotic plankton. We established the ability of the MGT approach to capture interspecies associations through the analysis of a nitrogen-fixing haptophyte-cyanobacterial symbiotic association. This MGT collection provides a valuable resource for analyses of eukaryotic plankton in the open ocean by giving access to the genomic content and functional potential of many ecologically relevant eukaryotic species.

Transcription of a chromatin template involves the concerted interaction of many different proteins and protein complexes. Analyses of specific factors showed that these interactions change during stress and upon developmental switches. However, how the binding of multiple factors at any given locus is coordinated has been technically challenging to investigate. Here we used Epi-Decoder in yeast to systematically decode, at one transcribed locus, the chromatin binding changes of hundreds of proteins in parallel upon perturbation of transcription. By taking advantage of improved Epi-Decoder libraries, we observed broad rewiring of local chromatin proteomes following chemical inhibition of RNA polymerase. Rapid reduction of RNA polymerase II binding was accompanied by reduced binding of many other core transcription proteins and gain of chromatin remodelers. In quiescent cells, where strong transcriptional repression is induced by physiological signals, eviction of the core transcriptional machinery was accompanied by the appearance of quiescent cell–specific repressors and rewiring of the interactions of protein-folding factors and metabolic enzymes. These results show that Epi-Decoder provides a powerful strategy for capturing the temporal binding dynamics of multiple chromatin proteins under varying conditions and cell states. The systematic and comprehensive delineation of dynamic local chromatin proteomes will greatly aid in uncovering protein–protein relationships and protein functions at the chromatin template.

A time course experiment is a widely used design in the study of cellular processes such as differentiation or response to stimuli. In this paper, we propose time course regulatory analysis (TimeReg) as a method for the analysis of gene regulatory networks based on paired gene expression and chromatin accessibility data from a time course. TimeReg can be used to prioritize regulatory elements, to extract core regulatory modules at each time point, to identify key regulators driving changes of the cellular state, and to causally connect the modules across different time points. We applied the method to analyze paired chromatin accessibility and gene expression data from a retinoic acid (RA)–induced mouse embryonic stem cells (mESCs) differentiation experiment. The analysis identified 57,048 novel regulatory elements regulating cerebellar development, synapse assembly, and hindbrain morphogenesis, which substantially extended our knowledge of cis-regulatory elements during differentiation. Using single-cell RNA-seq data, we showed that the core regulatory modules can reflect the properties of different subpopulations of cells. Finally, the driver regulators are shown to be important in clarifying the relations between modules across adjacent time points. As a second example, our method on Ascl1-induced direct reprogramming from fibroblast to neuron time course data identified Id1/2 as driver regulators of early stage of reprogramming.

Cellular heterogeneity in gene expression is driven by cellular processes, such as cell cycle and cell-type identity, and cellular environment such as spatial location. The cell cycle, in particular, is thought to be a key driver of cell-to-cell heterogeneity in gene expression, even in otherwise homogeneous cell populations. Recent advances in single-cell RNA-sequencing (scRNA-seq) facilitate detailed characterization of gene expression heterogeneity and can thus shed new light on the processes driving heterogeneity. Here, we combined fluorescence imaging with scRNA-seq to measure cell cycle phase and gene expression levels in human induced pluripotent stem cells (iPSCs). By using these data, we developed a novel approach to characterize cell cycle progression. Although standard methods assign cells to discrete cell cycle stages, our method goes beyond this and quantifies cell cycle progression on a continuum. We found that, on average, scRNA-seq data from only five genes predicted a cell's position on the cell cycle continuum to within 14% of the entire cycle and that using more genes did not improve this accuracy. Our data and predictor of cell cycle phase can directly help future studies to account for cell cycle–related heterogeneity in iPSCs. Our results and methods also provide a foundation for future work to characterize the effects of the cell cycle on expression heterogeneity in other cell types.

Retrospective lineage tracing harnesses naturally occurring mutations in cells to elucidate single cell development. Common single-cell phylogenetic fate mapping methods have utilized highly mutable microsatellite loci found within the human genome. Such methods were limited by the introduction of in vitro noise through polymerase slippage inherent in DNA amplification, which we characterized to be approximately 10–100x higher than the in vivo replication mutation rate. Here, we present RETrace, a method for simultaneously capturing both microsatellites and methylation-informative cytosines to characterize both lineage and cell type, respectively, from the same single cell. An important unique feature of RETrace was the introduction of linear amplification of microsatellites in order to reduce in vitro amplification noise. We further coupled microsatellite capture with single-cell reduced representation bisulfite sequencing (scRRBS), to measure the CpG methylation status on the same cell for cell type inference. When compared to existing retrospective lineage tracing methods, RETrace achieved higher accuracy (88% triplet accuracy from an ex vivo HCT116 tree) at a higher cell division resolution (lowering the required number of cell division difference between single cells by approximately 100 divisions). Simultaneously, RETrace demonstrated the ability to capture on average 150,000 unique CpGs per single cell in order to accurately determine cell type. We further formulated additional developments that would allow high-resolution mapping on microsatellite-stable cells or tissues with RETrace. Overall, we present RETrace as a foundation for multi-omics lineage mapping and cell typing of single cells.

The human immune system relies on highly complex and diverse transcripts and the proteins they encode. These include transcripts encoding human leukocyte antigen (HLA) receptors as well as B cell and T cell receptors (BCR and TCR). Determining which alleles an individual possesses for each HLA gene (high-resolution HLA typing) is essential to establish donor–recipient compatibility in organ and bone marrow transplantations. In turn, the repertoires of millions of unique BCR and TCR transcripts in each individual carry a vast amount of health-relevant information. Both short-read RNA-seq-based HLA typing and BCR/TCR repertoire sequencing (AIRR-seq) currently rely on our incomplete knowledge of the genetic diversity at HLA and BCR/TCR loci. Here, we generated over 10,000,000 full-length cDNA sequences at a median accuracy of 97.9% using our nanopore sequencing-based Rolling Circle Amplification to Concatemeric Consensus (R2C2) protocol. We used this data set to (1) show that deep and accurate full-length cDNA sequencing can be used to provide isoform-level transcriptome analysis for more than 9000 loci, (2) generate accurate sequences of HLA alleles, and (3) extract detailed AIRR data for the analysis of the adaptive immune system. The HLA and AIRR analysis approaches we introduce here are untargeted and therefore do not require prior knowledge of the composition or genetic diversity of HLA and BCR/TCR loci.