Hox Gene Variation Drives Morphological Specialization of Humpback Grouper Cromileptes altivelis.

BackgroundEntamoeba histolytica is a significant cause of disease worldwide. However, little is known about the genetic diversity of the parasite. We re-sequenced the genomes of ten laboratory cultured lines of the eukaryotic pathogen Entamoeba histolytica in order to develop a picture of genetic diversity across the genome.ResultsThe extreme nucleotide composition bias and repetitiveness of the E. histolytica genome provide a challenge for short-read mapping, yet we were able to define putative single nucleotide polymorphisms in a large portion of the genome. The results suggest a rather low level of single nucleotide diversity, although genes and gene families with putative roles in virulence are among the more polymorphic genes. We did observe large differences in coverage depth among genes, indicating differences in gene copy number between genomes. We found evidence indicating that recombination has occurred in the history of the sequenced genomes, suggesting that E. histolytica may reproduce sexually.ConclusionsE. histolytica displays a relatively low level of nucleotide diversity across its genome. However, large differences in gene family content and gene copy number are seen among the sequenced genomes. The pattern of polymorphism indicates that E. histolytica reproduces sexually, or has done so in the past, which has previously been suggested but not proven.

We have examined the hypothesis that MHC ancestral haplotypes have a specific content of genes regulating the extent of autoimmune reactions. Gene copy number was quantitated by objective densitometry after PFGE was used to separate heterozygous AHs of different lengths. Initially we analyzed examples of known gene copy number at the C4 and 21 hydroxylase loci and showed that the approach provides predictable results. We then studied heterozygotes containing one characterized and one uncharacterized AH with particular attention to the gene copy number at the C4, Cyp21, and DRB loci. Each AH studied has a characteristic gene copy number at each locus studied. The same may be true of TNF, but other possibilities must be considered. AHs are markers for extensive chromosomal segments including particular numbers of several functional genes. Since AHs mark susceptibility to autoimmune disease, differences in gene copy number may be implicated.

Organisms have evolved various mechanisms to cope with the differences in the gene copy numbers between sexes caused by degeneration of Y and W sex chromosomes. Complete dosage compensation or at least expression balance between sexes has been reported predominantly in XX/XY systems, but rarely in ZZ/ZW systems. However, this often-reported pattern is based on comparisons of lineages where sex chromosomes evolved from nonhomologous genomic regions, potentially differing in sensitivity to differences in gene copy numbers. Here we document that two reptilian lineages (XX/XY iguanas and ZZ/ZW softshell turtles), which independently co-opted the same ancestral genomic region for the function of sex chromosomes, evolved different gene dose regulatory mechanisms. The independent co-option of the same genomic region for the role of sex chromosomes as in the iguanas and the softshell turtles offers great opportunity for testing evolutionary scenarios on sex chromosome evolution under the explicit control of the genomic background and gene identity. We show that the parallel loss of functional genes from the Y chromosome of the green anole and the W chromosome of the Florida softshell turtle led to different dosage compensation mechanisms. Our approach controlling for genetic background thus does not support that the variability in the regulation of gene dose differences is a consequence of ancestral autosomal gene content.

Copy number variations (CNVs) in the human genome are conventionally detected using high-throughput scanning technologies, such as comparative genomic hybridization and high-density single nucleotide polymorphism (SNP) microarrays, or relatively low-throughput techniques, such as quantitative polymerase chain reaction (PCR). All these approaches are limited in resolution and can at best distinguish a twofold (or 50%) difference in copy number. We have developed a new technology to study copy numbers using a platform known as the digital array, a nanofluidic biochip capable of accurately quantitating genes of interest in DNA samples. We have evaluated the digital array's performance using a model system, to show that this technology is exquisitely sensitive, capable of differentiating as little as a 15% difference in gene copy number (or between 6 and 7 copies of a target gene). We have also analyzed commercial DNA samples for their CYP2D6 copy numbers and confirmed that our results were consistent with those obtained independently using conventional techniques. In a screening experiment with breast cancer and normal DNA samples, the ERBB2 gene was found to be amplified in about 35% of breast cancer samples. The use of the digital array enables accurate measurement of gene copy numbers and is of significant value in CNV studies.

With the advent of real-time polymerase chain reaction (PCR), it is now possible to measure nucleic acid concentrations with an accuracy that was not possible only a few years ago. Examples are the analysis of gene expression or gene duplications/losses, where twofold differences in nucleic acid concentration have routinely been determined with almost 100% accuracy. As our primary interest is in prenatal diagnosis, we have investigated whether real-time PCR could be used for the diagnosis of chromosomal anomalies, in particular the aneuploidies such as trisomy 21, where the difference in copy number is only 50%. The feasibility of such an approach was first tested in a pilot study, in which we were able to demonstrate that trisomy 21 samples could be detected with 100% specificity. We have recently modified this test in order to permit the simultaneous analysis of trisomies 18 and 21, and have in a large scale analysis demonstrated that our approach can be used for the highly reproducible and robust detection of only 1.5-fold differences in gene copy number.

BackgroundCopy number variations (CNVs) of the Plasmodium falciparum multidrug resistance 1 (pfmdr1), P. falciparum plasmepsin2 (pfplasmepsin2) and P. falciparum GTP cyclohydrolase 1 (pfgch1) genes are associated with anti-malarial drug resistance in P. falciparum malaria. Droplet digital PCR (ddPCR) assays have been developed for accurate assessment of CNVs in several human genes. The aim of the present study was to develop and validate ddPCR assays for detection of the CNVs of P. falciparum genes associated with resistance to anti-malarial drugs.MethodsA multiplex ddPCR assay was developed to detect the CNVs in the pfmdr1 and pfplasmepsin2 genes, while a duplex ddPCR assay was developed to detect CNV in the pfgch1 gene. The gene copy number (GCN) quantification limit, as well as the accuracy and precision of the ddPCR assays were determined and compared to conventional quantitative PCR (qPCR). In order to reduce the cost of testing, a multiplex ddPCR assay of two target genes, pfmdr1 and pfplasmepsin2, was validated. In addition, the CNVs of genes of field samples collected from Thailand from 2015 to 2019 (n = 84) were assessed by ddPCR and results were compared to qPCR as the reference assay.ResultsThere were no significant differences between the GCN results obtained from uniplex and multiplex ddPCR assays for detection of CNVs in the pfmdr1 and pfplasmepsin2 genes (p = 0.363 and 0.330, respectively). Based on the obtained gene copy number quantification limit, the accuracy and percent relative standard deviation (%RSD) value of the multiplex ddPCR assay were 95% and 5%, respectively, for detection of the CNV of the pfmdr1 gene, and 91% and 5% for detection of the CNV of the pfplasmepsin2 gene. There was no significant difference in gene copy numbers assessed by uniplex or duplex ddPCR assays regarding CNV in the pfgch1 gene (p = 0.276). The accuracy and %RSD value of the duplex ddPCR assay were 95% and 4%, respectively, regarding pfgch1 GCN. In the P. falciparum field samples, pfmdr1 and pfplasmepsin2 GCNs were amplified in 15% and 27% of samples from Ubon Ratchathani, Thailand, while pfgch1 GCN was amplified in 50% of samples from Yala, Thailand. There was 100% agreement between the GCN results obtained from the ddPCR and qPCR assays (κ = 1.00). The results suggested that multiplex ddPCR assay is the optional assay for the accurate detection of gene copy number without requiring calibration standards, while the cost and required time are reduced. Based on the results of this study, criteria for GCN detection by ddPCR analysis were generated.ConclusionsThe developed ddPCR assays are simple, accurate, precise and cost-effective tools for detection of the CNVs in the pfmdr1, pfplasmepsin2 and pfgch1 genes of P. falciparum. The ddPCR assay is a useful additional tool for the surveillance of anti-malarial drug resistance.

Most plant mitochondrial genomes exist as subgenomic-size fragments apparently due to recombination between repetitive sequences. This leads to the possibility that independently replicating subgenomic domains could result in mitochondrial gene copy number variation. We show, through Southern-blot analysis of both restricted and intact mtDNA, that there are gene-specific copy number differences in the monocot Zea mays. Comparison of two different maize genotypes, B37(N) and B37(T), a cytoplasmic male-sterile strain, reveal fewer gene copy number differences for B37(T) than for B37(N). In contrast to maize, significant gene copy number differences are not detected in the dicot Brassica hirta. We also demonstrate that mitochondrial transcriptional rates in both species are apparently dependent on gene copy number since relative rates determined by run-on analysis are proportional to relative gene copy numbers. Thus a direct relationship exists between plant mitochondrial gene copy number and transcriptional rate.

Protozoan parasites of the genus Leishmania adapt to environmental change through chromosome and gene copy number variations. Only little is known about external or intrinsic factors that govern Leishmania genomic adaptation. Here, by conducting longitudinal genome analyses of 10 new Leishmania clinical isolates, we uncovered important differences in gene copy number among genetically highly related strains and revealed gain and loss of gene copies as potential drivers of long-term environmental adaptation in the field. In contrast, chromosome rather than gene amplification was associated with short-term environmental adaptation to in vitro culture. Karyotypic solutions were highly reproducible but unique for a given strain, suggesting that chromosome amplification is under positive selection and dependent on species- and strain-specific intrinsic factors. We revealed a progressive increase in read depth towards the chromosome ends for various Leishmania isolates, which may represent a nonclassical mechanism of telomere maintenance that can preserve integrity of chromosome ends during selection for fast in vitro growth. Together our data draw a complex picture of Leishmania genomic adaptation in the field and in culture, which is driven by a combination of intrinsic genetic factors that generate strain-specific phenotypic variations, which are under environmental selection and allow for fitness gain.IMPORTANCE Protozoan parasites of the genus Leishmania cause severe human and veterinary diseases worldwide, termed leishmaniases. A hallmark of Leishmania biology is its capacity to adapt to a variety of unpredictable fluctuations inside its human host, notably pharmacological interventions, thus, causing drug resistance. Here we investigated mechanisms of environmental adaptation using a comparative genomics approach by sequencing 10 new clinical isolates of the L. donovani, L. major, and L. tropica complexes that were sampled across eight distinct geographical regions. Our data provide new evidence that parasites adapt to environmental change in the field and in culture through a combination of chromosome and gene amplification that likely causes phenotypic variation and drives parasite fitness gains in response to environmental constraints. This novel form of gene expression regulation through genomic change compensates for the absence of classical transcriptional control in these early-branching eukaryotes and opens new venues for biomarker discovery.

In order to maximize recombinant protein expression in mammalian cells many factors need to be considered such as transfection method, vector construction, screening techniques and culture conditions. In addition, the host cell line can have a profound effect on the protein expression. However, auditioning or directly comparing host cell lines for optimal protein expression may be difficult since most transfection methods are based on random integration of the gene of interest into the host cell genome. Thus it is not possible to determine whether differences in expression between various host cell lines are due to the phenotype of the host cell itself or genetic factors such as gene copy number or gene location. To improve cell line generation, the ACE System was developed based on pre-engineered artificial chromosomes with multiple recombination acceptor sites. This system allows for targeted transfection and has been effectively used to rapidly generate stable CHO cell lines expressing high levels of monoclonal antibody. A key feature of the ACE System is the ability to isolate and purify ACEs containing the gene(s) of interest and transfect the same ACEs into different host cell lines. This feature allows the direct auditioning of host cells since the host cells have been transfected with ACEs that contain the same number of gene copies in the same genetic environment. To investigate this audition feature, three CHO host cell lines (CHOK1SV, CHO-S and DG44) were transfected with the same ACE containing gene copies of a human monoclonal IgG1 antibody. Clonal cell lines were generated allowing a direct comparison of antibody expression and stability between the CHO host cells. Results showed that the CHOK1SV host cell line expressed antibody at levels of more than two to five times that for DG44 and CHO-S host cell lines, respectively. To confirm that the ACE itself was not responsible for the low antibody expression seen in the CHO-S based clones, the ACE was isolated and purified from these cells and transfected back into fresh CHOK1SV cells. The resulting expression of the antibody from the ACE newly transfected into CHOK1SV increased fivefold compared to its expression in CHO-S and confirmed that the differences in expression between the different CHO host cells was due to the cell phenotype rather than differences in gene copy number and/or location. These results demonstrate the utility of the ACE System in providing a rapid and direct technique for auditioning host cell lines for optimal recombinant protein expression.

Polar bear (Ursus maritimus) and brown bear (Ursus arctos) are recently diverged species that inhabit vastly differing habitats. Thus, analysis of the polar bear and brown bear genomes represents a unique opportunity to investigate the evolutionary mechanisms and genetic underpinnings of rapid ecological adaptation in mammals. Copy number (CN) differences in genomic regions between closely related species can underlie adaptive phenotypes and this form of genetic variation has not been explored in the context of polar bear evolution. Here, we analyzed the CN profiles of 17 polar bears, 9 brown bears, and 2 black bears (Ursus americanus). We identified an average of 318 genes per individual that showed evidence of CN variation (CNV). Nearly 200 genes displayed species-specific CN differences between polar bear and brown bear species. Principal component analysis of gene CN provides strong evidence that CNV evolved rapidly in the polar bear lineage and mainly resulted in CN loss. Olfactory receptors composed 47% of CN differentiated genes, with the majority of these genes being at lower CN in the polar bear. Additionally, we found significantly fewer copies of several genes involved in fatty acid metabolism as well as AMY1B, the salivary amylase-encoding gene in the polar bear. These results suggest that natural selection shaped patterns of CNV in response to the transition from an omnivorous to primarily carnivorous diet during polar bear evolution. Our analyses of CNV shed light on the genomic underpinnings of ecological adaptation during polar bear evolution.

Despite the fact, that monoclonal antibodies are the fastest growing group of biopharmaceuticals in development, this is not true for the IgM class, which remains as enigmatic as ever. While more examples of usefulness of IgMs for medical applications are emerging, their recombinant production is still not common. In our study, stable monoclonal IgM producing CHO DG44 and HEK 293 cell lines, expressing two model IgM molecules (IgM-617 and IgM-012) were established. Recombinant cell lines were compared in regard of specific productivity, specific growth rate, maximal achieved antibody titer, gene copy numbers and transcription levels of transgene. IgM-617 cell lines were identified as high while IgM-012 clones were low producers. Although differences in gene copy numbers as well as in transcription levels were observed, they did not seem to be a limitation. Levels of relevant endoplasmic reticulum-stress related proteins were analyzed and no indications of unfolded protein response were detected. This could indicate that the difference in the intrinsic protein stability of our model proteins (as was previously observed on purified samples) might cause lower yields of IgM-012. Transcriptomics and/or proteomics follow up studies might be necessary for identification of potential bottlenecks in IgM producing cell lines.

Allelic imbalance is a common feature of many malignancies. We have measured allelic imbalance in genomic DNA from the breast cancer cell lines T47D, MDA-MB-231, two antiestrogen sensitive (MCF7N and MCF7L) and two resistant MCF7 cell lines (MMU2 and LCC9) using single nucleotide polymorphism (SNP) oligonucleotide microarrays. DNA from MCF7(L) and MMU2 cells was also analysed by comparative genome hybridisation (CGH) to compare with SNP microarray data. Proteins previously determined to be involved in disease progression were quantified by Western blot and compared to array data. The SNP and CGH array both detected cytogenetic abnormalities commonly found in breast cancer: amplification of chromosomes 11q13-14.1, 17q and 20q containing cyclin D1, BCAS1 and 3 (Breast Cancer Amplified Sequence) and AIB1 (Amplified in Breast cancer) genes; losses at 6q, 9p and X chromosomes, which included ERalpha (Estrogen Receptor alpha) and p16 ( INK4A ) genes. However the SNP chip array data additionally identified regions of loss of heterozygosity (LOH) followed by duplication of the remaining allele-uniparental disomy (UPD). Good concordance between SNP arrays and CGH analyses was observed, however there was poor correlation between gene copy number and protein levels between the cell lines. There were reductions in ERalpha, cyclin D1 and p27 protein levels whilst p21 protein levels were elevated in antiestrogen resistant MCF7 cell lines. Although protein levels varied there was no difference in gene copy number. This study shows SNP and CGH array analysis are powerful tools for analysis of allelic imbalance in breast cancer. However, the antiestrogen resistant phenotype was likely to be due to changes in gene expression and protein degradation rather than in altered gene copy number.

Independent origins of sociality in bees and ants are associated with independent expansions of particular odorant receptor (OR) gene subfamilies. In ants, one clade within the OR gene family, the 9-exon subfamily, has dramatically expanded. These receptors detect cuticular hydrocarbons (CHCs), key social signaling molecules in insects. It is unclear to what extent 9-exon OR subfamily expansion is associated with the independent evolution of sociality across Hymenoptera, warranting studies of taxa with independently derived social behavior. Here, we describe OR gene family evolution in the northern paper wasp, Polistes fuscatus, and compare it to four additional paper wasp species spanning ∼40 million years of evolutionary divergence. We find 200 putatively functional OR genes in P. fuscatus, matching predictions from neuroanatomy, and more than half of these are in the 9-exon subfamily. Most OR gene expansions are tandemly arrayed at orthologous loci in Polistes genomes, and microsynteny analysis shows species-specific gain and loss of 9-exon ORs within tandem arrays. There is evidence of episodic positive diversifying selection shaping ORs in expanded subfamilies. Values of omega (dN/dS) are higher among 9-exon ORs compared to other OR subfamilies. Within the Polistes OR gene tree, branches in the 9-exon OR clade experience relaxed negative (relaxed purifying) selection relative to other branches in the tree. Patterns of OR evolution within Polistes are consistent with 9-exon OR function in CHC perception by combinatorial coding, with both natural selection and neutral drift contributing to interspecies differences in gene copy number and sequence.

Optimized hematopoietic gene therapy requires vectors with strong expression in the desired target cell population and the ability to select for the expressing transduced cells. In the context of drug resistance selection of repopulating hematopoietic stem cells in the mouse, we examined tissue expression after transduced marrow transplantation of the drug selection gene, G156A mutant O6-methylguanine-DNA methyltransferase (G156A MGMT). To gain more experience with the rigor of the impact of selection on tissue-specific gene expression, we also asked whether there are expression differences between three different onco-retroviral backbones--MPSV, SF, and MFG. MGMT expression was compared after O6-benzylguanine (BG) and 1,3-bis (2-chloroethyl)-1-nitrosourea (BCNU) drug selection in vivo. After mice were transplanted with cells transduced with MPSV, MFG, or SF retroviral vectors expressing G156A MGMT and drug treated, nearly complete replacement by transduced progenitors was observed in the marrow. Each backbone supported MGMT expression in all four hematopoietic lineages in vivo indicating that MGMT-mediated selection is indeed robust. Expression in marrow, spleen, and thymus was very similar between the vectors and differences were most likely due to differences in gene copy number per selected cell. In primary and secondary recipients, the highest expression was observed in MFG and this was the vector that transduced at the greatest proviral copy number per cell. These data indicate that strong selection pressure using the MGMT gene to protect primary and secondary repopulating murine stem cells from the toxicity of BCNU. Regardless of the vector backbone used, multiorgan expression was observed without evidence of gene silencing. These data help establish mutant, BG-resistant MGMT as a potent selection gene for stem cell selection in vivo.

The genetic basis of schizophrenia is obscure. In an XX male patient with schizophrenia we previously showed that one X;Y translocation breakpoint was in pseudoautosomal region 1 (PAR1) with the effect that the proximal segment of PAR1 from the PAR1 boundary to acetylserotonin N-methyl transferase (ASMT) distally was triplicated in this patient. This study determined whether dosage imbalances of X-Y homologous regions in general are associated with schizophrenia. A multiplex semi-quantitative polymerase chain reaction assay was developed to quantify MIC2 gene as a representative of PAR1 and compare it with the SYBL1 gene which maps in pseudoautosomal region 2 (PAR2) and protocadherin XY (PCDHXY), located at Xq21.3. Each of these three loci was co-amplified with the autosomal gene MSX2 using Cy5-labelled primers and the products separated by electrophoresis in polyacrylamide gels. Results were expressed as ratios of peak area of the target gene to MSX2 which served as an internal dosage control. Using genomes with sex chromosome aneuploidies, the method was found sensitive enough to detect a two-fold difference in gene copy number. We confirmed the MIC2 triplication in the XX male patient but found no significant difference in gene dosage of MIC2, PCDHXY and SYBL1 in a panel of 17 patients with schizophrenia compared to controls. No evidence was obtained for gene dosage imbalances in MIC2, PCDHXY and SYBL1 in patients with schizophrenia.