Erik Lavington

We used long-read sequencing to produce a telomere-to-telomere genome assembly for the heterotrophic stramenopile protist Aurantiochytrium limacinum MYA-1381. Its ∼62 Mbp nuclear genome comprises 26 linear chromosomes with a novel configuration: subtelomeric rDNAs are interspersed with long repeated sequence elements denoted as LOng REpeated - TElomere And Rdna Spacers (LORE-TEARS). These repeats may play a role in chromosome end maintenance. A ∼300 Kbp circular herpesvirus-like genomic element is present at a high copy number. A 269 Kbp related virus-like element was found to reside between two complete sets of rRNA and LORE-TEAR sequences on one end of chromosome 15, indicating recent recombination between the viral and nuclear genome. Our data reveal new types of giant endogenous viral elements originating from herpes-like viruses and existing as either ‘stand-alone’ or integrated elements.

Until recently, most viruses detected and characterized were of economic significance, associated with agricultural and medical diseases. This was certainly true for the eukaryote‐infecting circular Rep (replication‐associated protein)‐encoding single‐stranded DNA (CRESS DNA) viruses, which were thought to be a relatively small group of viruses. With the explosion of metagenomic sequencing over the past decade and increasing use of rolling‐circle replication for sequence amplification, scientists have identified and annotated copious numbers of novel CRESS DNA viruses – many without known hosts but which have been found in association with eukaryotes. Similar advances in cellular genomics have revealed that many eukaryotes have endogenous sequences homologous to viral Reps, which not only provide ‘fossil records’ to reconstruct the evolutionary history of CRESS DNA viruses but also reveal potential host species for viruses known by their sequences alone. The Rep protein is a conserved protein that all CRESS DNA viruses use to assist rolling‐circle replication that is known to be endogenized in a few eukaryotic species (notably tobacco and water yam). A systematic search for endogenous Rep‐like sequences in GenBank's non‐redundant eukaryotic database was performed using tBLASTn. We utilized relaxed search criteria for the capture of integrated Rep sequence within eukaryotic genomes, identifying 93 unique species with an endogenized fragment of Rep in their nuclear, plasmid (one species), mitochondrial (six species) or chloroplast (eight species) genomes. These species come from 19 different phyla, scattered across the eukaryotic tree of life. Exogenous and endogenous CRESS DNA viral Rep tree topology suggested potential hosts for one family of uncharacterized viruses and supports a primarily fungal host range for genomoviruses.

Cassava mosaic disease (CMD) represents a serious threat to cassava, a major root crop for more than 300 million Africans. CMD is caused by single-stranded DNA begomoviruses that evolve rapidly, making it challenging to develop durable disease resistance. In addition to the evolutionary forces of mutation, recombination and reassortment, factors such as climate, agriculture practices and the presence of DNA satellites may impact viral diversity. To gain insight into the factors that alter and shape viral diversity in planta, we used high-throughput sequencing to characterize the accumulation of nucleotide diversity after inoculation of infectious clones corresponding to African cassava mosaic virus (ACMV) and East African cassava mosaic Cameroon virus (EACMCV) in the susceptible cassava landrace Kibandameno. We found that vegetative propagation had a significant effect on viral nucleotide diversity, while temperature and a satellite DNA did not have measurable impacts in our study. EACMCV diversity increased linearly with the number of vegetative propagation passages, while ACMV diversity increased for a time and then decreased in later passages. We observed a substitution bias toward C→T and G→A for mutations in the viral genomes consistent with field isolates. Non-coding regions excluding the promoter regions of genes showed the highest levels of nucleotide diversity for each genome component. Changes in the 5′ intergenic region of DNA-A resembled the sequence of the cognate DNA-B sequence. The majority of nucleotide changes in coding regions were non-synonymous, most with predicted deleterious effects on protein structure, indicative of relaxed selection pressure over six vegetative passages. Overall, these results underscore the importance of knowing how cropping practices affect viral evolution and disease progression.

Cassava mosaic disease (CMD) represents a serious threat to cassava, a major root crop for more than 300 million Africans. CMD is caused by single-stranded DNA begomoviruses that evolve rapidly, making it challenging to develop durable disease resistance. In addition to the evolutionary forces of mutation, recombination, and reassortment, factors such as climate, agriculture practices, and the presence of DNA satellites may impact viral diversity. To gain insight into the factors that alter and shape viral diversity in planta , we used high-throughput sequencing to characterize the accumulation of nucleotide diversity after inoculation of infectious clones corresponding to African cassava mosaic virus (ACMV) and East African cassava mosaic Cameroon virus (EACMCV) in the susceptible cassava landrace Kibandameno. We found that vegetative propagation had a significant effect on viral nucleotide diversity, while temperature and a satellite DNA did not have measurable impacts in our study. EACMCV diversity increased linearly with the number of vegetative propagation passages, while ACMV diversity increased for a time and then decreased in later passages. We observed a substitution bias toward C→T and G→A for mutations in the viral genomes consistent with field isolates. Non-coding regions excluding the promoter regions of genes showed the highest levels of nucleotide diversity for each genome component. Changes in the 5’ intergenic region of DNA-A resembled the sequence of the cognate DNA-B sequence. The majority of nucleotide changes in coding regions were non-synonymous, most with predicted deleterious effects on protein structure, indicative of relaxed selection pressure over 6 vegetative passages. Overall, these results underscore the importance of knowing how cropping practices affect viral evolution and disease progression.

Chromosomal inversions are an ubiquitous feature of genetic variation. Theoretical models describe several mechanisms by which inversions can drive adaptation and be maintained as polymorphisms. While inversions have been shown previously to be under selection, or contain genetic variation under selection, the specific phenotypic consequences of inversions leading to their maintenance remain unclear. Here we use genomic sequence and expression data from the Drosophila Genetic Reference Panel to explore the effects of two cosmopolitan inversions, In(2L)t and In(3R)Mo, on patterns of transcriptional variation. We demonstrate that each inversion has a significant effect on transcript abundance for hundreds of genes across the genome. Inversion affected loci (IAL) appear both within inversions as well as on unlinked chromosomes. Importantly, IAL do not appear to be influenced by the previously reported genome-wide expression correlation structure. We found that five genes involved with sterol uptake, four of which are Niemann-Pick Type 2 orthologs, are upregulated in flies with In(3R)Mo but do not have SNPs in LD with the inversion. We speculate that this upregulation is driven by genetic variation in mod(mdg4) that is in LD with In(3R)Mo. We find that there is little evidence for regional or position effect of inversions on gene expression at the chromosomal level but do find evidence for the distal breakpoint of In(3R)Mo interrupting one gene and possibly disassociating the two flanking genes from regulatory elements.

Chromosomal inversions are an ubiquitous feature of genetic variation. Theoretical models describe several mechanisms by which inversions can drive adaptation and be maintained as polymorphisms. While inversions have been shown previously to be under selection, or contain genetic variation under selection, the specific phenotypic consequences of inversions leading to their maintenance remain unclear. Here we use genomic sequence and expression data from the Drosophila Genetic Reference Panel to explore the effects of two cosmopolitan inversions, In(2L)t and In(3R)Mo , on patterns of transcriptional variation. We demonstrate that each inversion has a significant effect on transcript abundance for hundreds of genes across the genome. Inversion affected loci (IAL) appear both within inversions as well as on unlinked chromosomes. Importantly, IAL do not appear to be influenced by the previously reported genome-wide expression correlation structure. We found that five genes involved with sterol uptake, four of which are Niemann-Pick Type 2 orthologs, are upregulated in flies with In(3R)Mo but do not have SNPs in LD with the inversion. We speculate that this upregulation is driven by genetic variation in mod(mdg4) that is in LD with In(3R)Mo . We find that there is little evidence for regional or position effect of inversions on gene expression at the chromosomal level but do find evidence for the distal breakpoint of In(3R)Mo interrupting one gene and possibly disassociating the two flanking genes from regulatory elements.

Very little information exists for long-term changes in genetic variation in natural populations. Here we take the unique opportunity to compare a set of data for SNPs in 15 metabolic genes from eastern US collections of Drosophila melanogaster that span a large latitudinal range and represent two collections separated by 12 to 13 years. We also expand this to a 22-year interval for the Adh gene and approximately 30 years for the G6pd and Pgd genes. During these intervals, five genes showed a statistically significant change in average SNP allele frequency corrected for latitude. While much remains unchanged, we see five genes where latitudinal clines have been lost or gained and two where the slope significantly changes. The long-term frequency shift towards a southern favored Adh S allele reported in Australia populations is not observed in the eastern US over a period of 21 years. There is no general pattern of southern-favored or northern-favored alleles increasing in frequency across the genes. This observation points to the fluid nature of some allelic variation over this time period and the action of selective responses or migration that may be more regional than uniformly imposed across the cline.

There is a connection between nutrient inputs, energy-sensing pathways, lifespan variation and aging. Despite the role of metabolic enzymes in energy homeostasis and their metabolites as nutrient signals, little is known about how their gene expression impacts lifespan. In this report, we use P-element mutagenesis in Drosophila to study the effect on lifespan of reductions in expression of seven central metabolic enzymes, and contrast the effects on normal diet and dietary restriction. The major observation is that for five of seven genes, the reduction of gene expression extends lifespan on one or both diets. Two genes are involved in redox balance, and we observe that lower activity genotypes significantly extend lifespan. The hexokinases also show extension of lifespan with reduced gene activity. Since both affect the ATP/ADP ratio, this connects with the role of AMP-activated protein kinase as an energy sensor in regulating lifespan and mediating caloric restriction. These genes possess significant expression variation in natural populations, and our experimental genotypes span this level of natural activity variation. Our studies link the readout of energy state with the perturbation of the genes of central metabolism and demonstrate their effect on lifespan.