Complementary CRISPR screen highlights the contrasting role of membrane-bound and soluble ICAM-1 in regulating antigen-specific tumor cell killing by cytotoxic T cells.

Cytotoxic CD8 +T lymphocytes (CTLs) are key players of adaptive anti-tumor immunity based on their ability to specifically recognize and destroy tumor cells. Many cancer immunotherapies rely on unleashing CTL function. However, tumors can evade killing through strategies which are not yet fully elucidated. To provide deeper insight into tumor evasion mechanisms in an antigen-dependent manner, we established a human co-culture system composed of tumor and primary immune cells. Using this system, we systematically investigated intrinsic regulators of tumor resistance by conducting a complementary CRISPR screen approach. By harnessing CRISPR activation (CRISPRa) and CRISPR knockout (KO) technology in parallel, we investigated gene gain-of-function as well as loss-of-function across genes with annotated function in a colon carcinoma cell line. CRISPRa and CRISPR KO screens uncovered 187 and 704 hits, respectively, with 60 gene hits overlapping between both. These data confirmed the role of interferon-γ (IFN-γ), tumor necrosis factor α (TNF-α) and autophagy pathways and uncovered novel genes implicated in tumor resistance to killing. Notably, we discovered that ILKAP encoding the integrin-linked kinase-associated serine/threonine phosphatase 2 C, a gene previously unknown to play a role in antigen specific CTL-mediated killing, mediate tumor resistance independently from regulating antigen presentation, IFN-γ or TNF-α responsiveness. Moreover, our work describes the contrasting role of soluble and membrane-bound ICAM-1 in regulating tumor cell killing. The deficiency of membrane-bound ICAM-1 (mICAM-1) or the overexpression of soluble ICAM-1 (sICAM-1) induced resistance to CTL killing, whereas PD-L1 overexpression had no impact. These results highlight the essential role of ICAM-1 at the immunological synapse between tumor and CTL and the antagonist function of sICAM-1.

Senescent cells suppress macrophage-mediated corpse removal via upregulation of the CD47-QPCT/L axis.

Progressive accrual of senescent cells in aging and chronic diseases is associated with detrimental effects in tissue homeostasis. We found that senescent fibroblasts and epithelia were not only refractory to macrophage-mediated engulfment and removal, but they also paralyzed the ability of macrophages to remove bystander apoptotic corpses. Senescent cell-mediated efferocytosis suppression (SCES) was independent of the senescence-associated secretory phenotype (SASP) but instead required direct contact between macrophages and senescent cells. SCES involved augmented senescent cell expression of CD47 coinciding with increased CD47-modifying enzymes QPCT/L. SCES was reversible by interfering with the SIRPα-CD47-SHP-1 axis or QPCT/L activity. While CD47 expression increased in human and mouse senescent cells in vitro and in vivo, another ITIM-containing protein, CD24, contributed to SCES specifically in human epithelial senescent cells where it compensated for genetic deficiency in CD47. Thus, CD47 and CD24 link the pathogenic effects of senescent cells to homeostatic macrophage functions, such as efferocytosis, which we hypothesize must occur efficiently to maintain tissue homeostasis.

Antifibrotic Drug Nintedanib Inhibits CSF1R to Promote IL-4-associated Tissue Repair Macrophages.

Profibrotic and prohomeostatic macrophage phenotypes remain ill-defined, both in vivo and in vitro, impeding the successful development of drugs that reprogram macrophages as an attractive therapeutic approach to manage fibrotic disease. The goal of this study was to reveal profibrotic and prohomeostatic macrophage phenotypes that could guide the design of new therapeutic approaches targeting macrophages to treat fibrotic disease. This study used nintedanib, a broad kinase inhibitor approved for idiopathic pulmonary fibrosis, to dissect lung macrophage phenotypes during fibrosis-linked inflammation by combining in vivo and in vitro bulk and single-cell RNA-sequencing approaches. In the bleomycin model, nintedanib drove the expression of IL-4/IL-13-associated genes important for tissue regeneration and repair at early and late time points in lung macrophages. These findings were replicated in vitro in mouse primary bone marrow-derived macrophages exposed to IL-4/IL-13 and nintedanib. In addition, nintedanib promoted the expression of IL-4/IL-13 pathway genes in human macrophages in vitro. The molecular mechanism was connected to inhibition of the colony stimulating factor 1 (CSF1) receptor in both human and mouse macrophages. Moreover, nintedanib counterbalanced the effects of TNF on IL-4/IL-13 in macrophages to promote expression of IL-4/IL-13-regulated tissue repair genes in fibrotic contexts in vivo and in vitro. This study demonstrates that one of nintedanib's antifibrotic mechanisms is to increase IL-4 signaling in macrophages through inhibition of the CSF1 receptor, resulting in the promotion of tissue repair phenotypes.

Automation enables high-throughput and reproducible single-cell transcriptomics library preparation.

Next-generation sequencing (NGS) has revolutionized genomics, decreasing sequencing costs and allowing researchers to draw correlations between diseases and DNA or RNA changes. Technical advances have enabled the analysis of RNA expression changes between single cells within a heterogeneous population, known as single-cell RNA-seq (scRNA-seq). Despite resolving transcriptomes of cellular subpopulations, scRNA-seq has not replaced RNA-seq, due to higher costs and longer hands-on time. Here, we developed an automated workflow to increase throughput (up to 48 reactions) and to reduce by 75% the hands-on time of scRNA-seq library preparation, using the 10X Genomics Single Cell 3' kit. After gel bead-in-emulsion (GEM) generation on the 10X Genomics Chromium Controller, cDNA amplification was performed, and the product was normalized and subjected to either the manual, standard library preparation method or a fully automated, walk-away method using a Biomek i7 Hybrid liquid handler. Control metrics showed that both quantity and quality of the single-cell gene expression libraries generated were equivalent in size and yield. Key scRNA-seq downstream quality metrics, such as unique molecular identifiers count, mitochondrial RNA content, and cell and gene counts, further showed high correlations between automated and manual workflows. Using the UMAP dimensionality reduction technique to visualize all cells, we were able to further correlate the results observed between the manual and automated methods (R=0.971). The method developed here allows for the fast, error-free, and reproducible multiplex generation of high-quality single-cell gene expression libraries.

Universal Gene Correction Approaches for ß-hemoglobinopathies Using CRISPR-Cas9 and Adeno-Associated Virus Serotype 6 Donor Templates.

Mutations in the human ß-globin gene are the cause of ß-hemoglobinopathies, one of the most common inherited single-gene blood disorders in the world. Novel therapeutic approaches are based on lentiviral vectors (LVs) or CRISPR-Cas9-mediated gene disruption to express adult hemoglobin (HbA), or to reactivate the completely functional fetal hemoglobin, respectively. Nonetheless, LVs present a risk of insertional mutagenesis, while gene-disrupting transcription factors (BCL11A, KLF1) involved in the fetal-to-adult hemoglobin switch might generate dysregulation of other cellular processes. Therefore, universal gene addition/correction approaches combining CRISPR-Cas9 and homology directed repair (HDR) by delivering a DNA repair template through adeno-associated virus could mitigate the limitations of both lentiviral gene transfer and gene disruption strategies, ensuring targeted integration and controlled transgene expression. In this study, we attained high rates of gene addition (up to 12%) and gene correction (up to 38%) in hematopoietic stem and progenitor cells from healthy donors without any cell sorting/enrichment or the application of HDR enhancers. Furthermore, these approaches were tested in heterozygous (ß0/ß+) and homozygous (ß0/ß0, ß+/ß+) ß-thalassemia patients, achieving a significant increase in HbA and demonstrating the universal therapeutic potential of this study for the treatment of ß-hemoglobinopathies.

Genetic evolution of in situ follicular neoplasia to aggressive B-cell lymphoma of germinal center subtype.

In situ follicular neoplasia (ISFN) is the earliest morphologically identifiable precursor of follicular lymphoma (FL). Although it is genetically less complex than FL and has low risk for progression, ISFN already harbors secondary genetic alterations, in addition to the defining t(14;18)(q32;q21) translocation. FL, in turn, frequently progresses to diffuse large B-cell lymphoma (DLBCL) or high-grade B-cell lymphoma (HGBL). By BCL2 staining of available reactive lymphoid tissue obtained at any time point in patients with aggressive B-cell lymphoma (BCL), we identified ten paired cases of ISFN and DLBCL/HGBL, including six de novo tumors and four tumors transformed from FL as an intermediate step, and investigated their clonal evolution using microdissection and next-generation sequencing. A clonal relationship between ISFN and aggressive BCL was established by immunoglobulin and/or BCL2 rearrangements and/or the demonstration of shared somatic mutations for all ten cases. Targeted sequencing revealed CREBBP, KMT2D, EZH2, TNFRSF14 and BCL2 as the genes most frequently mutated already in ISFN. Based on the distribution of private and shared mutations, two patterns of clonal evolution were evident. In most cases, the aggressive lymphoma, ISFN and, when present, FL revealed divergent evolution from a common progenitor, whereas linear evolution with sequential accumulation of mutations was less frequent. In conclusion, we demonstrate for the first time that t(14;18)+ aggressive BCL can arise from ISFN without clinically evident FL as an intermediate step and that during this progression, branched evolution is common.

Comparative targeting analysis of KLF1, BCL11A, and HBG1/2 in CD34+ HSPCs by CRISPR/Cas9 for the induction of fetal hemoglobin.

ß-hemoglobinopathies are caused by abnormal or absent production of hemoglobin in the blood due to mutations in the ß-globin gene (HBB). Imbalanced expression of adult hemoglobin (HbA) induces strong anemia in patients suffering from the disease. However, individuals with natural-occurring mutations in the HBB cluster or related genes, compensate this disparity through γ-globin expression and subsequent fetal hemoglobin (HbF) production. Several preclinical and clinical studies have been performed in order to induce HbF by knocking-down genes involved in HbF repression (KLF1 and BCL11A) or disrupting the binding sites of several transcription factors in the γ-globin gene (HBG1/2). In this study, we thoroughly compared the different CRISPR/Cas9 gene-disruption strategies by gene editing analysis and assessed their safety profile by RNA-seq and GUIDE-seq. All approaches reached therapeutic levels of HbF after gene editing and showed similar gene expression to the control sample, while no significant off-targets were detected by GUIDE-seq. Likewise, all three gene editing platforms were established in the GMP-grade CliniMACS Prodigy, achieving similar outcome to preclinical devices. Based on this gene editing comparative analysis, we concluded that BCL11A is the most clinically relevant approach while HBG1/2 could represent a promising alternative for the treatment of ß-hemoglobinopathies.

Whole genome characterization of autochthonous Bos taurus brachyceros and introduced Bos indicus indicus cattle breeds in Cameroon regarding their adaptive phenotypic traits and pathogen resistance.

BACKGROUND: African indigenous taurine cattle display unique adaptive traits shaped by husbandry management, regional climate and exposure to endemic pathogens. They are less productive with respect to milk and meat production which has been associated with amongst others, small size, traditional beliefs, husbandry practices, limited feed resources, disease burden and lack of sustained breeding for trait improvement. This resulted in the severe dwindling of their population size rendering them vulnerable to extinction. The Namchi taurine cattle breed is referred to as [Namchi (Doayo)] and shows resistance traits against trypanosome infection and exposure to tick infestation. Nonetheless, the historically later introduced Zebu cattle are the main cattle breeds in Africa today, even though they suffer more from locally prevailing pathogens. By using a whole genome sequencing approach, we sequenced with high depth for the first time the genomes of five cattle breeds from Cameroon in order to provide a valuable genetic resource for future African cattle breeding: the Namchi, an endangered trypano-tolerant taurine breed, the Kapsiki, an indigenous trypano-susceptible taurine breed, and three Zebu (Bos indicus indicus) breeds: Ngaoundere Gudali, White Fulani and Red Fulani. RESULTS: Approximately 167 Gigabases of raw sequencing data were generated for each breed and mapped to the cattle reference genomes ARS-UCD1.2 and UMD3.1.The coverage was 103 to 140-fold when aligning the reads to ARS-UCD1.2 with an average mapping rate of ~ 99%, and 22 to 30-fold when aligning the reads to UMD3.1 with an average mapping rate of ~ 64%. The single nucleotide polymorphisms (SNPs) obtained from analysis using the genome ARS-UCD1.2 were compared with reference genomes of European Bos taurus Holstein, the Asian Bos indicus Brahman, and the African trypanotolerant N'Dama breeds. A total of ~ 100 million (M) SNPs were identified and 7.7 M of those were breed-specific. An approximately 11.1 M constituted of small insertions and deletions. By using only breed-specific non-synonymous variants we identified genes as genetic signatures and associated Gene Ontology (GO) terms that could explain certain cattle-breed specific phenotypes such as increased tolerance against trypanosome parasites in the Namchi breed and heat tolerance in the Kapsiki breed. Phylogenetic analysis grouped, except for Namchi, the Bos taurus breeds Kapsiki, N'Dama and Holstein together while the B. indicus breeds White and Red Fulani, Gudali and Brahman clustered separately. The deviating result for Namchi indicates a hybrid status of the selected animal with a recent introgression of Zebu genes into its genome. CONCLUSIONS: The findings provide the first comprehensive set of genome-wide variant data of the most important Cameroonian cattle breeds. The genomic data shall constitute a foundation for breed amelioration whilst exploiting the heritable traits and support conservation efforts for the endangered local cattle breeds.

Gene correction of HBB mutations in CD34+ hematopoietic stem cells using Cas9 mRNA and ssODN donors.

BACKGROUND: ß-Thalassemia is an inherited hematological disorder caused by mutations in the human hemoglobin beta (HBB) gene that reduce or abrogate ß-globin expression. Although lentiviral-mediated expression of ß-globin and autologous transplantation is a promising therapeutic approach, the risk of insertional mutagenesis or low transgene expression is apparent. However, targeted gene correction of HBB mutations with programmable nucleases such as CRISPR/Cas9, TALENs, and ZFNs with non-viral repair templates ensures a higher safety profile and endogenous expression control. METHODS: We have compared three different gene-editing tools (CRISPR/Cas9, TALENs, and ZFNs) for their targeting efficiency of the HBB gene locus. As a proof of concept, we studied the personalized gene-correction therapy for a common ß-thalassemia splicing variant HBBIVS1-110 using Cas9 mRNA and several optimally designed single-stranded oligonucleotide (ssODN) donors in K562 and CD34+ hematopoietic stem cells (HSCs). RESULTS: Our results exhibited that indel frequency of CRISPR/Cas9 was superior to TALENs and ZFNs (P < 0.0001). Our designed sgRNA targeting the site of HBBIVS1-110 mutation showed indels in both K562 cells (up to 77%) and CD34+ hematopoietic stem cells-HSCs (up to 87%). The absolute quantification by next-generation sequencing showed that up to 8% site-specific insertion of the NheI tag was achieved using Cas9 mRNA and a chemically modified ssODN in CD34+ HSCs. CONCLUSION: Our approach provides guidance on non-viral gene correction in CD34+ HSCs using Cas9 mRNA and chemically modified ssODN. However, further optimization is needed to increase the homology directed repair (HDR) to attain a real clinical benefit for ß-thalassemia.

Two independent sulfation processes regulate mouth-form plasticity in the nematode Pristionchus pacificus.

Sulfation of biomolecules, like phosphorylation, is one of the most fundamental and ubiquitous biochemical modifications with important functions during detoxification. This process is reversible, involving two enzyme classes: a sulfotransferase, which adds a sulfo group to a substrate; and a sulfatase that removes the sulfo group. However, unlike phosphorylation, the role of sulfation in organismal development is poorly understood. In this study, we find that two independent sulfation events regulate the development of mouth morphology in the nematode Pristionchus pacificus. This nematode has the ability to form two alternative mouth morphologies depending on environmental cues, an example of phenotypic plasticity. We found that, in addition to a previously described sulfatase, a sulfotransferase is involved in regulating the mouth-form dimorphism in P. pacificus However, it is unlikely that both of these sulfation-associated enzymes act upon the same substrates, as they are expressed in different cell types. Furthermore, animals mutant in genes encoding both enzymes show condition-dependent epistatic interactions. Thus, our study highlights the role of sulfation-associated enzymes in phenotypic plasticity of mouth structures in Pristionchus.

Identification and Functional Characterization of Phosphorylation Sites of the Human Papillomavirus 31 E8

The papillomavirus E2 protein regulates transcription, replication, and nuclear retention of viral genomes. Phosphorylation of E2 in the hinge region has been suggested to modulate protein stability, DNA-binding activity, and chromosomal attachment. The papillomavirus E8^E2 protein shares the hinge domain with E2 and acts as a repressor of viral replication. Mass spectrometry analyses of human papillomavirus 31 (HPV31) E8^E2 and E2 proteins identify phosphorylated S78, S81, and S100 in E8^E2 and S266 and S269 in E2 in their hinge regions. Phos-tag analyses of wild-type and mutant proteins indicate that S78 is a major phosphorylation site in E8^E2, but the corresponding S266 in E2 is not. Phosphorylation at S78 regulates E8^E2's repression activity of reporter constructs, whereas the corresponding E2 mutants do not display a phenotype. Phosphorylation at S78 does not alter E8^E2's protein stability, nuclear localization, or binding to DNA or to cellular NCoR/SMRT complexes. Surprisingly, in the context of HPV31 genomes, mutation of E8^E2 S78 does not modulate viral replication or transcription in undifferentiated or differentiated cells. However, comparative transcriptome analyses of differentiated HPV31 E8^E2 S78A and S78E cell lines reveal that the expression of a small number of cellular genes is changed. Validation experiments suggest that the transcription of the cellular LYPD2 gene is altered in a phospho-S78 E8^E2-dependent manner. In summary, our data suggest that phosphorylation of S78 in E8^E2 regulates its repression activity by a novel mechanism, and this seems to be important for the modulation of host cell gene expression but not viral replication.IMPORTANCE Posttranslational modification of viral proteins is a common feature to modulate their activities. Phosphorylation of serine residues S298 and S301 in the hinge region of the bovine papillomavirus type 1 E2 protein has been shown to restrict viral replication. The papillomavirus E8^E2 protein shares the hinge domain with E2 and acts as a repressor of viral replication. A large fraction of HPV31 E8^E2 is phosphorylated at S78 in the hinge region, and this is important for E8^E2's repression activity. Surprisingly, phosphorylation at S78 in E8^E2 has no impact on viral replication in tissue culture but rather seems to modulate the expression of a small number of cellular genes. This may indicate that phosphorylation of viral transcription factors serves to broaden their target gene specificity.

Transcriptomic profile of cystic fibrosis patients identifies type I interferon response and ribosomal stalk proteins as potential modifiers of disease severity.

Cystic Fibrosis (CF) is the most common monogenic disease among people of Western European descent and caused by mutations in the CFTR gene. However, the disease severity is immensely variable even among patients with similar CFTR mutations due to the possible effect of 'modifier genes'. To identify genetic modifiers, we applied RNA-seq based transcriptomic analyses in CF patients with a mild and severe lung phenotype. Global gene expression and enrichment analyses revealed that genes of the type I interferon response and ribosomal stalk proteins are potential modifiers of CF related lung dysfunction. The results provide a new set of CF modifier genes with possible implications as new therapeutic targets for the treatment of CF.

Duplications and Positive Selection Drive the Evolution of Parasitism-Associated Gene Families in the Nematode Strongyloides papillosus.

Gene duplication is a major mechanism playing a role in the evolution of phenotypic complexity and in the generation of novel traits. By comparing parasitic and nonparasitic nematodes, a recent study found that the evolution of parasitism in Strongyloididae is associated with a large expansion in the Astacin and CAP gene families.To gain novel insights into the developmental processes in the sheep parasite Strongyloides papillosus, we sequenced transcriptomes of different developmental stages and sexes. Overall, we found that the majority of genes are developmentally regulated and have one-to-one orthologs in the diverged S. ratti genome. Together with the finding of similar expression profiles between S. papillosus and S. ratti, these results indicate a strong evolutionary constraint acting against change at sequence and expression levels. However, the comparison between parasitic and free-living females demonstrates a quite divergent pattern that is mostly due to the previously mentioned expansion in the Astacin and CAP gene families. More detailed phylogenetic analysis of both gene families shows that most members date back to single expansion events early in the Strongyloides lineage and have undergone subfunctionalization resulting in clusters that are highly expressed either in infective larvae or in parasitic females. Finally, we found increased evidence for positive selection in both gene families relative to the genome-wide expectation.In summary, our study reveals first insights into the developmental transcriptomes of S. papillosus and provides a detailed analysis of sequence and expression evolution in parasitism-associated gene families.

Succession and dynamics of Pristionchus nematodes and their microbiome during decomposition of Oryctes borbonicus on La Réunion Island.

Insects and nematodes represent the most species-rich animal taxa and they occur together in a variety of associations. Necromenic nematodes of the genus Pristionchus are found on scarab beetles with more than 30 species known from worldwide samplings. However, little is known about the dynamics and succession of nematodes and bacteria during the decomposition of beetle carcasses. Here, we study nematode and bacterial succession of the decomposing rhinoceros beetle Oryctes borbonicus on La Réunion Island. We show that Pristionchus pacificus exits the arrested dauer stage seven days after the beetles´ deaths. Surprisingly, new dauers are seen after 11 days, suggesting that some worms return to the dauer stage after one reproductive cycle. We used high-throughput sequencing of the 16S rRNA genes of decaying beetles, beetle guts and nematodes to study bacterial communities in comparison to soil. We find that soil environments have the most diverse bacterial communities. The bacterial community of living and decaying beetles are more stable but one single bacterial family dominates the microbiome of decaying beetles. In contrast, the microbiome of nematodes is relatively similar even across different families. This study represents the first characterization of the dynamics of nematode-bacterial interactions during the decomposition of insects.

Life History Responses and Gene Expression Profiles of the Nematode Pristionchus pacificus Cultured on Cryptococcus Yeasts.

Nematodes, the earth's most abundant metazoa are found in all ecosystems. In order to survive in diverse environments, they have evolved distinct feeding strategies and they can use different food sources. While some nematodes are specialists, including parasites of plants and animals, others such as Pristionchus pacificus are omnivorous feeders, which can live on a diet of bacteria, protozoans, fungi or yeast. In the wild, P. pacificus is often found in a necromenic association with beetles and is known to be able to feed on a variety of microbes as well as on nematode prey. However, in laboratory studies Escherichia coli OP50 has been used as standard food source, similar to investigations in Caenorhabditis elegans and it is unclear to what extent this biases the obtained results and how relevant findings are in real nature. To gain first insight into the variation in traits induced by a non-bacterial food source, we study Pristionchus-fungi interactions under laboratory conditions. After screening different yeast strains, we were able to maintain P. pacificus for at least 50-60 generations on Cryptococcus albidus and Cryptococcus curvatus. We describe life history traits of P. pacificus on both yeast strains, including developmental timing, survival and brood size. Despite a slight developmental delay and problems to digest yeast cells, which are both reflected at a transcriptomic level, all analyses support the potential of Cryptococcus strains as food source for P. pacificus. In summary, our work establishes two Cryptococcus strains as alternative food source for P. pacificus and shows change in various developmental, physiological and morphological traits, including the transcriptomic profiles.

First insights into the nature and evolution of antisense transcription in nematodes.

BACKGROUND: The development of multicellular organisms is coordinated by various gene regulatory mechanisms that ensure correct spatio-temporal patterns of gene expression. Recently, the role of antisense transcription in gene regulation has moved into focus of research. To characterize genome-wide patterns of antisense transcription and to study their evolutionary conservation, we sequenced a strand-specific RNA-seq library of the nematode Pristionchus pacificus. RESULTS: We identified 1112 antisense configurations of which the largest group represents 465 antisense transcripts (ASTs) that are fully embedded in introns of their host genes. We find that most ASTs show homology to protein-coding genes and are overrepresented in proteomic data. Together with the finding, that expression levels of ASTs and host genes are uncorrelated, this indicates that most ASTs in P. pacificus do not represent non-coding RNAs and do not exhibit regulatory functions on their host genes. We studied the evolution of antisense gene pairs across 20 nematode genomes, showing that the majority of pairs is lineage-specific and even the highly conserved vps-4, ddx-27, and sel-2 loci show abundant structural changes including duplications, deletions, intron gains and loss of antisense transcription. In contrast, host genes in general, are remarkably conserved and encode exceptionally long introns leading to unusually large blocks of conserved synteny. CONCLUSIONS: Our study has shown that in P. pacificus antisense transcription as such does not define non-coding RNAs but is rather a feature of highly conserved genes with long introns. We hypothesize that the presence of regulatory elements imposes evolutionary constraint on the intron length, but simultaneously, their large size makes them a likely target for translocation of genomic elements including protein-coding genes that eventually end up as ASTs.

Draft Genome of the Scarab Beetle Oryctes borbonicus on La Réunion Island.

Beetles represent the largest insect order and they display extreme morphological, ecological and behavioral diversity, which makes them ideal models for evolutionary studies. Here, we present the draft genome of the scarab beetle Oryctes borbonicus, which has a more basal phylogenetic position than the two previously sequenced pest species Tribolium castaneum and Dendroctonus ponderosae providing the potential for sequence polarization. Oryctes borbonicus is endemic to La Réunion, an island located in the Indian Ocean, and is the host of the nematode Pristionchus pacificus, a well-established model organism for integrative evolutionary biology. At 518 Mb, the O. borbonicus genome is substantially larger and encodes more genes than T. castaneum and D. ponderosae We found that only 25% of the predicted genes of O. borbonicus are conserved as single copy genes across the nine investigated insect genomes, suggesting substantial gene turnover within insects. Even within beetles, up to 21% of genes are restricted to only one species, whereas most other genes have undergone lineage-specific duplications and losses. We illustrate lineage-specific duplications using detailed phylogenetic analysis of two gene families. This study serves as a reference point for insect/coleopteran genomics, although its original motivation was to find evidence for potential horizontal gene transfer (HGT) between O. borbonicus and P. pacificus The latter was previously shown to be the recipient of multiple horizontally transferred genes including some genes from insect donors. However, our study failed to provide any clear evidence for additional HGTs between the two species.

Ancient gene duplications have shaped developmental stage-specific expression in Pristionchus pacificus.

BACKGROUND: The development of multicellular organisms is accompanied by gene expression changes in differentiating cells. Profiling stage-specific expression during development may reveal important insights into gene sets that contributed to the morphological diversity across the animal kingdom. RESULTS: We sequenced RNA-seq libraries throughout a developmental timecourse of the nematode Pristionchus pacificus. The transcriptomes reflect early larval stages, adult worms including late larvae, and growth-arrested dauer larvae and allowed the identification of developmentally regulated gene clusters. Our data reveals similar trends as previous transcriptome profiling of dauer worms and represents the first expression data for early larvae in P. pacificus. Gene expression clusters characterizing early larval stages show most significant enrichments of chaperones, while collagens are most significantly enriched in transcriptomes of late larvae and adult worms. By combining expression data with phylogenetic analysis, we found that developmentally regulated genes are found in paralogous clusters that have arisen through lineage-specific duplications after the split from the Caenorhabditis elegans branch. CONCLUSIONS: We propose that gene duplications of developmentally regulated genes represent a plausible evolutionary mechanism to increase the dosage of stage-specific expression. Consequently, this may contribute to the substantial divergence in expression profiles that has been observed across larger evolutionary time scales.

Microevolution of Duplications and Deletions and Their Impact on Gene Expression in the Nematode Pristionchus pacificus.

The evolution of diversity across the animal kingdom has been accompanied by tremendous gene loss and gain. While comparative genomics has been fruitful to characterize differences in gene content across highly diverged species, little is known about the microevolution of structural variations that cause these differences in the first place. In order to investigate the genomic impact of structural variations, we made use of genomic and transcriptomic data from the nematode Pristionchus pacificus, which has been established as a satellite model to Caenorhabditis elegans for comparative biology. We exploit the fact that P. pacificus is a highly diverse species for which various genomic data including the draft genome of a sister species P. exspectatus is available. Based on resequencing coverage data for two natural isolates we identified large (> 2 kb) deletions and duplications relative to the reference strain. By restriction to completely syntenic regions between P. pacificus and P. exspectatus, we were able to polarize the comparison and to assess the impact of structural variations on expression levels. We found that while loss of genes correlates with lack of expression, duplication of genes has virtually no effect on gene expression. Further investigating expression of individual copies at sites that segregate between the duplicates, we found in the majority of cases only one of the copies to be expressed. Nevertheless, we still find that certain gene classes are strongly depleted in deletions as well as duplications, suggesting evolutionary constraint acting on synteny. In summary, our results are consistent with a model, where most structural variations are either deleterious or neutral and provide first insights into the microevolution of structural variations in the P. pacificus genome.

Deoxyribonucleic acid methylation profiling of single human blastocysts by methylated CpG-island amplification coupled with CpG-island microarray.

OBJECTIVE: To study whether methylated CpG-island (CGI) amplification coupled with microarray (MCAM) can be used to generate DNA (deoxyribonucleic acid) methylation profiles from single human blastocysts. DESIGN: A pilot microarray study with methylated CpG-island amplification applied to human blastocyst genomic DNA and hybridized on CpG-island microarrays. SETTING: University research laboratory. PATIENT(S): Five cryopreserved sibling 2-pronuclear zygotes that were surplus to requirements for clinical treatment by in vitro fertilization were donated with informed consent from a patient attending Bourn Hall Clinic, Cambridge, United Kingdom. INTERVENTION(S): None. MAIN OUTCOME MEASURE(S): Successful generation of genome-wide DNA methylation profiles at CpG islands from individual human blastocysts, with common genomic regions of DNA methylation identified between embryos. RESULT(S): Between 472 and 734 CpG islands were methylated in each blastocyst, with 121 CpG islands being commonly methylated in all 5 blastocysts. A further 159 CGIs were commonly methylated in 4 of the 5 tested blastocysts. Methylation was observed at a number of CGIs within imprinted-gene, differentially methylated regions (DMRs), including placental and preimplantation-specific DMRs. CONCLUSION(S): The MCAM method is capable of providing comprehensive DNA methylation data in individual human blastocysts.