Copy number variations (CNVs) in the Neurexin 1 (NRXN1) gene, which encodes a presynaptic protein involved in neurotransmitter release, are some of the most frequently observed single-gene variants associated with autism spectrum disorder (ASD). To address the functional contribution of NRXN1 CNVs to behavioral phenotypes relevant to ASD, we carried out systematic behavioral phenotyping of an allelic series of Nrxn1 mouse models: one carrying promoter and exon 1 deletion abolishing Nrxn1α transcription, one carrying exon 9 deletion disrupting Nrxn1α protein translation, and one carrying an intronic deletion with no observable effect on Nrxn1α expression. We found that homozygous loss of Nrxn1α resulted in enhanced aggression in males, reduced affiliative social behaviors in females, and significantly altered circadian activities in both sexes. Heterozygous or homozygous loss of Nrxn1α affected the preference for social novelty in male mice, and notably, enhanced repetitive motor skills and motor coordination in both sexes. In contrast, mice bearing an intronic deletion of Nrxn1 did not display alterations in any of the behaviors assessed. These findings demonstrate the importance of Nrxn1α gene dosage in regulating social, circadian, and motor functions, and the variables of sex and genomic positioning of CNVs in the expression of autism-related phenotypes. Importantly, mice with heterozygous loss of Nrxn1, as found in numerous autistic individuals, show an elevated propensity to manifest autism-related phenotypes, supporting the use of models with this genomic architecture to study ASD etiology and assess additional genetic variants associated with autism.

The X-chromosome is among the largest human chromosomes. It differs from autosomes by a number of important features including hemizygosity in males, an almost complete inactivation of one copy in females, and unique patterns of recombination. We used data from the Catalog of Published Genome Wide Association Studies to compare densities of the GWAS-detected SNPs on the X-chromosome and autosomes. The density of GWAS-detected SNPs on the X-chromosome is 6-fold lower compared to the density of the GWAS-detected SNPs on autosomes. Differences between the X-chromosome and autosomes cannot be explained by differences in the overall SNP density, lower X-chromosome coverage by genotyping platforms or low call rate of X-chromosomal SNPs. Similar differences in the density of GWAS-detected SNPs were found in female-only GWASs (e.g. ovarian cancer GWASs). We hypothesized that the lower density of GWAS-detected SNPs on the X-chromosome compared to autosomes is not a result of a methodological bias, e.g. differences in coverage or call rates, but has a real underlying biological reason-a lower density of functional SNPs on the X-chromosome versus autosomes. This hypothesis is supported by the observation that (i) the overall SNP density of X-chromosome is lower compared to the SNP density on autosomes and that (ii) the density of genic SNPs on the X-chromosome is lower compared to autosomes while densities of intergenic SNPs are similar.

Photoreceptor cells (PRCs) are the light-detecting cells of the retina. Such cells can be non-invasively imaged using optical coherence tomography (OCT) which is used in clinical settings to diagnose and monitor ocular diseases. Here we present the largest genome-wide association study of PRC morphology to date utilising quantitative phenotypes extracted from OCT images within the UK Biobank. We discovered 111 loci associated with the thickness of one or more of the PRC layers, many of which had prior associations to ocular phenotypes and pathologies, and 27 with no prior associations. We further identified 10 genes associated with PRC thickness through gene burden testing using exome data. In both cases there was a significant enrichment for genes involved in rare eye pathologies, in particular retinitis pigmentosa. There was evidence for an interaction effect between common genetic variants, VSX2 involved in eye development and PRPH2 known to be involved in retinal dystrophies. We further identified a number of genetic variants with a differential effect across the macular spatial field. Our results suggest a continuum between common and rare variation which impacts retinal structure, sometimes leading to disease.

Hundreds of genetic variants implicated in Mendelian disease have been characterized in dogs and commercial screening is being offered for most of them worldwide. There is typically limited information available regarding the broader population frequency of variants and uncertainty regarding their functional and clinical impact in ancestry backgrounds beyond the discovery breed. Genetic panel screening of disease-associated variants, commercially offered directly to the consumer or via a veterinary clinician, provides an opportunity to establish large-scale cohorts with phenotype data available to address open questions related to variant prevalence and relevance. We screened the largest canine cohort examined in a single study to date (1,054,293 representative dogs from our existing cohort of 3.5 million; a total of 811,628 mixed breed dogs and 242,665 purebreds from more than 150 countries) to examine the prevalence and distribution of a total of 250 genetic disease-associated variants in the general population. Electronic medical records from veterinary clinics were available for 43.5% of the genotyped dogs, enabling the clinical impact of variants to be investigated. We provide detailed frequencies for all tested variants across breeds and find that 57% of dogs carry at least one copy of a studied Mendelian disease-associated variant. Focusing on a subset of variants, we provide evidence of full penetrance for 10 variants, and plausible evidence for clinical significance of 22 variants, on diverse breed backgrounds. Specifically, we report that inherited hypocatalasia is a notable oral health condition, confirm that factor VII deficiency presents as subclinical bleeding propensity and verify two genetic causes of reduced leg length. We further assess genome-wide heterozygosity levels in over 100 breeds, and show that a reduction in genome-wide heterozygosity is associated with an increased Mendelian disease variant load. The accumulated knowledge represents a resource to guide discussions on genetic test relevance by breed.

Adaptive immunity's success relies on the extraordinary diversity of protein receptors on B and T cell membranes. Despite this diversity, the existence of public receptors shared by many individuals gives hope for developing population-wide vaccines and therapeutics. Using probabilistic modeling, we show many of these public receptors are shared by chance in healthy individuals. This predictable overlap is driven not only by biases in the random generation process of receptors, as previously reported, but also by their common functional selection. However, the model underestimates sharing between repertoires of individuals infected with SARS-CoV-2, suggesting strong specific antigen-driven convergent selection. We exploit this discrepancy to identify COVID-associated receptors, which we validate against datasets of receptors with known viral specificity. We study their properties in terms of sequence features and network organization, and use them to design an accurate diagnostic tool for predicting SARS-CoV-2 status from repertoire data.

Genetic studies of disease progression can be used to identify factors that may influence survival or prognosis, which may differ from factors that influence on disease susceptibility. Studies of disease progression feed directly into therapeutics for disease, whereas studies of incidence inform prevention strategies. However, studies of disease progression are known to be affected by collider (also known as "index event") bias since the disease progression phenotype can only be observed for individuals who have the disease. This applies equally to observational and genetic studies, including genome-wide association studies and Mendelian randomisation (MR) analyses. In this paper, our aim is to review several statistical methods that can be used to detect and adjust for index event bias in studies of disease progression, and how they apply to genetic and MR studies using both individual- and summary-level data. Methods to detect the presence of index event bias include the use of negative controls, a comparison of associations between risk factors for incidence in individuals with and without the disease, and an inspection of Miami plots. Methods to adjust for the bias include inverse probability weighting (with individual-level data), or Slope-Hunter and Dudbridge et al.'s index event bias adjustment (when only summary-level data are available). We also outline two approaches for sensitivity analysis. We then illustrate how three methods to minimise bias can be used in practice with two applied examples. Our first example investigates the effects of blood lipid traits on mortality from coronary heart disease, while our second example investigates genetic associations with breast cancer mortality.

The worldwide sheep population comprises more than 1000 breeds. Together, these exhibit a considerable morphological diversity, which has not been extensively investigated at the molecular level. Here, we analyze whole-genome sequencing individuals of 1,098 domestic sheep from 154 breeds, and 69 wild sheep from seven Ovis species. On average, we detected 6.8%, 1.0% and 0.2% introgressed sequence in domestic sheep originating from Iranian mouflon, urial and argali, respectively, with rare introgressions from other wild species. Interestingly, several introgressed haplotypes contributed to the morphological differentiations across sheep breeds, such as a RXFP2 haplotype from Iranian mouflon conferring the spiral horn trait, a MSRB3 haplotype from argali strongly associated with ear morphology, and a VPS13B haplotype probably originating from urial and mouflon possibly associated with facial traits. Our results reveal that introgression events from wild Ovis species contributed to the high rate of morphological differentiation in sheep breeds, but also to individual variation within breeds. We propose that long divergent haplotypes are a ubiquitous source of phenotypic variation that allows adaptation to a variable environment, and that these remain intact in the receiving population probably due to reduced recombination.

Structural variations (SVs) are a key type of cancer genomic alterations, contributing to oncogenesis and progression of many cancers, including colorectal cancer (CRC). However, SVs in CRC remain difficult to be reliably detected due to limited SV-detection capacity of the commonly used short-read sequencing. This study investigated the somatic SVs in 21 pairs of CRC samples by Nanopore whole-genome long-read sequencing. 5200 novel somatic SVs from 21 CRC patients (494 SVs / patient) were identified. A 4.9-Mbp long inversion that silences APC expression (confirmed by RNA-seq) and an 11.2-kbp inversion that structurally alters CFTR were identified. Two novel gene fusions that might functionally impact the oncogene RNF38 and the tumor-suppressor SMAD3 were detected. RNF38 fusion possesses metastasis-promoting ability confirmed by in vitro migration and invasion assay, and in vivo metastasis experiments. This work highlighted the various applications of long-read sequencing in cancer genome analysis, and shed new light on how somatic SVs structurally alter critical genes in CRC. The investigation on somatic SVs via nanopore sequencing revealed the potential of this genomic approach in facilitating precise diagnosis and personalized treatment of CRC.

The molecular mechanism of tumor metastasis, especially how metastatic tumor cells colonize in a distant site, remains poorly understood. Here we reported that ARHGAP15, a Rho GTPase activating protein, enhanced gastric cancer (GC) metastatic colonization, which was quite different from its reported role as a tumor suppressor gene in other cancers. It was upregulated in metastatic lymph nodes and significantly associated with a poor prognosis. Ectopic expression of ARHGAP15 promoted metastatic colonization of gastric cancer cells in murine lungs and lymph nodes in vivo or protected cells from oxidative-related death in vitro. However, genetic downregulation of ARHGAP15 had the opposite effect. Mechanistically, ARHGAP15 inactivated RAC1 and then decreased intracellular accumulation of reactive oxygen species (ROS), thus enhancing the antioxidant capacity of colonizing tumor cells under oxidative stress. This phenotype could be phenocopied by inhibition of RAC1 or rescued by the introduction of constitutively active RAC1 into cells. Taken together, these findings suggested a novel role of ARHGAP15 in promoting gastric cancer metastasis by quenching ROS through inhibiting RAC1 and its potential value for prognosis estimation and targeted therapy.

Chromosome movements and licensing of synapsis must be tightly regulated during early meiosis to ensure accurate chromosome segregation and avoid aneuploidy, although how these steps are coordinated is not fully understood. Here we show that GRAS-1, the worm homolog of mammalian GRASP/Tamalin and CYTIP, coordinates early meiotic events with cytoskeletal forces outside the nucleus. GRAS-1 localizes close to the nuclear envelope (NE) in early prophase I and interacts with NE and cytoskeleton proteins. Delayed homologous chromosome pairing, synaptonemal complex (SC) assembly, and DNA double-strand break repair progression are partially rescued by the expression of human CYTIP in gras-1 mutants, supporting functional conservation. However, Tamalin, Cytip double knockout mice do not exhibit obvious fertility or meiotic defects, suggesting evolutionary differences between mammals. gras-1 mutants show accelerated chromosome movement during early prophase I, implicating GRAS-1 in regulating chromosome dynamics. GRAS-1-mediated regulation of chromosome movement is DHC-1-dependent, placing it acting within the LINC-controlled pathway, and depends on GRAS-1 phosphorylation at a C-terminal S/T cluster. We propose that GRAS-1 coordinates the early steps of homology search and licensing of SC assembly by regulating the pace of chromosome movement in early prophase I.

Mediation analysis is commonly used to identify mechanisms and intermediate factors between causes and outcomes. Studies drawing on polygenic scores (PGSs) can readily employ traditional regression-based procedures to assess whether trait M mediates the relationship between the genetic component of outcome Y and outcome Y itself. However, this approach suffers from attenuation bias, as PGSs capture only a (small) part of the genetic variance of a given trait. To overcome this limitation, we developed MA-GREML: a method for Mediation Analysis using Genome-based Restricted Maximum Likelihood (GREML) estimation. Using MA-GREML to assess mediation between genetic factors and traits comes with two main advantages. First, we circumvent the limited predictive accuracy of PGSs that regression-based mediation approaches suffer from. Second, compared to methods employing summary statistics from genome-wide association studies, the individual-level data approach of GREML allows to directly control for confounders of the association between M and Y. In addition to typical GREML parameters (e.g., the genetic correlation), MA-GREML estimates (i) the effect of M on Y, (ii) the direct effect (i.e., the genetic variance of Y that is not mediated by M), and (iii) the indirect effect (i.e., the genetic variance of Y that is mediated by M). MA-GREML also provides standard errors of these estimates and assesses the significance of the indirect effect. We use analytical derivations and simulations to show the validity of our approach under two main assumptions, viz., that M precedes Y and that environmental confounders of the association between M and Y are controlled for. We conclude that MA-GREML is an appropriate tool to assess the mediating role of trait M in the relationship between the genetic component of Y and outcome Y. Using data from the US Health and Retirement Study, we provide evidence that genetic effects on Body Mass Index (BMI), cognitive functioning and self-reported health in later life run partially through educational attainment. For mental health, we do not find significant evidence for an indirect effect through educational attainment. Further analyses show that the additive genetic factors of these four outcomes do partially (cognition and mental health) and fully (BMI and self-reported health) run through an earlier realization of these traits.

Transposable elements (TE) are selfish genetic elements that can cause harmful mutations. In Drosophila, it has been estimated that half of all spontaneous visible marker phenotypes are mutations caused by TE insertions. Several factors likely limit the accumulation of exponentially amplifying TEs within genomes. First, synergistic interactions between TEs that amplify their harm with increasing copy number are proposed to limit TE copy number. However, the nature of this synergy is poorly understood. Second, because of the harm posed by TEs, eukaryotes have evolved systems of small RNA-based genome defense to limit transposition. However, as in all immune systems, there is a cost of autoimmunity and small RNA-based systems that silence TEs can inadvertently silence genes flanking TE insertions. In a screen for essential meiotic genes in Drosophila melanogaster, a truncated Doc retrotransposon within a neighboring gene was found to trigger the germline silencing of ald, the Drosophila Mps1 homolog, a gene essential for proper chromosome segregation in meiosis. A subsequent screen for suppressors of this silencing identified a new insertion of a Hobo DNA transposon in the same neighboring gene. Here we describe how the original Doc insertion triggers flanking piRNA biogenesis and local gene silencing. We show that this local gene silencing occurs in cis and is dependent on deadlock, a component of the Rhino-Deadlock-Cutoff (RDC) complex, to trigger dual-strand piRNA biogenesis at TE insertions. We further show how the additional Hobo insertion leads to de-silencing by reducing flanking piRNA biogenesis triggered by the original Doc insertion. These results support a model of TE-mediated gene silencing by piRNA biogenesis in cis that depends on local determinants of transcription. This may explain complex patterns of off-target gene silencing triggered by TEs within populations and in the laboratory. It also provides a mechanism of sign epistasis among TE insertions, illuminates the complex nature of their interactions and supports a model in which off-target gene silencing shapes the evolution of the RDC complex.

X-chromosome inactivation (XCI) silences one X in female cells to balance sex-differences in X-dosage. A subset of X-linked genes escape XCI, but the extent to which this phenomenon occurs and how it varies across tissues and in a population is as yet unclear. To characterize incidence and variability of escape across individuals and tissues, we conducted a transcriptomic study of escape in adipose, skin, lymphoblastoid cell lines and immune cells in 248 healthy individuals exhibiting skewed XCI. We quantify XCI escape from a linear model of genes' allelic fold-change and XIST-based degree of XCI skewing. We identify 62 genes, including 19 lncRNAs, with previously unknown patterns of escape. We find a range of tissue-specificity, with 11% of genes escaping XCI constitutively across tissues and 23% demonstrating tissue-restricted escape, including cell type-specific escape across immune cells of the same individual. We also detect substantial inter-individual variability in escape. Monozygotic twins share more similar escape than dizygotic twins, indicating that genetic factors may underlie inter-individual differences in escape. However, discordant escape also occurs within monozygotic co-twins, suggesting environmental factors also influence escape. Altogether, these data indicate that XCI escape is an under-appreciated source of transcriptional differences, and an intricate phenotype impacting variable trait expressivity in females.

Circadian clock and chromatin-remodeling complexes are tightly intertwined systems that regulate rhythmic gene expression. The circadian clock promotes rhythmic expression, timely recruitment, and/or activation of chromatin remodelers, while chromatin remodelers regulate accessibility of clock transcription factors to the DNA to influence expression of clock genes. We previously reported that the BRAHMA (BRM) chromatin-remodeling complex promotes the repression of circadian gene expression in Drosophila. In this study, we investigated the mechanisms by which the circadian clock feeds back to modulate daily BRM activity. Using chromatin immunoprecipitation, we observed rhythmic BRM binding to clock gene promoters despite constitutive BRM protein expression, suggesting that factors other than protein abundance are responsible for rhythmic BRM occupancy at clock-controlled loci. Since we previously reported that BRM interacts with two key clock proteins, CLOCK (CLK) and TIMELESS (TIM), we examined their effect on BRM occupancy to the period (per) promoter. We observed reduced BRM binding to the DNA in clk null flies, suggesting that CLK is involved in enhancing BRM occupancy to initiate transcriptional repression at the conclusion of the activation phase. Additionally, we observed reduced BRM binding to the per promoter in flies overexpressing TIM, suggesting that TIM promotes BRM removal from DNA. These conclusions are further supported by elevated BRM binding to the per promoter in flies subjected to constant light and experiments in Drosophila tissue culture in which the levels of CLK and TIM are manipulated. In summary, this study provides new insights into the reciprocal regulation between the circadian clock and the BRM chromatin-remodeling complex.

Neonicotinoid insecticides target insect nicotinic acetylcholine receptors (nAChRs) and their adverse effects on non-target insects are of serious concern. We recently found that cofactor TMX3 enables robust functional expression of insect nAChRs in Xenopus laevis oocytes and showed that neonicotinoids (imidacloprid, thiacloprid, and clothianidin) exhibited agonist actions on some nAChRs of the fruit fly (Drosophila melanogaster), honeybee (Apis mellifera) and bumblebee (Bombus terrestris) with more potent actions on the pollinator nAChRs. However, other subunits from the nAChR family remain to be explored. We show that the Dα3 subunit co-exists with Dα1, Dα2, Dβ1, and Dβ2 subunits in the same neurons of adult D. melanogaster, thereby expanding the possible nAChR subtypes in these cells alone from 4 to 12. The presence of Dα1 and Dα2 subunits reduced the affinity of imidacloprid, thiacloprid, and clothianidin for nAChRs expressed in Xenopus laevis oocytes, whereas the Dα3 subunit enhanced it. RNAi targeting Dα1, Dα2 or Dα3 in adults reduced expression of targeted subunits but commonly enhanced Dβ3 expression. Also, Dα1 RNAi enhanced Dα7 expression, Dα2 RNAi reduced Dα1, Dα6, and Dα7 expression and Dα3 RNAi reduced Dα1 expression while enhancing Dα2 expression, respectively. In most cases, RNAi treatment of either Dα1 or Dα2 reduced neonicotinoid toxicity in larvae, but Dα2 RNAi enhanced neonicotinoid sensitivity in adults reflecting the affinity-reducing effect of Dα2. Substituting each of Dα1, Dα2, and Dα3 subunits by Dα4 or Dβ3 subunit mostly increased neonicotinoid affinity and reduced efficacy. These results are important because they indicate that neonicotinoid actions involve the integrated activity of multiple nAChR subunit combinations and counsel caution in interpreting neonicotinoid actions simply in terms of toxicity.

Animal traits develop through the expression and action of numerous regulatory and realizator genes that comprise a gene regulatory network (GRN). For each GRN, its underlying patterns of gene expression are controlled by cis-regulatory elements (CREs) that bind activating and repressing transcription factors. These interactions drive cell-type and developmental stage-specific transcriptional activation or repression. Most GRNs remain incompletely mapped, and a major barrier to this daunting task is CRE identification. Here, we used an in silico method to identify predicted CREs (pCREs) that comprise the GRN which governs sex-specific pigmentation of Drosophila melanogaster. Through in vivo assays, we demonstrate that many pCREs activate expression in the correct cell-type and developmental stage. We employed genome editing to demonstrate that two CREs control the pupal abdomen expression of trithorax, whose function is required for the dimorphic phenotype. Surprisingly, trithorax had no detectable effect on this GRN's key trans-regulators, but shapes the sex-specific expression of two realizator genes. Comparison of sequences orthologous to these CREs supports an evolutionary scenario where these trithorax CREs predated the origin of the dimorphic trait. Collectively, this study demonstrates how in silico approaches can shed novel insights on the GRN basis for a trait's development and evolution.

Specification of migratory cell fate from a stationary population is complex and indispensable both for metazoan development as well for the progression of the pathological condition like tumor metastasis. Though this cell fate transformation is widely prevalent, the molecular understanding of this phenomenon remains largely elusive. We have employed the model of border cells (BC) in Drosophila oogenesis and identified germline activity of an RNA binding protein, Cup that limits acquisition of migratory cell fate from the neighbouring follicle epithelial cells. As activation of JAK-STAT in the follicle cells is critical for BC specification, our data suggest that Cup, non-cell autonomously restricts the domain of JAK-STAT by activating Notch in the follicle cells. Employing genetics and Delta endocytosis assay, we demonstrate that Cup regulates Delta recycling in the nurse cells through Rab11GTPase thus facilitating Notch activation in the adjacent follicle cells. Since Notch and JAK-STAT are antagonistic, we propose that germline Cup functions through Notch and JAK-STAT to modulate BC fate specification from their static epithelial progenitors.

Yeast cells maintain an intricate network of nutrient signaling pathways enabling them to integrate information on the availability of different nutrients and adjust their metabolism and growth accordingly. Cells that are no longer capable of integrating this information, or that are unable to make the necessary adaptations, will cease growth and eventually die. Here, we studied the molecular basis underlying the synthetic lethality caused by loss of the protein kinase Sch9, a key player in amino acid signaling and proximal effector of the conserved growth-regulatory TORC1 complex, when combined with either loss of the cyclin-dependent kinase (CDK) Pho85 or loss of its inhibitor Pho81, which both have pivotal roles in phosphate sensing and cell cycle regulation. We demonstrate that it is specifically the CDK-cyclin pair Pho85-Pho80 or the partially redundant CDK-cyclin pairs Pho85-Pcl6/Pcl7 that become essential for growth when Sch9 is absent. Interestingly, the respective three CDK-cyclin pairs regulate the activity and distribution of the phosphatidylinositol-3 phosphate 5-kinase Fab1 on endosomes and vacuoles, where it generates phosphatidylinositol-3,5 bisphosphate that serves to recruit both TORC1 and its substrate Sch9. In addition, Pho85-Pho80 directly phosphorylates Sch9 at Ser726, and to a lesser extent at Thr723, thereby priming Sch9 for its subsequent phosphorylation and activation by TORC1. The TORC1-Sch9 signaling branch therefore integrates Pho85-mediated information at different levels. In this context, we also discovered that loss of the transcription factor Pho4 rescued the synthetic lethality caused by loss of Pho85 and Sch9, indicating that both signaling pathways also converge on Pho4, which appears to be wired to a feedback loop involving the high-affinity phosphate transporter Pho84 that fine-tunes Sch9-mediated responses.

Pharmacological vitamin C (VC) is a potential natural compound for cancer treatment. However, the mechanism underlying its antitumor effects remains unclear. In this study, we found that pharmacological VC significantly inhibits the mTOR (including mTORC1 and mTORC2) pathway activation and promotes GSK3-FBXW7-mediated Rictor ubiquitination and degradation by increasing the cellular ROS. Moreover, we identified that HMOX1 is a checkpoint for pharmacological-VC-mediated mTOR inactivation, and the deletion of FBXW7 or HMOX1 suppresses the regulation of pharmacological VC on mTOR activation, cell size, cell viability, and autophagy. More importantly, it was observed that the inhibition of mTOR by pharmacological VC supplementation in vivo produces positive therapeutic responses in tumor growth, while HMOX1 deficiency rescues the inhibitory effect of pharmacological VC on tumor growth. These results demonstrate that VC influences cellular activities and tumor growth by inhibiting the mTOR pathway through Rictor and HMOX1, which may have therapeutic potential for cancer treatment.

Recently, distinct mutational footprints observed in metastatic tumors, secondary malignancies and normal human tissues have been demonstrated to be caused by the exposure to several chemotherapeutic drugs. These characteristic mutations originate from specific lesions caused by these chemicals to the DNA of exposed cells. However, it is unknown whether the exposure to these chemotherapies leads to a specific footprint of larger chromosomal aberrations. Here, we address this question exploiting whole genome sequencing data of metastatic tumors obtained from patients exposed to different chemotherapeutic drugs. As a result, we discovered a specific copy number footprint across tumors from patients previously exposed to platinum-based therapies. This footprint is characterized by a significant increase in the number of chromosomal fragments of copy number 1-4 and size smaller than 10 Mb in exposed tumors with respect to their unexposed counterparts (median 14-387% greater across tumor types). The number of chromosomal fragments characteristic of the platinum-associated CN footprint increases significantly with the activity of the well known platinum-related footprint of single nucleotide variants across exposed tumors.

[This corrects the article DOI: 10.1371/journal.pgen.1010150.].

Admixture graphs are mathematical structures that describe the ancestry of populations in terms of divergence and merging (admixing) of ancestral populations as a graph. An admixture graph consists of a graph topology, branch lengths, and admixture proportions. The branch lengths and admixture proportions can be estimated using numerous numerical optimization methods, but inferring the topology involves a combinatorial search for which no polynomial algorithm is known. In this paper, we present a reversible jump MCMC algorithm for sampling high-probability admixture graphs and show that this approach works well both as a heuristic search for a single best-fitting graph and for summarizing shared features extracted from posterior samples of graphs. We apply the method to 11 Native American and Siberian populations and exploit the shared structure of high-probability graphs to characterize the relationship between Saqqaq, Inuit, Koryaks, and Athabascans. Our analyses show that the Saqqaq is not a good proxy for the previously identified gene flow from Arctic people into the Na-Dene speaking Athabascans.

Integrative and conjugative elements (ICEs) are major contributors to genome plasticity in bacteria. ICEs reside integrated in the chromosome of a host bacterium and are passively propagated during chromosome replication and cell division. When activated, ICEs excise from the chromosome and may be transferred through the ICE-encoded conjugation machinery into a recipient cell. Integration into the chromosome of the new host generates a stable transconjugant. Although integration into the chromosome of a new host is critical for the stable acquisition of ICEs, few studies have directly investigated the molecular events that occur in recipient cells during generation of a stable transconjugant. We found that integration of ICEBs1, an ICE of Bacillus subtilis, occurred several generations after initial transfer to a new host. Premature integration in new hosts led to cell death and hence decreased fitness of the ICE and transconjugants. Host lethality due to premature integration was caused by rolling circle replication that initiated in the integrated ICEBs1 and extended into the host chromosome, resulting in catastrophic genome instability. Our results demonstrate that the timing of integration of an ICE is linked to cessation of autonomous replication of the ICE, and that perturbing this linkage leads to a decrease in ICE and host fitness due to a loss of viability of transconjugants. Linking integration to cessation of autonomous replication appears to be a conserved regulatory scheme for mobile genetic elements that both replicate and integrate into the chromosome of their host.

Dietary nutrient composition is essential for shaping important fitness traits and behaviours. Many organisms are protein limited, and for Drosophila melanogaster this limitation manifests at the level of the single most limiting essential Amino Acid (AA) in the diet. The identity of this AA and its effects on female fecundity is readily predictable by a procedure called exome matching in which the sum of AAs encoded by a consumer's exome is used to predict the relative proportion of AAs required in its diet. However, the exome matching calculation does not weight AA contributions to the overall profile by protein size or expression. Here, we update the exome matching calculation to include these weightings. Surprisingly, although nearly half of the transcriptome is differentially expressed when comparing male and female flies, we found that creating transcriptome-weighted exome matched diets for each sex did not enhance their fecundity over that supported by exome matching alone. These data indicate that while organisms may require different amounts of dietary protein across conditions, the relative proportion of the constituent AAs remains constant. Interestingly, we also found that exome matched AA profiles are generally conserved across taxa and that the composition of these profiles might be explained by energetic and elemental limitations on microbial AA synthesis. Thus, it appears that ecological constraints amongst autotrophs shape the relative proportion of AAs that are available across trophic levels and that this constrains biomass composition.