Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5′-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

Isolation of DNA sequences with homology to a degenerate FaRP oligonucleotide from the crayfish Procambarus clarkii

The neuropeptide Th1RFamide with the sequence Phe-Met-Arg-Phe-amide was originally isolated in the clam Macrocallista nimbosa (price and Greenberg, 1977). Since its discovery, a large family ofFl\1RFamide-related peptides termed FaRPs have been found to be present in all major animal phyla with functions ranging from modulation of neuronal activity to alteration of muscular contractions. However, little is known about the genetics encoding these peptides, especially in invertebrates. As FaRP-encoding genes have yet to be investigated in the invertebrate Malacostracean subphylum, the isolation and characterization ofFaRP-encoding DNA and mRNA was pursued in this project. The immediate aims of this thesis were: (1) to amplify mRNA sequences of Procambarus clarkii using a degenerate oligonucleotide primer deduced from the common amino acid sequence ofisolated Procambarus FaRPS, (2) to determine if these amplification products encode FaRP gene sequences, and (3) to create a selective cDNA library of sequences recognized by the degenerate oligonucleotide primer. The polymerase chain reaction - rapid amplification of cDNA ends (PCR-RACE) is a procedure in which a single gene-specific primer is used in conjunction with a generalized 3' or 5' primer to amplify copies ofthe region between a single point in the transcript and the 3' or 5' end of cDNA of interest (Frohman et aI., 1988). PCRRACE reactions were optimized with respect to primers used, buffer composition, cycle number, nature ofgenetic substrate to be amplified, annealing, extension and denaturation temperatures and times, and use of reamplification procedures. Amplification products were cloned into plasmid vectors and recombinant products were isolated, as were the recombinant plaques formed in the selective cDNA library. Labeled amplification products were hybridized to recombinant bacteriophage to determine ligated amplification product presence. When sequenced, the five isolated PCR-RACE amplification products were determined not to possess FaRP-encoding sequences. The 200bp, 450bp, and 1500bp sequences showed homology to the Caenorhabditis elegans cosmid K09A11, which encodes for cytochrome P450; transfer-RNA; transposase; and tRNA-Tyr, while the 500bp and 750bp sequences showed homology with the complete genome of the Vaccinia virus. Under the employed amplification conditions the degenerate oligonucleotide primer was observed to bind to and to amplify sequences with either 9 or 10bp of 17bp identity. The selective cDNA library was obselVed to be of extremely low titre. When library titre was increased, white. plaques were isolated. Amplification analysis of eight isolated Agt11 sequences from these plaques indicated an absence of an insertion sequence. The degenerate 17 base oligonucleotide primer synthesized from the common amino acid sequence ofisolated Procambarus FaRPs was thus determined to be non-specific in its binding under the conditions required for its use, and to be insufficient for the isolation and identification ofFaRP-encoding sequences. A more specific primer oflonger sequence, lower degeneracy, and higher melting temperature (TJ is recommended for further investigation into the FaRP-encoding genes of Procambarlls clarkii.

Regulation of gene expression in the dinoflagellate Lingulodinium polyedrum

Les dinoflagellés sont des eucaryotes unicellulaires que l’on retrouve autant en eau douce qu’en milieu marin. Ils sont particulièrement connus pour causer des fleurs d’algues toxiques nommées ‘marée-rouge’, ainsi que pour leur symbiose avec les coraux et pour leur importante contribution à la fixation du carbone dans les océans. Au point de vue moléculaire, ils sont aussi connus pour leur caractéristiques nucléaires uniques, car on retrouve généralement une quantité immense d’ADN dans leurs chromosomes et ceux-ci sont empaquetés et condensés sous une forme cristalline liquide au lieu de nucléosomes. Les gènes encodés par le noyau sont souvent présents en multiples copies et arrangés en tandem et aucun élément de régulation transcriptionnelle, y compris la boite TATA, n’a encore été observé. L’organisation unique de la chromatine des dinoflagellés suggère que différentes stratégies sont nécessaires pour contrôler l’expression des gènes de ces organismes. Dans cette étude, j’ai abordé ce problème en utilisant le dinoflagellé photosynthétique Lingulodinium polyedrum comme modèle. L. polyedrum est d’un intérêt particulier, car il a plusieurs rythmes circadiens (journalier). À ce jour, toutes les études sur l’expression des gènes lors des changements circadiens ont démontrées une régulation à un niveau traductionnel. Pour mes recherches, j’ai utilisé les approches transcriptomique, protéomique et phosphoprotéomique ainsi que des études biochimiques pour donner un aperçu de la mécanique de la régulation des gènes des dinoflagellés, ceci en mettant l’accent sur l’importance de la phosphorylation du système circadien de L. polyedrum. L’absence des protéines histones et des nucléosomes est une particularité des dinoflagellés. En utilisant la technologie RNA-Seq, j’ai trouvé des séquences complètes encodant des histones et des enzymes modifiant les histones. L polyedrum exprime donc des séquences conservées codantes pour les histones, mais le niveau d’expression protéique est plus faible que les limites de détection par immunodétection de type Western. Les données de séquençage RNA-Seq ont également été utilisées pour générer un transcriptome, qui est une liste des gènes exprimés par L. polyedrum. Une recherche par homologie de séquences a d’abord été effectuée pour classifier les transcrits en diverses catégories (Gene Ontology; GO). Cette analyse a révélé une faible abondance des facteurs de transcription et une surprenante prédominance, parmi ceux-ci, des séquences à domaine Cold Shock. Chez L. polyedrum, plusieurs gènes sont répétés en tandem. Un alignement des séquences obtenues par RNA-Seq avec les copies génomiques de gènes organisés en tandem a été réalisé pour examiner la présence de transcrits polycistroniques, une hypothèse formulée pour expliquer le manque d’élément promoteur dans la région intergénique de la séquence de ces gènes. Cette analyse a également démontré une très haute conservation des séquences codantes des gènes organisés en tandem. Le transcriptome a également été utilisé pour aider à l’identification de protéines après leur séquençage par spectrométrie de masse, et une fraction enrichie en phosphoprotéines a été déterminée comme particulièrement bien adapté aux approches d’analyse à haut débit. La comparaison des phosphoprotéomes provenant de deux périodes différentes de la journée a révélée qu’une grande partie des protéines pour lesquelles l’état de phosphorylation varie avec le temps est reliées aux catégories de liaison à l’ARN et de la traduction. Le transcriptome a aussi été utilisé pour définir le spectre des kinases présentes chez L. polyedrum, qui a ensuite été utilisé pour classifier les différents peptides phosphorylés qui sont potentiellement les cibles de ces kinases. Plusieurs peptides identifiés comme étant phosphorylés par la Casein Kinase 2 (CK2), une kinase connue pour être impliquée dans l’horloge circadienne des eucaryotes, proviennent de diverses protéines de liaison à l’ARN. Pour évaluer la possibilité que quelques-unes des multiples protéines à domaine Cold Shock identifiées dans le transcriptome puissent moduler l’expression des gènes de L. polyedrum, tel qu’observé chez plusieurs autres systèmes procaryotiques et eucaryotiques, la réponse des cellules à des températures froides a été examinée. Les températures froides ont permis d’induire rapidement un enkystement, condition dans laquelle ces cellules deviennent métaboliquement inactives afin de résister aux conditions environnementales défavorables. Les changements dans le profil des phosphoprotéines seraient le facteur majeur causant la formation de kystes. Les phosphosites prédits pour être phosphorylés par la CK2 sont la classe la plus fortement réduite dans les kystes, une découverte intéressante, car le rythme de la bioluminescence confirme que l’horloge a été arrêtée dans le kyste.

The South American Plasmodium falciparum var gene repertoire is limited, highly shared and possibly lacks several antigenic types

The Plasmodium falciparum var gene family encodes large variant antigens, which are important virulence factors, and also targets of the humoral host response. The frequently observed mild outcomes of falciparum malaria in many places of the Amazon area prompted us to ask whether a globally restricted variant (var) gene repertoire is present in currently circulating and older isolates of this area. By exhaustive analysis of var gene tags from 89 isolates and clones taken during many years from all over the Brazilian Amazon, we estimate that there are probably no more than 350-430 distinct sequence types, less than for any similar sized area studied so far. Detailed analysis of the var tags from genetically distinct clones obtained from single isolates revealed restricted and redundant repertoires suggesting either a low incidence of infective bites or restricted variant gene diversity in inoculated parasites. Additionally, we found a structuring of var gene repertoires observed as a higher pairwise typing sharing in isolates from the same microregion compared to isolates from different regions. Fine analysis of translated var tags revealed that certain Distinct Sequence Identifiers (DSIDs) were differently represented in Brazilian/South American isolates when compared to datasets from other continents. By global alignment of worldwide var DBL alpha sequences and sorting in groups with more than 76% identity, 125 clusters were formed and more than half of all genes were found in nine clusters with 50 or more sequences. While Brazilian/South American sequences were represented only in 64 groups, African sequences were found in the majority of clusters. DSID type 1 related sequences accumulated almost completely in one single cluster, indicating that limited recombination occurs in these specific var gene types. These data demonstrate the so far highest pairwise type sharing values for the var gene family in isolates from all over an entire subcontinent. The apparent lack of specific sequences types suggests that the P. falciparum transmission dynamics in the whole Amazon are probably different from any other endemic region studied and possibly interfere with the parasite`s ability to efficiently diversify its variant gene repertoires. (C) 2010 Elsevier B.V. All rights reserved.

Development of Y-chromosome-Specific SCAR Markers Conserved in Taurine, Zebu and Bubaline Cattle

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genetic relatedness within the genus Campylobacter inferred from rpoB sequences

The genus Campylobacter comprises 17 species, some of which are important animal and human pathogens. To gain more insight into the genetic relatedness of this genus and to improve the molecular tools available for diagnosis, a universal sequencing approach was established for the gene encoding the beta-subunit of RNA polymerase (rpoB) for the genus Campylobacter. A total of 59 strains, including the type strains of currently recognized species as well as field isolates, were investigated in the study. A primer set specific for Campylobacter species enabled straightforward amplification and sequencing of a 530 bp fragment of the rpoB gene. The 16S rRNA gene sequences of all of the strains were determined in parallel. A good congruence was obtained between 16S rRNA and rpoB gene sequence-based trees within the genus Campylobacter. The branching of the rpoB tree was similar to that of the 16S rRNA gene tree, even though a few discrepancies were observed for certain species. The resolution of the rpoB gene within the genus Campylobacter was generally much higher than that of the 16S rRNA gene sequence, resulting in a clear separation of most species and even some subspecies. The universally applicable amplification and sequencing approach for partial rpoB gene sequence determination provides a powerful tool for DNA sequence-based discrimination of Campylobacter species.

Comparative phylogenies of the housekeeping genes atpD, infB and rpoB and the 16S rRNA gene within the Pasteurellaceae

Phylogenies of housekeeping gene and 16S rRNA gene sequences were compared to improve the classification of the bacterial family Pasteurellaceae and knowledge of the evolutionary relationships of its members. Deduced partial protein sequences of the housekeeping genes atpD, infB and rpoB were compared in 28, 36 and 28 representative taxa of the Pasteurellaceae, respectively. The monophyly of representatives of the genus Gallibacterium was recognized by analysis of all housekeeping genes, while members of Mannheimia, Actinobacillus sensu stricto and the core group of Pasteurella sensu stricto formed monophyletic groups with two out of three housekeeping genes. Representatives of Mannheimia, Actinobacillus sensu stricto, [Haemophilus] ducreyi and [Pasteurella] trehalosi formed a monophyletic unit by analysis of all three housekeeping genes, which was in contrast to the 16S rRNA gene-derived phylogeny, where these taxa occurred at separate positions in the phylogenetic tree. Representatives of the Rodent, Avian and Aphrophilus-Haemophilus 16S rRNA gene groups were weakly supported by phylogenetic analysis of housekeeping genes. Phylogenies derived by comparison of the housekeeping genes diverged significantly from the 16S rRNA gene-derived phylogeny as evaluated by the likelihood ratio test. A low degree of congruence was also observed between the individual housekeeping gene-derived phylogenies. Estimates on speciation derived from 16S rRNA and housekeeping gene sequence comparisons resulted in quite different evolutionary scenarios for members of the Pasteurellaceae. The phylogeny based on the housekeeping genes supported observed host associations between Mannheimia, Actinobacillus sensu stricto and [Pasteurella] trehalosi and animals with paired hooves.

Transcriptional regulation of the chicken MLC3f gene

The expression of the chicken fast skeletal myosin alkali light chain (MLC) 3f is subject to complex patterns of control by developmental and physiologic signals. Regulation over MLC3f gene expression is thought to be exerted primarily at the transcriptional level. The purpose of this dissertation was to identify cis-acting elements on the 5$\sp\prime$ flanking region of chicken MLC3f gene that are important for transcriptional regulation. The results show that the 5$\sp\prime$ flanking region of MLC3f gene contains multiple cis-acting elements. The nucleotide sequence of these elements demonstrates a high degree of conservation between different species and are also found in the 5$\sp\prime$ flanking regions of many muscle protein genes. The first regulatory region is located between $-$185 and $-$150 bp from the transcription start site and contains an AT-rich element. Linker scanner analyses have revealed that this element has a positive effect on transcription of the MLC3f promoter. Furthermore, when linked to a heterologous viral promoter, it can enhance reporter gene expression in a muscle-specific manner, independent of distance or orientation.^ The second regulatory region is located between $-$96 and $-$64 from the transcription start site. Sequences downstream of $-$96 have the capacity to drive muscle-specific reporter gene expression, although the region between $-$96 and $-$64 has no intrinsic enhancer-like activity. Linker scanner analyses have identified a GC-rich motif that required efficient transcription of the MLC3f promoter. Mutations to this region of DNA results in diminished capacity to drive reporter gene expression and is correlated with disruption of the ability to bind sequence-specific transcription factors. These sequence-specific DNA-binding proteins were detected in both muscle and non-muscle extracts. The results suggest that the mere presence or absence of transcription factors cannot be solely responsible for regulation of MLC3f expression and that tissue-specific expression may arise from complex interactions with muscle-specific, as well as more ubiquitous transcription factors with multiple regulatory elements on the gene. ^

Stable isotope ratios of foraminifera and sea surface temperature estimation of sediment core W2804A-14

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C (epsilon p) in a central equatorial Pacific sediment core that spans the last ~255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic compostion of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon p, derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 (c e) based on a new, site-specific calibration of the relationship between epsilon p and c e. The calibration was based on reassessment of existing epsilon p versus c e data, which support a physiologically based model in which epsilon p is inversely related to c e. Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henry's law and the temperature determined from the alkenone-unsaturation index UK 37. Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon p and 1/c e. These are discussed in detail and it is concluded that the observed record of epsilon p most probably reflects significant variations in Delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from ~110 µatm during glacial intervals (ocean > atmosphere) to ~60 µatm during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p < 0.01) inverse correlations of Delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.

Neogene benthic foraminiferal faunas of a southwest Pacific bathyal depth transect

Temporal changes in benthic foraminiferal assemblages were quantitatively examined (> 63 µm fraction) in four southwest Pacific deep-sea Neogene sequences in a depth transect between approximately 1300 and 3200 m to assist in evaluating paleoeeanographic history. The most conspicuous changes in benthic foraminiferal assemblages occurred in association with paleoclimatic changes defined at least in part by oxygen isotopic changes. The largest, centered at ~15 Ma (early Middle Miocene), is represented by an increase in the relative frequencies of Epistominella exigua, which underwent a major upward depth migration at that time. This was contemporaneous with the well-known positive oxygen isotopic shift in the early Middle Miocene. In Sites 588 and 590, most of the increase in relative abundances of E. exigua occurred during the middle to later part of the ~80 shift, following major growth of the east Antarctic ice sheet. Later assemblage changes occurred at 8.5 and 6.5 Ma. These associations indicate that the benthic foraminiferal assemblages in this depth transect largely adjusted to changes in deep waters related to Antarctic cryospheric evolution. In general, the Neogene benthic foraminiferal assemblages in this region underwent little change during the last 23 million years. This faunal conservatism suggests that deep-sea environments underwent relatively little change in the southwest Pacific during much of the Neogene. Although paleoceanographic changes did occur, partly in response to highlatitude cryospheric evolution, these were not of sufficient magnitude to create major deep-sea faunal changes in this part of the ocean. The benthic foraminiferal assemblages are dominated by individuals smaller than 150 µm. Most taxonomic turnover occurred in the larger (> 150 µm) size fractions.

Extensive purification of a putative RNA polymerase I holoenzyme from plants that accurately initiates rRNA gene transcription in vitro

RNA polymerase I (pol I) is a nuclear enzyme whose function is to transcribe the duplicated genes encoding the precursor of the three largest ribosomal RNAs. We report a cell-free system from broccoli (Brassica oleracea) inflorescence that supports promoter-dependent RNA pol I transcription in vitro. The transcription system was purified extensively by DEAE-Sepharose, Biorex 70, Sephacryl S300, and Mono Q chromatography. Activities required for pre-rRNA transcription copurified with the polymerase on all four columns, suggesting their association as a complex. Purified fractions programmed transcription initiation from the in vivo start site and utilized the same core promoter sequences required in vivo. The complex was not dissociated in 800 mM KCl and had a molecular mass of nearly 2 MDa based on gel filtration chromatography. The most highly purified fractions contain ≈30 polypeptides, two of which were identified immunologically as RNA polymerase subunits. These data suggest that the occurrence of a holoenzyme complex is probably not unique to the pol II system but may be a general feature of eukaryotic nuclear polymerases.

A gene spans the pseudoautosomal boundary in mice

The X and Y chromosomes of the mouse, like those of other mammals, are heteromorphic over most of their length, but at the distal ends of the chromosomes is a region of sequence identity, the pseudoautosomal region (PAR), where the chromosomes pair and recombine during male meiosis. The point at which the PAR diverges into X- and Y-specific sequences is called the pseudoautosomal boundary. We have completed a genomic walk from the X-specific Amelogenin gene to the PAR. Analysis of this region revealed that the pseudoautosomal boundary of mice is located within an intron of a transcribed gene that encodes a novel RING finger protein. The first three of the exons of the gene are located on the X chromosome whereas the 3′ exons of the gene are located on both X and Y chromosomes. This unusual arrangement may indicate that the gene is in a state of transition from pseudoautosomal to X-unique and provides evidence for a process of attrition of the pseudoautosomal region on the Y chromosome.

Specific molecular recognition of mixed nucleic acid sequences: An aromatic dication that binds in the DNA minor groove as a dimer

Phenylamidine cationic groups linked by a furan ring (furamidine) and related compounds bind as monomers to AT sequences of DNA. An unsymmetric derivative (DB293) with one of the phenyl rings of furamidine replaced with a benzimidazole has been found by quantitative footprinting analyses to bind to GC-containing sites on DNA more strongly than to pure AT sequences. NMR structural analysis and surface plasmon resonance binding results clearly demonstrate that DB293 binds in the minor groove at specific GC-containing sequences of DNA in a highly cooperative manner as a stacked dimer. Neither the symmetric bisphenyl nor bisbenzimidazole analogs of DB293 bind significantly to the GC containing sequences. DB293 provides a paradigm for design of compounds for specific recognition of mixed DNA sequences and extends the boundaries for small molecule-DNA recognition.

Efficient and specific ribozyme-mediated reduction of bovine alpha-lactalbumin expression in double transgenic mice.

Transgenic mice carrying a bovine alpha-lactalbumin (alpha-lac) specific ribozyme gene under the transcriptional control of the mouse mammary tumor virus long terminal repeat were generated and cross-bred with animals that highly express a bovine alpha-lac transgene (0.4 mg of alpha-lac/ml(-1) of milk). The ribozyme contains the hammerhead catalytic domain, flanked by 12-nt sequences complementary to the 3' untranslated region of bovine alpha-lac transcript. High-level expression of the ribozyme gene was detected by Northern blot analysis in the mammary gland of 7-8 day lactating transgenic mice, from 3 of 12 lines analyzed. Heterozygous expression of the ribozyme resulted in a reduction in the levels of the target mRNA to 78, 58, and 50% of that observed in the nonribozyme transgenic littermate controls for three independent lines. The ribozyme-mediated reduction in the levels of the bovine protein paralleled that observed for the mRNA, and was positively correlated with the level of expression of the ribozyme. In nonribozyme expressing transgenic mice, the level of bovine alpha-lac mRNA and protein was not affected. The specificity of this activity is demonstrated by the absence of a reduction in the levels of the endogenous murine alpha-lac mRNA or protein. These results demonstrate the feasibility of ribozyme-mediated down-regulation of highly-expressed transcripts in transgenic animals.

The genomic and physical organization of Ty1-copia-like sequences as a component of large genomes in Pinus elliottii var. elliottii and other gymnosperms.

A DNA sequence, TPE1, representing the internal domain of a Ty1-copia retroelement, was isolated from genomic DNA of Pinus elliottii Engelm. var. elliottii (slash pine). Genomic Southern analysis showed that this sequence, carrying partial reverse transcriptase and integrase gene sequences, is highly amplified within the genome of slash pine and part of a dispersed element >4.8 kbp. Fluorescent in situ hybridization to metaphase chromosomes shows that the element is relatively uniformly dispersed over all 12 chromosome pairs and is highly abundant in the genome. It is largely excluded from centromeric regions and intercalary chromosomal sites representing the 18S-5.8S-25S rRNA genes. Southern hybridization with specific DNA probes for the reverse transcriptase gene shows that TPE1 represents a large subgroup of heterogeneous Ty1-copia retrotransposons in Pinus species. Because no TPE1 transcription could be detected, it is most likely an inactive element--at least in needle tissue. Further evidence for inactivity was found in recombinant reverse transcriptase and integrase sequences. The distribution of TPE1 within different gymnosperms that contain Ty1-copia group retrotransposons, as shown by a PCR assay, was investigated by Southern hybridization. The TPE1 family is highly amplified and conserved in all Pinus species analyzed, showing a similar genomic organization in the three- and five-needle pine species investigated. It is also present in spruce, bald cypress (swamp cypress), and in gingko but in fewer copies and a different genomic organization.