Exome sequencing resolves apparent incidental findings and reveals further complexity of SH3TC2 variant alleles causing Charcot-Marie-Tooth neuropathy

BACKGROUND

The debate regarding the relative merits of whole genome sequencing (WGS) versus exome sequencing (ES) centers around comparative cost, average depth of coverage for each interrogated base, and their relative efficiency in the identification of medically actionable variants from the myriad of variants identified by each approach. Nevertheless, few genomes have been subjected to both WGS and ES, using multiple next generation sequencing platforms. In addition, no personal genome has been so extensively analyzed using DNA derived from peripheral blood as opposed to DNA from transformed cell lines that may either accumulate mutations during propagation or clonally expand mosaic variants during cell transformation and propagation.

METHODS

We investigated a genome that was studied previously by SOLiD chemistry using both ES and WGS, and now perform six independent ES assays (Illumina GAII (x2), Illumina HiSeq (x2), Life Technologies' Personal Genome Machine (PGM) and Proton), and one additional WGS (Illumina HiSeq).

RESULTS

We compared the variants identified by the different methods and provide insights into the differences among variants identified between ES runs in the same technology platform and among different sequencing technologies. We resolved the true genotypes of medically actionable variants identified in the proband through orthogonal experimental approaches. Furthermore, ES identified an additional SH3TC2 variant (p.M1?) that likely contributes to the phenotype in the proband.

CONCLUSIONS

ES identified additional medically actionable variant calls and helped resolve ambiguous single nucleotide variants (SNV) documenting the power of increased depth of coverage of the captured targeted regions. Comparative analyses of WGS and ES reveal that pseudogenes and segmental duplications may explain some instances of apparent disease mutations in unaffected individuals.

Rare complete knockouts in humans: population distribution and significant role in autism spectrum disorders

To characterize the role of rare complete human knockouts in autism spectrum disorders (ASDs), we identify genes with homozygous or compound heterozygous loss-of-function (LoF) variants (defined as nonsense and essential splice sites) from exome sequencing of 933 cases and 869 controls. We identify a 2-fold increase in complete knockouts of autosomal genes with low rates of LoF variation (≤ 5% frequency) in cases and estimate a 3% contribution to ASD risk by these events, confirming this observation in an independent set of 563 probands and 4,605 controls. Outside the pseudoautosomal regions on the X chromosome, we similarly observe a significant 1.5-fold increase in rare hemizygous knockouts in males, contributing to another 2% of ASDs in males. Taken together, these results provide compelling evidence that rare autosomal and X chromosome complete gene knockouts are important inherited risk factors for ASD.

Mind the gap: upgrading genomes with Pacific Biosciences RS long-read sequencing technology

Many genomes have been sequenced to high-quality draft status using Sanger capillary electrophoresis and/or newer short-read sequence data and whole genome assembly techniques. However, even the best draft genomes contain gaps and other imperfections due to limitations in the input data and the techniques used to build draft assemblies. Sequencing biases, repetitive genomic features, genomic polymorphism, and other complicating factors all come together to make some regions difficult or impossible to assemble. Traditionally, draft genomes were upgraded to “phase 3 finished” status using time-consuming and expensive Sanger-based manual finishing processes. For more facile assembly and automated finishing of draft genomes, we present here an automated approach to finishing using long-reads from the Pacific Biosciences RS (PacBio) platform. Our algorithm and associated software tool, PBJelly, (publicly available at https://sourceforge.net/projects/pb-jelly/) automates the finishing process using long sequence reads in a reference-guided assembly process. PBJelly also provides “lift-over” co-ordinate tables to easily port existing annotations to the upgraded assembly. Using PBJelly and long PacBio reads, we upgraded the draft genome sequences of a simulated Drosophila melanogaster, the version 2 draft Drosophila pseudoobscura, an assembly of the Assemblathon 2.0 budgerigar dataset, and a preliminary assembly of the Sooty mangabey. With 24× mapped coverage of PacBio long-reads, we addressed 99% of gaps and were able to close 69% and improve 12% of all gaps in D. pseudoobscura. With 4× mapped coverage of PacBio long-reads we saw reads address 63% of gaps in our budgerigar assembly, of which 32% were closed and 63% improved. With 6.8× mapped coverage of mangabey PacBio long-reads we addressed 97% of gaps and closed 66% of addressed gaps and improved 19%. The accuracy of gap closure was validated by comparison to Sanger sequencing on gaps from the original D. pseudoobscura draft assembly and shown to be dependent on initial reference quality.

Abstract 5113: Functional analysis of genomic variants identified through whole exome sequencing for susceptibility to lymphocytic leukemia

The underlying genetic basis of many childhood cancers remains largely unknown, which places a significant focus on the discovery and understanding of cancer susceptibility genes. Our research is focused on the identification of novel high risk childhood cancer susceptibility genes through next generation sequencing methods and functional assays. One kindred contained four individuals diagnosed during childhood with lymphocytic leukemia or lymphoma transmitted in an autosomal dominant pattern of inheritance. Whole exome sequencing was performed on the constitutional DNA of three of the patients diagnosed with childhood cancer and one non-transmitting parent as an internal control. This analysis resulted in over 1000 unique single nucleotide variant (SNV) sites shared among the three affected family members. Of these SNVs, greater than 95% were found in control databases and excluded. Also, any variant found in the non-transmitting parent was eliminated. Systematic bioinformatics analysis resulted in a short list of seven missense SNVs predicted to be functionally important. Within this list is the missense mutation of a conserved residue (L254P) of Human Cytosolic 5′ Nucleotidase (NT5C1A). NT5C1A has a role in nucleoside metabolism as this protein primarily catalyzes the production of adenosine through the dephosphorylation of AMP. Adenosine metabolism is essential for lymphocyte maturation and viability, making NT5C1A an attractive leukemia susceptibility candidate to pursue further through functional genomic assays. Prior studies have shown that overexpression of wild type NT5C1A in HEK293 cells improves cell survival after treatment with purine and pyrimidine analogs also used in leukemia treatment. We created multiple HEK293 cell lines that stably overexpress either the NT5C1A wild type or L254P mutant protein. These cells were tested for quantitative changes in NT5C1A mediated cytotoxicity to the nucleoside analogs, Cladribine and Gemcitabine. We have replicated the prior findings in cells overexpressing wild type NT5C1A and demonstrated that resistance is dependent on the amount of exogenous NT5C1A protein expressed. Importantly, overexpression of the L254P protein is deficient in mediating resistance to these analogs, which is consistent with this novel variant interfering with NT5C1A activity. This finding has led us to hypothesize that the L254P mutation leads to a reduction in protein function, most likely through reduced substrate binding, which may increase AMP levels in early lymphocytes and predispose to tumorigenesis. Future directions include experiments to determine the L254P protein's substrate binding affinity and activity and the introduction of the mutation into murine leukemia models as well as analysis of L254P as a pharmacogenetic variant.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 5113. doi:1538-7445.AM2012-5113

Patterns and rates of exonic de novo mutations in autism spectrum disorders

Autism spectrum disorders (ASD) are believed to have genetic and environmental origins, yet in only a modest fraction of individuals can specific causes be identified. To identify further genetic risk factors, here we assess the role of de novo mutations in ASD by sequencing the exomes of ASD cases and their parents (n = 175 trios). Fewer than half of the cases (46.3%) carry a missense or nonsense de novo variant, and the overall rate of mutation is only modestly higher than the expected rate. In contrast, the proteins encoded by genes that harboured de novo missense or nonsense mutations showed a higher degree of connectivity among themselves and to previous ASD genes as indexed by protein-protein interaction screens. The small increase in the rate of de novo events, when taken together with the protein interaction results, are consistent with an important but limited role for de novo point mutations in ASD, similar to that documented for de novo copy number variants. Genetic models incorporating these data indicate that most of the observed de novo events are unconnected to ASD; those that do confer risk are distributed across many genes and are incompletely penetrant (that is, not necessarily sufficient for disease). Our results support polygenic models in which spontaneous coding mutations in any of a large number of genes increases risk by 5- to 20-fold. Despite the challenge posed by such models, results from de novo events and a large parallel case-control study provide strong evidence in favour of CHD8 and KATNAL2 as genuine autism risk factors.

Whole-genome sequencing for optimized patient management

Whole-genome sequencing of patient DNA can facilitate diagnosis of a disease, but its potential for guiding treatment has been under-realized. We interrogated the complete genome sequences of a 14-year-old fraternal twin pair diagnosed with dopa (3,4-dihydroxyphenylalanine)-responsive dystonia (DRD; Mendelian Inheritance in Man #128230). DRD is a genetically heterogeneous and clinically complex movement disorder that is usually treated with l-dopa, a precursor of the neurotransmitter dopamine. Whole-genome sequencing identified compound heterozygous mutations in the SPR gene encoding sepiapterin reductase. Disruption of SPR causes a decrease in tetrahydrobiopterin, a cofactor required for the hydroxylase enzymes that synthesize the neurotransmitters dopamine and serotonin. Supplementation of l-dopa therapy with 5-hydroxytryptophan, a serotonin precursor, resulted in clinical improvements in both twins.

Characterization of single-nucleotide variation in Indian-origin rhesus macaques (Macaca mulatta)

BACKGROUND

Rhesus macaques are the most widely utilized nonhuman primate model in biomedical research. Previous efforts have validated fewer than 900 single nucleotide polymorphisms (SNPs) in this species, which limits opportunities for genetic studies related to health and disease. Extensive information about SNPs and other genetic variation in rhesus macaques would facilitate valuable genetic analyses, as well as provide markers for genome-wide linkage analysis and the genetic management of captive breeding colonies.

RESULTS

We used the available rhesus macaque draft genome sequence, new sequence data from unrelated individuals and existing published sequence data to create a genome-wide SNP resource for Indian-origin rhesus monkeys. The original reference animal and two additional Indian-origin individuals were resequenced to low coverage using SOLiD™ sequencing. We then used three strategies to validate SNPs: comparison of potential SNPs found in the same individual using two different sequencing chemistries, and comparison of potential SNPs in different individuals identified with either the same or different sequencing chemistries. Our approach validated approximately 3 million SNPs distributed across the genome. Preliminary analysis of SNP annotations suggests that a substantial number of these macaque SNPs may have functional effects. More than 700 non-synonymous SNPs were scored by Polyphen-2 as either possibly or probably damaging to protein function and these variants now constitute potential models for studying functional genetic variation relevant to human physiology and disease.

CONCLUSIONS

Resequencing of a small number of animals identified greater than 3 million SNPs. This provides a significant new information resource for rhesus macaques, an important research animal. The data also suggests that overall genetic variation is high in this species. We identified many potentially damaging non-synonymous coding SNPs, providing new opportunities to identify rhesus models for human disease.

Whole-genome sequencing in a patient with Charcot-Marie-Tooth neuropathy

BACKGROUND

Whole-genome sequencing may revolutionize medical diagnostics through rapid identification of alleles that cause disease. However, even in cases with simple patterns of inheritance and unambiguous diagnoses, the relationship between disease phenotypes and their corresponding genetic changes can be complicated. Comprehensive diagnostic assays must therefore identify all possible DNA changes in each haplotype and determine which are responsible for the underlying disorder. The high number of rare, heterogeneous mutations present in all humans and the paucity of known functional variants in more than 90% of annotated genes make this challenge particularly difficult. Thus, the identification of the molecular basis of a genetic disease by means of whole-genome sequencing has remained elusive. We therefore aimed to assess the usefulness of human whole-genome sequencing for genetic diagnosis in a patient with Charcot-Marie-Tooth disease.

METHODS

We identified a family with a recessive form of Charcot-Marie-Tooth disease for which the genetic basis had not been identified. We sequenced the whole genome of the proband, identified all potential functional variants in genes likely to be related to the disease, and genotyped these variants in the affected family members.

RESULTS

We identified and validated compound, heterozygous, causative alleles in SH3TC2 (the SH3 domain and tetratricopeptide repeats 2 gene), involving two mutations, in the proband and in family members affected by Charcot-Marie-Tooth disease. Separate subclinical phenotypes segregated independently with each of the two mutations; heterozygous mutations confer susceptibility to neuropathy, including the carpal tunnel syndrome.

CONCLUSIONS

As shown in this study of a family with Charcot-Marie-Tooth disease, whole-genome sequencing can identify clinically relevant variants and provide diagnostic information to inform the care of patients.

Discovery of novel microRNAs in female reproductive tract using next generation sequencing

MicroRNAs (miRNAs) are small non-coding RNAs that mediate post-transcriptional gene silencing. Over 700 human miRNAs have currently been identified, many of which are mutated or de-regulated in diseases. Here we report the identification of novel miRNAs through deep sequencing the small RNAome (<30 nt) of over 100 tissues or cell lines derived from human female reproductive organs in both normal and disease states. These specimens include ovarian epithelium and ovarian cancer, endometrium and endometriomas, and uterine myometrium and uterine smooth muscle tumors. Sequence reads not aligning with known miRNAs were each mapped to the genome to extract flanking sequences. These extended sequence regions were folded in silico to identify RNA hairpins. Sequences demonstrating the ability to form a stem loop structure with low minimum free energy (<-25 kcal) and predicted Drosha and Dicer cut sites yielding a mature miRNA sequence matching the actual sequence were considered putative novel miRNAs. Additional confidence was achieved when putative novel hairpins assembled a collection of sequences highly similar to the putative mature miRNA but with heterogeneous 3'-ends. A confirmed novel miRNA fulfilled these criteria and had its "star" sequence in our collection. We found 7 distinct confirmed novel miRNAs, and 51 additional novel miRNAs that represented highly confident predictions but without detectable star sequences. Our novel miRNAs were detectable in multiple samples, but expressed at low levels and not specific to any one tissue or cell type. To date, this study represents the largest set of samples analyzed together to identify novel miRNAs.

A link between mir-100 and FRAP1/mTOR in clear cell ovarian cancer

MicroRNAs (miRNAs) are small noncoding RNAs that direct gene regulation through translational repression and degradation of complementary mRNA. Although miRNAs have been implicated as oncogenes and tumor suppressors in a variety of human cancers, functional roles for individual miRNAs have not been described in clear cell ovarian carcinoma, an aggressive and chemoresistant subtype of ovarian cancer. We performed deep sequencing to comprehensively profile miRNA expression in 10 human clear cell ovarian cancer cell lines compared with normal ovarian surface epithelial cultures and discovered 54 miRNAs that were aberrantly expressed. Because of the critical roles of the phosphatidylinositol 3-kinase/v-akt murine thymoma viral oncogene homolog 1/mammalian target of rapamycin (mTOR) pathway in clear cell ovarian cancer, we focused on mir-100, a putative tumor suppressor that was the most down-regulated miRNA in our cancer cell lines, and its up-regulated target, FRAP1/mTOR. Overexpression of mir-100 inhibited mTOR signaling and enhanced sensitivity to the rapamycin analog RAD001 (everolimus), confirming the key relationship between mir-100 and the mTOR pathway. Furthermore, overexpression of the putative tumor suppressor mir-22 repressed the EVI1 oncogene, which is known to suppress apoptosis by stimulating phosphatidylinositol 3-kinase/v-akt murine thymoma viral oncogene homolog 1 signaling. In addition to these specific effects, reversing the expression of mir-22 and the putative oncogene mir-182 had widespread effects on target and nontarget gene populations that ultimately caused a global shift in the cancer gene signature toward a more normal state. Our experiments have revealed strong candidate miRNAs and their target genes that may contribute to the pathogenesis of clear cell ovarian cancer, thereby highlighting alternative therapeutic strategies for the treatment of this deadly cancer.

Expression profiling of microRNAs by deep sequencing

MicroRNAs are short non-coding RNAs that regulate the stability and translation of mRNAs. Profiling experiments, using microarray or deep sequencing technology, have identified microRNAs that are preferentially expressed in certain tissues, specific stages of development, or disease states such as cancer. Deep sequencing utilizes massively parallel sequencing, generating millions of small RNA sequence reads from a given sample. Profiling of microRNAs by deep sequencing measures absolute abundance and allows for the discovery of novel microRNAs that have eluded previous cloning and standard sequencing efforts. Public databases provide in silico predictions of microRNA gene targets by various algorithms. To better determine which of these predictions represent true positives, microRNA expression data can be integrated with gene expression data to identify putative microRNA:mRNA functional pairs. Here we discuss tools and methodologies for the analysis of microRNA expression data from deep sequencing.

Mouse let-7 miRNA populations exhibit RNA editing that is constrained in the 5'-seed/ cleavage/anchor regions and stabilize predicted mmu-let-7a:mRNA duplexes

Massively parallel sequencing of millions of < 30-nt RNAs expressed in mouse ovary, embryonic pancreas (E14.5), and insulin-secreting beta-cells (betaTC-3) reveals that approximately 50% of the mature miRNAs representing mostly the mmu-let-7 family display internal insertion/deletions and substitutions when compared to precursor miRNA and the mouse genome reference sequences. Approximately, 12%-20% of species associated with mmu-let-7 populations exhibit sequence discrepancies that are dramatically reduced in nucleotides 3-7 (5'-seed) and 10-15 (cleavage and anchor sites). This observation is inconsistent with sequencing error and leads us to propose that the changes arise predominantly from post-transcriptional RNA-editing activity operating on miRNA:target mRNA complexes. Internal nucleotide modifications are most enriched at the ninth nucleotide position. A common ninth base edit of U-to-G results in a significant increase in stability of down-regulated let-7a targets in inhibin-deficient mice (Inha-/-). An excess of U-insertions (14.8%) over U-deletions (1.5%) and the presence of cleaved intermediates suggest that a mammalian TUTase (terminal uridylyl transferase) mediated dUTP-dependent U-insertion/U-deletion cycle may be a possible mechanism. We speculate that mRNA target site-directed editing of mmu-let-7a duplex-bulges stabilizes "loose" miRNA:mRNA target associations and functions to expand the target repertoire and/or enhance mRNA decay over translational repression. Our results also demonstrate that the systematic study of sequence variation within specific RNA classes in a given cell type from millions of sequences generated by next-generation sequencing (NGS) technologies ("intranomics") can be used broadly to infer functional constraints on specific parts of completely uncharacterized RNAs.

Novel microRNA candidates and miRNA-mRNA pairs in embryonic stem (ES) cells

BACKGROUND

MicroRNAS (miRNAS: a class of short non-coding RNAs) are emerging as important agents of post transcriptional gene regulation and integral components of gene networks. MiRNAs have been strongly linked to stem cells, which have a remarkable dual role in development. They can either continuously replenish themselves (self-renewal), or differentiate into cells that execute a limited number of specific actions (pluripotence).

METHODOLOGY/PRINCIPAL FINDINGS

In order to identify novel miRNAs from narrow windows of development we carried out an in silico search for micro-conserved elements (MCE) in adult tissue progenitor transcript sequences. A plethora of previously unknown miRNA candidates were revealed including 545 small RNAs that are enriched in embryonic stem (ES) cells over adult cells. Approximately 20% of these novel candidates are down-regulated in ES (Dicer(-/-)) ES cells that are impaired in miRNA maturation. The ES-enriched miRNA candidates exhibit distinct and opposite expression trends from mmu-mirs (an abundant class in adult tissues) during retinoic acid (RA)-induced ES cell differentiation. Significant perturbation of trends is found in both miRNAs and novel candidates in ES (GCNF(-/-)) cells, which display loss of repression of pluripotence genes upon differentiation.

CONCLUSION/SIGNIFICANCE

Combining expression profile information with miRNA target prediction, we identified miRNA-mRNA pairs that correlate with ES cell pluripotence and differentiation. Perturbation of these pairs in the ES (GCNF(-/-)) mutant suggests a role for miRNAs in the core regulatory networks underlying ES cell self-renewal, pluripotence and differentiation.

Computational and transcriptional evidence for microRNAs in the honey bee genome

BACKGROUND

Non-coding microRNAs (miRNAs) are key regulators of gene expression in eukaryotes. Insect miRNAs help regulate the levels of proteins involved with development, metabolism, and other life history traits. The recently sequenced honey bee genome provides an opportunity to detect novel miRNAs in both this species and others, and to begin to infer the roles of miRNAs in honey bee development.

RESULTS

Three independent computational surveys of the assembled honey bee genome identified a total of 65 non-redundant candidate miRNAs, several of which appear to have previously unrecognized orthologs in the Drosophila genome. A subset of these candidate miRNAs were screened for expression by quantitative RT-PCR and/or genome tiling arrays and most predicted miRNAs were confirmed as being expressed in at least one honey bee tissue. Interestingly, the transcript abundance for several known and novel miRNAs displayed caste or age-related differences in honey bees. Genes in proximity to miRNAs in the bee genome are disproportionately associated with the Gene Ontology terms 'physiological process', 'nucleus' and 'response to stress'.

CONCLUSION

Computational approaches successfully identified miRNAs in the honey bee and indicated previously unrecognized miRNAs in the well-studied Drosophila melanogaster genome despite the 280 million year distance between these insects. Differentially transcribed miRNAs are likely to be involved in regulating honey bee development, and arguably in the extreme developmental switch between sterile worker bees and highly fertile queens.

Forward neutral pion production in p + p and d + Au collisions at square root sNN=200 GeV

Measurements of the production of forward pi0 mesons from p + p and d + Au collisions at square root sNN=200 GeV are reported. The p + p yield generally agrees with next-to-leading order perturbative QCD calculations. The d + Au yield per binary collision is suppressed as eta increases, decreasing to approximately 30% of the p + p yield at eta =4.00, well below shadowing expectations. Exploratory measurements of azimuthal correlations of the forward pi0 with charged hadrons at eta approximately 0 show a recoil peak in p + p that is suppressed in d + Au at low pion energy. These observations are qualitatively consistent with a saturation picture of the low-x gluon structure of heavy nuclei.

Distributions of charged hadrons associated with high transverse momentum particles in pp and Au+Au collisions at sqrt[sNN]=200 GeV

Charged hadrons in [EQUATION: SEE TEXT] associated with particles of [EQUATION: SEE TEXT] are reconstructed in pp and Au+Au collisions at sqrt[sNN]=200 GeV. The associated multiplicity and p magnitude sum are found to increase from pp to central Au+Au collisions. The associated p distributions, while similar in shape on the nearside, are significantly softened on the awayside in central Au+Au relative to pp and not much harder than that of inclusive hadrons. The results, consistent with jet quenching, suggest that the awayside fragments approach equilibration with the medium traversed.

Multistrange Baryon elliptic flow in Au+Au collisions at square root of sNN=200 GeV

We report on the first measurement of elliptic flow v2(pT) of multistrange baryons Xi- +Xi+ and Omega- + Omega+ in heavy-ion collisions. In minimum-bias Au+Au collisions at square root of s(NN)=200 GeV, a significant amount of elliptic flow, comparable to other nonstrange baryons, is observed for multistrange baryons which are expected to be particularly sensitive to the dynamics of the partonic stage of heavy-ion collisions. The pT dependence of v2 of the multistrange baryons confirms the number of constituent quark scaling previously observed for lighter hadrons. These results support the idea that a substantial fraction of the observed collective motion is developed at the early partonic stage in ultrarelativistic nuclear collisions at the Relativistic Heavy Ion Collider.

Multiplicity and pseudorapidity distributions of photons in Au+Au collisions at square root of (S(NN)) = 62.4 GeV

We present the first measurement of pseudorapidity distribution of photons in the region 2.3 < or = eta < or = 3.7 for different centralities in Au+Au collisions at square root of (S(NN)) = 62.4 GeV. We find that the photon yield scales with the number of participating nucleons at all collision centralities studied. The pseudorapidity distribution of photons, dominated by pi0 decays, has been compared to those of charged pions, photons, and inclusive charged particles from heavy-ion and nucleon-nucleon collisions at various energies. The photon production has been shown to be consistent with the energy and centrality independent limiting fragmentation scenario.

Omega- and Omega+ production in central Pb + Pb collisions at 40 and 158A GeV

Results are presented on Omega production in central Pb+Pb collisions at 40 and 158A GeV beam energy. For the first time in heavy ion reactions, rapidity distributions and total yields were measured for the sum Omega(-) + Omega(+) at 40A GeV and for Omega(-) and Omega(+) separately at 158A GeV. The yields are strongly underpredicted by the string-hadronic UrQMD model but agree better with predictions from hadron gas models.

Open charm yields in d+Au collisions at squareroot[sNN]=200 GeV

Midrapidity open charm spectra from direct reconstruction of D0(D0)-->K-/+pi+/- in d+Au collisions and indirect electron-positron measurements via charm semileptonic decays in p+p and d+Au collisions at squareroot[sNN]=200 GeV are reported. The D0(D0) spectrum covers a transverse momentum (pT) range of 0.1<pT<3 GeV/c, whereas the electron spectra cover a range of 1<pT<4 GeV/c. The electron spectra show approximate binary collision scaling between p+p and d+Au collisions. From these two independent analyses, the differential cross section per nucleon-nucleon binary interaction at midrapidity for open charm production from d+Au collisions at BNL RHIC is dsigma(NN)cc/dy=0.30+/-0.04(stat)+/-0.09(syst) mb. The results are compared to theoretical calculations. Implications for charmonium results in A+A collisions are discussed.

System-size dependence of strangeness production in nucleus-nucleus collisions at squareroot[sNN]=17.3 GeV

Emission of pi+/-, K+/-, phi, and Lambda was measured in near-central C+C and Si+Si collisions at 158 AGeV beam energy. Together with earlier data for p+p, S+S, and Pb+Pb, the system-size dependence of relative strangeness production in nucleus-nucleus collisions is obtained. Its fast rise and the saturation observed at about 60 participating nucleons can be understood as the onset of the formation of coherent systems of increasing size.

Azimuthal anisotropy and correlations at large transverse momenta in p + p and Au + Au collisions at square root sNN=200 GeV

Results on high transverse momentum charged particle emission with respect to the reaction plane are presented for Au + Au collisions at square root s(NN)=200 GeV. Two- and four-particle correlations results are presented as well as a comparison of azimuthal correlations in Au + Au collisions to those in p + p at the same energy. The elliptic anisotropy v(2) is found to reach its maximum at p(t) approximately 3 GeV/c, then decrease slowly and remain significant up to p(t) approximately 7-10 GeV/c. Stronger suppression is found in the back-to-back high-p(t) particle correlations for particles emitted out of plane compared to those emitted in plane. The centrality dependence of v(2) at intermediate p(t) is compared to simple models based on jet quenching.