Gene duplication: a drive for phenotypic diversity and cause of human disease

Gene duplication is one of the key factors driving genetic innovation, i.e., producing novel genetic variants. Although the contribution of whole-genome and segmental duplications to phenotypic diversity across species is widely appreciated, the phenotypic spectrum and potential pathogenicity of small-scale duplications in individual genomes are less well explored. This review discusses the nature of small-scale duplications and the phenotypes produced by such duplications. Phenotypic variation and disease phenotypes induced by duplications are more diverse and widespread than previously anticipated, and duplications are a major class of disease-related genomic variation. Pathogenic duplications particularly involve dosage-sensitive genes with both similar and dissimilar over- and underexpression phenotypes, and genes encoding proteins with a propensity to aggregate. Phenotypes related to human-specific copy number variation in genes regulating environmental responses and immunity are increasingly recognized. Small genomic duplications containing defense-related genes also contribute to complex common phenotypes.

Fatal adenovirus hepatitis during maintenance therapy for childhood acute lymphoblastic leukemia

Disseminated adenoviral infection with hepatitis is rare in children undergoing standard chemotherapy. We report on a 3(1/2)-year-old male with fatal adenovirus hepatitis receiving maintenance chemotherapy for acute lymphoblastic leukemia (ALL). Adenoviral hepatitis was proven by histology, viral culture, and PCR in a liver biopsy. Quantitative real-time PCR in the peripheral blood showed adenoviral DNA copy number >10(9)/ml. Despite aggressive supportive care and antiviral treatment with cidofovir, the patient died rapidly due to fulminant liver failure. Diagnostic and treatment options for adenovirus infection remain unsatisfactory for these patients. We propose suggestions for diagnosis and therapy.

Enhancement of heterologous expression of alkaline phytase in Pichia pastors

Phytic acid is the major storage form of phosphorus and inositol in seeds and legumes. It forms insoluble phytate salts by chelating with positively charged mineral ions. Non-ruminant animals are not able to digest phytate due to the lack of phytases in their GI tracks, thus the undigested phytate is excreted leading to environmental contamination. Supplementation with phytases in animal feed has proven to be an effective strategy to alleviate nutritional and environmental issues. The unique catalytic and thermal stability properties of alkaline phytase from lily pollen (LlALP) suggest that it has the potential to be useful as a feed supplement. Our goal is to develop a method for the production of substantial amounts of rLlALP for animal feed and structural studies. rLlALP2 has been successfully expressed in the yeast, Pichia pastoris. However, expression yield was modest (8-10 mg/L). Gene copy number has been identified as an important parameter in enhancing protein yields. Multicopy clones were selected using Zeocin-resistance-based vectors and challenging transformants to high Zeocin levels under different conditions. Data indicate that increasing selection pressure led to the generation of clones with amplification of both rLlAlp2 and Zeor genes and the two genes were not equally amplified. Additionally, clones generated by step-wise methods led to clones with greater amplification. The effects of transgene copy number and gene sequence optimization on expression levels of rLlALP2 were examined. The data indicate that increasing the copy number of rLlAlp2 in transformed clones was detrimental to expression level. The use of a sequence-optimized rLlAlp2 (op-rLlAlp2) increased expression yield of the active enzyme by 25-50%, suggesting that transcription and translation efficiency are not major bottlenecks in the production of rLlALP2. Lowering induction temperature to 20 oC led to an increase in enzyme activity of 1.2 to 20-fold, suggesting that protein folding or post-translational processes may be limiting factors for rLlALP2 production. Cumulatively, optimization of copy number, gene sequence optimization and reduced temperature led to increase of rLlALP2 enzyme activity by three-fold (25-30 mg/L). In an effort to simplify the purification process of rLlALP2, extracellular expression of phytase was investigated. Extracellular expression is dependent on the presence of an appropriate secretion signal upstream of the transgene native signal peptide(s) present in the transgene may also influence secretion efficiency. The data suggest that deletion of both N- and C-terminal signal peptides of rLlALP2 enhanced α-mating factor (α-MF)-driven secretion of LlALP2 by four-fold. The secretion signal peptide of chicken egg white lysozyme was ineffective in secretion rLlALP2 in P. pastoris. To enhance rLlALP2 secretion, effectiveness of the strong inducible promoter (PAOX1) was compared with the constitutive promoter (PGAP). The intracellular yield of rLlALP2 was about four-fold greater under the control of PGAP compared to PAOX1 and extracellular expression level of rLlALP2 was around eight-fold (75-100 mg/L) greater. The successful production of active rLlALP2 in P. pastoris will allow us to conduct the animal feed supplementation studies and structural studies.

PHYLOGENETIC ANALYSES, GENE DOSAGE EFFECT AND SELECTION IN MEMBERS OF THE LATERAL BOUNDARY DOMAIN (LBD) GENES TRANSCRIPTION FACTOR FAMILY

Wood formation is an economically and environmentally important process and has played a significant role in the evolution of terrestrial plants. Despite its significance, the molecular underpinnings of the process are still poorly understood. We have previously shown that four Lateral Boundary Domain (LBD) transcription factors have important roles in the regulation of wood formation with two (LBD1 and LBD4) involved in secondary phloem and ray cell development and two (LBD15 and LBD18) in secondary xylem formation. Here, we used comparative phylogenetic analyses to test potential roles of the four LBD genes in the evolution of woodiness. We studied the copy number and variation in DNA and amino acid sequences of the four LBDs in a wide range of woody and herbaceous plant taxa with fully sequenced and annotated genomes. LBD1 showed the highest gene copy number across the studied species, and LBD1 gene copy number was strongly and significantly correlated with the level of ray seriation. The lianas, cucumber and grape, with multiseriate ray cells showed the highest gene copy number (12 and 11, respectively). Because lianas’ growth habit requires significant twisting and bending, the less lignified ray parenchyma cells likely facilitate stem flexibility and maintenance of xylem conductivity. We further demonstrate conservation of amino acids in the LBD18 protein sequences that are specific to woody taxa. Neutrality tests showed evidence for strong purifying selection on these gene regions across various orders, indicating adaptive convergent evolution of LBD18. Structural modeling demonstrates that the conserved amino acids have a significant impact on the tertiary protein structure and thus are likely of significant functional importance.

Regulation of phosphoinositide 3-kinase expression in health and disease

Both the biology and the therapeutic potential of the phosphoinositide 3-kinase (PI3K) signalling axis have been the subject of intense investigation; however, little is known about the regulation of PI3K expression. Emerging evidence indicates that PI3K levels change in response to cellular stimulation with insulin and nuclear receptor ligands, and during various physiological and pathological processes including differentiation, regeneration, hypertension and cancer. Recently identified mechanisms that control PI3K production include increased gene copy number in cancer, and transcriptional regulation of the p110alpha PI3K gene by FOXO3a, NF-kappaB and p53, and of the PI3K regulatory subunits by STAT3, EBNA-2 and SREBP. In most instances, however, the impact of alterations in PI3K expression on PI3K signalling and disease remains to be established.

Clustering subspecies of Aeromonas salmonicida using IS630 typing.

BACKGROUND The insertion element IS630 found in Aeromonas salmonicida belongs to the IS630-Tc1-mariner superfamily of transposons. It is present in multiple copies and represents approximately half of the IS present in the genome of A. salmonicida subsp. salmonicida A449. RESULTS By using High Copy Number IS630 Restriction Fragment Length Polymorphism (HCN-IS630-RFLP), strains of various subspecies of Aeromonas salmonicida showed conserved or clustering patterns, thus allowing their differentiation from each other. Fingerprints of A. salmonicida subsp. salmonicida showed the highest homogeneity while 'atypical' A. salmonicida strains were more heterogeneous. IS630 typing also differentiated A. salmonicida from other Aeromonas species. The copy number of IS630 in Aeromonas salmonicida ranges from 8 to 35 and is much lower in other Aeromonas species. CONCLUSIONS HCN-IS630-RFLP is a powerful tool for subtyping of A. salmonicida. The high stability of IS630 insertions in A. salmonicida subsp. salmonicida indicates that it might have played a role in pathoadaptation of A. salmonicida which has reached an optimal configuration in the highly virulent and specific fish pathogen A. salmonicida subsp. salmonicida.

A non-coding genomic duplication at the HMX1 locus is associated with crop ears in highland cattle.

Highland cattle with congenital crop ears have notches of variable size on the tips of both ears. In some cases, cartilage deformation can be seen and occasionally the external ears are shortened. We collected 40 cases and 80 controls across Switzerland. Pedigree data analysis confirmed a monogenic autosomal dominant mode of inheritance with variable expressivity. All affected animals could be traced back to a single common ancestor. A genome-wide association study was performed and the causative mutation was mapped to a 4 Mb interval on bovine chromosome 6. The H6 family homeobox 1 (HMX1) gene was selected as a positional and functional candidate gene. By whole genome re-sequencing of an affected Highland cattle, we detected 6 non-synonymous coding sequence variants and two variants in an ultra-conserved element at the HMX1 locus with respect to the reference genome. Of these 8 variants, only a non-coding 76 bp genomic duplication (g.106720058_106720133dup) located in the conserved region was perfectly associated with crop ears. The identified copy number variation probably results in HMX1 misregulation and possible gain-of-function. Our findings confirm the role of HMX1 during the development of the external ear. As it is sometimes difficult to phenotypically diagnose Highland cattle with slight ear notches, genetic testing can now be used to improve selection against this undesired trait.

RNA-dependent genome processing during nuclear differentiation: the model systems of stichotrichous ciliates

We introduce ciliated protozoa, and more specifically the stichotrichous ciliates Oxytricha and Stylonychia, as biological model systems for the analysis of programmed DNA-reorganization processes during nuclear differentiation. These include DNA excision, DNA elimination, reordering of gene segments and specific gene amplification. We show that small nuclear RNAs specify DNA sequences to be excised or retained, but also discuss the need for a RNA template molecule derived from the parental nucleus for these processes. This RNA template guides reordering of gene segments to become functional genes and determines gene copy number in the differentiated nucleus. Since the template is derived from the parental macronucleus, gene reordering and DNA amplification are inherited in a non-Mendelian epigenetic manner.

Estrogen receptor β expression and androgen receptor phosphorylation correlate with a poor clinical outcome in hormone-naive prostate cancer and are elevated in castration-resistant disease

Patients with advanced prostate cancer (PC) are usually treated with androgen withdrawal. While this therapy is initially effective, nearly all PCs become refractory to it. As hormone receptors play a crucial role in this process, we constructed a tissue microarray consisting of PC samples from 107 hormone-naïve (HN) and 101 castration-resistant (CR) PC patients and analyzed the androgen receptor (AR) gene copy number and the protein expression profiles of AR, Serin210-phosphorylated AR (pAR(210)), estrogen receptor (ER)β, ERα and the proliferation marker Ki67. The amplification of the AR gene was virtually restricted to CR PC and was significantly associated with increased AR protein expression (P<0.0001) and higher tumor cell proliferation (P=0.001). Strong AR expression was observed in a subgroup of HN PC patients with an adverse prognosis. In contrast, the absence of AR expression in CR PC was significantly associated with a poor overall survival. While pAR(210) was predominantly found in CR PC patients (P<0.0001), pAR(210) positivity was observed in a subgroup of HN PC patients with a poor survival (P<0.05). Epithelial ERα expression was restricted to CR PC cells (9%). ERβ protein expression was found in 38% of both HN and CR PCs, but was elevated in matched CR PC specimens. Similar to pAR(210), the presence of ERβ in HN patients was significantly associated with an adverse prognosis (P<0.005). Our results strongly suggest a major role for pAR(210) and ERβ in HN PC. The expression of these markers might be directly involved in CR tumor growth.

Prediction of outcome in patients with low-grade squamous intraepithelial lesions by fluorescence in situ hybridization analysis of human papillomavirus, TERC, and MYC

BACKGROUND Cytology is an excellent method with which to diagnose preinvasive lesions of the uterine cervix, but it suffers from limited specificity for clinically significant lesions. Supplementary methods might predict the natural course of the detected lesions. The objective of the current study was to test whether a multicolor fluorescence in situ hybridization (FISH) assay might help to stratify abnormal results of Papanicolaou tests. METHODS A total of 219 liquid-based cytology specimens of low-grade squamous intraepithelial lesions (LSIL), 49 atypical squamous cells of undetermined significance (ASCUS) specimens, 52 high-grade squamous intraepithelial lesion (HSIL) specimens, and 50 normal samples were assessed by FISH with probes for the human papillomavirus (HPV), MYC, and telomerase RNA component (TERC). Subtyping of HPV by polymerase chain reaction (PCR) was performed in a subset of cases (n=206). RESULTS There was a significant correlation found between HPV detection by FISH and PCR (P<.0001). In patients with LSILs, the presence of HPV detected by FISH was significantly associated with disease progression (P<.0001). An increased MYC and/or TERC gene copy number (>2 signals in>10% of cells) prevailed in 43% of ASCUS specimens and was more frequent in HSIL (85%) than in LSIL (33%) (HSIL vs LSIL: P<.0001). Increased TERC gene copy number was significantly correlated with progression of LSIL (P<.01; odds ratio, 7.44; area under the receiver operating characteristic curve, 0.73; positive predictive value, 0.30; negative predictive value, 0.94) CONCLUSIONS: The detection of HPV by FISH analysis is feasible in liquid-based cytology and is significantly correlated with HPV analysis by PCR. The analysis of TERC gene copy number may be useful for risk stratification in patients with LSIL.

A genome-wide association study of resistance to HIV infection in highly exposed uninfected individuals with hemophilia A.

Human genetic variation contributes to differences in susceptibility to HIV-1 infection. To search for novel host resistance factors, we performed a genome-wide association study (GWAS) in hemophilia patients highly exposed to potentially contaminated factor VIII infusions. Individuals with hemophilia A and a documented history of factor VIII infusions before the introduction of viral inactivation procedures (1979-1984) were recruited from 36 hemophilia treatment centers (HTCs), and their genome-wide genetic variants were compared with those from matched HIV-infected individuals. Homozygous carriers of known CCR5 resistance mutations were excluded. Single nucleotide polymorphisms (SNPs) and inferred copy number variants (CNVs) were tested using logistic regression. In addition, we performed a pathway enrichment analysis, a heritability analysis, and a search for epistatic interactions with CCR5 Δ32 heterozygosity. A total of 560 HIV-uninfected cases were recruited: 36 (6.4%) were homozygous for CCR5 Δ32 or m303. After quality control and SNP imputation, we tested 1 081 435 SNPs and 3686 CNVs for association with HIV-1 serostatus in 431 cases and 765 HIV-infected controls. No SNP or CNV reached genome-wide significance. The additional analyses did not reveal any strong genetic effect. Highly exposed, yet uninfected hemophiliacs form an ideal study group to investigate host resistance factors. Using a genome-wide approach, we did not detect any significant associations between SNPs and HIV-1 susceptibility, indicating that common genetic variants of major effect are unlikely to explain the observed resistance phenotype in this population.

Rtr1 is the Saccharomyces cerevisiae homolog of a novel family of RNA polymerase II-binding proteins.

Cells must rapidly sense and respond to a wide variety of potentially cytotoxic external stressors to survive in a constantly changing environment. In a search for novel genes required for stress tolerance in Saccharomyces cerevisiae, we identified the uncharacterized open reading frame YER139C as a gene required for growth at 37 degrees C in the presence of the heat shock mimetic formamide. YER139C encodes the closest yeast homolog of the human RPAP2 protein, recently identified as a novel RNA polymerase II (RNAPII)-associated factor. Multiple lines of evidence support a role for this gene family in transcription, prompting us to rename YER139C RTR1 (regulator of transcription). The core RNAPII subunits RPB5, RPB7, and RPB9 were isolated as potent high-copy-number suppressors of the rtr1Delta temperature-sensitive growth phenotype, and deletion of the nonessential subunits RPB4 and RPB9 hypersensitized cells to RTR1 overexpression. Disruption of RTR1 resulted in mycophenolic acid sensitivity and synthetic genetic interactions with a number of genes involved in multiple phases of transcription. Consistently, rtr1Delta cells are defective in inducible transcription from the GAL1 promoter. Rtr1 constitutively shuttles between the cytoplasm and nucleus, where it physically associates with an active RNAPII transcriptional complex. Taken together, our data reveal a role for members of the RTR1/RPAP2 family as regulators of core RNAPII function.

The Role of Acidic Fibroblast Growth Factor and Fibroblast Growth Factor Receptor 4 in High Grade Serous Carcinoma

Background: High grade serous carcinoma whether ovarian, tubal or primary peritoneal, continues to be the most lethal gynecologic malignancy in the USA. Although combination chemotherapy and aggressive surgical resection has improved survival in the past decade the majority of patients still succumb to chemo-resistant disease recurrence. It has recently been reported that amplification of 5q31-5q35.3 is associated with poor prognosis in patients with high grade serous ovarian carcinoma. Although the amplicon contains over 50 genes, it is notable for the presence of several members of the fibroblast growth factor signaling axis. In particular acidic fibroblast growth factor (FGF1) has been demonstrated to be one of the driving genes in mediating the observed prognostic effect of the amplicon in ovarian cancer patients. This study seeks to further validate the prognostic value of fibroblast growth receptor 4 (FGFR4), another candidate gene of the FGF/FGFR axis located in the same amplicon. The emphasis will be delineating the role the FGF1/FGFR4 signaling axis plays in high grade serous ovarian carcinoma; and test the feasibility of targeting the FGF1/FGFR4 axis therapeutically. Materials and Methods: Spearman and Pearson correlation studies on data generated from array CGH and transcriptome profiling analyses on 51 microdissected tumor samples were used to identify genes located on chromosome 5q31-35.3 that showed significant correlation between DNA and mRNA copy numbers. Significant correlation between FGF1 and FGFR4 DNA copy numbers was further validated by qPCR analysis on DNA isolated from 51 microdissected tumor samples. Immunolocalization and quantification of FGFR4 expression were performed on paraffin embedded tissue samples from 183 cases of high-grade serous ovarian carcinoma. The expression was then correlated with clinical data to assess impact on survival. The expression of FGF1 and FGFR4 in vitro was quantified by real-time PCR and western blotting in six high-grade serous ovarian carcinoma cell lines and compared to those in human ovarian surface epithelial cells to identify overexpression. The effect of FGF1 on these cell lines after serum starvation was quantified for in vitro cellular proliferation, migration/invasion, chemoresistance and survival utilizing a combination of commercially available colorimetric, fluorometric and electrical impedance assays. FGFR4 expression was then transiently silenced via siRNA transfection and the effects on response to FGF1, cellular proliferation, and migration were quantified. To identify relevant cellular pathways involved, responsive cell lines were transduced with different transcription response elements using the Cignal-Lenti reporter system and treated with FGF1 with and without transient FGFR4 knock down. This was followed by western blot confirmation for the relevant phosphoproteins. Anti-FGF1 antibodies and FGFR trap proteins were used to attempt inhibition of FGF mediated phenotypic changes and relevant signaling in vitro. Orthotopic intraperitoneal tumors were established in nude mice using serous cell lines that have been previously transfected with luciferase expressing constructs. The mice were then treated with FGFR trap protein. Tumor progression was then followed via bioluminescent imaging. The FGFR4 gene from 52 clinical samples was sequenced to screen for mutations. Results: FGFR4 DNA and mRNA copy numbers were significantly correlated and FGFR4 DNA copy number was significantly correlated with that of FGF1. Survival of patients with high FGFR4 expressing tumors was significantly shorter that those with low expression(median survival 28 vs 55 month p< 0.001) In a multivariate cox regression model FGFR expression significantly increased risk of death (HR 2.1, p<0.001). FGFR4 expression was significantly higher in all cell lines tested compared to HOSE, OVCA432 cell line in particular had very high expression suggesting amplification. FGF1 was also particularly overexpressed in OVCA432. FGF1 significantly increased cell survival after serum deprivation in all cell lines. Transient knock down of FGFR4 caused significant reduction in cell migration and proliferation in vitro and significantly decreased the proliferative effects of FGF1 in vitro. FGFR1, FGFR4 traps and anti-FGF1 antibodies did not show activity in vitro. OVCA432 transfected with the cignal lenti reporter system revealed significant activation of MAPK, NFkB and WNT pathways, western blotting confirmed the results. Reverse phase protein array (RPPA) analysis also showed activation of MAPK, AKT, WNT pathways and down regulation of E Cadherin. FGFR trap protein significantly reduced tumor growth in vivo in an orthotopic mouse model. Conclusions: Overexpression and amplification of several members of the FGF signaling axis present on the amplicon 5q31-35.3 is a negative prognostic indicator in high grade serous ovarian carcinoma and may drive poor survival associated with that amplicon. Activation of The FGF signaling pathway leads to downstream activation of MAPK, AKT, WNT and NFkB pathways leading to a more aggressive cancer phenotype with increased tumor growth, evasion of apoptosis and increased migration and invasion. Inhibition of FGF pathway in vivo via FGFR trap protein leads to significantly decreased tumor growth in an orthotopic mouse model.

DUAL ENERGY COMPUTED TOMOGRAPHY FOR PROTON THERAPY TREATMENT PLANNING

Proton radiation therapy is gaining popularity because of the unique characteristics of its dose distribution, e.g., high dose-gradient at the distal end of the percentage-depth-dose curve (known as the Bragg peak). The high dose-gradient offers the possibility of delivering high dose to the target while still sparing critical organs distal to the target. However, the high dose-gradient is a double-edged sword: a small shift of the highly conformal high-dose area can cause the target to be substantially under-dosed or the critical organs to be substantially over-dosed. Because of that, large margins are required in treatment planning to ensure adequate dose coverage of the target, which prevents us from realizing the full potential of proton beams. Therefore, it is critical to reduce uncertainties in the proton radiation therapy. One major uncertainty in a proton treatment is the range uncertainty related to the estimation of proton stopping power ratio (SPR) distribution inside a patient. The SPR distribution inside a patient is required to account for tissue heterogeneities when calculating dose distribution inside the patient. In current clinical practice, the SPR distribution inside a patient is estimated from the patient’s treatment planning computed tomography (CT) images based on the CT number-to-SPR calibration curve. The SPR derived from a single CT number carries large uncertainties in the presence of human tissue composition variations, which is the major drawback of the current SPR estimation method. We propose to solve this problem by using dual energy CT (DECT) and hypothesize that the range uncertainty can be reduced by a factor of two from currently used value of 3.5%. A MATLAB program was developed to calculate the electron density ratio (EDR) and effective atomic number (EAN) from two CT measurements of the same object. An empirical relationship was discovered between mean excitation energies and EANs existing in human body tissues. With the MATLAB program and the empirical relationship, a DECT-based method was successfully developed to derive SPRs for human body tissues (the DECT method). The DECT method is more robust against the uncertainties in human tissues compositions than the current single-CT-based method, because the DECT method incorporated both density and elemental composition information in the SPR estimation. Furthermore, we studied practical limitations of the DECT method. We found that the accuracy of the DECT method using conventional kV-kV x-ray pair is susceptible to CT number variations, which compromises the theoretical advantage of the DECT method. Our solution to this problem is to use a different x-ray pair for the DECT. The accuracy of the DECT method using different combinations of x-ray energies, i.e., the kV-kV, kV-MV and MV-MV pair, was compared using the measured imaging uncertainties for each case. The kV-MV DECT was found to be the most robust against CT number variations. In addition, we studied how uncertainties propagate through the DECT calculation, and found general principles of selecting x-ray pairs for the DECT method to minimize its sensitivity to CT number variations. The uncertainties in SPRs estimated using the kV-MV DECT were analyzed further and compared to those using the stoichiometric method. The uncertainties in SPR estimation can be divided into five categories according to their origins: the inherent uncertainty, the DECT modeling uncertainty, the CT imaging uncertainty, the uncertainty in the mean excitation energy, and SPR variation with proton energy. Additionally, human body tissues were divided into three tissue groups – low density (lung) tissues, soft tissues and bone tissues. The uncertainties were estimated separately because their uncertainties were different under each condition. An estimate of the composite range uncertainty (2s) was determined for three tumor sites – prostate, lung, and head-and-neck, by combining the uncertainty estimates of all three tissue groups, weighted by their proportions along typical beam path for each treatment site. In conclusion, the DECT method holds theoretical advantages in estimating SPRs for human tissues over the current single-CT-based method. Using existing imaging techniques, the kV-MV DECT approach was capable of reducing the range uncertainty from the currently used value of 3.5% to 1.9%-2.3%, but it is short to reach our original goal of reducing the range uncertainty by a factor of two. The dominant source of uncertainties in the kV-MV DECT was the uncertainties in CT imaging, especially in MV CT imaging. Further reduction in beam hardening effect, the impact of scatter, out-of-field object etc. would reduce the Hounsfeld Unit variations in CT imaging. The kV-MV DECT still has the potential to reduce the range uncertainty further.

The Draft Assembly of the Radically Organized Stylonychia lemnae Macronuclear Genome.

Stylonychia lemnae is a classical model single-celled eukaryote, and a quintessential ciliate typified by dimorphic nuclei: A small, germline micronucleus and a massive, vegetative macronucleus. The genome within Stylonychia's macronucleus has a very unusual architecture, comprised variably and highly amplified "nanochromosomes," each usually encoding a single gene with a minimal amount of surrounding noncoding DNA. As only a tiny fraction of the Stylonychia genes has been sequenced, and to promote research using this organism, we sequenced its macronuclear genome. We report the analysis of the 50.2-Mb draft S. lemnae macronuclear genome assembly, containing in excess of 16,000 complete nanochromosomes, assembled as less than 20,000 contigs. We found considerable conservation of fundamental genomic properties between S. lemnae and its close relative, Oxytricha trifallax, including nanochromosomal gene synteny, alternative fragmentation, and copy number. Protein domain searches in Stylonychia revealed two new telomere-binding protein homologs and the presence of linker histones. Among the diverse histone variants of S. lemnae and O. trifallax, we found divergent, coexpressed variants corresponding to four of the five core nucleosomal proteins (H1.2, H2A.6, H2B.4, and H3.7) suggesting that these ciliates may possess specialized nucleosomes involved in genome processing during nuclear differentiation. The assembly of the S. lemnae macronuclear genome demonstrates that largely complete, well-assembled highly fragmented genomes of similar size and complexity may be produced from one library and lane of Illumina HiSeq 2000 shotgun sequencing. The provision of the S. lemnae macronuclear genome sets the stage for future detailed experimental studies of chromatin-mediated, RNA-guided developmental genome rearrangements.