A vaccine against congenital human cytomegalovirus (CMV) infection is a major public health priority. Congenital CMV causes substantial long-term morbidity, particularly sensorineural hearing loss (SNHL), in newborns, and the public health impact of this infection on maternal and child health is underrecognized. Although progress toward development of a vaccine has been limited by an incomplete understanding of the correlates of protective immunity for the fetus, knowledge about some of the key components of the maternal immune response necessary for preventing transplacental transmission is accumulating. Moreover, although there have been concerns raised about observations indicating that maternal seropositivity does not fully prevent recurrent maternal CMV infections during pregnancy, it is becoming increasing clear that preconception immunity does confer some measure of protection against both CMV transmission and CMV disease (if transmission occurs) in the newborn infant. Although the immunity to CMV conferred by both infection and vaccination is imperfect, there are encouraging data emerging from clinical trials demonstrating the immunogenicity and potential efficacy of candidate CMV vaccines. In the face of the knowledge that between 20,000 and 30,000 infants are born with congenital CMV in the United States every year, there is an urgent and compelling need to accelerate the pace of vaccine trials. In this minireview, we summarize the status of CMV vaccines in clinical trials and provide a perspective on what would be required for a CMV immunization program to become incorporated into clinical practice.

Reactive immunization with a single-dose cholera vaccine that could rapidly (within days) protect immunologically naive individuals during virgin soil epidemics, when cholera reaches immunologically naive populations that have not experienced cholera for decades, would facilitate cholera control. One dose of attenuated Vibrio cholerae O1 classical Inaba vaccine CVD 103-HgR (Vaxchora) containing ≥2 x 10⁸ CFU induces vibriocidal antibody seroconversion (a correlate of protection) in >90% of U.S. adults. A previous CVD 103-HgR commercial formulation required ≥2 x 10⁹ CFU to elicit high levels of seroconversion in populations in developing countries. We compared the vibriocidal responses of Malians (individuals 18 to 45 years old) randomized to ingest a single ≥2 x 10⁸-CFU standard dose (n = 50) or a ≥2 x 10⁹-CFU high dose (n = 50) of PaxVax CVD 103-HgR with buffer or two doses (n = 50) of Shanchol inactivated cholera vaccine (the immunologic comparator). To maintain blinding, participants were dosed twice 2 weeks apart; CVD 103-HgR recipients ingested placebo 2 weeks before or after ingesting vaccine. Seroconversion (a ≥4-fold vibriocidal titer rise) between the baseline and 14 days after CVD 103-HgR ingestion and following the first and second doses of Shanchol were the main outcomes measured. By day 14 postvaccination, the rates of seroconversion after ingestion of a single standard dose and a high dose of CVD 103-HgR were 71.7% (33/46 participants) and 83.3% (40/48 participants), respectively. The rate of seroconversion following the first dose of Shanchol, 56.0% (28/50 participants), was significantly lower than that following the high dose of CVD 103-HgR (P = 0.003). The vibriocidal geometric mean titer (GMT) of the high dose of CVD 103-HgR exceeded the GMT of the standard dose at day 14 (214 versus 95, P = 0.045) and was ~2-fold higher than the GMT on day 7 and day 14 following the first Shanchol dose (P > 0.05). High-dose CVD 103-HgR is recommended for accelerated evaluation in developing countries to assess its efficacy and practicality in field situations. (This study has been registered at ClinicalTrials.gov under registration no. NCT02145377.)

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 December 2019, a novel coronavirus causing an infectious respiratory disease (COVID-19) was identified, which since then has developed into a pandemic with higher rates of mortality in older individuals and those with underlying medical conditions.¹ Multiple sclerosis (MS) is an immune-mediated neurologic disease which requires long-term treatment with immunotherapies that have been shown to increase the risk of infections.² As a result, there is significant anxiety among patients and neurologists during the pandemic regarding the infection outcome in this patient population. We present a patient with MS treated with fingolimod who was diagnosed with COVID-19 and had a favorable outcome.

Objective

To identify and characterize myeloid cell populations within the CSF of patients with MS and anti-myelin oligodendrocyte glycoprotein (MOG) disorder by high-resolution single-cell gene expression analysis.

Objective

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

Objective

To study the immunomodulatory effect of dimethyl fumarate (DF) on granulocyte macrophage colony-stimulating factor (GM-CSF) production in CD4⁺ T cells in experimental autoimmune encephalomyelitis (EAE) and human peripheral blood mononuclear cells (PBMCs).

Objective

To investigate the immunologic impact of a single cycle of rituximab (RTX) in children and adolescents with immune-mediated disorders, we evaluated B cells and immunoglobulin levels of 20 patients with neuroimmunologic, nephrologic, dermatologic, and rheumatologic disorders treated under recommended guidelines.

Objective

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

Alemtuzumab is approved for the treatment of relapsing-remitting MS and is used off-label for patients with chronic lymphocytic leukemia and as induction and antirejection therapy in kidney transplant recipients.¹ Guillain-Barré syndrome (GBS) or chronic inflammatory demyelinating polyradiculoneuropathy (CIDP) complicating alemtuzumab treatment was reported in 9 patients with hematologic malignancy or MS.^1–3 The risk of GBS or CIDP in solid organ transplant recipients treated with alemtuzumab is unknown.

Stiff person syndrome (SPS) is a rare neurologic disorder characterized by progressively worsening rigidity and spasms of the axial and limb muscles. Dyspnea has been recently recognized as a common symptom in SPS,¹ and life-threatening respiratory crises have been occasionally reported and suspected to be responsible for sudden death in these patients.^2,3 The pathophysiologic mechanisms of these respiratory manifestations remain unclear. Some authors have hypothesized that rigidity and/or spasm of the muscles of the trunk could prevent normal rib cage movements and excursion of the diaphragm.¹

Objective

To determine whether serum neurofilament light chain (sNfL) levels are associated with recent MRI activity in patients with relapsing-remitting MS (RRMS).

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

Rock is a material that is affected by wear, and the curvature of the asperities on a rock joint surface increases with the degree of wear after shearing. Based on the Greenwood and Williamson (GW) model, a new model considering the change of asperity curvature is proposed to explain the wear behaviour of rock joints. First, the shear stiffness formula for a joint surface is derived when the asperity curvature is constant, which shows that the shear stiffness increases with increase of asperity curvature. According to the Mohr–Coulomb criterion, the yield position of a single asperity under normal force and tangential friction force is discussed. Then, the critical normal force for a single asperity at a specific friction coefficient is obtained, which shows that the normal force corresponds to the curvature radius of the asperity. A rough surface model with multi-level curvature radius is proposed. With increase of normal force, the higher-order asperities gradually fail and the curvature radius become larger. A specific pressure value excites a specific radius of curvature, and the larger the pressure, the larger the radius of curvature. The relation between the normal force and the curvature radius is proposed and a shear stiffness formula considering the change of curvature radius of the asperity is derived. The proposed model is verified on the basis of the published experimental results. The calculation results of the proposed model can reflect the test results well: for a given joint surface, with increase in normal force the joint surface gradually becomes smooth; for different joint surfaces, with increase in roughness, the joint surface is more easily smoothed.

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.

The decline or drying up of groundwater sources near a tunnel route is damaging to groundwater users. Therefore, forecasting the impact of a tunnel on nearby groundwater sources is a challenging task in tunnel design. In this study, numerical and analytical approaches were applied to the Qomroud water conveyance tunnel (located in Lorestan province, Iran) to assess the impact of tunnelling on the nearby extraction water wells. Using simulation of groundwater-level fluctuation owing to tunnelling, the drawdown at the well locations was determined. From the drawdowns and using Dupuit's equation, the depletion of well flow rates after tunnelling was estimated. To evaluate the results, observed well flow rates before and after tunnelling were compared with the predicted flow rates. The observed and estimated water well flows (before and after tunnelling) showed a regression factor of 0.64, pointing to satisfactory results

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

The regulation of transposable element (TE) activity by small RNAs is a ubiquitous feature of germlines. However, despite the obvious benefits to the host in terms of ensuring the production of viable gametes and maintaining the integrity of the genomes they carry, it remains controversial whether TE regulation evolves adaptively. We examined the emergence and evolutionary dynamics of repressor alleles after P-elements invaded the Drosophila melanogaster genome in the mid-twentieth century. In many animals including Drosophila, repressor alleles are produced by transpositional insertions into piRNA clusters, genomic regions encoding the Piwi-interacting RNAs (piRNAs) that regulate TEs. We discovered that ~94% of recently collected isofemale lines in the Drosophila melanogaster Genetic Reference Panel (DGRP) contain at least one P-element insertion in a piRNA cluster, indicating that repressor alleles are produced by de novo insertion at an exceptional rate. Furthermore, in our sample of approximately 200 genomes, we uncovered no fewer than 80 unique P-element insertion alleles in at least 15 different piRNA clusters. Finally, we observe no footprint of positive selection on P-element insertions in piRNA clusters, suggesting that the rapid evolution of piRNA-mediated repression in D. melanogaster was driven primarily by mutation. Our results reveal for the first time how the unique genetic architecture of piRNA production, in which numerous piRNA clusters can encode regulatory small RNAs upon transpositional insertion, facilitates the nonadaptive rapid evolution of repression.

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.