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

Transposon-mediated genic rearrangements underlie variation in small RNA pathways

Transposable elements (TEs) can alter host gene structure and expression, whereas host organisms develop mechanisms to repress TE activities. In the nematode Caenorhabditis elegans, a small interfering RNA pathway dependent on the helicase ERI-6/7 primarily silences retrotransposons and recent genes of likely viral origin. By studying gene expression variation among wild C. elegans strains, we found that structural variants and transposon remnants likely underlie expression variation in eri-6/7 and the pathway targets. We further found that multiple insertions of the DNA transposons, Polintons, reshuffled the eri-6/7 locus and induced inversion of eri-6 in some wild strains. In the inverted configuration, gene function was previously shown to be repaired by unusual trans-splicing mediated by direct repeats. We identified that these direct repeats originated from terminal inverted repeats of Polintons. Our findings highlight the role of host-transposon interactions in driving rapid host genome diversification among natural populations and shed light on evolutionary novelty in genes and splicing mechanisms.

Atlas of telomeric repeat diversity in Arabidopsis thaliana

BACKGROUND

Telomeric repeat arrays at the ends of chromosomes are highly dynamic in composition, but their repetitive nature and technological limitations have made it difficult to assess their true variation in genome diversity surveys.

RESULTS

We have comprehensively characterized the sequence variation immediately adjacent to the canonical telomeric repeat arrays at the very ends of chromosomes in 74 genetically diverse Arabidopsis thaliana accessions. We first describe several types of distinct telomeric repeat units and then identify evolutionary processes such as local homogenization and higher-order repeat formation that shape diversity of chromosome ends. By comparing largely isogenic samples, we also determine repeat number variation of the degenerate and variant telomeric repeat array at both the germline and somatic levels. Finally, our analysis of haplotype structure uncovers chromosome end-specific patterns in the distribution of variant telomeric repeats, and their linkage to the more proximal non-coding region.

CONCLUSIONS

Our findings illustrate the spectrum of telomeric repeat variation at multiple levels in A. thaliana-in germline and soma, across all chromosome ends, and across genetic groups-thereby expanding our knowledge of the evolution of chromosome ends.

Identifying transgene insertions in Caenorhabditis elegans genomes with Oxford Nanopore sequencing

Genetically modified organisms are commonly used in disease research and agriculture but the precise genomic alterations underlying transgenic mutations are often unknown. The position and characteristics of transgenes, including the number of independent insertions, influences the expression of both transgenic and wild-type sequences. We used long-read, Oxford Nanopore Technologies (ONT) to sequence and assemble two transgenic strains of Caenorhabditis elegans commonly used in the research of neurodegenerative diseases: BY250 (pPdat-1::GFP) and UA44 (GFP and human α-synuclein), a model for Parkinson's research. After scaffolding to the reference, the final assembled sequences were ∼102 Mb with N50s of 17.9 Mb and 18.0 Mb, respectively, and L90s of six contiguous sequences, representing chromosome-level assemblies. Each of the assembled sequences contained more than 99.2% of the Nematoda BUSCO genes found in the C. elegans reference and 99.5% of the annotated C. elegans reference protein-coding genes. We identified the locations of the transgene insertions and confirmed that all transgene sequences were inserted in intergenic regions, leaving the organismal gene content intact. The transgenic C. elegans genomes presented here will be a valuable resource for Parkinson's research as well as other neurodegenerative diseases. Our work demonstrates that long-read sequencing is a fast, cost-effective way to assemble genome sequences and characterize mutant lines and strains.

Multiplex DNA fluorescence in situ hybridization to analyze maternal vs. paternal C. elegans chromosomes

Recent advances in microscopy have enabled studying chromosome organization at the single-molecule level, yet little is known about inherited chromosome organization. Here we adapt single-molecule chromosome tracing to distinguish two C. elegans strains (N2 and HI) and find that while their organization is similar, the N2 chromosome influences the folding parameters of the HI chromosome, in particular the step size, across generations. Furthermore, homologous chromosomes overlap frequently, but alignment between homologous regions is rare, suggesting that transvection is unlikely. We present a powerful tool to investigate chromosome architecture and to track the parent of origin.

Subtelomeric plasticity contributes to gene family expansion in the human parasitic flatworm Schistosoma mansoni

BACKGROUND

The genomic region that lies between the telomere and chromosome body, termed the subtelomere, is heterochromatic, repeat-rich, and frequently undergoes rearrangement. Within this region, large-scale structural changes enable gene diversification, and, as such, large multicopy gene families are often found at the subtelomere. In some parasites, genes associated with proliferation, invasion, and survival are often found in these regions, where they benefit from the subtelomere's highly plastic, rapidly changing nature. The increasing availability of complete (or near complete) parasite genomes provides an opportunity to investigate these typically poorly defined and overlooked genomic regions and potentially reveal relevant gene families necessary for the parasite's lifestyle.

RESULTS

Using the latest chromosome-scale genome assembly and hallmark repeat richness observed at chromosome termini, we have identified and characterised the subtelomeres of Schistosoma mansoni, a metazoan parasitic flatworm that infects over 250 million people worldwide. Approximately 12% of the S. mansoni genome is classified as subtelomeric, and, in line with other organisms, we find these regions to be gene-poor but rich in transposable elements. We find that S. mansoni subtelomeres have undergone extensive interchromosomal recombination and that these sites disproportionately contribute to the 2.3% of the genome derived from segmental duplications. This recombination has led to the expansion of subtelomeric gene clusters containing 103 genes, including the immunomodulatory annexins and other gene families with unknown roles. The largest of these is a 49-copy plexin domain-containing protein cluster, exclusively expressed in the tegument-the tissue located at the host-parasite physical interface-of intramolluscan life stages.

CONCLUSIONS

We propose that subtelomeric regions act as a genomic playground for trial-and-error of gene duplication and subsequent divergence. Owing to the importance of subtelomeric genes in other parasites, gene families implicated in this subtelomeric expansion within S. mansoni warrant further characterisation for a potential role in parasitism.

A nematode model to evaluate microdeletion phenotype expression

Microdeletion syndromes are genetic diseases caused by multilocus chromosomal deletions too small to be detected by karyotyping. They are typified by complex pleiotropic developmental phenotypes that depend both on the extent of the deletion and variations in genetic background. Microdeletion alleles cause a wide array of consequences involving multiple pathways. How simultaneous haploinsufficiency of numerous adjacent genes leads to complex and variable pleiotropic phenotypes is not well understood. CRISPR/Cas9 genome editing has been shown to induce microdeletion-like alleles at a meaningful rate. Here, we describe a microdeletion allele in Caenorhabditis elegans recovered during a CRISPR/Cas9 genome editing experiment. We mapped the allele to chromosome V, balanced it with a reciprocal translocation crossover suppressor, and precisely defined the breakpoint junction. The allele simultaneously removes 32 protein-coding genes, yet animals homozygous for this mutation are viable as adults. Homozygous animals display a complex phenotype including maternal effect lethality, producing polynucleated embryos that grow into uterine tumors, vulva morphogenesis defects, body wall distensions, uncoordinated movement, and a shortened life span typified by death by bursting. Our work provides an opportunity to explore the complexity and penetrance of microdeletion phenotypes in a simple genetic model system.

Telomeric repeat evolution in the phylum Nematoda revealed by high-quality genome assemblies and subtelomere structures

Telomeres are composed of tandem arrays of telomeric-repeat motifs (TRMs) and telomere-binding proteins (TBPs), which are responsible for ensuring end-protection and end-replication of chromosomes. TRMs are highly conserved owing to the sequence specificity of TBPs, although significant alterations in TRM have been observed in several taxa, except Nematoda. We used public whole-genome sequencing data sets to analyze putative TRMs of 100 nematode species and determined that three distinct branches included specific novel TRMs, suggesting that evolutionary alterations in TRMs occurred in Nematoda. We focused on one of the three branches, the Panagrolaimidae family, and performed a de novo assembly of four high-quality draft genomes of the canonical (TTAGGC) and novel TRM (TTAGAC) isolates; the latter genomes revealed densely clustered arrays of the novel TRM. We then comprehensively analyzed the subtelomeric regions of the genomes to infer how the novel TRM evolved. We identified DNA damage-repair signatures in subtelomeric sequences that were representative of consequences of telomere maintenance mechanisms by alternative lengthening of telomeres. We propose a hypothetical scenario in which TTAGAC-containing units are clustered in subtelomeric regions and pre-existing TBPs capable of binding both canonical and novel TRMs aided the evolution of the novel TRM in the Panagrolaimidae family.

Comprehensive detection of structural variation and transposable element differences between wild type laboratory lineages of C. elegans

Genomic structural variations (SVs) and transposable elements (TEs) can be significant contributors to genome evolution, altered gene expression, and risk of genetic diseases. Recent advancements in long-read sequencing have greatly improved the quality of de novo genome assemblies and enhanced the detection of sequence variants at the scale of hundreds or thousands of bases. Comparisons between two diverged wild isolates of Caenorhabditis elegans, the Bristol and Hawaiian strains, have been widely utilized in the analysis of small genetic variations. Genetic drift, including SVs and rearrangements of repeated sequences such as TEs, can occur over time from long-term maintenance of wild type isolates within the laboratory. To comprehensively detect both large and small structural variations as well as TEs due to genetic drift, we generated de novo genome assemblies and annotations for each strain from our lab collection using both long- and short-read sequencing and compared our assemblies and annotations with that of other lab wild type strains. Within our lab assemblies, we annotate over 3.1Mb of sequence divergence between the Bristol and Hawaiian isolates: 337,584 SNPs, 94,503 small insertion-deletions (<50bp), and 4,334 structural variations (>50bp). Further, we define the location and movement of specific DNA TEs between N2 Bristol and CB4856 Hawaiian wild type isolates. Specifically, we find the N2 Bristol genome has 20.6% more TEs from the Tc1/mariner family than the CB4856 Hawaiian genome. Moreover, we identified Zator elements as the most abundant and mobile TE family in the genome. Using specific TE sequences with unique SNPs, we also identify 38 TEs that moved intrachromosomally and 9 TEs that moved interchromosomally between the N2 Bristol and CB4856 Hawaiian genomes. By comparing the de novo genome assembly of our lab collection Bristol isolate to the VC2010 Bristol assembly, we also reveal that lab lineages display over 2 Mb of total variation: 1,162 SNPs, 1,528 indels, and 897 SVs with 95% of the variation due to SVs. Overall, our work demonstrates the unique contribution of SVs and TEs to variation and genetic drift between wild type laboratory strains assumed to be isogenic despite growing evidence of genetic drift and phenotypic variation.

Identification of a spontaneously arising variant affecting thermotaxis behavior in a recombinant inbred Caenorhabditis elegans line

Analyses of the contributions of genetic variants in wild strains to phenotypic differences have led to a more complete description of the pathways underlying cellular functions. Causal loci are typically identified via interbreeding of strains with distinct phenotypes in order to establish recombinant inbred lines (RILs). Since the generation of RILs requires growth for multiple generations, their genomes may contain not only different combinations of parental alleles but also genetic changes that arose de novo during the establishment of these lines. Here, we report that in the course of generating RILs between Caenorhabditis elegans strains that exhibit distinct thermotaxis behavioral phenotypes, we identified spontaneously arising variants in the ttx-1 locus. ttx-1 encodes the terminal selector factor for the AFD thermosensory neurons, and loss-of-function mutations in ttx-1 abolish thermotaxis behaviors. The identified genetic changes in ttx-1 in the RIL are predicted to decrease ttx-1 function in part via specifically affecting a subset of AFD-expressed ttx-1 isoforms. Introduction of the relevant missense mutation in the laboratory C. elegans strain via gene editing recapitulates the thermotaxis behavioral defects of the RIL. Our results suggest that spontaneously occurring genomic changes in RILs may complicate identification of loci contributing to phenotypic variation, but that these mutations may nevertheless lead to the identification of important causal molecules and mechanisms.

Effect of recombination on genetic diversity of Caenorhabditis elegans

Greater molecular divergence and genetic diversity are present in regions of high recombination in many species. Studies describing the correlation between variant abundance and recombination rate have long focused on recombination in the context of linked selection models, whereby interference between linked sites under positive or negative selection reduces genetic diversity in regions of low recombination. Here, we show that indels, especially those of intermediate sizes, are enriched relative to single nucleotide polymorphisms in regions of high recombination in C. elegans. To explain this phenomenon, we reintroduce an alternative model that emphasizes the mutagenic effect of recombination. To extend the analysis, we examine the variants with a phylogenetic context and discuss how different models could be examined together. The number of variants generated by recombination in natural populations could be substantial including possibly the majority of some indel subtypes. Our work highlights the potential importance of a mutagenic effect of recombination, which could have a significant role in the shaping of natural genetic diversity.

Telomerase-independent survival leads to a mosaic of complex subtelomere rearrangements in Chlamydomonas reinhardtii

Telomeres and subtelomeres, the genomic regions located at chromosome extremities, are essential for genome stability in eukaryotes. In the absence of the canonical maintenance mechanism provided by telomerase, telomere shortening induces genome instability. The landscape of the ensuing genome rearrangements is not accessible by short-read sequencing. Here, we leverage Oxford Nanopore Technologies long-read sequencing to survey the extensive repertoire of genome rearrangements in telomerase mutants of the model green microalga Chlamydomonas reinhardtii In telomerase-mutant strains grown for hundreds of generations, most chromosome extremities were capped by short telomere sequences that were either recruited de novo from other loci or maintained in a telomerase-independent manner. Other extremities did not end with telomeres but only with repeated subtelomeric sequences. The subtelomeric elements, including rDNA, were massively rearranged and involved in breakage-fusion-bridge cycles, translocations, recombinations, and chromosome circularization. These events were established progressively over time and displayed heterogeneity at the subpopulation level. New telomere-capped extremities composed of sequences originating from more internal genomic regions were associated with high DNA methylation, suggesting that de novo heterochromatin formation contributes to the restoration of chromosome end stability in C. reinhardtii The diversity of alternative strategies present in the same organism to maintain chromosome integrity and the variety of rearrangements found in telomerase mutants are remarkable, and illustrate genome plasticity at short timescales.

Beyond the reference: gene expression variation and transcriptional response to RNA interference in Caenorhabditis elegans

Though natural systems harbor genetic and phenotypic variation, research in model organisms is often restricted to a reference strain. Focusing on a reference strain yields a great depth of knowledge but potentially at the cost of breadth of understanding. Furthermore, tools developed in the reference context may introduce bias when applied to other strains, posing challenges to defining the scope of variation within model systems. Here, we evaluate how genetic differences among 5 wild Caenorhabditis elegans strains affect gene expression and its quantification, in general and after induction of the RNA interference (RNAi) response. Across strains, 34% of genes were differentially expressed in the control condition, including 411 genes that were not expressed at all in at least 1 strain; 49 of these were unexpressed in reference strain N2. Reference genome mapping bias caused limited concern: despite hyperdiverse hotspots throughout the genome, 92% of variably expressed genes were robust to mapping issues. The transcriptional response to RNAi was highly strain- and target-gene-specific and did not correlate with RNAi efficiency, as the 2 RNAi-insensitive strains showed more differentially expressed genes following RNAi treatment than the RNAi-sensitive reference strain. We conclude that gene expression, generally and in response to RNAi, differs across C. elegans strains such that the choice of strain may meaningfully influence scientific inferences. Finally, we introduce a resource for querying gene expression variation in this dataset at https://wildworm.biosci.gatech.edu/rnai/.

Benchmarking datasets for assembly-based variant calling using high-fidelity long reads

BACKGROUND

Recent advances in long-read sequencing technologies have enabled accurate identification of all genetic variants in individuals or cells; this procedure is known as variant calling. However, benchmarking studies on variant calling using different long-read sequencing technologies are still lacking.

RESULTS

We used two Caenorhabditis elegans strains to measure several variant calling metrics. These two strains shared true-positive genetic variants that were introduced during strain generation. In addition, both strains contained common and distinguishable variants induced by DNA damage, possibly leading to false-positive estimation. We obtained accurate and noisy long reads from both strains using high-fidelity (HiFi) and continuous long-read (CLR) sequencing platforms, and compared the variant calling performance of the two platforms. HiFi identified a 1.65-fold higher number of true-positive variants on average, with 60% fewer false-positive variants, than CLR did. We also compared read-based and assembly-based variant calling methods in combination with subsampling of various sequencing depths and demonstrated that variant calling after genome assembly was particularly effective for detection of large insertions, even with 10 × sequencing depth of accurate long-read sequencing data.

CONCLUSIONS

By directly comparing the two long-read sequencing technologies, we demonstrated that variant calling after genome assembly with 10 × or more depth of accurate long-read sequencing data allowed reliable detection of true-positive variants. Considering the high cost of HiFi sequencing, we herein propose appropriate methodologies for performing cost-effective and high-quality variant calling: 10 × assembly-based variant calling. The results of the present study may facilitate the development of methods for identifying all genetic variants at the population level.

Caenorhabditis elegans as a model system to evaluate neuroprotective potential of nano formulations

The impact of neurodegenerative illnesses on society is significant, but the mechanisms leading to neuronal malfunction and death in these conditions remain largely unknown despite identifying essential disease genes. To pinpoint the mechanisms behind the pathophysiology of neurodegenerative diseases, several researchers have turned to nematode C. elegans instead of using mammals. Since C. elegans is transparent, free-living, and amenable to culture, it has several benefits. As a result, all the neurons in C. elegans can be easily identified, and their connections are understood. Human proteins linked to Neurodegeneration can be made to express in them. It is also possible to analyze how C. elegans orthologs of the genes responsible for human neurodegenerative diseases function. In this article, we focused at some of the most important C. elegans neurodegeneration models that accurately represent many elements of human neurodegenerative illness. It has been observed that studies using the adaptable C. elegans have helped us in better understanding of human diseases. These studies have used it to replicate several aspects of human neurodegeneration. A nanotech approach involves engineering materials or equipments interacting with biological systems at the molecular level to trigger physiological responses by increasing stimulation, responding, and interacting with target sites while minimizing side effects, thus revolutionizing the treatment and diagnosis of neurodegenerative diseases. Nanotechnologies are being used to treat neurological disorders and deliver nanoscale drugs. This review explores the current and future uses of these nanotechnologies as innovative therapeutic modalities in treatment of neurodegenerative diseases using C elegans as an experimental model.

A beginner's guide to assembling a draft genome and analyzing structural variants with long-read sequencing technologies

Advances in long-read DNA sequencing technologies have enabled researchers to obtain high-quality genomes and finely resolve structural variants (SVs) in many species, even from small laboratories. The hands-on protocol presented here will guide you through the process of analyzing three different types of publicly available Drosophila melanogaster datasets obtained using current long-read sequencing technologies. We hope that this protocol will help in guiding researchers who are new to the process of long-read sequencing analysis.

Multiple Pristionchus pacificus genomes reveal distinct evolutionary dynamics between de novo candidates and duplicated genes

The birth of new genes is a major molecular innovation driving phenotypic diversity across all domains of life. Although repurposing of existing protein-coding material by duplication is considered the main process of new gene formation, recent studies have discovered thousands of transcriptionally active sequences as a rich source of new genes. However, differential loss rates have to be assumed to reconcile the high birth rates of these incipient de novo genes with the dominance of ancient gene families in individual genomes. Here, we test this rapid turnover hypothesis in the context of the nematode model organism Pristionchus pacificus We extended the existing species-level phylogenomic framework by sequencing the genomes of six divergent P. pacificus strains. We used these data to study the evolutionary dynamics of different age classes and categories of origin at a population level. Contrasting de novo candidates with new families that arose by duplication and divergence from known genes, we find that de novo candidates are typically shorter, show less expression, and are overrepresented on the sex chromosome. Although the contribution of de novo candidates increases toward young age classes, multiple comparisons within the same age class showed significantly higher attrition in de novo candidates than in known genes. Similarly, young genes remain under weak evolutionary constraints with de novo candidates representing the fastest evolving subcategory. Altogether, this study provides empirical evidence for the rapid turnover hypothesis and highlights the importance of the evolutionary timescale when quantifying the contribution of different mechanisms toward new gene formation.

Intraspecific de novo gene birth revealed by presence-absence variant genes in Caenorhabditis elegans

Genes embed their evolutionary history in the form of various alleles. Presence-absence variants (PAVs) are extreme cases of such alleles, where a gene present in one haplotype does not exist in another. Because PAVs may result from either birth or death of a gene, PAV genes and their alternative alleles, if available, can represent a basis for rapid intraspecific gene evolution. Using long-read sequencing technologies, this study traced the possible evolution of PAV genes in the PD1074 and CB4856 C. elegans strains as well as their alternative alleles in 14 other wild strains. We updated the CB4856 genome by filling 18 gaps and identified 46 genes and 7,460 isoforms from both strains not annotated previously. We verified 328 PAV genes, out of which 46 were C. elegans-specific. Among these possible newly born genes, 12 had alternative alleles in other wild strains; in particular, the alternative alleles of three genes showed signatures of active transposons. Alternative alleles of three other genes showed another type of signature reflected in accumulation of small insertions or deletions. Research on gene evolution using both species-specific PAV genes and their alternative alleles may provide new insights into the process of gene evolution.

Long-read sequencing infers a mechanism for copy number variation of template for alternative lengthening of telomeres in a wild C. elegans strain

Template for alternative lengthening of telomeres 1 (TALT1) is a specific sequence used to protect chromosomal ends from telomere damage first identified in the CB4856 strain of Caenorhabditis elegans . Here we assembled the genome of DL226, a wild strain with one more copy of TALT1-like sequences in its genome compared to those of CB4856, using long-read DNA sequencing. We found that a five-copy array of short telomeric repeats and TALT1s present in CB4856 were changed to a six-copy array due to the duplication of the third copy; there was an additional damage-repair trace in the new short telomeric repeat near the newly replicated TALT1.

Chromosome-level reference genomes for two strains of Caenorhabditis briggsae: an improved platform for comparative genomics

The publication of the Caenorhabditis briggsae reference genome in 2003 enabled the first comparative genomics studies between C. elegans and C. briggsae, shedding light on the evolution of genome content and structure in the Caenorhabditis genus. However, despite being widely used, the currently available C. briggsae reference genome is substantially less complete and structurally accurate than the C. elegans reference genome. Here, we used high-coverage Oxford Nanopore long-read and chromosome conformation capture data to generate chromosome-level reference genomes for two C. briggsae strains: QX1410, a new reference strain closely related to the laboratory AF16 strain, and VX34, a highly divergent strain isolated in China. We also sequenced 99 recombinant inbred lines (RILs) generated from reciprocal crosses between QX1410 and VX34 to create a recombination map and identify chromosomal domains. Additionally, we used both short- and long-read RNA sequencing (RNA-seq) data to generate high-quality gene annotations. By comparing these new reference genomes to the current reference, we reveal that hyper-divergent haplotypes cover large portions of the C. briggsae genome, similar to recent reports in C. elegans and C. tropicalis. We also show that the genomes of selfing Caenorhabditis species have undergone more rearrangement than their outcrossing relatives, which has biased previous estimates of rearrangement rate in Caenorhabditis. These new genomes provide a substantially improved platform for comparative genomics in Caenorhabditis and narrow the gap between the quality of genomic resources available for C. elegans and C. briggsae.

Genomic architecture of 5S rDNA cluster and its variations within and between species

BACKGROUND

Ribosomal DNAs (rDNAs) are arranged in purely tandem repeats, preventing them from being reliably assembled onto chromosomes during generation of genome assembly. The uncertainty of rDNA genomic structure presents a significant barrier for studying their function and evolution.

RESULTS

Here we generate ultra-long Oxford Nanopore Technologies (ONT) and short NGS reads to delineate the architecture and variation of the 5S rDNA cluster in the different strains of C. elegans and C. briggsae. We classify the individual rDNA's repeating units into 25 types based on the unique sequence variations in each unit of C. elegans (N2). We next perform assembly of the cluster by taking advantage of the long reads that carry these units, which led to an assembly of 5S rDNA cluster consisting of up to 167 consecutive 5S rDNA units in the N2 strain. The ordering and copy number of various rDNA units are consistent with the separation time between strains. Surprisingly, we observed a drastically reduced level of variation in the unit composition in the 5S rDNA cluster in the C. elegans CB4856 and C. briggsae AF16 strains than in the C. elegans N2 strain, suggesting that N2, a widely used reference strain, is likely to be defective in maintaining the 5S rDNA cluster stability compared with other wild isolates of C. elegans or C. briggsae.

CONCLUSIONS

The results demonstrate that Nanopore DNA sequencing reads are capable of generating assembly of highly repetitive sequences, and rDNA units are highly dynamic both within and between population(s) of the same species in terms of sequence and copy number. The detailed structure and variation of the 5S rDNA units within the rDNA cluster pave the way for functional and evolutionary studies.

Nematode chromosomes

The nematode Caenorhabditis elegans has shed light on many aspects of eukaryotic biology, including genetics, development, cell biology, and genomics. A major factor in the success of C. elegans as a model organism has been the availability, since the late 1990s, of an essentially gap-free and well-annotated nuclear genome sequence, divided among 6 chromosomes. In this review, we discuss the structure, function, and biology of C. elegans chromosomes and then provide a general perspective on chromosome biology in other diverse nematode species. We highlight malleable chromosome features including centromeres, telomeres, and repetitive elements, as well as the remarkable process of programmed DNA elimination (historically described as chromatin diminution) that induces loss of portions of the genome in somatic cells of a handful of nematode species. An exciting future prospect is that nematode species may enable experimental approaches to study chromosome features and to test models of chromosome evolution. In the long term, fundamental insights regarding how speciation is integrated with chromosome biology may be revealed.

Comparison of structural variants in the whole genome sequences of two Medicago truncatula ecotypes: Jemalong A17 and R108

BACKGROUND

Structural variants (SVs) constitute a large proportion of the genomic variation that results in phenotypic variation in plants. However, they are still a largely unexplored feature in most plant genomes. Here, we present the whole-genome landscape of SVs between two model legume Medicago truncatula ecotypes-Jemalong A17 and R108- that have been extensively used in various legume biology studies.

RESULTS

To catalogue SVs, we first resolved the previously published R108 genome assembly (R108 v1.0) to chromosome-scale using 124 × Hi-C data, resulting in a high-quality genome assembly. The inter-chromosomal reciprocal translocations between chromosomes 4 and 8 were confirmed by performing syntenic analysis between the two genomes. Combined with the Hi-C data, it appears that these translocation events had a significant effect on chromatin organization. Using both whole-genome and short-read alignments, we identified the genomic landscape of SVs between the two genomes, some of which may account for several phenotypic differences, including their differential responses to aluminum toxicity and iron deficiency, and the development of different anthocyanin leaf markings. We also found extensive SVs within the nodule-specific cysteine-rich gene family which encodes antimicrobial peptides essential for terminal bacteroid differentiation during nitrogen-fixing symbiosis.

CONCLUSIONS

Our results provide a near-complete R108 genome assembly and the first genomic landscape of SVs obtained by comparing two M. truncatula ecotypes. This may provide valuable genomic resources for the functional and molecular research of legume biology in the future.

Multi-omics analysis identifies essential regulators of mitochondrial stress response in two wild-type C. elegans strains

The mitochondrial unfolded protein response (UPRmt) is a promising pharmacological target for aging and age-related diseases. However, the integrative analysis of the impact of UPRmt activation on different signaling layers in animals with different genetic backgrounds is lacking. Here, we applied systems approaches to investigate the effect of UPRmt induced by doxycycline (Dox) on transcriptome, proteome, and lipidome in two genetically divergent worm strains, named N2 and CB4856. From the integrated omics datasets, we found that Dox prolongs lifespan of both worm strains through shared and strain-specific mechanisms. Specifically, Dox strongly impacts mitochondria, upregulates defense response, and lipid metabolism, while decreasing triglycerides. We further validated that lipid genes acs-2/20 and fat-7/6 were required for Dox-induced UPRmt and longevity in N2 and CB4856 worms, respectively. Our data have translational value as they indicate that the beneficial effects of Dox-induced UPRmt on lifespan are consistent across different genetic backgrounds through different regulators.

•

Dox extends lifespan of N2 and CB4856 via shared and strain-specific mechanisms

•

Dox controls mitochondria, defense responses, and lipid metabolism in both strains

•

Dox-mediated longevity requires acs-2/20 in N2 and fat-7/6 in CB4856 worms

Interface of Human/Wildlife Interactions: An Example of a Bold Coyote (Canis latrans) in Atlanta, GA, USA

There is arguably no other North American species that better illustrates the complexities of the human-wildlife interface than the coyote. In this study, a melanistic coyote in metropolitan Atlanta, Georgia was exhibiting unusually bold behaviors that included encounters with humans, domestic dogs, and attempts to enter homes. After tracking this coyote (nicknamed Carmine) across a highly urbanized landscape with participatory science, including at least 80 publicly reported sightings, he was captured and relocated to a wildlife sanctuary. Genome-wide analyses revealed 92.8% coyote ancestry, 1.7% gray wolf ancestry, and 5.5% domestic dog ancestry. The dog alleles in Carmine’s genome were estimated to have been acquired by his ancestors 14–29 years ago. Despite his bold behavior, Carmine did not carry any mutations known to shape hypersociability in canines. He did, however, carry a single copy of the dominant mutation responsible for his melanistic coat color. This detailed study of Carmine dispels common assumptions about the reticent coyote personality and the origins of behavior. His unusual bold behavior created a higher level of human-coyote interaction. He now serves as a public ambassador for human-wildlife coexistence, urging the global community to reconsider mythologies about wildlife and promote coexistence with them in landscapes significantly altered by human activity in our rapidly changing world.

Architecture and evolution of subtelomeres in the unicellular green alga Chlamydomonas reinhardtii

In most eukaryotes, subtelomeres are dynamic genomic regions populated by multi-copy sequences of different origins, which can promote segmental duplications and chromosomal rearrangements. However, their repetitive nature has complicated the efforts to sequence them, analyse their structure and infer how they evolved. Here, we use recent genome assemblies of Chlamydomonas reinhardtii based on long-read sequencing to comprehensively describe the subtelomere architecture of the 17 chromosomes of this model unicellular green alga. We identify three main repeated elements present at subtelomeres, which we call Sultan, Subtile and Suber, alongside three chromosome extremities with ribosomal DNA as the only identified component of their subtelomeres. The most common architecture, present in 27 out of 34 subtelomeres, is a heterochromatic array of Sultan elements adjacent to the telomere, followed by a transcribed Spacer sequence, a G-rich microsatellite and transposable elements. Sequence similarity analyses suggest that Sultan elements underwent segmental duplications within each subtelomere and rearranged between subtelomeres at a much lower frequency. Analysis of other green algae reveals species-specific repeated elements that are shared across subtelomeres, with an overall organization similar to C. reinhardtii. This work uncovers the complexity and evolution of subtelomere architecture in green algae.

Architecture and evolution of subtelomeres in the unicellular green alga Chlamydomonas reinhardtii

In most eukaryotes, subtelomeres are dynamic genomic regions populated by multi-copy sequences of different origins, which can promote segmental duplications and chromosomal rearrangements. However, their repetitive nature has complicated the efforts to sequence them, analyse their structure and infer how they evolved. Here, we use recent genome assemblies of Chlamydomonas reinhardtii based on long-read sequencing to comprehensively describe the subtelomere architecture of the 17 chromosomes of this model unicellular green alga. We identify three main repeated elements present at subtelomeres, which we call Sultan, Subtile and Suber, alongside three chromosome extremities with ribosomal DNA as the only identified component of their subtelomeres. The most common architecture, present in 27 out of 34 subtelomeres, is a heterochromatic array of Sultan elements adjacent to the telomere, followed by a transcribed Spacer sequence, a G-rich microsatellite and transposable elements. Sequence similarity analyses suggest that Sultan elements underwent segmental duplications within each subtelomere and rearranged between subtelomeres at a much lower frequency. Analysis of other green algae reveals species-specific repeated elements that are shared across subtelomeres, with an overall organization similar to C. reinhardtii. This work uncovers the complexity and evolution of subtelomere architecture in green algae.

Balancing selection maintains hyper-divergent haplotypes in Caenorhabditis elegans

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.

Punctuated Loci on Chromosome IV Determine Natural Variation in Orsay Virus Susceptibility of Caenorhabditis elegans Strains Bristol N2 and Hawaiian CB4856

Host-pathogen interactions play a major role in evolutionary selection and shape natural genetic variation. The genetically distinct Caenorhabditis elegans strains, Bristol N2 and Hawaiian CB4856, are differentially susceptible to the Orsay virus (OrV). Here, we report the dissection of the genetic architecture of susceptibility to OrV infection. We compare OrV infection in the relatively resistant wild-type CB4856 strain to the more susceptible canonical N2 strain. To gain insight into the genetic architecture of viral susceptibility, 52 fully sequenced recombinant inbred lines (CB4856 × N2 RILs) were exposed to OrV. This led to the identification of two loci on chromosome IV associated with OrV resistance. To verify the two loci and gain additional insight into the genetic architecture controlling virus infection, introgression lines (ILs) that together cover chromosome IV, were exposed to OrV. Of the 27 ILs used, 17 had an CB4856 introgression in an N2 background, and 10 had an N2 introgression in a CB4856 background. Infection of the ILs confirmed and fine-mapped the locus underlying variation in OrV susceptibility, and we found that a single nucleotide polymorphism in cul-6 may contribute to the difference in OrV susceptibility between N2 and CB4856. An allele swap experiment showed the strain CB4856 became as susceptible as the N2 strain by having an N2 cul-6 allele, although having the CB4856 cul-6 allele did not increase resistance in N2. In addition, we found that multiple strains with nonoverlapping introgressions showed a distinct infection phenotype from the parental strain, indicating that there are punctuated locations on chromosome IV determining OrV susceptibility. Thus, our findings reveal the genetic complexity of OrV susceptibility in C. elegans and suggest that viral susceptibility is governed by multiple genes.IMPORTANCE Genetic variation determines the viral susceptibility of hosts. Yet, pinpointing which genetic variants determine viral susceptibility remains challenging. Here, we have exploited the genetic tractability of the model organism Caenorhabditis elegans to dissect the genetic architecture of Orsay virus infection. Our results provide novel insight into natural determinants of Orsay virus infection.

Long-read sequencing and de novo genome assemblies reveal complex chromosome end structures caused by telomere dysfunction at the single nucleotide level

Karyotype change and subsequent evolution is triggered by chromosome fusion and rearrangement events, which often occur when telomeres become dysfunctional. Telomeres protect linear chromosome ends from DNA damage responses (DDRs), and telomere dysfunction may result in genome instability. However, the complex chromosome end structures and the other possible consequences of telomere dysfunction have rarely been resolved at the nucleotide level due to the lack of the high-throughput methods needed to analyse these highly repetitive regions. Here we applied long-read sequencing technology to Caenorhabditis elegans survivor lines that emerged after telomere dysfunction. The survivors have preserved traces of DDRs in their genomes and our data revealed that variants generated by telomere dysfunction are accumulated along all chromosomes. The reconstruction of the chromosome end structures through de novo genome assemblies revealed diverse types of telomere damage processing at the nucleotide level. When telomeric repeats were totally eroded by telomere dysfunction, DDRs were mostly terminated by chromosome fusion events. We also partially reconstructed the most complex end structure and its DDR signatures, which would have been accumulated via multiple cell divisions. These finely resolved chromosome end structures suggest possible mechanisms regarding the repair processes after telomere dysfunction, providing insights into chromosome evolution in nature.

Evolution: Drivers in a Traffic Jam

Genetic loci coding for a toxin and its antidote behave like selfish elements. Two new studies find an accumulation of such elements in one species.

DNA repair, recombination, and damage signaling

DNA must be accurately copied and propagated from one cell division to the next, and from one generation to the next. To ensure the faithful transmission of the genome, a plethora of distinct as well as overlapping DNA repair and recombination pathways have evolved. These pathways repair a large variety of lesions, including alterations to single nucleotides and DNA single and double-strand breaks, that are generated as a consequence of normal cellular function or by external DNA damaging agents. In addition to the proteins that mediate DNA repair, checkpoint pathways have also evolved to monitor the genome and coordinate the action of various repair pathways. Checkpoints facilitate repair by mediating a transient cell cycle arrest, or through initiation of cell suicide if DNA damage has overwhelmed repair capacity. In this chapter, we describe the attributes of Caenorhabditis elegans that facilitate analyses of DNA repair, recombination, and checkpoint signaling in the context of a whole animal. We review the current knowledge of C. elegans DNA repair, recombination, and DNA damage response pathways, and their role during development, growth, and in the germ line. We also discuss how the analysis of mutational signatures in C. elegans is helping to inform cancer mutational signatures in humans.

Natural variation in the sequestosome-related gene, sqst-5, underlies zinc homeostasis in Caenorhabditis elegans

Zinc is an essential trace element that acts as a co-factor for many enzymes and transcription factors required for cellular growth and development. Altering intracellular zinc levels can produce dramatic effects ranging from cell proliferation to cell death. To avoid such fates, cells have evolved mechanisms to handle both an excess and a deficiency of zinc. Zinc homeostasis is largely maintained via zinc transporters, permeable channels, and other zinc-binding proteins. Variation in these proteins might affect their ability to interact with zinc, leading to either increased sensitivity or resistance to natural zinc fluctuations in the environment. We can leverage the power of the roundworm nematode Caenorhabditis elegans as a tractable metazoan model for quantitative genetics to identify genes that could underlie variation in responses to zinc. We found that the laboratory-adapted strain (N2) is resistant and a natural isolate from Hawaii (CB4856) is sensitive to micromolar amounts of exogenous zinc supplementation. Using a panel of recombinant inbred lines, we identified two large-effect quantitative trait loci (QTL) on the left arm of chromosome III and the center of chromosome V that are associated with zinc responses. We validated and refined both QTL using near-isogenic lines (NILs) and identified a naturally occurring deletion in sqst-5, a sequestosome-related gene, that is associated with resistance to high exogenous zinc. We found that this deletion is relatively common across strains within the species and that variation in sqst-5 is associated with zinc resistance. Our results offer a possible mechanism for how organisms can respond to naturally high levels of zinc in the environment and how zinc homeostasis varies among individuals.

Draft Genomes of Amaranthus tuberculatus, Amaranthus hybridus, and Amaranthus palmeri

Amaranthus tuberculatus, Amaranthus hybridus, and Amaranthus palmeri are agronomically important weed species. Here, we present the most contiguous draft assemblies of these three species to date. We utilized a combination of Pacific Biosciences long-read sequencing and chromatin contact mapping information to assemble and order sequences of A. palmeri to near-chromosome-level resolution, with scaffold N50 of 20.1 Mb. To resolve the issues of heterozygosity and coassembly of alleles in diploid species, we adapted the trio binning approach to produce haplotype assemblies of A. tuberculatus and A. hybridus. This approach resulted in an improved assembly of A. tuberculatus, and the first genome assembly for A. hybridus, with contig N50s of 2.58 and 2.26 Mb, respectively. Species-specific transcriptomes and information from related species were used to predict transcripts within each assembly. Syntenic comparisons of these species and Amaranthus hypochondriacus identified sites of genomic rearrangement, including duplication and translocation, whereas genetic map construction within A. tuberculatus highlighted the need for further ordering of the A. hybridus and A. tuberculatus contigs. These multiple reference genomes will accelerate genomic studies in these species to further our understanding of weedy evolution within Amaranthus.

Stoichiometric interactions explain spindle dynamics and scaling across 100 million years of nematode evolution

The spindle shows remarkable diversity, and changes in an integrated fashion, as cells vary over evolution. Here, we provide a mechanistic explanation for variations in the first mitotic spindle in nematodes. We used a combination of quantitative genetics and biophysics to rule out broad classes of models of the regulation of spindle length and dynamics, and to establish the importance of a balance of cortical pulling forces acting in different directions. These experiments led us to construct a model of cortical pulling forces in which the stoichiometric interactions of microtubules and force generators (each force generator can bind only one microtubule), is key to explaining the dynamics of spindle positioning and elongation, and spindle final length and scaling with cell size. This model accounts for variations in all the spindle traits we studied here, both within species and across nematode species spanning over 100 million years of evolution.

Caenorhabditis elegans dauers vary recovery in response to bacteria from natural habitat

Many species use dormant stages for habitat selection by tying recovery to informative external cues. Other species have an undiscerning strategy in which they recover randomly despite having advanced sensory systems. We investigated whether elements of a species' habitat structure and life history can bar it from developing a discerning recovery strategy. The nematode Caenorhabditis elegans has a dormant stage called the dauer larva that disperses between habitat patches. On one hand, C. elegans colonization success is profoundly influenced by the bacteria found in its habitat patches, so we might expect this to select for a discerning strategy. On the other hand, C. elegans' habitat structure and life history suggest that there is no fitness benefit to varying recovery, which might select for an undiscerning strategy. We exposed dauers of three genotypes to a range of bacteria acquired from the worms' natural habitat. We found that C. elegans dauers recover in all conditions but increase recovery on certain bacteria depending on the worm's genotype, suggesting a combination of undiscerning and discerning strategies. Additionally, the worms' responses did not match the bacteria's objective quality, suggesting that their decision is based on other characteristics.

Comparative Analysis of Clinical and Environmental Strains of Exophiala spinifera by Long-Reads Sequencing and RNAseq Reveal Adaptive Strategies

Exophiala spinifera, a capsule-producing black yeast, is overrepresented as agent of disseminated infection in humans with inherited dysfunction of the CARD9 gene. In a review of published caspase recruitment domain-containing protein 9 (CARD9) deficiency cases, black fungi were linked to mutations other than those prevalent in yeast and dermatophyte cases, and were found to respond to a larger panel of cytokines. Here, we sequenced and annotated the genomes of BMU 08022 from a patient with CARD9 deficiency and two environmental strains, BMU 00051 and BMU 00047. We performed genomic and transcriptomic analysis for these isolates including published black yeasts genomes, using a combination of long-read (PACBIO) and short-read (Illumina) sequencing technologies with a hybrid assembly strategy. We identified the virulence factors, fitness, and the major genetic and gene expression differences between the strains with RNAseq technology. Genome assembly reached sub-chromosome level with between 12,043 and 12,130 predicted genes. The number of indels identified in the clinical strain was higher than observed in environmental strains. We identify a relatively large core genome of 9,887 genes. Moreover, substantial syntenic rearrangements of scaffolds I and III in the CARD9-related isolate were detected. Seventeen gene clusters were involved in the production of secondary metabolites. PKS-cluster 17 was consistently found to be absent in the clinical strain. Comparative transcriptome analysis demonstrated that 16 single-copy genes were significantly differentially expressed upon incubation in brain-heart infusion broth vs. Sabouraud glucose broth. Most of the single-copy genes upregulated with Brain Heart Infusion (BHI) were transporters. There were 48 unique genes differentially expressed exclusively to the clinical strain in two different media, including genes from various metabolic processes and transcriptional regulation. Up-regulated genes in the clinical strain with Gene Ontology (GO) enrichment are mainly involved in transmembrane transport, biosynthetic process and metabolic process. This study has provided novel insights into understanding of strain-differences in intrinsic virulence of the species and indicated that intraspecific variability may be related to habitat choice. This indicates that strains of E. spinifera are differentially prone to cause infection in susceptible patient populations, and provides clues for future studies exploring the mechanisms of pathogenic and adaptive strategies of black yeasts in immunodeficient patients.

Repair and Reconstruction of Telomeric and Subtelomeric Regions and Genesis of New Telomeres: Implications for Chromosome Evolution

DNA damage repair within telomeres are suppressed to maintain the integrity of linear chromosomes, but the accidental activation of repairs can lead to genome instability. This review develops the concept that mechanisms to repair DNA damage in telomeres contribute to genetic variability and karyotype evolution, rather than catastrophe. Spontaneous breaks in telomeres can be repaired by telomerase, but in some cases DNA repair pathways are activated, and can cause chromosomal rearrangements or fusions. The resultant changes can also affect subtelomeric regions that are adjacent to telomeres. Subtelomeres are actively involved in such chromosomal changes, and are therefore the most variable regions in the genome. The case of Caenorhabditis elegans in the context of changes of subtelomeric structures revealed by long-read sequencing is also discussed. Theoretical and methodological issues covered in this review will help to explore the mechanism of chromosome evolution by reconstruction of chromosomal ends in nature.

A spontaneous complex structural variant in rcan-1 increases exploratory behavior and laboratory fitness of Caenorhabditis elegans

Over long evolutionary timescales, major changes to the copy number, function, and genomic organization of genes occur, however, our understanding of the individual mutational events responsible for these changes is lacking. In this report, we study the genetic basis of adaptation of two strains of C. elegans to laboratory food sources using competition experiments on a panel of 89 recombinant inbred lines (RIL). Unexpectedly, we identified a single RIL with higher relative fitness than either of the parental strains. This strain also displayed a novel behavioral phenotype, resulting in higher propensity to explore bacterial lawns. Using bulk-segregant analysis and short-read resequencing of this RIL, we mapped the change in exploration behavior to a spontaneous, complex rearrangement of the rcan-1 gene that occurred during construction of the RIL panel. We resolved this rearrangement into five unique tandem inversion/duplications using Oxford Nanopore long-read sequencing. rcan-1 encodes an ortholog to human RCAN1/DSCR1 calcipressin gene, which has been implicated as a causal gene for Down syndrome. The genomic rearrangement in rcan-1 creates two complete and two truncated versions of the rcan-1 coding region, with a variety of modified 5’ and 3’ non-coding regions. While most copy-number variations (CNVs) are thought to act by increasing expression of duplicated genes, these changes to rcan-1 ultimately result in the reduction of its whole-body expression due to changes in the upstream regions. By backcrossing this rearrangement into a common genetic background to create a near isogenic line (NIL), we demonstrate that both the competitive advantage and exploration behavioral changes are linked to this complex genetic variant. This NIL strain does not phenocopy a strain containing an rcan-1 loss-of-function allele, which suggests that the residual expression of rcan-1 is necessary for its fitness effects. Our results demonstrate how colonization of new environments, such as those encountered in the laboratory, can create evolutionary pressure to modify gene function. This evolutionary mismatch can be resolved by an unexpectedly complex genetic change that simultaneously duplicates and diversifies a gene into two uniquely regulated genes. Our work shows how complex rearrangements can act to modify gene expression in ways besides increased gene dosage.

Evolution acts on genetic variants that modify phenotypes that increase the likelihood of staying alive and passing on these genetic changes to subsequent generations (i.e. fitness). There is general interest in understanding the types of genetic variants that can increase fitness in specific environments. One route that fitness can be increased is through changes in behavior, such as finding new food sources. Here, we identify a spontaneous genetic change that increases exploration behavior and fitness of animals in laboratory environments. Interestingly, this genetic change is not a simple genetic change that deletes or changes the sequence of a protein product, but rather a complex structural variant that simultaneously duplicates the rcan-1 gene and also modifies its expression in a number of tissues. Our work demonstrates how a complex structural change can duplicate a gene, modify the DNA control regions that determine its cellular sites of action, and confer a fitness advantage that could lead to its spread in a population.

The Telomere Paradox: Stable Genome Preservation with Rapidly Evolving Proteins

Telomeres ensure chromosome length homeostasis and protection from catastrophic end-to-end chromosome fusions. All eukaryotes require this essential, strictly conserved telomere-dependent genome preservation. However, recent evolutionary analyses of mammals, plants, and flies report pervasive rapid evolution of telomere proteins. The causes of this paradoxical observation - that unconserved machinery underlies an essential, conserved function - remain enigmatic. Indeed, these fast-evolving telomere proteins bind, extend, and protect telomeric DNA, which itself evolves slowly in most systems. We hypothesize that the universally fast-evolving subtelomere - the telomere-adjacent, repetitive sequence - is a primary driver of the 'telomere paradox'. Under this model, radical sequence changes in the subtelomere perturb subtelomere-dependent, telomere functions. Compromised telomere function then spurs adaptation of telomere proteins to maintain telomere length homeostasis and protection. We propose an experimental framework that leverages both protein divergence and subtelomeric sequence divergence to test the hypothesis that subtelomere sequence evolution shapes recurrent innovation of telomere machinery.

Telomere-Specialized Retroelements in Drosophila: Adaptive Symbionts of the Genome, Neutral, or in Conflict?

Linear chromosomes shorten in every round of replication. In Drosophila, telomere-specialized long interspersed retrotransposable elements (LINEs) belonging to the jockey clade offset this shortening by forming head-to-tail arrays at Drosophila telomere ends. As such, these telomeric LINEs have been considered adaptive symbionts of the genome, protecting it from premature decay, particularly as Drosophila lacks a conventional telomerase holoenzyme. However, as reviewed here, recent work reveals a high degree of variation and turnover in the telomere-specialized LINE lineages across Drosophila. There appears to be no absolute requirement for LINE activity to maintain telomeres in flies, hence the suggestion that the telomere-specialized LINEs may instead be neutral or in conflict with the host, rather than adaptive.

Tightly linked antagonistic-effect loci underlie polygenic phenotypic variation in C. elegans

Recent work has provided strong empirical support for the classic polygenic model for trait variation. Population-based findings suggest that most regions of genome harbor variation affecting most traits. Here, we use the approach of experimental genetics to show that, indeed, most genomic regions carry variants with detectable effects on growth and reproduction in Caenorhabditis elegans populations sensitized by nickel stress. Nine of 15 adjacent intervals on the X chromosome, each encompassing ∼0.001 of the genome, have significant effects when tested individually in near-isogenic lines (NILs). These intervals have effects that are similar in magnitude to those of genome-wide significant loci that we mapped in a panel of recombinant inbred advanced intercross lines (RIAILs). If NIL-like effects were randomly distributed across the genome, the RIAILs would exhibit phenotypic variance that far exceeds the observed variance. However, the NIL intervals are arranged in a pattern that significantly reduces phenotypic variance relative to a random arrangement; adjacent intervals antagonize one another, cancelling each other's effects. Contrary to the expectation of small additive effects, our findings point to large-effect variants whose effects are masked by epistasis or linkage disequilibrium between alleles of opposing effect.