Nonrandom Genomic Distribution Research Articles

Functional characterization of archaic-specific variants in mitonuclear genes: insights from comparative analysis in S. cerevisiae.

Neanderthal and Denisovan hybridisation with modern humans has generated a non-random genomic distribution of introgressed regions, the result of drift and selection dynamics. Cross-species genomic incompatibility and more efficient removal of slightly deleterious archaic variants have been proposed as selection-based processes involved in the post-hybridisation purge of archaic introgressed regions. Both scenarios require the presence of functionally different alleles across Homo species onto which selection operated differently according to which populations hosted them, but only a few of these variants have been pinpointed so far. In order to identify functionally divergent archaic variants removed in humans, we focused on mitonuclear genes, which are underrepresented in the genomic landscape of archaic humans. We searched for non-synonymous, fixed, archaic-derived variants present in mitonuclear genes, rare or absent in human populations. We then compared the functional impact of archaic and human variants in the model organism Saccharomyces cerevisiae. Notably, a variant within the mitochondrial tyrosyl-tRNA synthetase 2 (YARS2) gene exhibited a significant decrease in respiratory activity and a substantial reduction of Cox2 levels, a proxy for mitochondrial protein biosynthesis, coupled with the accumulation of the YARS2 protein precursor and a lower amount of mature enzyme. Our work suggests that this variant is associated with mitochondrial functionality impairment, thus contributing to the purging of archaic introgression in YARS2. While different molecular mechanisms may have impacted other mitonuclear genes, our approach can be extended to the functional screening of mitonuclear genetic variants present across species and populations.

Functionally Related Genes Cluster into Genomic Regions That Coordinate Transcription at a Distance in Saccharomyces cerevisiae.

Balancing gene expression is a fundamental challenge of all cell types. To properly regulate transcription on a genome-wide level, there are myriad mechanisms employed by the cell. One layer to this regulation is through spatial positioning, with particular chromosomal loci exerting an influence on transcription throughout a region. Many coregulated gene families utilize spatial positioning to coordinate transcription, with functionally related genes clustering together which can allow coordinated expression via adjacent gene coregulation. The mechanisms underlying this process have not been elucidated, though there are many coregulated gene families that exhibit this genomic distribution. In the present study, we tested for a role for the enhancer-promoter (EP) hypothesis, which demonstrates that regulatory elements can exert transcriptional effects over a broad distance, in coordinating transcriptional coregulation using budding yeast, Saccharomyces cerevisiae We empirically validated the EP model, finding that the genomic distance a promoter can affect varies by locus, which can profoundly affect levels of transcription, phenotype, and the extent of transcriptional disruption throughout a genomic region. Using the nitrogen metabolism, ribosomal protein, toxin response, and heat shock gene families as our test case, we report functionally clustered genes localize to genomic loci that are more conducive to transcriptional regulation at a distance compared to the unpaired members of the same families. Furthermore, we report that the coregulation of functional clusters is dependent, in part, on chromatin maintenance and remodeling, providing one mechanism underlying adjacent gene coregulation.IMPORTANCE The two-dimensional, physical positioning of genes along a chromosome can impact proper transcriptional regulation throughout a genomic region. The transcription of neighboring genes is correlated in a genome-wide manner, which is a characteristic of eukaryotes. Many coregulated gene families can be found clustered with another member of the same set-which can result in adjacent gene coregulation of the pair. Due to the myriad gene families that exhibit a nonrandom genomic distribution, there are likely multiple mechanisms working in concert to properly regulate transcriptional coordination of functionally clustered genes. In this study, we utilized budding yeast in an attempt to elucidate mechanisms that underlie this coregulation: testing and empirically validating the enhancer-promoter hypothesis in this species and reporting that functionally related genes cluster to genomic regions that are more conducive to transcriptional regulation at a distance. These clusters rely, in part, on chromatin maintenance and remodelers to maintain proper transcriptional coordination. Our work provides insight into the mechanisms underlying adjacent gene coregulation.

Systematic Identification, Characterization, and Conservation of Adjacent-Gene Coregulation in the Budding Yeast Saccharomyces cerevisiae.

It is essential that cells orchestrate gene expression for the specific niche that they occupy, and this often requires coordination of the expression of large sets of genes. There are multiple regulatory systems that exist for modulation of gene expression, including the adjacent-gene coregulation of the rRNA and ribosome biogenesis and ribosomal protein families. Both gene families exhibit a nonrandom genomic distribution, often clustered directly adjacent to another member of the same family, which results in a tighter transcriptional coordination among adjacent paired genes than that of the unpaired genes within each regulon and can result in a shared promoter that coordinates expression of the pairs. This nonrandom genomic distribution has been seen in a few functionally related gene families, and many of these functional pairings are conserved across divergent fungal lineages. To date, the significance of these observations has not been extended in a systematic way to characterize how prevalent the role of adjacent-gene coregulation is in transcriptional regulation. In the present study, we systematically analyzed the transcriptional coherence of the functional pairs compared to the singletons within all gene families defined by the Gene Ontology Slim designation, using Saccharomyces cerevisiae as a model system, finding that clusters exhibit a tighter transcriptional correlation under specific contexts. We found that the longer a functional pairing is conserved the tighter its response to broad stress and nutritional responses, that roughly 25% of gene families exhibit a nonrandom genomic distribution, and that many of these clusters are conserved. This suggests that adjacent-gene coregulation is a widespread, yet underappreciated, transcriptional mechanism.IMPORTANCE The spatial positioning of genes throughout the genome arrangement can alter their expression in many eukaryotic organisms. Often this results in a genomic context-specific effect on transcription. One example of this is through the clustering of functionally related genes, which results in adjacent-gene coregulation in the budding yeast Saccharomyces cerevisiae In the present study, we set out to systematically characterize the prevalence of this phenomenon, finding the genomic organization of functionally related genes into clusters is a characteristic of myriad gene families. These arrangements are found in many evolutionarily divergent fungi and thus represent a widespread, yet underappreciated, layer of transcriptional regulation.

Obtaining retrotransposon sequences, analysis of their genomic distribution and use of retrotransposon-derived genetic markers in lentil (Lens culinaris Medik.).

Retrotransposons with long terminal repeats (LTR-RTs) are widespread mobile elements in eukaryotic genomes. We obtained a total of 81 partial LTR-RT sequences from lentil corresponding to internal retrotransposon components and LTRs. Sequences were obtained by PCR from genomic DNA. Approximately 37% of the LTR-RT internal sequences presented premature stop codons, pointing out that these elements must be non-autonomous. LTR sequences were obtained using the iPBS technique which amplifies sequences between LTR-RTs. A total of 193 retrotransposon-derived genetic markers, mainly iPBS, were used to obtain a genetic linkage map from 94 F7 inbred recombinant lines derived from the cross between the cultivar Lupa and the wild ancestor L. culinaris subsp. orientalis. The genetic map included 136 markers located in eight linkage groups. Clusters of tightly linked retrotransposon-derived markers were detected in linkage groups LG1, LG2, and LG6, hence denoting a non-random genomic distribution. Phylogenetic analyses identified the LTR-RT families in which internal and LTR sequences are included. Ty3-gypsy elements were more frequent than Ty1-copia, mainly due to the high Ogre element frequency in lentil, as also occurs in other species of the tribe Vicieae. LTR and internal sequences were used to analyze in silico their distribution among the contigs of the lentil draft genome. Up to 8.8% of the lentil contigs evidenced the presence of at least one LTR-RT similar sequence. A statistical analysis suggested a non-random distribution of these elements within of the lentil genome. In most cases (between 97% and 72%, depending on the LTR-RT type) none of the internal sequences flanked by the LTR sequence pair was detected, suggesting that defective and non-autonomous LTR-RTs are very frequent in lentil. Results support that LTR-RTs are abundant and widespread throughout of the lentil genome and that they are a suitable source of genetic markers useful to carry out further genetic analyses.

The Drosophila early ovarian transcriptome provides insight to the molecular causes of recombination rate variation across genomes

BackgroundEvidence in yeast indicates that gene expression is correlated with recombination activity and double-strand break (DSB) formation in some hotspots. Studies of nucleosome occupancy in yeast and mice also suggest that open chromatin influences the formation of DSBs. In Drosophila melanogaster, high-resolution recombination maps show an excess of DSBs within annotated transcripts relative to intergenic sequences. The impact of active transcription on recombination landscapes, however, remains unexplored in a multicellular organism. We then investigated the transcription profile during early meiosis in D. melanogaster females to obtain a glimpse at the relevant transcriptional dynamics during DSB formation, and test the specific hypothesis that DSBs preferentially target transcriptionally active genomic regions.ResultsOur study of transcript profiles of early- and late-meiosis using mRNA-seq revealed, 1) significant differences in gene expression, 2) new genes and exons, 3) parent-of-origin effects on transcription in early-meiosis stages, and 4) a nonrandom genomic distribution of transcribed genes. Importantly, genomic regions that are more actively transcribed during early meiosis show higher rates of recombination, and we ruled out DSB preference for genic regions that are not transcribed.ConclusionsOur results provide evidence in a multicellular organism that transcription during the initial phases of meiosis increases the likelihood of DSB and give insight into the molecular determinants of recombination rate variation across the D. melanogaster genome. We propose that a model where variation in gene expression plays a role altering the recombination landscape across the genome could provide a molecular, heritable and plastic mechanism to observed patterns of recombination variation, from the high level of intra-specific variation to the known influence of environmental factors and stress conditions.

Non-Random Lentiviral Vector Insertions in Bone Marrow Progenitors from Fanconi Anemia Patients

Improved protocols using lentiviral vectors have been established with minimal cytokine exposure and short transduction times proving more suitable for overcoming the disease-specific challenge in correcting functionally defective hematopoietic stem cells (HSCs) of Fanconi Anemia (FA) patients. Bone marrow (BM) cells from FA patients were transduced ex vivo with lentiviral vectors (LVs) expressing FANCA and/or EGFP using optimized conditions to preserve the repopulating properties of the primitive hematopoietic stem cells (manuscript submitted). In a forward preclinical screening of possible LV-induced side effects we analyzed the insertional inventory in colonies generated by FA BM cells previously transduced with the LVs. We have established and optimized DNA and RNA isolation procedures for minimal cell numbers, suitable for large scale screening of colony forming cell (CFC) derived colonies by linear amplification-mediated PCR (LAM-PCR) and massive parallel pyrosequencing (454 GS Flx system; Roche). This approach is applicable for detecting early indicators of clonal selection, and is based on the analysis of common integration sites (CIS) and non-random distribution of vector insertions in particular genomic loci. From a total of 180 CFC-derived colonies expressing the EGFP LV marker gene, 298 vector insertions could be sequenced and mapped to the human genome. The analysis of vector targeted gene coding regions showed a non-random genomic distribution of LV insertions, with a significant overrepresentation of RefSeq genes that are part of distinct functional categories. Accordingly vector associated genes are predominantly involved in cellular signal cascades regulated by the MAP Kinase family known to be involved in a wide variety of cellular processes such as proliferation, differentiation, transcription regulation and development. Apart from the observed high integration frequency in genes (>80%), partial loss of vector LTR nucleotides was detected in >10% of the integrants (3–25bp). Notably, >20% of the lentiviral insertions were found to be located in CIS of predominantly 2nd order. Further screening assays of LV transduced CFC-derived colonies will allow a deeper investigation in the functional consequences of such CIS targeting in gene therapy protocols of FA. However our results suggest that the LV transduction of FA BM progenitors leads to a relatively high frequency of insertions in CIS which may be indicative of an insertion based (specific) selection mechanism. We herby show that the ex vivo large scale integration site analyses of CFC-derived colonies from patients considered to undergo gene therapeutic treatments constitutes a robust approach, which combined with mouse preclinical biosafety studies will help to improve the safety of clinical gene therapy protocols. The non-random distribution of LV integrations in CIS associated genes and in genes involved in particular cellular pathways may be indicative for the altered biochemical pathways characteristic of FA stem cells, with reported defects in DNA repair and self-renewal.

The Subtelomere of Oryza sativa Chromosome 3 Short Arm as a Hot Bed of New Gene Origination in Rice

Despite general observations of non-random genomic distribution of new genes, it is unclear whether or not new genes preferentially occur in certain genomic regions driven by related molecular mechanisms. Using 1.5 Mb of genomic sequences from short arms of chromosome 3 of Oryza glaberrima and O. punctata, we conducted a comparative genomic analysis with the reference O. sativa ssp. japonica genome. We identified a 60-kb segment located in the middle of the subtelomeric region of chromosome 3, which is unique to the species O. sativa. The region contained gene duplicates that occurred in Asian cultivated rice species that diverged from the ancestor of Asian and African cultivated rice one million years ago (MYA). For the 12 genes and one complete retrotransposon identified in this segment in O. sativa ssp. japonica, we searched for their parental genes. The high similarity between duplicated paralogs further supports the recent origination of these genes. We found that this segment was recently generated through multiple independent gene recombination and transposon insertion events. Among the 12 genes, we found that five had chimeric gene structures derived from multiple parental genes. Nine out of the 12 new genes seem to be functional, as suggested by Ka/Ks analysis and the presence of cDNA and/or MPSS data. Furthermore, for the eight transcribed genes, at least two genes could be classified as defense or stress response-related genes. Given these findings, and the fact that subtelomeres are associated with high rates of recombination and transcription, it is likely that subtelomeres may facilitate gene recombination and transposon insertions and serve as hot spots for new gene origination in rice genomes.

Nonrandom Representation of Sex-Biased Genes on Chicken Z Chromosome

Several lines of evidence suggest that the X chromosome of various animal species has an unusual complement of genes with sex-biased or sex-specific expression. However, the study of the X chromosome gene content in different organisms provided conflicting results. The most striking contrast concerns the male-biased genes, which were reported to be almost depleted from the X chromosome in Drosophila but overrepresented on the X chromosome in mammals. To elucidate the reason for these discrepancies, we analysed the gene content of the Z chromosome in chicken. Our analysis of the publicly available expressed sequence tags (EST) data and genome draft sequence revealed a significant underrepresentation of ovary-specific genes on the chicken Z chromosome. For the brain-expressed genes, we found a significant enrichment of male-biased genes but an indication of underrepresentation of female-biased genes on the Z chromosome. This is the first report on the nonrandom gene content in a homogametic sex chromosome of a species with heterogametic female individuals. Further comparison of gene contents of the independently evolved X and Z sex chromosomes may offer new insight into the evolutionary processes leading to the nonrandom genomic distribution of sex-biased and sex-specific genes.

NONRANDOM DISTRIBUTION OF GENES WITH SEX-BIASED EXPRESSION IN THE CHICKEN GENOME

Evolutionary theory predicts that sexually antagonistic genes should show a nonrandom genomic distribution with sex chromosomes usually being enriched for such genes. However, empirical observations from model organisms (Drosophila melanogaster, Caenorhabditis elegans, mammals) on the genomic location of genes with sex-biased expression have provided conflicting data and are not easily explained by a unified framework based on standard models of the evolution of sexually antagonistic genes. Previous studies have been confined to organisms with male heterogamety, meaning that effects related to homo- or heterozygosity of sex chromosomes cannot be separated from effects related to sex-specific characteristics. We therefore studied the genomic distribution of genes with sex-biased expression in the chicken, that is, in an organism with female heterogamety (males ZZ, females ZW). From the abundance of transcripts in expressed sequence tag libraries, we found an underrepresentation of female-specific genes (germ line and somatic tissue) and an overrepresentation of male-specific genes (somatic) on the Z chromosome. This is consistent with theoretical predictions only if mutations beneficial to one sex generally tend to be at least partly dominant (h > 0.5). We also note that sexual selection for a male-biased trait is facilitated by Z-linkage, because sons in organisms with female heterogamety will always inherit a Z chromosome from their fathers.

NONRANDOM DISTRIBUTION OF GENES WITH SEX-BIASED EXPRESSION IN THE CHICKEN GENOME

Evolutionary theory predicts that sexually antagonistic genes should show a nonrandom genomic distribution with sex chromosomes usually being enriched for such genes. However, empirical observations from model organisms (Drosophila melanogaster, Caenorhabditis elegans, mammals) on the genomic location of genes with sex-biased expression have provided conflicting data and are not easily explained by a unified framework based on standard models of the evolution of sexually antagonistic genes. Previous studies have been confined to organisms with male heterogamety, meaning that effects related to homo- or heterozygosity of sex chromosomes cannot be separated from effects related to sex-specific characteristics. We therefore studied the genomic distribution of genes with sex-biased expression in the chicken, that is, in an organism with female heterogamety (males ZZ, females ZW). From the abundance of transcripts in expressed sequence tag libraries, we found an underrepresentation of female-specific genes (germ line and somatic tissue) and an overrepresentation of male-specific genes (somatic) on the Z chromosome. This is consistent with theoretical predictions only if mutations beneficial to one sex generally tend to be at least partly dominant (h > 0.5). We also note that sexual selection for a male-biased trait is facilitated by Z-linkage, because sons in organisms with female heterogamety will always inherit a Z chromosome from their fathers.

Insertion bias and purifying selection of retrotransposons in the Arabidopsis thaliana genome

An analysis of the distribution and age of LTR retrotransposons in the Arabidopsis genome revealed that Pseudoviridae insert randomly along the chromosome and have been recently active whereas Metaviridae were more active in the past and preferentially target heterochromatin.

Transposable elements and the evolution of genome organization in mammals.

All mammalian transposable elements characterized to date appear to be nonrandomly distributed in the mammalian genome. While no element has been found to be exclusively restricted in its chromosomal location, LINE elements and some retrovirus-like elements are preferentially accumulated in G-banding regions of the chromosomes, and in some cases in the sex chromosomes, while SINE elements occur preferentially in R-banding regions. Four mechanisms are presented which may explain the nonrandom genomic distribution of mammalian transposons: i) sequence-specific insertion, ii) S-phase insertion, iii) ectopic excision, and iv) recombinational editing. Some of the available data are consistent with each of these four models, but no single model is sufficient to explain all of the existing data.