Variation in Position Effect Variegation Within a Natural Population

The Structure of Simple Satellite Variation in the Human Genome and Its Correlation With Centromere Ancestry

Although repetitive DNA forms much of the human genome, its study is challenging due to limitations in assembly and alignment of repetitive short-reads. We have deployed k-Seek, software that detects tandem repeats embedded in single reads, on 2,504 human genomes from the 1,000 Genomes Project to quantify the variation and abundance of simple satellites (repeat units <20 bp). We find that the ancestral monomer of Human Satellite 3 makes up the largest portion of simple satellite content in humans (mean of ∼8 Mb). We discovered ∼50,000 rare tandem repeats that are not detected in the T2T-CHM13v2.0 assembly, including undescribed variants of telomericand pericentromeric repeats. We find broad homogeneity of the most abundant repeats across populations, except for AG-rich repeats which are more abundant in African individuals. We also find cliques of highly similar AG- and AT-rich satellites that are interspersed and form higher-order structures that covary in copy number across individuals, likely through concerted amplification via unequal exchange. Finally, we use pericentromeric polymorphisms to estimate centromeric genetic relatedness between individuals and find a strong predictive relationship between centromeric lineages and pericentromeric simple satellite abundances. In particular, ancestral monomers of Human Satellite 2 and Human Satellite 3 abundances correlate with clusters of centromeric ancestry on chromosome 16 and chromosome 9, with some clusters structured by population. These results provide new descriptions of the population dynamics that underlie the evolution of simple satellites in humans.

Natural genetic variation in Drosophila melanogaster reveals genes associated with Coxiella burnetii infection

The gram-negative bacterium Coxiella burnetii is the causative agent of Query (Q) fever in humans and coxiellosis in livestock. Host genetics are associated with C. burnetii pathogenesis both in humans and animals; however, it remains unknown if specific genes are associated with severity of infection. We employed the Drosophila Genetics Reference Panel to perform a genome-wide association study to identify host genetic variants that affect host survival to C. burnetii infection. The genome-wide association study identified 64 unique variants (P < 10⁻⁵) associated with 25 candidate genes. We examined the role each candidate gene contributes to host survival during C. burnetii infection using flies carrying a null mutation or RNAi knockdown of each candidate. We validated 15 of the 25 candidate genes using at least one method. This is the first report establishing involvement of many of these genes or their homologs with C. burnetii susceptibility in any system. Among the validated genes, FER and tara play roles in the JAK/STAT, JNK, and decapentaplegic/TGF-β signaling pathways which are components of known innate immune responses to C. burnetii infection. CG42673 and DIP-ε play roles in bacterial infection and synaptic signaling but have no previous association with C. burnetii pathogenesis. Furthermore, since the mammalian ortholog of CG13404 (PLGRKT) is an important regulator of macrophage function, CG13404 could play a role in host susceptibility to C. burnetii through hemocyte regulation. These insights provide a foundation for further investigation regarding the genetics of C. burnetii susceptibility across a wide variety of hosts.

Integrating GWAS and Transcriptomics to Identify the Molecular Underpinnings of Thermal Stress Responses in Drosophila melanogaster

Thermal tolerance of an organism depends on both the ability to dynamically adjust to a thermal stress and preparatory developmental processes that enhance thermal resistance. However, the extent to which standing genetic variation in thermal tolerance alleles influence dynamic stress responses vs. preparatory processes is unknown. Here, using the model species Drosophila melanogaster, we used a combination of Genome Wide Association mapping (GWAS) and transcriptomic profiling to characterize whether genes associated with thermal tolerance are primarily involved in dynamic stress responses or preparatory processes that influence physiological condition at the time of thermal stress. To test our hypotheses, we measured the critical thermal minimum (CTmin) and critical thermal maximum (CTmax) of 100 lines of the Drosophila Genetic Reference Panel (DGRP) and used GWAS to identify loci that explain variation in thermal limits. We observed greater variation in lower thermal limits, with CTmin ranging from 1.81 to 8.60°C, while CTmax ranged from 38.74 to 40.64°C. We identified 151 and 99 distinct genes associated with CTmin and CTmax, respectively, and there was strong support that these genes are involved in both dynamic responses to thermal stress and preparatory processes that increase thermal resistance. Many of the genes identified by GWAS were involved in the direct transcriptional response to thermal stress (72/151 for cold; 59/99 for heat), and overall GWAS candidates were more likely to be differentially expressed than other genes. Further, several GWAS candidates were regulatory genes that may participate in the regulation of stress responses, and gene ontologies related to development and morphogenesis were enriched, suggesting many of these genes influence thermal tolerance through effects on development and physiological status. Overall, our results suggest that thermal tolerance alleles can influence both dynamic plastic responses to thermal stress and preparatory processes that improve thermal resistance. These results also have utility for directly comparing GWAS and transcriptomic approaches for identifying candidate genes associated with thermal tolerance.

The impact of genetic background and cell lineage on the level and pattern of gene expression in position effect variegation

Background

Chromatin-based transcriptional silencing is often described as a stochastic process, largely because of the mosaic expression observed in position effect variegation (PEV), where a euchromatic reporter gene is silenced in some cells as a consequence of juxtaposition with heterochromatin. High levels of variation in PEV phenotypes are commonly observed in reporter stocks. To ascertain whether background mutations are the major contributors to this variation, we asked how much of the variation is determined by genetic variants segregating in the population, examining both the level and pattern of expression using the fruit fly, Drosophila melanogaster, as the model.

Results

Using selective breeding of a fourth chromosome PEV reporter line, 39C-12, we isolated two inbred lines exhibiting contrasting degrees of variegation (A1: low expression, D1: high expression). Within each inbred population, remarkable similarity is observed in the degree of variegation: 90% of the variation between the two inbred lines in the degree of silencing can be explained by genotype. Further analyses suggest that this result reflects the combined effect of multiple independent trans-acting loci. While the initial observations are based on a PEV phenotype scored in the fly eye (hsp70-white reporter), similar degrees of silencing were observed using a beta-gal reporter scored across the whole fly. Further, the pattern of variegation becomes almost identical within each inbred line; significant pigment enrichment in the same quadrant of the eye was found for both A1 and D1 lines despite different degrees of expression.

Conclusions

The results indicate that background genetic variants play the major role in determining the variable degrees of PEV commonly observed in laboratory stocks. Interestingly, not only does the degree of variegation become consistent in inbred lines, the patterns of variegation also appear similar. Combining these observations with the spreading model for local heterochromatin formation, we propose an augmented stochastic model to describe PEV in which the genetic background drives the overall level of silencing, working with the cell lineage-specific regulatory environment to determine the on/off probability at the reporter locus in each cell. This model acknowledges cell type-specific events in the context of broader genetic impacts on heterochromatin formation.

Sex-specific adaptation and genomic responses to Y chromosome presence in female reproductive and neural tissues

Y chromosomes typically harbour a small number of genes and an abundance of repetitive sequences. In Drosophila, the Y chromosome comprises multimegabase long segments of repetitive DNA and a handful of protein-coding genes. In mammals, the Y chromosome also harbours a disproportionally high abundance of repeats. Here, we built on a Drosophila melanogaster model in which the Y chromosome is decoupled from sexual determination. Genotypes were genetically identical for the autosomes, X chromosome, and mitochondria, but differ by the presence or dose of the Y chromosome. Addition of an extra Y chromosome had limited impact in males. However, the presence of a Y chromosome in females induced a disproportionate response in genes expressed in the ovaries as well as genes encoded by the mitochondrial genome. Furthermore, the data revealed significant consequences of Y chromosome presence in larvae neuronal tissue. This included the repression of genes implicated in reproductive behaviour, courtship, mating and synaptic function. Our findings exhibit the Y chromosome as a hotspot for sex-specific adaptation. They suggest roles for natural selection on Y-linked genetic elements exerting impact on sex-specific tissues as well as somatic tissues shared by males and females.

The Y Chromosome Modulates Splicing and Sex-Biased Intron Retention Rates in Drosophila

The Drosophila Y chromosome is a 40-Mb segment of mostly repetitive DNA; it harbors a handful of protein-coding genes and a disproportionate amount of satellite repeats, transposable elements, and multicopy DNA arrays. Intron retention (IR) is a type of alternative splicing (AS) event by which one or more introns remain within the mature transcript. IR recently emerged as a deliberate cellular mechanism to modulate gene expression levels and has been implicated in multiple biological processes. However, the extent of sex differences in IR and the contribution of the Y chromosome to the modulation of AS and IR rates has not been addressed. Here we showed pervasive IR in the fruit fly Drosophila melanogaster with thousands of novel IR events, hundreds of which displayed extensive sex bias. The data also revealed an unsuspected role for the Y chromosome in the modulation of AS and IR. The majority of sex-biased IR events introduced premature termination codons and the magnitude of sex bias was associated with gene expression differences between the sexes. Surprisingly, an extra Y chromosome in males (X^YY genotype) or the presence of a Y chromosome in females (X^XY genotype) significantly modulated IR and recapitulated natural differences in IR between the sexes. Our results highlight the significance of sex-biased IR in tuning sex differences and the role of the Y chromosome as a source of variable IR rates between the sexes. Modulation of splicing and IR rates across the genome represent new and unexpected outcomes of the Drosophila Y chromosome.