Tandem Gene Arrays Research Articles

Short Synthetic Terminators for Assembly of Transcription Units in Vitro and Stable Chromosomal Integration in Yeast S. cerevisiae.

Assembly of synthetic genetic circuits is central to synthetic biology. Yeast S. cerevisiae, in particular, has proven to be an ideal chassis for synthetic genome assemblies by exploiting its efficient homologous recombination. However, this property of efficient homologous recombination poses a problem for multigene assemblies in yeast, since repeated usage of standard parts, such as transcriptional terminators, can lead to rearrangements of the repeats in assembled DNA constructs in vivo. To address this issue in developing a library of orthogonal genetic components for yeast, we designed a set of short synthetic terminators based on a consensus sequence with random linkers to avoid repetitive sequences. We constructed a series of expression vectors with these synthetic terminators for efficient assembly of synthetic genes using Gateway recombination reactions. We also constructed two BAC (bacterial artificial chromosome) vectors for assembling multiple transcription units with the synthetic terminators in vitro and their integration in the yeast genome. The tandem array of synthetic genes integrated in the genome by this method is highly stable because there are few homologous segments in the synthetic constructs. Using this system of assembly and genomic integration of transcription units, we tested the synthetic terminators and their influence on the proximal transcription units. Although all the synthetic terminators have the common consensus with the identical length, they showed different activities and impacts on the neighboring transcription units.

Dinoflagellate Gene Structure and Intron Splice Sites in a Genomic Tandem Array.

Dinoflagellates are one of the last major lineages of eukaryotes for which little is known about genome structure and organization. We report here the sequence and gene structure of a clone isolated from a cosmid library which, to our knowledge, represents the largest contiguously sequenced, dinoflagellate genomic, tandem gene array. These data, combined with information from a large transcriptomic library, allowed a high level of confidence of every base pair call. This degree of confidence is not possible with PCR‐based contigs. The sequence contains an intron‐rich set of five highly expressed gene repeats arranged in tandem. One of the tandem repeat gene members contains an intron 26,372 bp long. This study characterizes a splice site consensus sequence for dinoflagellate introns. Two to nine base pairs around the 3′ splice site are repeated by an identical two to nine base pairs around the 5′ splice site. The 5′ and 3′ splice sites are in the same locations within each repeat so that the repeat is found only once in the mature mRNA. This identically repeated intron boundary sequence might be useful in gene modeling and annotation of genomes.

Comparative genomics identifies the Magnaporthe oryzae avirulence effector AvrPi9 that triggers Pi9‐mediated blast resistance in rice

We identified the Magnaporthe oryzae avirulence effector AvrPi9 cognate to rice blast resistance gene Pi9 by comparative genomics of requisite strains derived from a sequential planting method. AvrPi9 encodes a small secreted protein that appears to localize in the biotrophic interfacial complex and is translocated to the host cell during rice infection. AvrPi9 forms a tandem gene array with its paralogue proximal to centromeric region of chromosome 7. AvrPi9 is expressed highly at early stages during initiation of blast disease. Virulent isolate strains contain Mg-SINE within the AvrPi9 coding sequence. Loss of AvrPi9 did not lead to any discernible defects during growth or pathogenesis in M. oryzae. This study reiterates the role of diverse transposable elements as off-switch agents in acquisition of gain-of-virulence in the rice blast fungus. The prevalence of AvrPi9 correlates well with the avirulence pathotype in diverse blast isolates from the Philippines and China, thus supporting the broad-spectrum resistance conferred by Pi9 in different rice growing areas. Our results revealed that Pi9 and Piz-t at the Pi2/9 locus activate race specific resistance by recognizing sequence-unrelated AvrPi9 and AvrPiz-t genes, respectively.

Mono-allelic VSG expression by RNA polymerase I in Trypanosoma brucei: Expression site control from both ends?

Trypanosoma brucei is a vector borne, lethal protistan parasite of humans and livestock in sub-Saharan Africa. Antigenic variation of its cell surface coat enables the parasite to evade adaptive immune responses and to live freely in the blood of its mammalian hosts. The coat consists of ten million copies of variant surface glycoprotein (VSG) that is expressed from a single VSG gene, drawn from a large repertoire and located near the telomere at one of fifteen so-called bloodstream expression sites (BESs). Thus, antigenic variation is achieved by switching to the expression of a different VSG gene. A BES is a tandem array of expression site-associated genes and a terminal VSG gene. It is polycistronically transcribed by a multifunctional RNA polymerase I (RNAPI) from a short promoter that is located 45–60kb upstream of the VSG gene. The mechanism(s) restricting VSG expression to a single BES are not well understood. There is convincing evidence that epigenetic silencing and transcription attenuation play important roles. Furthermore, recent data indicated that there is regulation at the level of transcription initiation and that, surprisingly, the VSG mRNA appears to have a role in restricting VSG expression to a single gene. Here, we review BES expression regulation and propose a model in which telomere-directed, epigenetic BES silencing is opposed by BES promoter-directed, activated RNAPI transcription.

The complete genome of Blastobotrys (Arxula) adeninivorans LS3 - a yeast of biotechnological interest.

BackgroundThe industrially important yeast Blastobotrys (Arxula) adeninivorans is an asexual hemiascomycete phylogenetically very distant from Saccharomyces cerevisiae. Its unusual metabolic flexibility allows it to use a wide range of carbon and nitrogen sources, while being thermotolerant, xerotolerant and osmotolerant.ResultsThe sequencing of strain LS3 revealed that the nuclear genome of A. adeninivorans is 11.8 Mb long and consists of four chromosomes with regional centromeres. Its closest sequenced relative is Yarrowia lipolytica, although mean conservation of orthologs is low. With 914 introns within 6116 genes, A. adeninivorans is one of the most intron-rich hemiascomycetes sequenced to date. Several large species-specific families appear to result from multiple rounds of segmental duplications of tandem gene arrays, a novel mechanism not yet described in yeasts. An analysis of the genome and its transcriptome revealed enzymes with biotechnological potential, such as two extracellular tannases (Atan1p and Atan2p) of the tannic-acid catabolic route, and a new pathway for the assimilation of n-butanol via butyric aldehyde and butyric acid.ConclusionsThe high-quality genome of this species that diverged early in Saccharomycotina will allow further fundamental studies on comparative genomics, evolution and phylogenetics. Protein components of different pathways for carbon and nitrogen source utilization were identified, which so far has remained unexplored in yeast, offering clues for further biotechnological developments. In the course of identifying alternative microorganisms for biotechnological interest, A. adeninivorans has already proved its strengthened competitiveness as a promising cell factory for many more applications.

Measuring chromatin interaction dynamics on the second time scale at single-copy genes.

The chromatin immunoprecipitation (ChIP) assay is widely used to capture interactions between chromatin and regulatory proteins, but it is unknown how stable most native interactions are. Although live-cell imaging suggests short-lived interactions at tandem gene arrays, current methods cannot measure rapid binding dynamics at single-copy genes. We show, by using a modified ChIP assay with subsecond temporal resolution, that the time dependence of formaldehyde cross-linking can be used to extract in vivo on and off rates for site-specific chromatin interactions varying over a ~100-fold dynamic range. By using the method, we show that a regulatory process can shift weakly bound TATA-binding protein to stable promoter interactions, thereby facilitating transcription complex formation. This assay provides an approach for systematic, quantitative analyses of chromatin binding dynamics in vivo.

Quantifying the transcriptional output of single alleles in single living mammalian cells

Transcription kinetics of actively transcribing genes in vivo have generally been measured using tandem gene arrays. However, tandem arrays do not reflect the endogenous state of genome organization in which genes appear as single alleles. Here we present a robust technique for the quantification of mRNA synthesis from a single allele in real time in single living mammalian cells. The protocol describes how to generate cell clones harboring an MS2-tagged allele and how to detect in vivo transcription from this tagged allele at high spatial and temporal resolution throughout the cell cycle. Quantification of nascent mRNAs produced from the single tagged allele is performed using RNA fluorescence in situ hybridization (FISH) and live-cell imaging. Subsequent analyses and data modeling detailed in the protocol include measurements of transcription rates of RNA polymerase II, determination of the number of polymerases recruited to the tagged allele and measurement of the spacing between polymerases. Generation of the cells containing the single tagged alleles should take up to 1 month; RNA FISH or live-cell imaging will require an additional week.

Dinoflagellate tandem array gene transcripts are highly conserved and not polycistronic

Dinoflagellates are an important component of the marine biota, but a large genome with high-copy number (up to 5,000) tandem gene arrays has made genomic sequencing problematic. More importantly, little is known about the expression and conservation of these unusual gene arrays. We assembled de novo a gene catalog of 74,655 contigs for the dinoflagellate Lingulodinium polyedrum from RNA-Seq (Illumina) reads. The catalog contains 93% of a Lingulodinium EST dataset deposited in GenBank and 94% of the enzymes in 16 primary metabolic KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways, indicating it is a good representation of the transcriptome. Analysis of the catalog shows a marked underrepresentation of DNA-binding proteins and DNA-binding domains compared with other algae. Despite this, we found no evidence to support the proposal of polycistronic transcription, including a marked underrepresentation of sequences corresponding to the intergenic spacers of two tandem array genes. We also have used RNA-Seq to assess the degree of sequence conservation in tandem array genes and found their transcripts to be highly conserved. Interestingly, some of the sequences in the catalog have only bacterial homologs and are potential candidates for horizontal gene transfer. These presumably were transferred as single-copy genes, and because they are now all GC-rich, any derived from AT-rich contexts must have experienced extensive mutation. Our study not only has provided the most complete dinoflagellate gene catalog known to date, it has also exploited RNA-Seq to address fundamental issues in basic transcription mechanisms and sequence conservation in these algae.

The Impact of Genome Triplication on Tandem Gene Evolution in Brassica rapa

Whole genome duplication (WGD) and tandem duplication (TD) are both important modes of gene expansion. However, how WGD influences tandemly duplicated genes is not well studied. We used Brassica rapa, which has undergone an additional genome triplication (WGT) and shares a common ancestor with Arabidopsis thaliana, Arabidopsis lyrata, and Thellungiella parvula, to investigate the impact of genome triplication on tandem gene evolution. We identified 2,137, 1,569, 1,751, and 1,135 tandem gene arrays in B. rapa, A. thaliana, A. lyrata, and T. parvula respectively. Among them, 414 conserved tandem arrays are shared by the three species without WGT, which were also considered as existing in the diploid ancestor of B. rapa. Thus, after genome triplication, B. rapa should have 1,242 tandem arrays according to the 414 conserved tandems. Here, we found 400 out of the 414 tandems had at least one syntenic ortholog in the genome of B. rapa. Furthermore, 294 out of the 400 shared syntenic orthologs maintain tandem arrays (more than one gene for each syntenic hit) in B. rapa. For the 294 tandem arrays, we obtained 426 copies of syntenic paralogous tandems in the triplicated genome of B. rapa. In this study, we demonstrated that tandem arrays in B. rapa were dramatically fractionated after WGT when compared either to non-tandem genes in the B. rapa genome or to the tandem arrays in closely related species that have not experienced a recent whole genome polyploidization event.

Visualizing the Transcription Cycle of Endogenous RNA Polymerase II in Single Living Cells

How efficient is transcription in vivo? Recent fluorescence microscopy studies have begun to confront this question by measuring the binding times of RNA polymerase II (pol II) and other components of the transcription machinery in single living cells. This has revealed an extremely dynamic transcription machine that appears to operate with components in continual flux. However, it remains debatable whether or not these rapid dynamics lead to efficient transcription, mainly because it has been difficult to distinguish different phases of the transcription cycle in vivo. For example, using GFP it is possible to visualize the movement of pol II, but directly distinguishing the uninitiated fraction from the initiated or elongating fractions has not yet been possible. To overcome this difficulty, we have loaded fluorescent antibody fragments (Fab) against unphosphorylated and phosphorylated forms of pol II into living cells containing an inducible tandem gene array. Using this unique system, we are quantifying the accumulation of uninitiated (unphosphorylated), initiated (ser 5 phosphorylated), and elongating (ser 2 phosphorylated) forms of pol II to the gene array after induction. Our observations suggest transcription is quite efficient, with pol II being recruited to activated genes within ∼3 minutes time, after which ∼70% are initiated within 30 seconds and ∼35% proceed to elongation within another 30 seconds. We are currently investigating how these transcription dynamics correlate with histone modification dynamics at the gene array.

The evolution and consequences of snaR family transposition in primates

The small NF90 associated RNA (snaR) family of small noncoding RNAs (ncRNA) appears to have evolved from retrotransposon ancestors at or soon after pivotal stages in primate evolution. snaRs are thought to be derived from a FLAM C-like (free left Alu monomer) element through multiple short insertion/deletion (indel) and nucleotide (nt) substitution events. Tracing snaR’s complex evolutionary history through primate genomes led to the recent discovery of two novel retrotransposons: the Alu/snaR related (ASR) and catarrhine ancestor of snaR (CAS) elements. ASR elements are present in the genomes of Simiiformes, CAS elements are present in Old World Monkeys and apes, and snaRs are restricted to the African Great Apes (Homininae, including human, gorilla, chimpanzee and bonobo). Unlike their ancestors, snaRs have disseminated by multiple rounds of segmental duplication of a larger encompassing element. This process has produced large tandem gene arrays in humans and possibly precipitated the accelerated evolution of snaR. Furthermore, snaR segmental duplication created a new form of chorionic gonadotropin β subunit (CGβ) gene, recently classified as Type II CGβ, which has altered mRNA tissue expression and can generate a novel short peptide.

Chromosome and gene copy number variation allow major structural change between species and strains of Leishmania.

Leishmania parasites cause a spectrum of clinical pathology in humans ranging from disfiguring cutaneous lesions to fatal visceral leishmaniasis. We have generated a reference genome for Leishmania mexicana and refined the reference genomes for Leishmania major, Leishmania infantum, and Leishmania braziliensis. This has allowed the identification of a remarkably low number of genes or paralog groups (2, 14, 19, and 67, respectively) unique to one species. These were found to be conserved in additional isolates of the same species. We have predicted allelic variation and find that in these isolates, L. major and L. infantum have a surprisingly low number of predicted heterozygous SNPs compared with L. braziliensis and L. mexicana. We used short read coverage to infer ploidy and gene copy numbers, identifying large copy number variations between species, with 200 tandem gene arrays in L. major and 132 in L. mexicana. Chromosome copy number also varied significantly between species, with nine supernumerary chromosomes in L. infantum, four in L. mexicana, two in L. braziliensis, and one in L. major. A significant bias against gene arrays on supernumerary chromosomes was shown to exist, indicating that duplication events occur more frequently on disomic chromosomes. Taken together, our data demonstrate that there is little variation in unique gene content across Leishmania species, but large-scale genetic heterogeneity can result through gene amplification on disomic chromosomes and variation in chromosome number. Increased gene copy number due to chromosome amplification may contribute to alterations in gene expression in response to environmental conditions in the host, providing a genetic basis for disease tropism.

Tandem gene arrays, plastic chromosomal organizations

This short article presents an overview of tandem gene arrays (TGAs) in hemiascomycete yeasts. In silico and in vivo analyses are combined to address structural, functional and evolutionary aspects of these particular chromosomal structures. Genomic instability of TGAs is discussed. We conclude that TGAs are generally dynamic regions of the genome in that they are the seats of chromosomal rearrangement events. In addition, they are often breeding grounds of new genes for a rapid adaptation of cells to demands of the environment.

An Exceptional Gene: Evolution of the TSPY Gene Family in Humans and Other Great Apes

The TSPY gene stands out from all other human protein-coding genes because of its high copy number and tandemly-repeated organization. Here, we review its evolutionary history in great apes in order to assess whether these unusual properties are more likely to result from a relaxation of constraint or an unusual functional role. Detailed comparisons with chimpanzee are possible because a finished sequence of the chimpanzee Y chromosome is available, together with more limited data from other apes. These comparisons suggest that the human-chimpanzee ancestral Y chromosome carried a tandem array of TSPY genes which expanded on the human lineage while undergoing multiple duplication events followed by pseudogene formation on the chimpanzee lineage. The protein coding region is the most highly conserved of the multi-copy Y genes in human-chimpanzee comparisons, and the analysis of the dN/dS ratio indicates that TSPY is evolutionarily highly constrained, but may have experienced positive selection after the human-chimpanzee split. We therefore conclude that the exceptionally high copy number in humans is most likely due to a human-specific but unknown functional role, possibly involving rapid production of a large amount of TSPY protein at some stage during spermatogenesis.

The amplification and evolution of orthologous 22-kDa α-prolamin tandemly arrayed genes in coix, sorghum and maize genomes

Tandemly arrayed genes (TAGs) account for about one-third of the duplicated genes in eukaryotic genomes. They provide raw genetic material for biological evolution, and play important roles in genome evolution. The 22-kDa prolamin genes in cereal genomes represent typical TAG organization, and provide the good material to investigate gene amplification of TAGs in closely related grass genomes. Here, we isolated and sequenced the Coix 22-kDa prolamin (coixin) gene cluster (283 kb), and carried out a comparative analysis with orthologous 22-kDa prolamin gene clusters from maize and sorghum. The 22-kDa prolamin gene clusters descended from orthologous ancestor genes, but underwent independent gene amplification paths after the separation of these species, therefore varied dramatically in sequence and organization. Our analysis indicated that the gene amplification model of 22-kDa prolamin gene clusters can be divided into three major stages. In the first stage, rare gene duplications occurred from the ancestor gene copy accidentally. In the second stage, rounds of gene amplification occurred by unequal crossing over to form tandem gene array(s). In the third stage, gene array was further diverged by other genomic activities, such as transposon insertions, segmental rearrangements, etc. Unlike their highly conserved sequences, the amplified 22-kDa prolamin genes diverged rapidly at their expression capacities and expression levels. Such processes had no apparent correlation to age or order of amplified genes within TAG cluster, suggesting a fast evolving nature of TAGs after gene amplification. These results provided insights into the amplification and evolution of TAG families in grasses.

Natural Diversity in Flowering Responses of Arabidopsis thaliana Caused by Variation in a Tandem Gene Array

Tandemly arrayed genes that belong to gene families characterize genomes of many organisms. Gene duplication and subsequent relaxation of selection can lead to the establishment of paralogous cluster members that may evolve along different trajectories. Here, we report on the structural variation in MADS AFFECTING FLOWERING 2 (MAF2) gene, one member of the tandemly duplicated cluster of MADS-box-containing transcription factors in Arabidopsis thaliana. The altered gene structure at the MAF2 locus is present as a moderate-frequency polymorphism in Arabidopsis and leads to the extensive diversity in transcript patterns due to alternative splicing. Rearrangements at the MAF2 locus are associated with an early flowering phenotype in BC(5) lines. The lack of suppression of flowering time in a MAF2-insertion line expressing the MAF2-specific artificial miRNA suggests that these MAF2 variants are behaving as loss-of-function alleles. The variation in gene architecture is also associated with segregation distortion, which may have facilitated the spread and the establishment of the corresponding alleles throughout the Eurasian range of the A. thaliana population.

Radical Redesign of a Tandem Array of Four R67 Dihydrofolate Reductase Genes Yields a Functional, Folded Protein Possessing 45 Substitutions

R67 dihydrofolate reductase (DHFR) is a plasmid-encoded, type II enzyme. Four monomers (78 amino acids long) assemble into a homotetramer possessing 222 symmetry. In previous studies, a tandem array of four R67 DHFR gene copies was fused in frame to generate a functional monomer named Quad1. This protein possessed the essential tertiary structure of the R67 "parent". To facilitate mutagenesis reactions, restriction enzyme sites were introduced in the tandem gene array. S59A and H362L mutations were also added to minimize possible folding topologies; this protein product, named Quad3, possesses 10 substitutions and is functional. Since R67 DHFR possesses a stable scaffold, a large jump in sequence space was performed by the further addition of 45 amino acid substitutions. The mutational design utilized alternate sequences from other type II DHFRs. In addition, most of the mutations were positioned on the surface of the protein as well as in the disordered N-terminal sequence, which serves as the linker between the fused domains. The resulting Quad4 protein is quite functional; however, it is less stable than Quad1, suffering a DeltaDeltaG loss of 5 kcal/mol at pH 5. One unexpected result was formation of Quad4 dimers and higher order oligomers at pH 8. R67 DHFR, and its derivative Quad proteins, possesses a robust scaffold, capable of withstanding introduction of >or=55 substitutions.

The Genome Sequence of Trypanosoma brucei gambiense, Causative Agent of Chronic Human African Trypanosomiasis

Background Trypanosoma brucei gambiense is the causative agent of chronic Human African Trypanosomiasis or sleeping sickness, a disease endemic across often poor and rural areas of Western and Central Africa. We have previously published the genome sequence of a T. b. brucei isolate, and have now employed a comparative genomics approach to understand the scale of genomic variation between T. b. gambiense and the reference genome. We sought to identify features that were uniquely associated with T. b. gambiense and its ability to infect humans.Methods and FindingsAn improved high-quality draft genome sequence for the group 1 T. b. gambiense DAL 972 isolate was produced using a whole-genome shotgun strategy. Comparison with T. b. brucei showed that sequence identity averages 99.2% in coding regions, and gene order is largely collinear. However, variation associated with segmental duplications and tandem gene arrays suggests some reduction of functional repertoire in T. b. gambiense DAL 972. A comparison of the variant surface glycoproteins (VSG) in T. b. brucei with all T. b. gambiense sequence reads showed that the essential structural repertoire of VSG domains is conserved across T. brucei.ConclusionsThis study provides the first estimate of intraspecific genomic variation within T. brucei, and so has important consequences for future population genomics studies. We have shown that the T. b. gambiense genome corresponds closely with the reference, which should therefore be an effective scaffold for any T. brucei genome sequence data. As VSG repertoire is also well conserved, it may be feasible to describe the total diversity of variant antigens. While we describe several as yet uncharacterized gene families with predicted cell surface roles that were expanded in number in T. b. brucei, no T. b. gambiense-specific gene was identified outside of the subtelomeres that could explain the ability to infect humans.

Genome-wide computational prediction of tandem gene arrays: application in yeasts

BackgroundThis paper describes an efficient in silico method for detecting tandem gene arrays (TGAs) in fully sequenced and compact genomes such as those of prokaryotes or unicellular eukaryotes. The originality of this method lies in the search of protein sequence similarities in the vicinity of each coding sequence, which allows the prediction of tandem duplicated gene copies independently of their functionality.ResultsApplied to nine hemiascomycete yeast genomes, this method predicts that 2% of the genes are involved in TGAs and gene relics are present in 11% of TGAs. The frequency of TGAs with degenerated gene copies means that a significant fraction of tandem duplicated genes follows the birth-and-death model of evolution. A comparison of sequence identity distributions between sets of homologous gene pairs shows that the different copies of tandem arrayed paralogs are less divergent than copies of dispersed paralogs in yeast genomes. It suggests that paralogs included in tandem structures are more recent or more subject to the gene conversion mechanism than other paralogs.ConclusionThe method reported here is a useful computational tool to provide a database of TGAs composed of functional or nonfunctional gene copies. Such a database has obvious applications in the fields of structural and comparative genomics. Notably, a detailed study of the TGA catalog will make it possible to tackle the fundamental questions of the origin and evolution of tandem gene clusters.

Combined phylogeny and neighborhood analysis of the evolution of the ABC transporters conferring multiple drug resistance in hemiascomycete yeasts

BackgroundPleiotropic Drug Resistant transporters (PDR) are members of the ATP-Binding Cassette (ABC) subfamily which export antifungals and other xenobiotics in fungi and plants. This subfamily of transmembrane transporters has nine known members in Saccharomyces cerevisiae.We have analyzed the complex evolution of the pleiotropic drug resistance proteins (Pdrp) subfamily where gene duplications and deletions occur independently in individual genomes.This study was carried out on 62 Pdrp from nine hemiascomycetous species, seven of which span 6 of the 14 clades of the Saccharomyces complex while the two others species, Debaryomyces hansenii and Yarrowia lipolytica, are further apart from an evolutive point of view.ResultsCombined phylogenetic and neighborhood analyses enabled us to identify five Pdrp clusters in the Saccharomyces complex. Three of them comprise orthologs of the Pdrp sensu stricto, Pdr5p, Pdr10p, Pdr12p, Pdr15p, Snq2p and YNR070wp. The evolutive pathway of the orthologs of Snq2 and YNR070w is particularly complex due to a tandem gene array in Eremothecium gossypii, Kluyveromyces lactis and Saccharomyces (Lachancea) kluyveri. This pathway and different cases of duplications and deletions were clarified by using a neighborhood analysis based on synteny.For the two distant species, Yarrowia lipolytica and Debaryomyces hansenii, no neighborhood evidence is available for these clusters and many homologs of Pdr5 and Pdr15 are phylogenetically assigned to species-based clusters.Two other clusters comprise the orthologs of the sensu lato Pdrp, Aus1p/Pdr11p and YOL075cp respectively. The evolutionary pathway of these clusters is simpler. Nevertheless, orthologs of these genes are missing in some species.ConclusionNumerous duplications were traced among the Hemiascomycetous Pdrp studied. The role of the Whole Genome Duplication (WGD) is sorted out and our analyses confirm the common ancestrality of Pdr5p and Pdr15p. A tandem gene array is observed in Eremothecium gossypii. One of the copies is the ortholog of Snq2 while the other one is lost in the post-WGD species.The neighborhood analysis provides an efficient method to trace the history of genes and disentangle the orthology and paralogy relationships.