The complete sequence and comparative analysis of ape sex chromosomes.

Apes possess two sex chromosomes-the male-specific Y chromosome and the X chromosome, which is present in both males and females. The Y chromosome is crucial for male reproduction, with deletions being linked to infertility1. The X chromosome is vital for reproduction and cognition2. Variation in mating patterns and brain function among apes suggests corresponding differences in their sex chromosomes. However, owing to their repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the methodology developed for the telomere-to-telomere (T2T) human genome, we produced gapless assemblies of the X and Y chromosomes for five great apes (bonobo (Pan paniscus), chimpanzee (Pan troglodytes), western lowland gorilla (Gorilla gorilla gorilla), Bornean orangutan (Pongo pygmaeus) and Sumatran orangutan (Pongo abelii)) and a lesser ape (the siamang gibbon (Symphalangus syndactylus)), and untangled the intricacies of their evolution. Compared with the X chromosomes, the ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements-owing to the accumulation of lineage-specific ampliconic regions, palindromes, transposable elements and satellites. Many Y chromosome genes expand in multi-copy families and some evolve under purifying selection. Thus, the Y chromosome exhibits dynamic evolution, whereas the X chromosome is more stable. Mapping short-read sequencing data to these assemblies revealed diversity and selection patterns on sex chromosomes of more than 100 individual great apes. These reference assemblies are expected to inform human evolution and conservation genetics of non-human apes, all of which are endangered species.

Genomic-transcriptomic evolution in lung cancer and metastasis.

Intratumour heterogeneity (ITH) fuels lung cancer evolution, which leads to immune evasion and resistance to therapy1. Here, using paired whole-exome and RNA sequencing data, we investigate intratumour transcriptomic diversity in 354 non-small cell lung cancer tumours from 347 out of the first 421 patients prospectively recruited into the TRACERx study2,3. Analyses of 947 tumour regions, representing both primary and metastatic disease, alongside 96 tumour-adjacent normal tissue samples implicate the transcriptome as a major source of phenotypic variation. Gene expression levels and ITH relate to patterns of positive and negative selection during tumour evolution. We observe frequent copy number-independent allele-specific expression that is linked to epigenomic dysfunction. Allele-specific expression can also result in genomic-transcriptomic parallel evolution, which converges on cancer gene disruption. We extract signatures of RNA single-base substitutions and link their aetiology to the activity of the RNA-editing enzymes ADAR and APOBEC3A, thereby revealing otherwise undetected ongoing APOBEC activity in tumours. Characterizing the transcriptomes of primary-metastatic tumour pairs, we combine multiple machine-learning approaches that leverage genomic and transcriptomic variables to link metastasis-seeding potential to the evolutionary context of mutations and increased proliferation within primary tumour regions. These results highlight the interplay between the genome and transcriptome in influencing ITH, lung cancer evolution and metastasis.

The evolution of non-small cell lung cancer metastases in TRACERx.

Metastatic disease is responsible for the majority of cancer-related deaths1. We report the longitudinal evolutionary analysis of 126 non-small cell lung cancer (NSCLC) tumours from 421 prospectively recruited patients in TRACERx who developed metastatic disease, compared with a control cohort of 144 non-metastatic tumours. In 25% of cases, metastases diverged early, before the last clonal sweep in the primary tumour, and early divergence was enriched for patients who were smokers at the time of initial diagnosis. Simulations suggested that early metastatic divergence more frequently occurred at smaller tumour diameters (less than 8 mm). Single-region primary tumour sampling resulted in 83% of late divergence cases being misclassified as early, highlighting the importance of extensive primary tumour sampling. Polyclonal dissemination, which was associated with extrathoracic disease recurrence, was found in 32% of cases. Primary lymph node disease contributed to metastatic relapse in less than 20% of cases, representing a hallmark of metastatic potential rather than a route to subsequent recurrences/disease progression. Metastasis-seeding subclones exhibited subclonal expansions within primary tumours, probably reflecting positive selection. Our findings highlight the importance of selection in metastatic clone evolution within untreated primary tumours, the distinction between monoclonal versus polyclonal seeding in dictating site of recurrence, the limitations of current radiological screening approaches for early diverging tumours and the need to develop strategies to target metastasis-seeding subclones before relapse.

The evolution of lung cancer and impact of subclonal selection in TRACERx.

Lung cancer is the leading cause of cancer-associated mortality worldwide1. Here we analysed 1,644 tumour regions sampled at surgery or during follow-up from the first 421 patients with non-small cell lung cancer prospectively enrolled into the TRACERx study. This project aims to decipher lung cancer evolution and address the primary study endpoint: determining the relationship between intratumour heterogeneity and clinical outcome. In lung adenocarcinoma, mutations in 22 out of 40 common cancer genes were under significant subclonal selection, including classical tumour initiators such as TP53 and KRAS. We defined evolutionary dependencies between drivers, mutational processes and whole genome doubling (WGD) events. Despite patients having a history of smoking, 8% of lung adenocarcinomas lacked evidence of tobacco-induced mutagenesis. These tumours also had similar detection rates for EGFR mutations and for RET, ROS1, ALK and MET oncogenic isoforms compared with tumours in never-smokers, which suggests that they have a similar aetiology and pathogenesis. Large subclonal expansions were associated with positive subclonal selection. Patients with tumours harbouring recent subclonal expansions, on the terminus of a phylogenetic branch, had significantly shorter disease-free survival. Subclonal WGD was detected in 19% of tumours, and 10% of tumours harboured multiple subclonal WGDs in parallel. Subclonal, but not truncal, WGD was associated with shorter disease-free survival. Copy number heterogeneity was associated with extrathoracic relapse within 1 year after surgery. These data demonstrate the importance of clonal expansion, WGD and copy number instability in determining the timing and patterns of relapse in non-small cell lung cancer and provide a comprehensive clinical cancer evolutionary data resource.

Antibodies against endogenous retroviruses promote lung cancer immunotherapy.

B cells are frequently found in the margins of solid tumours as organized follicles in ectopic lymphoid organs called tertiary lymphoid structures (TLS)1,2. Although TLS have been found to correlate with improved patient survival and response to immune checkpoint blockade (ICB), the underlying mechanisms of this association remain elusive1,2. Here we investigate lung-resident B cell responses in patients from the TRACERx 421 (Tracking Non-Small-Cell Lung Cancer Evolution Through Therapy) and other lung cancer cohorts, and in a recently established immunogenic mouse model for lung adenocarcinoma3. We find that both human and mouse lung adenocarcinomas elicit local germinal centre responses and tumour-binding antibodies, and further identify endogenous retrovirus (ERV) envelope glycoproteins as a dominant anti-tumour antibody target. ERV-targeting B cell responses are amplified by ICB in both humans and mice, and by targeted inhibition of KRAS(G12C) in the mouse model. ERV-reactive antibodies exert anti-tumour activity that extends survival in the mouse model, and ERV expression predicts the outcome of ICB in human lung adenocarcinoma. Finally, we find that effective immunotherapy in the mouse model requires CXCL13-dependent TLS formation. Conversely, therapeutic CXCL13 treatment potentiates anti-tumour immunity and synergizes with ICB. Our findings provide a possible mechanistic basis for the association of TLS with immunotherapy response.

Tracking early lung cancer metastatic dissemination in TRACERx using ctDNA.

Circulating tumour DNA (ctDNA) can be used to detect and profile residual tumour cells persisting after curative intent therapy1. The study of large patient cohorts incorporating longitudinal plasma sampling and extended follow-up is required to determine the role of ctDNA as a phylogenetic biomarker of relapse in early-stage non-small-cell lung cancer (NSCLC). Here we developed ctDNA methods tracking a median of 200 mutations identified in resected NSCLC tissue across 1,069 plasma samples collected from 197 patients enrolled in the TRACERx study2. A lack of preoperative ctDNA detection distinguished biologically indolent lung adenocarcinoma with good clinical outcome. Postoperative plasma analyses were interpreted within the context of standard-of-care radiological surveillance and administration of cytotoxic adjuvant therapy. Landmark analyses of plasma samples collected within 120 days after surgery revealed ctDNA detection in 25% of patients, including 49% of all patients who experienced clinical relapse; 3 to 6 monthly ctDNA surveillance identified impending disease relapse in an additional 20% of landmark-negative patients. We developed a bioinformatic tool (ECLIPSE) for non-invasive tracking of subclonal architecture at low ctDNA levels. ECLIPSE identified patients with polyclonal metastatic dissemination, which was associated with a poor clinical outcome. By measuring subclone cancer cell fractions in preoperative plasma, we found that subclones seeding future metastases were significantly more expanded compared with non-metastatic subclones. Our findings will support (neo)adjuvant trial advances and provide insights into the process of metastatic dissemination using low-ctDNA-level liquid biopsy.

Glial expression of a steroidogenic enzyme underlies natural variation in hitchhiking behavior.

Phoresy is an interspecies interaction that facilitates spatial dispersal by attaching to a more mobile species. Hitchhiking species have evolved specific traits for physical contact and successful phoresy, but the regulatory mechanisms involved in such traits and their evolution are largely unexplored. The nematode Caenorhabditis elegans displays a hitchhiking behavior known as nictation during its stress-induced developmental stage. Dauer-specific nictation behavior has an important role in natural C. elegans populations, which experience boom-and-bust population dynamics. In this study, we investigated the nictation behavior of 137 wild C. elegans strains sampled throughout the world. We identified species-wide natural variation in nictation and performed a genome-wide association mapping. We show that the variants in the promoter of nta-1, encoding a putative steroidogenic enzyme, underlie differences in nictation. This difference is due to the changes in nta-1 expression in glial cells, which implies that glial steroid metabolism regulates phoretic behavior. Population genetic analysis and geographic distribution patterns suggest that balancing selection maintained two nta-1 haplotypes that existed in ancestral C. elegans populations. Our findings contribute to further understanding of the molecular mechanism of species interaction and the maintenance of genetic diversity within natural populations.

Metastatic triple negative breast cancer adapts its metabolism to destination tissues while retaining key metabolic signatures.

Triple negative breast cancer (TNBC) metastases are assumed to exhibit similar functions in different organs as in the original primary tumor. However, studies of metastasis are often limited to a comparison of metastatic tumors with primary tumors of their origin, and little is known about the adaptation to the local environment of the metastatic sites. We therefore used transcriptomic data and metabolic network analyses to investigate whether metastatic tumors adapt their metabolism to the metastatic site and found that metastatic tumors adopt a metabolic signature with some similarity to primary tumors of their destinations. The extent of adaptation, however, varies across different organs, and metastatic tumors retain metabolic signatures associated with TNBC. Our findings suggest that a combination of anti-metastatic approaches and metabolic inhibitors selected specifically for different metastatic sites, rather than solely targeting TNBC primary tumors, may constitute a more effective treatment approach.

Vaccinating humanitarian workers against COVID-19.

Objective: To describe the United Nations' (UN's) coronavirus disease 2019 (COVID-19) vaccination programme and its efforts to vaccinate frontline humanitarian personnel stationed in locations where access to COVID-19 vaccine was limited or absent. Methods: The vaccination programme was structured as a two-level operation: a global vaccine deployment support team and local vaccine deployment teams in each participating country, territory or administrative area. The central group, led by a global vaccine coordinator, oversaw medical, legal, financial, logistical, data, technological and communication aspects. Local vaccine deployment teams were led by coordinators who managed registration, logistics, communication and vaccine administration. The programme used World Health Organization-approved COVID-19 vaccines and developed prioritization criteria for distributing vaccine supplies. The programme ensured that vaccines for the UN personnel were not diverted from the populations they were intended to serve. Findings: The programme successfully formed 120 deployment teams across 152 eligible countries, territories and administrative areas, targeting approximately 673 000 individuals. By April 2023, 72 countries, territories and administrative areas had received over 470 000 doses, of which 337 072 doses were administered. Almost half of the doses administered (167 616) were to individuals in five UN hardship countries. Ninety-five severe adverse events were reported, but none led to any reported medical evacuation, permanent disability or death. Conclusion: The programme demonstrated effective global coordination and local implementation, adapting to diverse contexts and operational challenges. The model can serve as a guide for global actors for future health emergencies, or for deploying health aid at a regional or global scale.

Iterative Regression of Corrective Baselines (IRCB): A New Model for Quantitative Spectroscopy.

In this work, a new model with broad utility for quantitative spectroscopy development is reported. A primary objective of this work is to create a novel modeling procedure that may allow for higher automation of the model development process. The fundamental concept is simple yet powerful even for complex spectra and is employed with no additional preprocessing. This approach is applicable for several types of spectroscopic data to develop regression models that have similar or greater quality than the current methods. The key modeling steps are a matrix transformation and subsequent feature selection process that are collectively referred to as iterative regression of corrective baselines (IRCB). The transformed matrix (Xtransform) is a linearized form of the original X data set. Features from Xtransform that are predictive of Y can be ranked and selected by ordinary least-squares regression. The best features (rows of Xtransform) are linear depictions of Y that can be utilized to develop regression models with several machine learning models. The IRCB workflow is first detailed by using a case study of Fourier transform infrared (FTIR) spectroscopy for prepared solutions of a three-component mixture. Next, IRCB is applied and compared to benchmark results for the 2006 "Chimiométrie" near-infrared spectroscopy (NIR) soil composition challenge and Raman measurements of a simulated nuclear waste slurry.

Extensive Local Geographic Variation in Locoweed Toxin Produced by a Fungal Endophyte.

Legumes are notorious for coevolutionary arms races where chemical defenses are employed to ward off herbivores-particularly insect seed predators. Locoweeds are legumes containing the toxic alkaloid swainsonine which can poison livestock, but its role as a deterrent for insects is unknown. Swainsonine is produced by the fungal endophyte Alternaria section Undifilum, and the chemical composition of the toxin has been well characterized. Despite this knowledge, the ecological roles and evolutionary drivers of swainsonine toxins in locoweeds remain uncertain. Here, we quantify swainsonine concentrations and herbivory levels in the hyper-diverse locoweed Astragalus lentiginosus to evaluate its role as an evolved chemical defense. We found that A. lentiginosus shows considerable variation in swainsonine concentrations according to variety, in particular showing presence/absence variation at both population and local geographic scales. Surprisingly, herbivory levels from presumed generalist insects emerging from fruits showed no correlation with swainsonine concentrations. Conversely, seed and fruit herbivory levels linked to specialist Acanthoscelides seed beetles increased with concentrations of swainsonine-suggesting a possible coevolutionary arms race. Our results highlight that variation in endophyte-produced toxin systems may not follow classical expectations for geographic variation and ecological roles of plant chemicals. We discuss the implications of these results on plant-endophytic toxin systems and coevolutionary dynamics more broadly, highlighting a considerable need for more research in these systems.

Cystic Fibrosis Reprograms Airway Epithelial IL-33 Release and Licenses IL-33-Dependent Inflammation.

Rationale: Type 2 inflammation has been described in people with cystic fibrosis (CF). Whether loss of CFTR (cystic fibrosis transmembrane conductance regulator) function contributes directly to a type 2 inflammatory response has not been fully defined. Objectives: The potent alarmin IL-33 has emerged as a critical regulator of type 2 inflammation. We tested the hypothesis that CFTR deficiency increases IL-33 expression and/or release and deletion of IL-33 reduces allergen-induced inflammation in the CF lung. Methods: Human airway epithelial cells (AECs) grown from non-CF and CF cell lines and Cftr+/+ and Cftr-/- mice were used in this study. Pulmonary inflammation in Cftr+/+ and Cftr-/- mice with and without IL-33 or ST2 (IL-1 receptor-like 1) germline deletion was determined by histological analysis, BAL, and cytokine analysis. Measurements and Main Results: After allergen challenge, both CF human AECs and Cftr-/- mice had increased IL-33 expression compared with control AECs and Cftr+/+ mice, respectively. DUOX1 (dual oxidase 1) expression was increased in CF human AECs and Cftr-/- mouse lungs compared with control AECs and lungs from Cftr+/+ mice and was necessary for the increased IL-33 release in Cftr-/- mice compared with Cftr+/+ mice. IL-33 stimulation of Cftr-/- CD4+ T cells resulted in increased type 2 cytokine production compared with Cftr+/+ CD4+ T cells. Deletion of IL-33 or ST2 decreased both type 2 inflammation and neutrophil recruitment in Cftr-/- mice compared with Cftr+/+ mice. Conclusions: Absence of CFTR reprograms airway epithelial IL-33 release and licenses IL-33-dependent inflammation. Modulation of the IL-33/ST2 axis represents a novel therapeutic target in CF type 2-high and neutrophilic inflammation.

Are Patient-Centered Care, Healthcare Consumerism, and Trust in Physicians Compatible?: Positioning Analysis of the Patient-Provider Relationship.

Patient-centered care and healthcare consumerism are dominant models of the patient-provider relationship. Positioning theory was applied to examine consumers' and physicians' positions on patient-centered care and healthcare consumerism, along with their attitudes toward direct-to-consumer healthcare service advertising and trust in the medical profession. Surveys were conducted with a convenience sample of consumers and physicians respectively. Patient-centered care was the only theoretical construct that both consumers and physicians unequivocally embraced. Both groups were either ambivalent or skeptical of the other three concepts. Between the two groups, physicians exhibited a stronger endorsement of patient-centered care and more negative attitudes toward advertising than consumers. When the relationships among the theoretical constructs were examined, a negative correlation between patient-centered care and consumerism was found among consumers. Also, patient-centered care and trust were negatively correlated in both groups. Implications of these findings are discussed for strategic communication, consumer and physician education, and future research.

Higher Dose Noninvasive Transcutaneous Auricular Vagus Nerve Stimulation Increases Feeding Volumes and White Matter Microstructural Complexity in Open-Label Study of Infants Slated for Gastrostomy Tube.

OBJECTIVE: To determine whether transcutaneous auricular vagus nerve stimulation (taVNS) paired with twice daily bottle feeding increases the volume of oral feeds and white matter neuroplasticity in term-age-equivalent infants failing oral feeds and determined to need a gastrostomy tube. STUDY DESIGN: In this prospective, open-label study, 21 infants received taVNS paired with 2 bottle feeds for 2 - 3 weeks (2x). We compared 1) increase oral feeding volumes with 2x taVNS and previously reported once daily taVNS (1x) to determine a dose response, 2) number of infants who attained full oral feeding volumes, and 3) diffusional kurtosis imaging and magnetic resonance spectroscopy before and after treatment by paired t tests. RESULTS: All 2x taVNS treated infants significantly increased their feeding volumes compared with 10 days before treatment. Over 50% of 2x taVNS infants achieved full oral feeds but in a shorter time than 1x cohort (median 7 days [2x], 12.5 days [1x], P < .05). Infants attaining full oral feeds showed greater increase in radial kurtosis in the right corticospinal tract at the cerebellar peduncle and external capsule. Notably, 75% of infants of diabetic mothers failed full oral feeds, and their glutathione concentrations in the basal ganglia, a measure of central nervous system oxidative stress, were significantly associated with feeding outcome. CONCLUSIONS: In infants with feeding difficulty, increasing the number of daily taVNS-paired feeding sessions to twice-daily significantly accelerates response time but not the overall response rate of treatment. taVNS was associated with white matter motor tract plasticity in infants able to attain full oral feeds. TRIAL REGISTRATION: Clinicaltrials.gov (NCT04643808).

Phylogenetic patterns of bioactive secondary metabolites produced by fungal endosymbionts in morning glories (Ipomoeeae, Convolvulaceae).

Heritable fungal endosymbiosis is underinvestigated in plant biology and documented in only three plant families (Convolvulaceae, Fabaceae, and Poaceae). An estimated 40% of morning glory species in the tribe Ipomoeeae (Convolvulaceae) have associations with one of two distinct heritable, endosymbiotic fungi (Periglandula and Chaetothyriales) that produce the bioactive metabolites ergot alkaloids, indole diterpene alkaloids, and swainsonine, which have been of interest for their toxic effects on animals and potential medical applications. Here, we report the occurrence of ergot alkaloids, indole diterpene alkaloids, and swainsonine in the Convolvulaceae; and the fungi that produce them based on synthesis of previous studies and new indole diterpene alkaloid data from 27 additional species in a phylogenetic, geographic, and life-history context. We find that individual morning glory species host no more than one metabolite-producing fungal endosymbiont (with one possible exception), possibly due to costs to the host and overlapping functions of the alkaloids. The symbiotic morning glory lineages occur in distinct phylogenetic clades, and host species have significantly larger seed size than nonsymbiotic species. The distinct and widely distributed endosymbiotic relationships in the morning glory family and their alkaloids provide an accessible study system for understanding heritable plant-fungal symbiosis evolution and their potential functions for host plants.

Long-read sequencing reveals intra-species tolerance of substantial structural variations and new subtelomere formation in C. elegans.

Long-read sequencing technologies have contributed greatly to comparative genomics among species and can also be applied to study genomics within a species. In this study, to determine how substantial genomic changes are generated and tolerated within a species, we sequenced a C. elegans strain, CB4856, which is one of the most genetically divergent strains compared to the N2 reference strain. For this comparison, we used the Pacific Biosciences (PacBio) RSII platform (80×, N50 read length 11.8 kb) and generated de novo genome assembly to the level of pseudochromosomes containing 76 contigs (N50 contig = 2.8 Mb). We identified structural variations that affected as many as 2694 genes, most of which are at chromosome arms. Subtelomeric regions contained the most extensive genomic rearrangements, which even created new subtelomeres in some cases. The subtelomere structure of Chromosome VR implies that ancestral telomere damage was repaired by alternative lengthening of telomeres even in the presence of a functional telomerase gene and that a new subtelomere was formed by break-induced replication. Our study demonstrates that substantial genomic changes including structural variations and new subtelomeres can be tolerated within a species, and that these changes may accumulate genetic diversity within a species.

Local adaptation and spatiotemporal patterns of genetic diversity revealed by repeated sampling of Caenorhabditis elegans across the Hawaiian Islands.

The nematode Caenorhabditis elegans is among the most widely studied organisms, but relatively little is known about its natural ecology. Genetic diversity is low across much of the globe but high in the Hawaiian Islands and across the Pacific Rim. To characterize the niche and genetic diversity of C. elegans on the Hawaiian Islands and to explore how genetic diversity might be influenced by local adaptation, we repeatedly sampled nematodes over a three-year period, measured various environmental parameters at each sampling site, and whole-genome sequenced the C. elegans isolates that we identified. We found that the typical Hawaiian C. elegans niche comprises moderately moist native forests at high elevations (500-1,500 m) where ambient air temperatures are cool (15-20°C). Compared to other Caenorhabditis species found on the Hawaiian Islands (e.g., Caenorhabditis briggsae and Caenorhabditis tropicalis), we found that C. elegans were enriched in native habitats. We measured levels of genetic diversity and differentiation among Hawaiian C. elegans and found evidence of seven genetically distinct groups distributed across the islands. Then, we scanned these genomes for signatures of local adaptation and identified 18 distinct regions that overlap with hyper-divergent regions, which may be maintained by balancing selection and are enriched for genes related to environmental sensing, xenobiotic detoxification, and pathogen resistance. These results provide strong evidence of local adaptation among Hawaiian C. elegans and contribute to our understanding of the forces that shape genetic diversity on the most remote volcanic archipelago in the world.

Early Lung Disease Exhibits Bacteria-Dependent and -Independent Abnormalities in Cystic Fibrosis Pigs.

Rationale: Although it is clear that cystic fibrosis (CF) airway disease begins at a very young age, the early and subsequent steps in disease pathogenesis and the relative contribution of infection, mucus, and inflammation are not well understood. Objectives: As one approach to assessing the early contribution of infection, we tested the hypothesis that early and continuous antibiotics would decrease the airway bacterial burden. We believed that, if they do, this might reveal aspects of the disease that are more or less sensitive to decreasing infection. Methods: Three groups of pigs were studied from birth until ∼3 weeks of age: 1) wild-type, 2) CF, and 3) CF pigs treated continuously with broad-spectrum antibiotics from birth until study completion. Disease was assessed with chest computed tomography, histopathology, microbiology, and BAL. Measurements and Main Results: Disease was present by 3 weeks of age in CF pigs. Continuous antibiotics from birth improved chest computed tomography imaging abnormalities and airway mucus accumulation but not airway inflammation in the CF pig model. However, reducing bacterial infection did not improve two disease features already present at birth in CF pigs: air trapping and submucosal gland duct plugging. In the CF sinuses, antibiotics did not prevent the development of infection or disease or the number of bacteria but did alter the bacterial species. Conclusions: These findings suggest that CF airway disease begins immediately after birth and that early and continuous antibiotics impact some, but not all, aspects of CF lung disease development.

Harmonizing semantic annotations for computational models in biology.

Life science researchers use computational models to articulate and test hypotheses about the behavior of biological systems. Semantic annotation is a critical component for enhancing the interoperability and reusability of such models as well as for the integration of the data needed for model parameterization and validation. Encoded as machine-readable links to knowledge resource terms, semantic annotations describe the computational or biological meaning of what models and data represent. These annotations help researchers find and repurpose models, accelerate model composition and enable knowledge integration across model repositories and experimental data stores. However, realizing the potential benefits of semantic annotation requires the development of model annotation standards that adhere to a community-based annotation protocol. Without such standards, tool developers must account for a variety of annotation formats and approaches, a situation that can become prohibitively cumbersome and which can defeat the purpose of linking model elements to controlled knowledge resource terms. Currently, no consensus protocol for semantic annotation exists among the larger biological modeling community. Here, we report on the landscape of current annotation practices among the COmputational Modeling in BIology NEtwork community and provide a set of recommendations for building a consensus approach to semantic annotation.