Comparison of whole genomes has revealed large and frequent changes in the size of gene families. These changes occur because of high rates of both gene gain (via duplication) and loss (via deletion or pseudogenization), as well as the evolution of entirely new genes. Here we use the genomes of 12 fully sequenced Drosophila species to study the gain and loss of genes at unprecedented resolution. We find large numbers of both gains and losses, with over 40% of all gene families differing in size among the Drosophila. Approximately 17 genes are estimated to be duplicated and fixed in a genome every million years, a rate on par with that previously found in both yeast and mammals. We find many instances of extreme expansions or contractions in the size of gene families, including the expansion of several sex- and spermatogenesis-related families in D. melanogaster that also evolve under positive selection at the nucleotide level. Newly evolved gene families in our dataset are associated with a class of testes-expressed genes known to have evolved de novo in a number of cases. Gene family comparisons also allow us to identify a number of annotated D. melanogaster genes that are unlikely to encode functional proteins, as well as to identify dozens of previously unannotated D. melanogaster genes with conserved homologs in the other Drosophila. Taken together, our results demonstrate that the apparent stasis in total gene number among species has masked rapid turnover in individual gene gain and loss. It is likely that this genomic revolving door has played a large role in shaping the morphological, physiological, and metabolic differences among species.

The rapid digitization of genealogical and medical records enables the assembly of extremely large pedigree records spanning millions of individuals and trillions of pairs of relatives. Such pedigrees provide the opportunity to investigate the sociological and epidemiological history of human populations in scales much larger than previously possible. Linear mixed models (LMMs) are routinely used to analyze extremely large animal and plant pedigrees for the purposes of selective breeding. However, LMMs have not been previously applied to analyze population-scale human family trees. Here, we present Sparse Cholesky factorIzation LMM (Sci-LMM), a modeling framework for studying population-scale family trees that combines techniques from the animal and plant breeding literature and from human genetics literature. The proposed framework can construct a matrix of relationships between trillions of pairs of individuals and fit the corresponding LMM in several hours. We demonstrate the capabilities of Sci-LMM via simulation studies and by estimating the heritability of longevity and of reproductive fitness (quantified via number of children) in a large pedigree spanning millions of individuals and over five centuries of human history. Sci-LMM provides a unified framework for investigating the epidemiological history of human populations via genealogical records.

Endocytosis, the process by which cells internalize plasma membrane and associated cargo, is regulated extensively by posttranslational modifications. Previous studies suggested the potential involvement of scores of protein kinases in endocytic control, of which only a few have been validated in vivo. Here we show that the conserved NIMA-related kinases NEKL-2/NEK8/9 and NEKL-3/NEK6/7 (the NEKLs) control clathrin-mediated endocytosis in C. elegans. Loss of NEKL-2 or NEKL-3 activities leads to penetrant larval molting defects and to the abnormal localization of trafficking markers in arrested larvae. Using an auxin-based degron system, we also find that depletion of NEKLs in adult-stage C. elegans leads to gross clathrin mislocalization and to a dramatic reduction in clathrin mobility at the apical membrane. Using a non-biased genetic screen to identify suppressors of nekl molting defects, we identified several components and regulators of AP2, the major clathrin adapter complex acting at the plasma membrane. Strikingly, reduced AP2 activity rescues both nekl mutant molting defects as well as associated trafficking phenotypes, whereas increased levels of active AP2 exacerbate nekl defects. Moreover, in a unique example of mutual suppression, NEKL inhibition alleviates defects associated with reduced AP2 activity, attesting to the tight link between NEKL and AP2 functions. We also show that NEKLs are required for the clustering and internalization of membrane cargo required for molting. Notably, we find that human NEKs can rescue molting and trafficking defects in nekl mutant worms, suggesting that the control of intracellular trafficking is an evolutionarily conserved function of NEK family kinases.

Transcription of the human mitochondrial genome and correct processing of the two long polycistronic transcripts are crucial for oxidative phosphorylation. According to the tRNA punctuation model, nucleolytic processing of these large precursor transcripts occurs mainly through the excision of the tRNAs that flank most rRNAs and mRNAs. However, some mRNAs are not punctuated by tRNAs, and it remains largely unknown how these non-canonical junctions are resolved. The FASTK family proteins are emerging as key players in non-canonical RNA processing. Here, we have generated human cell lines carrying single or combined knockouts of several FASTK family members to investigate their roles in non-canonical RNA processing. The most striking phenotypes were obtained with loss of FASTKD4 and FASTKD5 and with their combined double knockout. Comprehensive mitochondrial transcriptome analyses of these cell lines revealed a defect in processing at several canonical and non-canonical RNA junctions, accompanied by an increase in specific antisense transcripts. Loss of FASTKD5 led to the most severe phenotype with marked defects in mitochondrial translation of key components of the electron transport chain complexes and in oxidative phosphorylation. We reveal that the FASTK protein family members are crucial regulators of non-canonical junction and non-coding mitochondrial RNA processing.

Archaeal pleomorphic viruses belonging to the Pleolipoviridae family represent an enigmatic group as they exhibit unique genomic features and are thought to have evolved through recombination with different archaeal plasmids. However, most of our understanding of the diversity and evolutionary trajectories of this clade comes from a handful of isolated representatives. Here we present 164 new genomes of pleolipoviruses obtained from metagenomic data of Australian hypersaline lakes and publicly available metagenomic data. We perform a comprehensive analysis on the diversity and evolutionary relationships of the newly discovered viruses and previously described pleolipoviruses. We propose to classify the viruses into five genera within the Pleolipoviridae family, with one new genus represented only by virus genomes retrieved in this study. Our data support the current hypothesis that pleolipoviruses reshaped their genomes through recombining with multiple different groups of plasmids, which is reflected in the diversity of their predicted replication strategies. We show that the proposed genus Epsilonpleolipovirus has evolutionary ties to pRN1-like plasmids from Sulfolobus, suggesting that this group could be infecting other archaeal phyla. Interestingly, we observed that the genome size of pleolipoviruses is correlated to the presence or absence of an integrase. Analyses of the host range revealed that all but one virus exhibit an extremely narrow range, and we show that the predicted tertiary structure of the spike protein is strongly associated with the host family, suggesting a specific adaptation to the host S-layer glycoprotein organization.

Antibiotic resistance is a serious global problem. Antibiotic resistance genes (ARG), which are widespread in environmental bacteria, can be transferred to pathogenic bacteria via horizontal gene transfer (HGT). Gut microbiomes are especially apt for the emergence and dissemination of ARG. Conjugation is the HGT route that is predominantly responsible for the spread of ARG. Little is known about conjugative elements of Gram-positive bacteria, including those of the phylum Firmicutes, which are abundantly present in gut microbiomes. A critical step in the conjugation process is the relaxase-mediated site- and strand-specific nick in the oriT region of the conjugative element. This generates a single-stranded DNA molecule that is transferred from the donor to the recipient cell via a connecting channel. Here we identified and characterized the relaxosome components oriT and the relaxase of the conjugative plasmid pLS20 of the Firmicute Bacillus subtilis. We show that the relaxase gene, named relLS20, is essential for conjugation, that it can function in trans and provide evidence that Tyr26 constitutes the active site residue. In vivo and in vitro analyses revealed that the oriT is located far upstream of the relaxase gene and that the nick site within oriT is located on the template strand of the conjugation genes. Surprisingly, the RelLS20 shows very limited similarity to known relaxases. However, more than 800 genes to which no function had been attributed so far are predicted to encode proteins showing significant similarity to RelLS20. Interestingly, these putative relaxases are encoded almost exclusively in Firmicutes bacteria. Thus, RelLS20 constitutes the prototype of a new family of relaxases. The identification of this novel relaxase family will have an important impact in different aspects of future research in the field of HGT in Gram-positive bacteria in general, and specifically in the phylum of Firmicutes, and in gut microbiome research.

PAX5, one of nine members of the mammalian paired box (PAX) family of transcription factors, plays an important role in B cell development. Approximately one-third of individuals with pre-B acute lymphoblastic leukemia (ALL) acquire heterozygous inactivating mutations of PAX5 in malignant cells, and heterozygous germline loss-of-function PAX5 mutations cause autosomal dominant predisposition to ALL. At least in mice, Pax5 is required for pre-B cell maturation, and leukemic remission occurs when Pax5 expression is restored in a Pax5-deficient mouse model of ALL. Together, these observations indicate that PAX5 deficiency reversibly drives leukemogenesis. PAX5 and its two most closely related paralogs, PAX2 and PAX8, which are not mutated in ALL, exhibit overlapping expression and function redundantly during embryonic development. However, PAX5 alone is expressed in lymphocytes, while PAX2 and PAX8 are predominantly specific to kidney and thyroid, respectively. We show that forced expression of PAX2 or PAX8 complements PAX5 loss-of-function mutation in ALL cells as determined by modulation of PAX5 target genes, restoration of immunophenotypic and morphological differentiation, and, ultimately, reduction of replicative potential. Activation of PAX5 paralogs, PAX2 or PAX8, ordinarily silenced in lymphocytes, may therefore represent a novel approach for treating PAX5-deficient ALL. In pursuit of this strategy, we took advantage of the fact that, in kidney, PAX2 is upregulated by extracellular hyperosmolarity. We found that hyperosmolarity, at potentially clinically achievable levels, transcriptionally activates endogenous PAX2 in ALL cells via a mechanism dependent on NFAT5, a transcription factor coordinating response to hyperosmolarity. We also found that hyperosmolarity upregulates residual wild type PAX5 expression in ALL cells and modulates gene expression, including in PAX5-mutant primary ALL cells. These findings specifically demonstrate that osmosensing pathways may represent a new therapeutic target for ALL and more broadly point toward the possibility of using gene paralogs to rescue mutations driving cancer and other diseases.

The Anaplastic Lymphoma Kinase (Alk) receptor tyrosine kinase (RTK) plays a critical role in the specification of founder cells (FCs) in the Drosophila visceral mesoderm (VM) during embryogenesis. Reporter gene and CRISPR/Cas9 deletion analysis reveals enhancer regions in and upstream of the Alk locus that influence tissue-specific expression in the amnioserosa (AS), the VM and the epidermis. By performing high throughput yeast one-hybrid screens (Y1H) with a library of Drosophila transcription factors (TFs) we identify Odd-paired (Opa), the Drosophila homologue of the vertebrate Zic family of TFs, as a novel regulator of embryonic Alk expression. Further characterization identifies evolutionarily conserved Opa-binding cis-regulatory motifs in one of the Alk associated enhancer elements. Employing Alk reporter lines as well as CRISPR/Cas9-mediated removal of regulatory elements in the Alk locus, we show modulation of Alk expression by Opa in the embryonic AS, epidermis and VM. In addition, we identify enhancer elements that integrate input from additional TFs, such as Binou (Bin) and Bagpipe (Bap), to regulate VM expression of Alk in a combinatorial manner. Taken together, our data show that the Opa zinc finger TF is a novel regulator of embryonic Alk expression.

Subtelomeres are duplication-rich, structurally variable regions of the human genome situated just proximal of telomeres. We report here that the most terminally located human subtelomeric genes encode a previously unrecognized third subclass of the Wiskott-Aldrich Syndrome Protein family, whose known members reorganize the actin cytoskeleton in response to extracellular stimuli. This new subclass, which we call WASH, is evolutionarily conserved in species as diverged as Entamoeba. We demonstrate that WASH is essential in Drosophila. WASH is widely expressed in human tissues, and human WASH protein colocalizes with actin in filopodia and lamellipodia. The VCA domain of human WASH promotes actin polymerization by the Arp2/3 complex in vitro. WASH duplicated to multiple chromosomal ends during primate evolution, with highest copy number reached in humans, whose WASH repertoires vary. Thus, human subtelomeres are not genetic junkyards, and WASH's location in these dynamic regions could have advantageous as well as pathologic consequences.

Sequencing family DNA samples provides an attractive alternative to population based designs to identify rare variants associated with human disease due to the enrichment of causal variants in pedigrees. Previous studies showed that genotype calling accuracy can be improved by modeling family relatedness compared to standard calling algorithms. Current family-based variant calling methods use sequencing data on single variants and ignore the identity-by-descent (IBD) sharing along the genome. In this study we describe a new computational framework to accurately estimate the IBD sharing from the sequencing data, and to utilize the inferred IBD among family members to jointly call genotypes in pedigrees. Through simulations and application to real data, we showed that IBD can be reliably estimated across the genome, even at very low coverage (e.g. 2X), and genotype accuracy can be dramatically improved. Moreover, the improvement is more pronounced for variants with low frequencies, especially at low to intermediate coverage (e.g. 10X to 20X), making our approach effective in studying rare variants in cost-effective whole genome sequencing in pedigrees. We hope that our tool is useful to the research community for identifying rare variants for human disease through family-based sequencing.

Maintenance of normal lipid homeostasis is crucial to heart function. On the other hand, the heart is now recognized to serve an important role in regulating systemic lipid metabolism; however, the molecular basis remains unclear. In this study, we identify the Drosophila Snail family of transcription factors (herein termed Sna TFs) as new mediators of the heart control of systemic lipid metabolism. Overexpression of Sna TF genes specifically in the heart promotes whole-body leanness whereas their knockdown in the heart promotes obesity. In addition, flies that are heterozygous for a snail deficiency chromosome also exhibit systemic obesity, and that cardiac-specific overexpression of Sna substantially reverses systemic obesity in these flies. We further show that genetically manipulating Sna TF levels in the fat body and intestine do not affect systemic lipid levels. Mechanistically, we find that flies bearing the overexpression or inhibition of Sna TFs in the postnatal heart only exhibit systemic lipid metabolic defects but not heart abnormalities. Cardiac-specific alterations of Sna TF levels also do not perturb cardiac morphology, viability, lipid metabolism or fly food intake. On the other hand, cardiac-specific manipulations of Sna TF levels alter lipogenesis and lipolysis gene expression, mitochondrial biogenesis and respiration, and lipid storage droplet 1 and 2 (Lsd-1 and Lsd-2) levels in the fat body. Together, our results reveal a novel and specific role of Sna TFs in the heart on systemic lipid homeostasis maintenance that is independent of cardiac development and function and involves the governance of triglyceride synthesis and breakdown, energy utilization, and lipid droplet dynamics in the fat body.

Chromatin regulation of eukaryotic genomes depends on the formation of nucleosome complexes between histone proteins and DNA. Histone variants, which are diversified by sequence or expression pattern, can profoundly alter chromatin properties. While variants in histone H2A and H3 families are well characterized, the extent of diversification of histone H2B proteins is less understood. Here, we report a systematic analysis of the histone H2B family in plants, which have undergone substantial divergence during the evolution of each major group in the plant kingdom. By characterising Arabidopsis H2Bs, we substantiate this diversification and reveal potential functional specialization that parallels the phylogenetic structure of emergent clades in eudicots. In addition, we identify a new class of highly divergent H2B variants, H2B.S, that specifically accumulate during chromatin compaction of dry seed embryos in multiple species of flowering plants. Our findings thus identify unsuspected diverse properties among histone H2B proteins in plants that has manifested into potentially novel groups of histone variants.

Regulation of the plant immune system is important for controlling the specificity and amplitude of responses to pathogens and in preventing growth-inhibiting autoimmunity that leads to reductions in plant fitness. In previous work, we reported that SRFR1, a negative regulator of effector-triggered immunity, interacts with SNC1 and EDS1. When SRFR1 is non-functional in the Arabidopsis accession Col-0, SNC1 levels increase, causing a cascade of events that lead to autoimmunity phenotypes. Previous work showed that some members of the transcriptional co-repressor family TOPLESS interact with SNC1 to repress negative regulators of immunity. Therefore, to explore potential connections between SRFR1 and TOPLESS family members, we took a genetic approach that examined the effect of each TOPLESS member in the srfr1 mutant background. The data indicated that an additive genetic interaction exists between SRFR1 and two members of the TOPLESS family, TPR2 and TPR3, as demonstrated by increased stunting and elevated PR2 expression in srfr1 tpr2 and srfr1 tpr2 tpr3 mutants. Furthermore, the tpr2 mutation intensifies autoimmunity in the auto-active snc1-1 mutant, indicating a novel role of these TOPLESS family members in negatively regulating SNC1-dependent phenotypes. This negative regulation can also be reversed by overexpressing TPR2 in the srfr1 tpr2 background. Similar to TPR1 that positively regulates snc1-1 phenotypes by interacting with SNC1, we show here that TPR2 directly binds the N-terminal domain of SNC1. In addition, TPR2 interacts with TPR1 in vivo, suggesting that the opposite functions of TPR2 and TPR1 are based on titration of SNC1-TPR1 complexes by TPR2 or altered functions of a SNC1-TPR1-TPR2 complex. Thus, this work uncovers diverse functions of individual members of the TOPLESS family in Arabidopsis and provides evidence for the additive effect of transcriptional and post-transcriptional regulation of SNC1.

The clinical utility of family history and genetic tests is generally well understood for simple Mendelian disorders and rare subforms of complex diseases that are directly attributable to highly penetrant genetic variants. However, little is presently known regarding the performance of these methods in situations where disease susceptibility depends on the cumulative contribution of multiple genetic factors of moderate or low penetrance. Using quantitative genetic theory, we develop a model for studying the predictive ability of family history and single nucleotide polymorphism (SNP)-based methods for assessing risk of polygenic disorders. We show that family history is most useful for highly common, heritable conditions (e.g., coronary artery disease), where it explains roughly 20%-30% of disease heritability, on par with the most successful SNP models based on associations discovered to date. In contrast, we find that for diseases of moderate or low frequency (e.g., Crohn disease) family history accounts for less than 4% of disease heritability, substantially lagging behind SNPs in almost all cases. These results indicate that, for a broad range of diseases, already identified SNP associations may be better predictors of risk than their family history-based counterparts, despite the large fraction of missing heritability that remains to be explained. Our model illustrates the difficulty of using either family history or SNPs for standalone disease prediction. On the other hand, we show that, unlike family history, SNP-based tests can reveal extreme likelihood ratios for a relatively large percentage of individuals, thus providing potentially valuable adjunctive evidence in a differential diagnosis.

The Drosophila Alp/Enigma family protein Zasp52 localizes to myotendinous junctions and Z-discs. It is required for terminal muscle differentiation and muscle attachment. Its vertebrate ortholog ZASP/Cypher also localizes to Z-discs, interacts with α-actinin through its PDZ domain, and is involved in Z-disc maintenance. Human mutations in ZASP cause myopathies and cardiomyopathies. Here we show that Drosophila Zasp52 is one of the earliest markers of Z-disc assembly, and we use a Zasp52-GFP fusion to document myofibril assembly by live imaging. We demonstrate that Zasp52 is required for adult Z-disc stability and pupal myofibril assembly. In addition, we show that two closely related proteins, Zasp66 and the newly identified Zasp67, are also required for adult Z-disc stability and are participating with Zasp52 in Z-disc assembly resulting in more severe, synergistic myofibril defects in double mutants. Zasp52 and Zasp66 directly bind to α-actinin, and they can also form a ternary complex. Our results indicate that Alp/Enigma family members cooperate in Z-disc assembly and myofibril formation; and we propose, based on sequence analysis, a novel class of PDZ domain likely involved in α-actinin binding.

Gene duplication facilitates functional diversification and provides greater phenotypic flexibility to an organism. Expanded gene families arise through repeated gene duplication but the extent of functional divergence that accompanies each paralogous gene is generally unexplored because of the difficulty in isolating the effects of single family members. The telomere-associated (TLO) gene family is a remarkable example of gene family expansion, with 14 members in the more pathogenic Candida albicans relative to two TLO genes in the closely-related species C. dubliniensis. TLO genes encode interchangeable Med2 subunits of the major transcriptional regulatory complex Mediator. To identify biological functions associated with each C. albicans TLO, expression of individual family members was regulated using a Tet-ON system and the strains were assessed across a range of phenotypes involved in growth and virulence traits. All TLOs affected multiple phenotypes and a single phenotype was often affected by multiple TLOs, including simple phenotypes such as cell aggregation and complex phenotypes such as virulence in a Galleria mellonella model of infection. No phenotype was regulated by all TLOs, suggesting neofunctionalization or subfunctionalization of ancestral properties among different family members. Importantly, regulation of three phenotypes could be mapped to individual polymorphic sites among the TLO genes, including an indel correlated with two phenotypes, growth in sucrose and macrophage killing. Different selective pressures have operated on the TLO sequence, with the 5' conserved Med2 domain experiencing purifying selection and the gene/clade-specific 3' end undergoing extensive positive selection that may contribute to the impact of individual TLOs on phenotypic variability. Therefore, expansion of the TLO gene family has conferred unique regulatory properties to each paralog such that it influences a range of phenotypes. We posit that the genetic diversity associated with this expansion contributed to C. albicans success as a commensal and opportunistic pathogen.

H3K9 methylation (H3K9me) is a conserved marker of heterochromatin, a transcriptionally silent chromatin structure. Knowledge of the mechanisms for regulating heterochromatin distribution is limited. The fission yeast JmjC domain-containing protein Epe1 localizes to heterochromatin mainly through its interaction with Swi6, a homologue of heterochromatin protein 1 (HP1), and directs JmjC-mediated H3K9me demethylation in vivo. Here, we found that loss of epe1 (epe1Δ) induced a red-white variegated phenotype in a red-pigment accumulation background that generated uniform red colonies. Analysis of isolated red and white colonies revealed that silencing of genes involved in pigment accumulation by stochastic ectopic heterochromatin formation led to white colony formation. In addition, genome-wide analysis of red- and white-isolated clones revealed that epe1Δ resulted in a heterogeneous heterochromatin distribution among clones. We found that Epe1 had an N-terminal domain distinct from its JmjC domain, which activated transcription in both fission and budding yeasts. The N-terminal transcriptional activation (NTA) domain was involved in suppression of ectopic heterochromatin-mediated red-white variegation. We introduced a single copy of Epe1 into epe1Δ clones harboring ectopic heterochromatin, and found that Epe1 could reduce H3K9me from ectopic heterochromatin but some of the heterochromatin persisted. This persistence was due to a latent H3K9me source embedded in ectopic heterochromatin. Epe1H297A, a canonical JmjC mutant, suppressed red-white variegation, but entirely failed to remove already-established ectopic heterochromatin, suggesting that Epe1 prevented stochastic de novo deposition of ectopic H3K9me in an NTA-dependent but JmjC-independent manner, while its JmjC domain mediated removal of H3K9me from established ectopic heterochromatin. Our results suggest that Epe1 not only limits the distribution of heterochromatin but also controls the balance between suppression and retention of heterochromatin-mediated epigenetic diversification.

Our recent study on the functional analysis of the Knickkopf protein from T. castaneum (TcKnk), indicated a novel role for this protein in protection of chitin from degradation by chitinases. Knk is also required for the laminar organization of chitin in the procuticle. During a bioinformatics search using this protein sequence as the query, we discovered the existence of a small family of three Knk-like genes (including the prototypical TcKnk) in the T. castaneum genome as well as in all insects with completed genome assemblies. The two additional Knk-like genes have been named TcKnk2 and TcKnk3. Further complexity arises as a result of alternative splicing and alternative polyadenylation of transcripts of TcKnk3, leading to the production of three transcripts (and by inference, three proteins) from this gene. These transcripts are named TcKnk3-Full Length (TcKnk3-FL), TcKnk3-5' and TcKnk3-3'. All three Knk-family genes appear to have essential and non-redundant functions. RNAi for TcKnk led to developmental arrest at every molt, while down-regulation of either TcKnk2 or one of the three TcKnk3 transcripts (TcKnk3-3') resulted in specific molting arrest only at the pharate adult stage. All three Knk genes appear to influence the total chitin content at the pharate adult stage, but to variable extents. While TcKnk contributes mostly to the stability and laminar organization of chitin in the elytral and body wall procuticles, proteins encoded by TcKnk2 and TcKnk3-3' transcripts appear to be required for the integrity of the body wall denticles and tracheal taenidia, but not the elytral and body wall procuticles. Thus, the three members of the Knk-family of proteins perform different essential functions in cuticle formation at different developmental stages and in different parts of the insect anatomy.

RNA interference (RNAi) utilizes small interfering RNAs (siRNAs) to direct silencing of specific genes through transcriptional and post-transcriptional mechanisms. The siRNA guides can originate from exogenous (exo-RNAi) or natural endogenous (endo-RNAi) sources of double-stranded RNA (dsRNA). In Caenorhabditis elegans, inactivation of genes that function in the endo-RNAi pathway can result in enhanced silencing of genes targeted by siRNAs from exogenous sources, indicating cross-regulation between the pathways. Here we show that members of another small RNA pathway, the mir-35-41 cluster of microRNAs (miRNAs) can regulate RNAi. In worms lacking miR-35-41, there is reduced expression of lin-35/Rb, the C. elegans homolog of the tumor suppressor Retinoblastoma gene, previously shown to regulate RNAi responsiveness. Genome-wide microarray analyses show that targets of endo-siRNAs are up-regulated in mir-35-41 mutants, a phenotype also displayed by lin-35/Rb mutants. Furthermore, overexpression of lin-35/Rb specifically rescues the RNAi hypersensitivity of mir-35-41 mutants. Although the mir-35-41 miRNAs appear to be exclusively expressed in germline and embryos, their effect on RNAi sensitivity is transmitted to multiple tissues and stages of development. Additionally, we demonstrate that maternal contribution of miR-35-41 or lin-35/Rb is sufficient to reduce RNAi effectiveness in progeny worms. Our results reveal that miRNAs can broadly regulate other small RNA pathways and, thus, have far reaching effects on gene expression beyond directly targeting specific mRNAs.

The Ewing sarcoma family of tumors (EFT) is a group of highly malignant small round blue cell tumors occurring in children and young adults. We report here the largest genomic survey to date of 101 EFT (65 tumors and 36 cell lines). Using a combination of whole genome sequencing and targeted sequencing approaches, we discover that EFT has a very low mutational burden (0.15 mutations/Mb) but frequent deleterious mutations in the cohesin complex subunit STAG2 (21.5% tumors, 44.4% cell lines), homozygous deletion of CDKN2A (13.8% and 50%) and mutations of TP53 (6.2% and 71.9%). We additionally note an increased prevalence of the BRCA2 K3326X polymorphism in EFT patient samples (7.3%) compared to population data (OR 7.1, p = 0.006). Using whole transcriptome sequencing, we find that 11% of tumors pathologically diagnosed as EFT lack a typical EWSR1 fusion oncogene and that these tumors do not have a characteristic Ewing sarcoma gene expression signature. We identify samples harboring novel fusion genes including FUS-NCATc2 and CIC-FOXO4 that may represent distinct small round blue cell tumor variants. In an independent EFT tissue microarray cohort, we show that STAG2 loss as detected by immunohistochemistry may be associated with more advanced disease (p = 0.15) and a modest decrease in overall survival (p = 0.10). These results significantly advance our understanding of the genomic and molecular underpinnings of Ewing sarcoma and provide a foundation towards further efforts to improve diagnosis, prognosis, and precision therapeutics testing.