Control of chromosomal rearrangements in {\it Escherichia coli\/}

A complete physical map of Escherichia coli K-12 strain MG1655 was constructed by digesting chromosomal DNA with the infrequently cutting restriction enzymes NotI, SfiI and XbaI and separating the fragments by pulsed field gel electrophoresis. The map was used to compare six K-12 strains of E. coli. Although several differences were noted and localized, the map of MG1655 was representative of all the K-12 strains tested. The maps were also used to analyze chromosomal rearrangements in the E. coli strain MG1655. The spontaneous and UV induced frequencies of tandem duplication formation were measured at several loci distributed around the chromosome. The spontaneous duplication frequency varied from 10$\sp{-5}$ to 10$\sp{-3}$ and increased at least ten-fold following mild UV irradiation treatment. Duplications of several regions of the chromosome, including the serA region and the metE region, were mapped using pulsed field gel electrophoresis. Duplications of serA were found to be large, ranging in size from 600 kb to 2100 kb. Several of the duplications isolated at serA were caused by ectopic recombination between IS5 elements and between IS186 elements. Duplications of the metE region, however, were almost exclusively the result of ectopic recombination between ribosomal RNA cistrons. Duplication frequencies were determined at both serA and metE in wild type and mismatch repair mutant strains (mutL, mutS, uvrD and recF). Even though all of the mismatch repair mutations increased duplication frequency of metE, the largest increases were observed in the mutL and mutS strains. Duplication frequency of serA was increased less dramatically by mutations in mismatch repair. Several duplications of metE isolated in a wild type and a mismatch repair mutant were mapped. The results showed that the same repeated sequences were used for duplication formation in the mismatch repair mutant as were used in the wild type strain. Several isolates showed evidence of multiple rearrangements indicating that mismatch repair may play a role in stabilizing the genome by controlling chromosomal rearrangement. ^

Origins and evolution of VNTR loci: The apolipoprotein B 3$\sp\prime$ VNTR

I studied the apolipoprotein (apo) B 3$\sp\prime$ variable number tandem repeat (VNTR) and did computer simulations of the stepwise mutation model to address four questions: (1) How did the apo B VNTR originate? (2) What is the mutational mechanism of repeat number change at the apo B VNTR? (3) To what extent are population and molecular level events responsible for the determination of the contemporary apo B allele frequency distribution? (4) Can VNTR allele frequency distributions be explained by a simple and conservative mutation-drift model? I used three general approaches to address these questions: (1) I characterized the apo B VNTR region in non-human primate species; (2) I constructed haplotypes of polymorphic markers flanking the apo B VNTR in a sample of individuals from Lorrain, France and studied the associations between the flanking-marker haplotypes and apo B VNTR size; (3) I did computer simulations of the one-step stepwise mutation model and compared the results to real data in terms of four allele frequency distribution characteristics.^ The results of this work have allowed me to conclude that the apo B VNTR originated after an initial duplication of a sequence which is still present as a single copy sequence in New World monkey species. I conclude that this locus did not originate by the transposition of an array of repeats from somewhere else in the genome. It is unlikely that recombination is the primary mutational mechanism. Furthermore, the clustered nature of these associations implicates a stepwise mutational mechanism. From the high frequencies of certain haplotype-allele size combinations, it is evident that population level events have also been important in the determination of the apo B VNTR allele frequency distribution. Results from computer simulations of the one-step stepwise mutation model have allowed me to conclude that bimodal and multimodal allele frequency distributions are not unexpected at loci evolving via stepwise mutation mechanisms. Short tandem repeat loci fit the stepwise mutation model best, followed by microsatellite loci. I therefore conclude that there are differences in the mutational mechanisms of VNTR loci as classed by repeat unit size. (Abstract shortened by UMI.) ^

CAOR: Context-aware Adaptive Opportunistic Routing in Mobile Ad-hoc Networks

Opportunistic routing (OR) employs a list of candidates to improve wireless transmission reliability. However, conventional list-based OR restricts the freedom of opportunism, since only the listed nodes are allowed to compete for packet forwarding. Additionally, the list is generated statically based on a single network metric prior to data transmission, which is not appropriate for mobile ad-hoc networks (MANETs). In this paper, we propose a novel OR protocol - Context-aware Adaptive Opportunistic Routing (CAOR) for MANETs. CAOR abandons the idea of candidate list and it allows all qualified nodes to participate in packet transmission. CAOR forwards packets by simultaneously exploiting multiple cross-layer context information, such as link quality, geographic progress, energy, and mobility.With the help of the Analytic Hierarchy Process theory, CAOR adjusts the weights of context information based on their instantaneous values to adapt the protocol behavior at run-time. Moreover, CAOR uses an active suppression mechanism to reduce packet duplication. Simulation results show that CAOR can provide efficient routing in highly mobile environments. The adaptivity feature of CAOR is also validated.

Independent polled mutations leading to complex gene expression differences in cattle

The molecular regulation of horn growth in ruminants is still poorly understood. To investigate this process, we collected 1019 hornless (polled) animals from different cattle breeds. High-density SNP genotyping confirmed the presence of two different polled associated haplotypes in Simmental and Holstein cattle co-localized on BTA 1. We refined the critical region of the Simmental polled mutation to 212 kb and identified an overlapping region of 932 kb containing the Holstein polled mutation. Subsequently, whole genome sequencing of polled Simmental and Holstein cows was used to determine polled associated genomic variants. By genotyping larger cohorts of animals with known horn status we found a single perfectly associated insertion/deletion variant in Simmental and other beef cattle confirming the recently published possible Celtic polled mutation. We identified a total of 182 sequence variants as candidate mutations for polledness in Holstein cattle, including an 80 kb genomic duplication and three SNPs reported before. For the first time we showed that hornless cattle with scurs are obligate heterozygous for one of the polled mutations. This is in contrast to published complex inheritance models for the bovine scurs phenotype. Studying differential expression of the annotated genes and loci within the mapped region on BTA 1 revealed a locus (LOC100848215), known in cow and buffalo only, which is higher expressed in fetal tissue of wildtype horn buds compared to tissue of polled fetuses. This implicates that the presence of this long noncoding RNA is a prerequisite for horn bud formation. In addition, both transcripts associated with polledness in goat and sheep (FOXL2 and RXFP2), show an overexpression in horn buds confirming their importance during horn development in cattle.

Delayed Haematological recovery after autologous stem cell transplantation is associated with favourable outcome in acute myeloid leukaemia.

Autologous stem cell transplantation (ASCT) is applied to consolidate first remission in patients with acute myeloid leukaemia (AML). However, outcome after ASCT widely varies among AML patients. We analyzed the prognostic significance of haematological recovery for neutrophils [absolute neutrophil count (ANC) >1·0 × 10(9) /l] and platelets (platelet count >20·0 × 10(9) /l), stratifying at day 20 after ASCT in 88 consecutive and homogeneously treated AML patients in first remission. We observed that patients with delayed recovery had better overall survival (OS; ANC: P < 0·0001 and platelets: P = 0·0062) and time to progression (TTP; ANC: P = 0·0003 and platelets: P = 0·0125). Delayed recovery was an independent marker for better OS and TTP in a multivariate analysis including age, gender, number of transfused CD34+ cells, cytogenetics, FLT3-internal tandem duplication and NPM1 mutation. Our results suggest that delayed neutrophil and platelet recovery is associated with longer OS and TTP in AML patients consolidated with ASCT in first remission.

Congenital aural atresia associated with agenesis of internal carotid artery in a girl with a FOXI3 deletion.

We report on the molecular characterization of a microdeletion of approximately 2.5 Mb at 2p11.2 in a female baby with left congenital aural atresia, microtia, and ipsilateral internal carotid artery agenesis. The deletion was characterized by fluorescence in situ hybridization, array comparative genomic hybridization, and whole genome re-sequencing. Among the genes present in the deleted region, we focused our attention on the FOXI3 gene. Foxi3 is a member of the Foxi class of Forkhead transcription factors. In mouse, chicken and zebrafish Foxi3 homologues are expressed in the ectoderm and endoderm giving rise to elements of the jaw as well as external, middle and inner ear. Homozygous Foxi3-/- mice have recently been generated and show a complete absence of the inner, middle, and external ears as well as severe defects in the jaw and palate. Recently, a 7-bp duplication within exon 1 of FOXI3 that produces a frameshift and a premature stop codon was found in hairless dogs. Mild malformations of the outer auditory canal (closed ear canal) and ear lobe have also been noted in a fraction of FOXI3 heterozygote Peruvian hairless dogs. Based on the phenotypes of Foxi3 mutant animals, we propose that FOXI3 may be responsible for the phenotypic features of our patient. Further characterization of the genomic region and the analysis of similar patients may help to demonstrate this point. © 2015 Wiley Periodicals, Inc.

16p11.2 600 kb Duplications confer risk for typical and atypical Rolandic epilepsy

Rolandic epilepsy (RE) is the most common idiopathic focal childhood epilepsy. Its molecular basis is largely unknown and a complex genetic etiology is assumed in the majority of affected individuals. The present study tested whether six large recurrent copy number variants at 1q21, 15q11.2, 15q13.3, 16p11.2, 16p13.11 and 22q11.2 previously associated with neurodevelopmental disorders also increase risk of RE. Our association analyses revealed a significant excess of the 600 kb genomic duplication at the 16p11.2 locus (chr16: 29.5-30.1 Mb) in 393 unrelated patients with typical (n = 339) and atypical (ARE; n = 54) RE compared with the prevalence in 65,046 European population controls (5/393 cases versus 32/65,046 controls; Fisher's exact test P = 2.83 × 10(-6), odds ratio = 26.2, 95% confidence interval: 7.9-68.2). In contrast, the 16p11.2 duplication was not detected in 1738 European epilepsy patients with either temporal lobe epilepsy (n = 330) and genetic generalized epilepsies (n = 1408), suggesting a selective enrichment of the 16p11.2 duplication in idiopathic focal childhood epilepsies (Fisher's exact test P = 2.1 × 10(-4)). In a subsequent screen among children carrying the 16p11.2 600 kb rearrangement we identified three patients with RE-spectrum epilepsies in 117 duplication carriers (2.6%) but none in 202 carriers of the reciprocal deletion. Our results suggest that the 16p11.2 duplication represents a significant genetic risk factor for typical and atypical RE.

Rapid and accurate G M1-gangliosidosis diagnosis using a parentage testing microsatellite

The G M1-gangliosidosis is an autosomal recessive lysosomal storage disease caused by structural defects of the beta-galactosidase gene (GLB1) which lead to a severe phenotypical impairment in homozygous individuals, whereas heterozygous carriers remain clinically normal. Currently employed DNA parentage tests include the analysis of microsatellites, which also have a diagnostic predictive value. The aim of this study was to provide a reliable tool for genotyping the canine GLB1 which can be effectively integrated in parentage testing investigations. For this purpose the association between the GLB1 gene and the AHT K253 microsatellite was analyzed in 30 Alaskan huskies (11 GLB1+/+, 17 GLB1+/- and 2 GLB1-/- dogs). The 143 bp AHT K253 microsatellite allele was identified only in GLB1+/- and GLB1-/- animals and was in strong linkage disequilibrium with the causative mutation for G M1-gangliosidosis, a 19 bp duplication within exon 15 of the GLB1 gene. The results of the present study revealed a 100% concordance between the previous established genotypes and those obtained after the analysis of the AHT K253 microsatellite. Thus, the genotype of the AHT K253 microsatellite, which is routinely determined during dog parentage testing, has a high predictive value for the G M1-gangliosidosis carrier status.

Global diversity, population stratification and selection of human copy number variation

In order to explore the diversity and selective signatures of duplication and deletion human copy number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single nucleotide variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load.

EXPRESSION OF THE HUMAN PLACENTAL LACTOGEN GENES

The placenta is the site of synthesis of various peptide and steroid hormones related to pregnancy. Human placental lactogen (hPL) is the predominant peptide hormone secreted by term placenta and its synthesis is tissue-specific and coupled to placenta development. The objective of this work was to study the structure and expression of the hPL.^ Poly(A('+))RNA from human term placenta was translated in a mouse-derived cell-free system. A major band corresponding to pre-hPL and a minor band comigrating with mature hPL, represent (TURN)15% of the total radioactively labeled proteins. Analysis of the poly(A('+))RNA showed a prominent band at approximately 860 nucleotides. A corresponding band was observed in Northern blots of total RNA, hybridized with {('32)P}-labeled recombinant plasmid containing a portion of hPL cDNA. Similar analyses of nuclear RNA showed at least four additional bands at 990, 1200, 1460 and 1760 nucleotides, respectively, which are likely precursors of hPL mRNA. Poly(A('+))RNA was used to construct a cDNA library, of which approximately 5% of the clones were found to hybridize to hPL DNA sequences. Heteroduplexes constructed between a clone containing a 815 bp hPL cDNA insert and a hPL genomic DNA clone revealed four small intervening sequences which can account for the lengths observed in hnRNA molecules.^ Recombinant plasmid HCS-pBR322 containing a 550 bp insert of a cDNA transcript of human placental lactogen (hPL) mRNA was ('3)H-labeled an hybridized in situ to human chromosome preparations. These experiments allowed assignment of the hPL and growth hormone (hGH) genes, which have over 90% nucleotide homology in their coding sequences, to band q22-24 of chromosome 17. A gene copy number experiment showed that both genes are present in (TURN)3 copies per haploid genome.^ Experiments were designed to determine if all members of the hPL gene cluster, consisting of four non-allelic genes, are transcribed in term placenta. Advantage was taken of differences in restriction endonuclease sites in the coding portions of the different hPL genes, to distinguish the putative cDNAs of the transcriptionally active genes. Two genes were found to be represented in the cDNA library and their cDNA transcripts were isolated and characterized. Three independent methods showed that their corresponding mRNAs are about equally represented in the hPL mRNA population. The two cDNAs code for prehPL proteins which differ at a single amino acid position. However the secreted hPLs have identical amino acid sequences. A tetramer insertion duplication was found in a palindrome area of the 3' untranslated region of one of the hPL mRNAs. ^

The p120-catenin/Kaiso signaling pathway

The canonical and non-canonical Wnt signaling pathways appear to interact with one another as a network in development, or when hyper-activated, in the progression of disease. A much studied key mediator of the canonical Wnt pathway, β-catenin, is characterized by a central armadillo-repeat domain that engages in multiple protein-protein interactions, such as those with cadherins functioning at cell-cell contact regions. In the nucleus, β-catenin forms a complex with the repressor TCF/LEF, promoting the activation of genes participating in processes such as proliferation, differentiation and stem cell survival. Somewhat similarly, the p120-catenin binds the distinct transcriptional repressor Kaiso, relieving Kaiso-mediated repression to promote gene activation. Here, employing Xenopus laevis, I report upon both downstream and upstream aspects of the p120-catenin/Kaiso pathway which was previously poorly understood. I first show that Kaiso, a BTB/POZ zinc-finger family member, directly represses canonical Wnt gene targets (Siamois, c-Fos, Cyclin-D1 and c-Myc) in conjunction with TCF. Depletion or dominant-negative inhibition of xKaiso results in Siamois de-repression, while xKaiso over-expression induces additional Siamois repression through recruitment of N-CoR co-repressor and chromatin modifications. Functional interdependencies are further corroborated by the capacity of Kaiso to suppress β-catenin-induced axis duplication. Thus, my work inter-relates the p120-catenin/Kaiso and β-catenin/TCF pathways at the level of specific gene promoters important in development and cancer progression. Regarding upstream aspects of the p120-catenin/Kaiso pathway, I collaboratively identified p120 in association with Frodo, a protein previously identified as a component of the canonical (β-catenin dependent) Wnt pathway. I determined that canonical Wnt signals result in Frodo-mediated stabilization of p120-catenin, resulting in the sequestration of Kaiso to the cytoplasm and thereby the activation (relief of repression) of gene targets. Developmental evidence supporting this view included findings that Frodo has the capacity to partially rescue Kaiso over-expression phenotypes in early Xenopus embryos. Taken together, my studies point to the convergence of p120-catenin/Kaiso and β-catenin/TCF signaling pathways at the level of gene transcription as well as at more upstream points during vertebrate development. ^

Coordination of murine limb muscle, skeleton and connective tissue development by Lmx1b

The underlying genetic defects of a congenital disease Nail-Patella Syndrome are loss-of-function mutations in the LMX1B gene. Lmx1b encodes a LIM-homeodomain transcription factor that is expressed specifically in the dorsal limb bud mesenchyme. Gain- and loss-of-function experiments suggest that Lmx1b is both necessary and sufficient to specify dorsal limb patterning. However, how Lmx1b coordinates patterning of the dorsal tissues in the limb, including muscle, skeleton and connective tissues, remains unknown. One possibility is that each tissue specifies its own pattern cell-autonomously, i.e., Lmx1b is expressed in tissues in which it functions and different tissues do not communicate with each other. Another possibility is that tissues that express Lmx1b interact with adjacent tissues and provide patterning information thereby directing the development of tissues non-cell-autonomously. Previous results showed that Lmx1b is expressed in limb connective tissue and skeleton, but is not expressed in muscle tissue. Moreover, muscles and muscle connective tissue are closely associated during development. Therefore, we hypothesize that Lmx1b controls limb muscle dorsal-ventral (DV) patterning through muscle connective tissue, but regulates skeleton and tendon/ligament development cell-autonomously. ^ To test this hypothesis, we first examined when and where the limb dorsal-ventral asymmetry is established during development. Subsequently, conditional knockout and overexpression experiments were performed to delete or activate Lmx1b in different tissues within the limb. Our results show that deletion of Lmx1b from whole limb mesenchyme results in all dorsal tissues, including muscle, tendon/ligament and skeleton, transforming into ventral structures. Skeleton-specific knockout of Lmx1b led to the dorsal duplication of distal sesamoid and metacarpal bones, but did not affect the pattern formation of other tissues, suggesting that Lmx1b controls skeleton development cell-autonomously. In addition, this skeleton-specific pattern alteration only occurs in distal limb tissues, not proximal limb tissues, indicating different regulatory mechanisms operate along the limb proximal-distal axis. Moreover, skeleton-specific ectopic expression of Lmx1b reveals a complementary skeletal-specific dorsalized phenotype. This result supports a cell-autonomous role for Lmx1b in dorsal-ventral skeletal patterning. This study enriched our understanding of limb development, and the insights from this research may also be applicable for the development of other organs. ^

Effects of congressional modernization laws on public access to innovative medical device technologies

The Federal Food and Drug Administration (FDA) and the Centers for Medicare and Medicaid (CMS) play key roles in making Class III, medical devices available to the public, and they are required by law to meet statutory deadlines for applications under review. Historically, both agencies have failed to meet their respective statutory requirements. Since these failures affect patient access and may adversely impact public health, Congress has enacted several “modernization” laws. However, the effectiveness of these modernization laws has not been adequately studied or established for Class III medical devices. ^ The aim of this research study was, therefore, to analyze how these modernization laws may have affected public access to medical devices. Two questions were addressed: (1) How have the FDA modernization laws affected the time to approval for medical device premarket approval applications (PMAs)? (2) How has the CMS modernization law affected the time to approval for national coverage decisions (NCDs)? The data for this research study were collected from publicly available databases for the period January 1, 1995, through December 31, 2008. These dates were selected to ensure that a sufficient period of time was captured to measure pre- and post-modernization effects on time to approval. All records containing original PMAs were obtained from the FDA database, and all records containing NCDs were obtained from the CMS database. Source documents, including FDA premarket approval letters and CMS national coverage decision memoranda, were reviewed to obtain additional data not found in the search results. Analyses were conducted to determine the effects of the pre- and post-modernization laws on time to approval. Secondary analyses of FDA subcategories were conducted to uncover any causal factors that might explain differences in time to approval and to compare with the primary trends. The primary analysis showed that the FDA modernization laws of 1997 and 2002 initially reduced PMA time to approval; after the 2002 modernization law, the time to approval began increasing and continued to increase through December 2008. The non-combined, subcategory approval trends were similar to the primary analysis trends. The combined, subcategory analysis showed no clear trends with the exception of non-implantable devices, for which time to approval trended down after 1997. The CMS modernization law of 2003 reduced NCD time to approval, a trend that continued through December 2008. This study also showed that approximately 86% of PMA devices do not receive NCDs. ^ As a result of this research study, recommendations are offered to help resolve statutory non-compliance and access issues, as follows: (1) Authorities should examine underlying causal factors for the observed trends; (2) Process improvements should be made to better coordinate FDA and CMS activities to include sharing data, reducing duplication, and establishing clear criteria for “safe and effective” and “reasonable and necessary”; (3) A common identifier should be established to allow tracking and trending of applications between FDA and CMS databases; (4) Statutory requirements may need to be revised; and (5) An investigation should be undertaken to determine why NCDs are not issued for the majority of PMAs. Any process improvements should be made without creating additional safety risks and adversely impacting public health. Finally, additional studies are needed to fully characterize and better understand the trends identified in this research study.^

Molecular Mechanisms of Vascular Disease in Patients with Rare Variants in MYH11

Thoracic aortic aneurysms and dissections (TAAD) are the primary disease affecting the thoracic ascending aorta, with an incidence rate of 10.4/100,000. Although about 20% of patients carry a mutation in a single gene that causes their disease, the remaining 80% of patients may also have genetic factors that increase their risk for developing TAAD. Many of the genes that predispose to TAAD encode proteins involved in smooth muscle cell (SMC) contraction and the disease-causing mutations are predicted to disrupt contractile function. SMCs are the predominant cell type in the ascending aortic wall. Mutations in MYH11, encoding the smooth muscle specific myosin heavy chain, are a rare cause of inherited TAAD. However, rare but recurrent non-synonymous variants in MYH11 are present in the general population but do not cause inherited TAAD. The goal of this study was to assess the potential role of these rare variants in vascular diseases. Two distinct variants were selected: the most commonly seen rare variant, MYH11 R247C, and a duplication of the chromosomal region spanning the MYH11 locus at 16p13.1. Genetic analyses indicated that both of these variants were significantly enriched in patients with TAAD compared with controls. A knock-in mouse model of the Myh11 R247C rare variant was generated, and these mice survive and reproduce normally. They have no structural abnormalities of the aorta or signs of aortic disease, but do have decreased aortic contractility. Myh11R247C/R247C mice also have increased proliferative response to vascular injury in vivo and increased proliferation of SMCs in vitro. Myh11R247C/R247C SMCs have decreased contractile gene and protein expression and are dedifferentiated. In fibroblasts, myosin force generation is required for maturation of focal adhesions, and enhancers of RhoA activity replace enhancers of Rac1 activity as maturation occurs. Consistent with these previous findings, focal adhesions are smaller in Myh11R247C/R247C SMCs, and there is decreased RhoA activation. A RhoA activator (CN03) rescues the dedifferentiated phenotype of Myh11R247C/R247C SMCs. Myh11R247C/R247C mice were bred with an existing murine model of aneurysm formation, the Acta2-/- mouse. Over time, mice carrying the R247C allele in conjunction with heterozygous or homozygous loss of Acta2 had significantly increased aortic diameter, and a more rapid accumulation of pathologic markers. These results suggest that the Myh11 R247C rare variant acts as a modifier gene increasing the risk for and severity of TAAD in mice. In patients with 16p13.1 duplications, aortic MYH11 expression is increased, but there is no corresponding increase in smooth muscle myosin heavy chain protein. Using SMCs that overexpress Myh11, we identified alterations in SMC phenotype leading to excessive protein turnover. All contractile proteins, not just myosin, are affected, and the proteins are turned over by autophagic degradation. Surprisingly, these cells are also more contractile compared with wild-type SMCs. The results described in this dissertation firmly establish that rare variants in MYH11 significantly affect the phenotype of SMCs. Further, the data suggests that these rare variants do increase the risk of TAAD via pathways involving altered SMC phenotype and contraction. Therefore, this study validates that these rare genetic variants alter vascular SMCs and provides model systems to explore the contribution of rare variants to disease.

Molecular and evolutionary genetics of the X -linked visual pigment genes in humans and New World monkeys

Normal humans have one red and at least one green visual pigment genes. These genes are tightly linked as tandem repeats on the X chromosome and each of them has six exons. There is only one X-linked visual pigment gene in New World monkeys (NWMs) but the locus has three polymorphic alleles encoding red, yellow and green visual pigments, respectively. The spectral properties of the squirrel monkey and the marmoset (both NWMs) have been studied and partial sequences of the three alleles are available. To study the evolutionary history of these X-linked opsin genes in humans and NWMs, coding and intron sequences of the three squirrel monkey alleles and the three marmoset alleles were amplified by PCR followed by subcloning and sequencing. Introns 2 and 4 of the human red and green pigment genes were also sequenced. The results obtained are as follows: (1) The sequences of introns 2 and 4 of the human red and green opsin genes are significantly more similar between the two genes than are coding sequences, contrary to the usual situation where coding regions are better conserved in evolution than are introns. The high similarities in the two introns are probably due to recent gene conversion events during evolution of the human lineage. (2) Phylogenetic analysis of both intron and exon sequences indicates that the phylogenetic tree of the available primate opsin genes is the same as the species tree. The two human genes were derived from a gene duplication event after the divergence of the human and NWM lineages. The three alleles in each of the two NWM species diverged after the split of the two NWMs but have persisted in the population for at least 5 million years. (3) Allelic gene conversion might have occurred between the three squirrel monkey alleles. (4) A model of additive effect of hydroxyl-bearing amino acids on spectral tuning is proposed by treating some unknown variables as groups. Under the assumption that some residues have no effect, it is found that at least five amino acid residues, at positions 178 (3 nm), 180 (5 nm), 230 ($-$4 nm), 277 (9 nm) and 285 (13 nm), have linear spectral tuning effects. (5) Adaptive evolution of the opsin genes to different spectral peaks was observed at four residues that are important for spectral tuning. ^