Search

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., Competing Interests: Competing interests statement:The authors declare no competing interest.

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 infertility 1 . The X chromosome is vital for reproduction and cognition 2 . 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., (© 2024. The Author(s).)

Patient-derived xenograft (PDX) models are widely used in cancer research. To investigate the genomic fidelity of non-small cell lung cancer PDX models, we established 48 PDX models from 22 patients enrolled in the TRACERx study. Multi-region tumor sampling increased successful PDX engraftment and most models were histologically similar to their parent tumor. Whole-exome sequencing enabled comparison of tumors and PDX models and we provide an adapted mouse reference genome for improved removal of NOD scid gamma (NSG) mouse-derived reads from sequencing data. PDX model establishment caused a genomic bottleneck, with models often representing a single tumor subclone. While distinct tumor subclones were represented in independent models from the same tumor, individual PDX models did not fully recapitulate intratumor heterogeneity. On-going genomic evolution in mice contributed modestly to the genomic distance between tumors and PDX models. Our study highlights the importance of considering primary tumor heterogeneity when using PDX models and emphasizes the benefit of comprehensive tumor sampling., (© 2024. The Author(s).)

This comprehensive narrative review delves into the intricate interplay between diet and inflammatory bowel disease (IBD), shedding light on the potential impact of dietary interventions in disease management. By analyzing nutritional interventions, risks, challenges, and future perspectives, this review serves as a vital resource for clinicians, researchers, and patients alike. The amalgamation of evidence underscores the significance of customizing dietary strategies for individual patients, considering disease phenotype and cultural factors. Through an exploration of dietary components' effects on IBD, including exclusive enteral nutrition and omega-3 fatty acids, this review offers pragmatic implementation advice and outlines avenues for further research. Bridging the gap between research findings and clinical applications, the review facilitates informed decision-making and patient-centric care. In the face of escalating IBD prevalence, this review emerges as an indispensable guide for healthcare professionals, empowering them to navigate the complexities of dietary management while enabling patients to actively participate in their care trajectory. Ultimately, this narrative review advances the understanding of diet's pivotal role in IBD management, fostering a more integrated approach to patient care and paving the way for improved research and policy initiatives in the field of inflammatory bowel diseases., Competing Interests: The authors have declared that no competing interests exist., (Copyright © 2023, Saeed et al.)

Summary: RNA sequencing (RNA-seq) can be applied to diverse tasks including quantifying gene expression, discovering quantitative trait loci and identifying gene fusion events. Although RNA-seq can detect germline variants, the complexities of variable transcript abundance, target capture and amplification introduce challenging sources of error. Here, we extend DeepVariant, a deep-learning-based variant caller, to learn and account for the unique challenges presented by RNA-seq data. Our DeepVariant RNA-seq model produces highly accurate variant calls from RNA-sequencing data, and outperforms existing approaches such as Platypus and GATK. We examine factors that influence accuracy, how our model addresses RNA editing events and how additional thresholding can be used to facilitate our models' use in a production pipeline., Supplementary Information: Supplementary data are available at Bioinformatics Advances online., (© The Author(s) 2023. Published by Oxford University Press.)

Here the Human Pangenome Reference Consortium presents a first draft of the human pangenome reference. The pangenome contains 47 phased, diploid assemblies from a cohort of genetically diverse individuals 1 . These assemblies cover more than 99% of the expected sequence in each genome and are more than 99% accurate at the structural and base pair levels. Based on alignments of the assemblies, we generate a draft pangenome that captures known variants and haplotypes and reveals new alleles at structurally complex loci. We also add 119 million base pairs of euchromatic polymorphic sequences and 1,115 gene duplications relative to the existing reference GRCh38. Roughly 90 million of the additional base pairs are derived from structural variation. Using our draft pangenome to analyse short-read data reduced small variant discovery errors by 34% and increased the number of structural variants detected per haplotype by 104% compared with GRCh38-based workflows, which enabled the typing of the vast majority of structural variant alleles per sample., (© 2023. The Author(s).)

Metastatic disease is responsible for the majority of cancer-related deaths 1 . 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., (© 2023. The Author(s).)

Intratumour heterogeneity (ITH) fuels lung cancer evolution, which leads to immune evasion and resistance to therapy 1 . 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 study 2,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., (© 2023. The Author(s).)

Circulating tumour DNA (ctDNA) can be used to detect and profile residual tumour cells persisting after curative intent therapy 1 . 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 study 2 . 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., (© 2023. The Author(s).)

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 elusive 1,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 adenocarcinoma 3 . 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., (© 2023. The Author(s).)

Lung cancer is the leading cause of cancer-associated mortality worldwide 1 . 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., (© 2023. The Author(s).)

Circular consensus sequencing with Pacific Biosciences (PacBio) technology generates long (10-25 kilobases), accurate 'HiFi' reads by combining serial observations of a DNA molecule into a consensus sequence. The standard approach to consensus generation, pbccs, uses a hidden Markov model. We introduce DeepConsensus, which uses an alignment-based loss to train a gap-aware transformer-encoder for sequence correction. Compared to pbccs, DeepConsensus reduces read errors by 42%. This increases the yield of PacBio HiFi reads at Q20 by 9%, at Q30 by 27% and at Q40 by 90%. With two SMRT Cells of HG003, reads from DeepConsensus improve hifiasm assembly contiguity (NG50 4.9 megabases (Mb) to 17.2 Mb), increase gene completeness (94% to 97%), reduce the false gene duplication rate (1.1% to 0.5%), improve assembly base accuracy (Q43 to Q45) and reduce variant-calling errors by 24%. DeepConsensus models could be trained to the general problem of analyzing the alignment of other types of sequences, such as unique molecular identifiers or genome assemblies., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Chemotaxis assays are used extensively to study behavioral responses of Caenorhabditis nematodes to environmental cues. These assays result in a chemotaxis index (CI) that denotes the behavioral response of a population of nematodes to a particular compound and can range from 1 (maximum attraction) to -1 (maximum avoidance). Traditional chemotaxis assays have low throughput because researchers must manually setup experimental populations and score CIs. Here, we describe an automated methodology that increases throughput by using liquid-handling robots to setup experimental populations and a custom image analysis package, ct, to automate the scoring of CIs from plate images., (Copyright: © 2022 by the authors.)

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., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

The distribution of fitness effects (DFE) for new mutations is one of the most theoretically important but difficult to estimate properties in population genetics. A crucial challenge to inferring the DFE from natural genetic variation is the sensitivity of the site frequency spectrum to factors like population size change, population substructure, genome structure, and nonrandom mating. Although inference methods aim to control for population size changes, the influence of nonrandom mating remains incompletely understood, despite being a common feature of many species. We report the DFE estimated from 326 genomes of Caenorhabditis elegans, a nematode roundworm with a high rate of self-fertilization. We evaluate the robustness of DFE inferences using simulated data that mimics the genomic structure and reproductive life history of C. elegans. Our observations demonstrate how the combined influence of self-fertilization, genome structure, and natural selection on linked sites can conspire to compromise estimates of the DFE from extant polymorphisms with existing methods. These factors together tend to bias inferences toward weakly deleterious mutations, making it challenging to have full confidence in the inferred DFE of new mutations as deduced from standing genetic variation in species like C. elegans. Improved methods for inferring the DFE are needed to appropriately handle strong linked selection and selfing. These results highlight the importance of understanding the combined effects of processes that can bias our interpretations of evolution in natural populations., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Large-scale ecological sampling can be difficult and costly, especially for organisms that are too small to be easily identified in a natural environment by eye. Typically, these microscopic floral and fauna are sampled by collecting substrates from nature and then separating organisms from substrates in the laboratory. In many cases, diverse organisms can be identified to the species-level using molecular barcodes. To facilitate large-scale ecological sampling of microscopic organisms, we used a geographic data-collection platform for mobile devices called Fulcrum that streamlines the organization of geospatial sampling data, substrate photographs, and environmental data at natural sampling sites. These sampling data are then linked to organism isolation data from the laboratory. Here, we describe the easyFulcrum R package, which can be used to clean, process, and visualize ecological field sampling and isolation data exported from the Fulcrum mobile application. We developed this package for wild nematode sampling, but it can be used with other organisms. The advantages of using Fulcrum combined with easyFulcrum are (1) the elimination of transcription errors by replacing manual data entry and/or spreadsheets with a mobile application, (2) the ability to clean, process, and visualize sampling data using a standardized set of functions in the R software environment, and (3) the ability to join disparate data to each other, including environmental data from the field and the molecularly defined identities of individual specimens isolated from samples., Competing Interests: The authors have declared that no competing interests exist.

Across diverse taxa, selfing species have evolved independently from outcrossing species thousands of times. The transition from outcrossing to selfing decreases the effective population size, effective recombination rate and heterozygosity within a species. These changes lead to a reduction in genetic diversity, and therefore adaptive potential, by intensifying the effects of random genetic drift and linked selection. Within the nematode genus Caenorhabditis, selfing has evolved at least three times, and all three species, including the model organism Caenorhabditis elegans, show substantially reduced genetic diversity relative to outcrossing species. Selfing and outcrossing Caenorhabditis species are often found in the same niches, but we still do not know how selfing species with limited genetic diversity can adapt to these environments. Here, we examine the whole-genome sequences from 609 wild C. elegans strains isolated worldwide and show that genetic variation is concentrated in punctuated hyper-divergent regions that cover 20% of the C. elegans reference genome. These regions are enriched in environmental response genes that mediate sensory perception, pathogen response and xenobiotic stress response. Population genomic evidence suggests that genetic diversity in these regions has been maintained by long-term balancing selection. Using long-read genome assemblies for 15 wild strains, we show that hyper-divergent haplotypes contain unique sets of genes and show levels of divergence comparable to levels found between Caenorhabditis species that diverged millions of years ago. These results provide an example of how species can avoid the evolutionary dead end associated with selfing.

Hawaiian isolates of the nematode species Caenorhabditis elegans have long been known to harbor genetic diversity greater than the rest of the worldwide population, but this observation was supported by only a small number of wild strains. To better characterize the niche and genetic diversity of Hawaiian C. elegans and other Caenorhabditis species, we sampled different substrates and niches across the Hawaiian islands. We identified hundreds of new Caenorhabditis strains from known species and a new species, Caenorhabditis oiwi . Hawaiian C. elegans are found in cooler climates at high elevations but are not associated with any specific substrate, as compared to other Caenorhabditis species. Surprisingly, admixture analysis revealed evidence of shared ancestry between some Hawaiian and non-Hawaiian C. elegans strains. We suggest that the deep diversity we observed in Hawaii might represent patterns of ancestral genetic diversity in the C. elegans species before human influence., Competing Interests: TC, SZ, DC, RT, SB, YW, KE, SH, DL, BR, GZ, Jv, KK, EA No competing interests declared, (© 2019, Crombie et al.)

To understand the genetic basis of complex traits, it is important to be able to efficiently phenotype many genetically distinct individuals. In the nematode Caenorhabditis elegans , individuals have been isolated from diverse populations around the globe and whole-genome sequenced. As a result, hundreds of wild strains with known genome sequences can be used for genome-wide association studies (GWAS). However, phenotypic analysis of these strains can be laborious, particularly for quantitative traits requiring multiple measurements per strain. Starvation resistance is likely a fitness-proximal trait for nematodes, and it is related to metabolic disease risk in humans. However, natural variation in C. elegans starvation resistance has not been systematically characterized, and precise measurement of the trait is time-intensive. Here, we developed a population-selection-and-sequencing-based approach to phenotype starvation resistance in a pool of 96 wild strains. We used restriction site-associated DNA sequencing (RAD-seq) to infer the frequency of each strain among survivors in a mixed culture over time during starvation. We used manual starvation survival assays to validate the trait data, confirming that strains that increased in frequency over time are starvation-resistant relative to strains that decreased in frequency. Further, we found that variation in starvation resistance is significantly associated with variation at a region on chromosome III. Using a near-isogenic line (NIL), we showed the importance of this genomic interval for starvation resistance. This study demonstrates the feasibility of using population selection and sequencing in an animal model for phenotypic analysis of quantitative traits, documents natural variation of starvation resistance in C. elegans , and identifies a genomic region that contributes to such variation., (Copyright © 2019 Webster et al.)

From quorum sensing in bacteria to pheromone signalling in social insects, chemical communication mediates interactions among individuals in local populations. In Caenorhabditis elegans, ascaroside pheromones can dictate local population density; high levels of pheromones inhibit the reproductive maturation of individuals. Little is known about how natural genetic diversity affects the pheromone responses of individuals from diverse habitats. Here, we show that a niche-associated variation in pheromone receptor genes contributes to natural differences in pheromone responses. We identified putative loss-of-function deletions that impair duplicated pheromone receptor genes (srg-36 and srg-37), which were previously shown to be lost in population-dense laboratory cultures. A common natural deletion in srg-37 arose recently from a single ancestral population that spread throughout the world; this deletion underlies reduced pheromone sensitivity across the global C. elegans population. We found that many local populations harbour individuals with a wild-type or a deletion allele of srg-37, suggesting that balancing selection has maintained the recent variation in this pheromone receptor gene. The two srg-37 genotypes are associated with niche diversity underlying boom-and-bust population dynamics. We hypothesize that human activities likely contributed to the gene flow and balancing selection of srg-37 variation through facilitating the migration of species and providing a favourable niche for the recently arisen srg-37 deletion.

Bleomycin is a powerful chemotherapeutic drug used to treat a variety of cancers. However, individual patients vary in their responses to bleomycin. The identification of genetic differences that underlie this response variation could improve treatment outcomes by tailoring bleomycin dosages to each patient. We used the model organism Caenorhabditis elegans to identify genetic determinants of bleomycin-response differences by performing linkage mapping on recombinants derived from a cross between the laboratory strain (N2) and a wild strain (CB4856). This approach identified a small genomic region on chromosome V that underlies bleomycin-response variation. Using near-isogenic lines, and strains with CRISPR-Cas9 mediated deletions and allele replacements, we discovered that a novel nematode-specific gene ( scb-1 ) is required for bleomycin resistance. Although the mechanism by which this gene causes variation in bleomycin responses is unknown, we suggest that a rare variant present in the CB4856 strain might cause differences in the potential stress-response function of scb-1 between the N2 and CB4856 strains, thereby leading to differences in bleomycin resistance., (Copyright © 2019 Brady et al.)

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., (© 2019 Kim et al.; Published by Cold Spring Harbor Laboratory Press.)

Phenotypic complexity is caused by the contributions of environmental factors and multiple genetic loci, interacting or acting independently. Studies of yeast and Arabidopsis often find that the majority of natural variation across phenotypes is attributable to independent additive quantitative trait loci (QTL). Detected loci in these organisms explain most of the estimated heritable variation. By contrast, many heritable components underlying phenotypic variation in metazoan models remain undetected. Before the relative impacts of additive and interactive variance components on metazoan phenotypic variation can be dissected, high replication and precise phenotypic measurements are required to obtain sufficient statistical power to detect loci contributing to this missing heritability. Here, we used a panel of 296 recombinant inbred advanced intercross lines of Caenorhabditis elegans and a high-throughput fitness assay to detect loci underlying responses to 16 different toxins, including heavy metals, chemotherapeutic drugs, pesticides, and neuropharmaceuticals. Using linkage mapping, we identified 82 QTL that underlie variation in responses to these toxins, and predicted the relative contributions of additive loci and genetic interactions across various growth parameters. Additionally, we identified three genomic regions that impact responses to multiple classes of toxins. These QTL hotspots could represent common factors impacting toxin responses. We went further to generate near-isogenic lines and chromosome substitution strains, and then experimentally validated these QTL hotspots, implicating additive and interactive loci that underlie toxin-response variation., (Copyright © 2018 Evans et al.)

Parasitic nematodes impose a debilitating health and economic burden across much of the world. Nematode resistance to anthelmintic drugs threatens parasite control efforts in both human and veterinary medicine. Despite this threat, the genetic landscape of potential resistance mechanisms to these critical drugs remains largely unexplored. Here, we exploit natural variation in the model nematodes Caenorhabditis elegans and Caenorhabditis briggsae to discover quantitative trait loci (QTL) that control sensitivity to benzimidazoles widely used in human and animal medicine. High-throughput phenotyping of albendazole, fenbendazole, mebendazole, and thiabendazole responses in panels of recombinant lines led to the discovery of over 15 QTL in C. elegans and four QTL in C. briggsae associated with divergent responses to these anthelmintics. Many of these QTL are conserved across benzimidazole derivatives, but others show drug and dose specificity. We used near-isogenic lines to recapitulate and narrow the C. elegans albendazole QTL of largest effect and identified candidate variants correlated with the resistance phenotype. These QTL do not overlap with known benzimidazole target resistance genes from parasitic nematodes and present specific new leads for the discovery of novel mechanisms of nematode benzimidazole resistance. Analyses of orthologous genes reveal conservation of candidate benzimidazole resistance genes in medically important parasitic nematodes. These data provide a basis for extending these approaches to other anthelmintic drug classes and a pathway towards validating new markers for anthelmintic resistance that can be deployed to improve parasite disease control.

Many chemotherapeutic drugs are differentially effective from one patient to the next. Understanding the causes of this variability is a critical step towards the development of personalized treatments and improvements to existing medications. Here, we investigate sensitivity to a group of anti-neoplastic drugs that target topoisomerase II using the model organism Caenorhabditis elegans. We show that wild strains of C. elegans vary in their sensitivity to these drugs, and we use an unbiased genetic approach to demonstrate that this natural variation is explained by a methionine-to-glutamine substitution in topoisomerase II (TOP-2). The presence of a non-polar methionine at this residue increases hydrophobic interactions between TOP-2 and its poison etoposide, as compared to a polar glutamine. We hypothesize that this stabilizing interaction results in increased genomic instability in strains that contain a methionine residue. The residue affected by this substitution is conserved from yeast to humans and is one of the few differences between the two human topoisomerase II isoforms (methionine in hTOPIIα and glutamine in hTOPIIβ). We go on to show that this amino acid difference between the two human topoisomerase isoforms influences cytotoxicity of topoisomerase II poisons in human cell lines. These results explain why hTOPIIα and hTOPIIβ are differentially affected by various poisons and demonstrate the utility of C. elegans in understanding the genetics of drug responses.

Summary: The variant call format (VCF) is a popular standard for storing genetic variation data. As a result, a large collection of tools has been developed that perform diverse analyses using VCF files. However, some tasks common to statistical and population geneticists have not been created yet. To streamline these types of analyses, we created novel tools that analyze or annotate VCF files and organized these tools into a command-line based utility named VCF-kit. VCF-kit adds essential utilities to process and analyze VCF files, including primer generation for variant validation, dendrogram production, genotype imputation from sequence data in linkage studies, and additional tools., Availability and Implementation: https://github.com/AndersenLab/VCF-kit., Contact: erik.andersen@northwestern.edu., (© The Author 2017. Published by Oxford University Press.)

Studies in model organisms have yielded considerable insights into the etiology of disease and our understanding of evolutionary processes. Caenorhabditis elegans is among the most powerful model organisms used to understand biology. However, C. elegans is not used as extensively as other model organisms to investigate how natural variation shapes traits, especially through the use of genome-wide association (GWA) analyses. Here, we introduce a new platform, the C. elegans Natural Diversity Resource (CeNDR) to enable statistical genetics and genomics studies of C. elegans and to connect the results to human disease. CeNDR provides the research community with wild strains, genome-wide sequence and variant data for every strain, and a GWA mapping portal for studying natural variation in C. elegans Additionally, researchers outside of the C. elegans community can benefit from public mappings and integrated tools for comparative analyses. CeNDR uses several databases that are continually updated through the addition of new strains, sequencing data, and association mapping results. The CeNDR data are accessible through a freely available web portal located at http://www.elegansvariation.org or through an application programming interface., (© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Telomeres are involved in the maintenance of chromosomes and the prevention of genome instability. Despite this central importance, significant variation in telomere length has been observed in a variety of organisms. The genetic determinants of telomere-length variation and their effects on organismal fitness are largely unexplored. Here, we describe natural variation in telomere length across the Caenorhabditis elegans species. We identify a large-effect variant that contributes to differences in telomere length. The variant alters the conserved oligonucleotide/oligosaccharide-binding fold of protection of telomeres 2 (POT-2), a homolog of a human telomere-capping shelterin complex subunit. Mutations within this domain likely reduce the ability of POT-2 to bind telomeric DNA, thereby increasing telomere length. We find that telomere-length variation does not correlate with offspring production or longevity in C. elegans wild isolates, suggesting that naturally long telomeres play a limited role in modifying fitness phenotypes in C. elegans., (Copyright © 2016 by the Genetics Society of America.)

Pre-B and pre-T lymphocytes must orchestrate a transition from a highly proliferative state to a quiescent one during development. Cyclin D3 is essential for these cells' proliferation, but little is known about its posttranslational regulation at this stage. Here, we show that the dual specificity tyrosine-regulated kinase 1A (DYRK1A) restrains Cyclin D3 protein levels by phosphorylating T283 to induce its degradation. Loss of DYRK1A activity, via genetic inactivation or pharmacologic inhibition in mice, caused accumulation of Cyclin D3 protein, incomplete repression of E2F-mediated gene transcription, and failure to properly couple cell cycle exit with differentiation. Expression of a nonphosphorylatable Cyclin D3 T283A mutant recapitulated these defects, whereas inhibition of Cyclin D:CDK4/6 mitigated the effects of DYRK1A inhibition or loss. These data uncover a previously unknown role for DYRK1A in lymphopoiesis, and demonstrate how Cyclin D3 protein stability is negatively regulated during exit from the proliferative phases of B and T cell development., (© 2015 Thompson et al.)

Objective: Measurements of neonatal metabolites are commonly used in newborn screening (NBS) programs to detect inborn errors of metabolism. Variation in these metabolites, particularly in infants born preterm (<37 weeks gestation), can result from multiple etiologies. We sought to evaluate the impact of maternal complications of pregnancy and environmental stressors on NBS metabolites., Methods: We examined 49 metabolic biomarkers obtained from routine NBS in 452 infants born preterm for association with maternal environmental stressors and complications of pregnancy., Results: Neonatal free carnitine (C0, p = 1.4 × 10(-7)), acetylcarnitine (C2, p = 2.7 × 10(-7)), octenoylcarnitine (C8:1, p = 5.2 × 10(-11)) and linoleoylcarnitine (C18:2, p = 9.1 × 10(-7)) were elevated in infants born to preeclamptic mothers. Similar elevations were observed in small for gestational age infants and in infants where labor was not initiated prior to delivery. When accounting for all three factors, associations remained strongest between acylcarnitines and preeclampsia., Conclusion: We observed that maternal conditions, particularly preeclampsia, influence NBS biomarkers. This is important for identifying maternal conditions that influence metabolites measured during routine NBS that are also markers of fetal growth and overall health.

Objectives: Identifying common clinical and environmental factors that influence newborn metabolic biomarkers will improve the utilization of metabolite panels for clinical diagnostic medicine., Design and Methods: Environmental effects including gender, season of birth, gestational age, birth weight, feeding method and age at time of collection were evaluated for over 50 metabolites collected by the Iowa Neonatal Metabolic Screening Program on 221,788 newborns over a six year period., Results: We replicated well known observations that low birth weight and preterm infants have higher essential amino acids and lower medium and long chain acylcarnitine levels than their term counterparts. Smaller, but still significant, differences were observed for gender and timing of sample collection, specifically the season in which the infant was born. Most intriguing were our findings of higher thyroid stimulating hormone in the winter months (P<1×10(-40)) which correlated with an increased false positive rate of congenital hypothyroidism in the winter (0.9%) compared to summer (0.6%). Previous studies, conducted globally, have identified an increased prevalence of suspected and confirmed cases of congenital hypothyroidism in the winter months. We found that the percentage of unresolved suspected cases were slightly higher in the winter (0.3% vs. 0.2%)., Conclusions: We identified differences in metabolites by gestational age, birth weight, gender and season. Some are widely reported such as gestational age and birth weight, while others such as the effect of seasonality are not as well studied., (Copyright © 2012 The Canadian Society of Clinical Chemists. All rights reserved.)

Nationally, newborn screening programs use 17-hydroxyprogesterone (17-OHP) as the biomarker to detect the rare but potentially fatal inherited disease, congenital adrenal hyperplasia. However, this biomarker is highly variable, with a high false-positive rate of detection, particularly in neonates born preterm. Several studies have examined various clinical and genetic factors to explain the variability of 17-OHP in preterm infants. The purpose of this study was to replicate previous clinical and genetic associations with 17-OHP in a well-characterized cohort of 762 preterm infants. We replicated previous findings that respiratory distress syndrome (p = 2 x 10(-3)) is associated with higher 17-OHP. Higher 17-OHP and false positives were significantly associated with lower gestational age and birth weight, as previously reported. Incorporating gestational age and birth weight together decreases the false-positive rate.