The genesis and evolution of homeobox gene clusters

Once called the 'Rosetta stone' of developmental biology, the homeobox continues to fascinate both evolutionary and developmental biologists. The birth of the homeotic, or Hox, gene cluster, and its subsequent evolution, has been crucial in mediating the major transitions in metazoan body plan. Comparative genomics studies indicate that the more recently discovered ParaHox and NK clusters were linked to the Hox cluster early in evolution, and that together they constituted a 'megacluster' of homeobox genes that conspicuously contributed to body-plan evolution.     

Conservation of gene linkage in dispersed vertebrate NK homeobox clusters

Nk homeobox genes are important regulators of many different developmental processes including muscle, heart, central nervous system and sensory organ development. They are thought to have arisen as part of the ANTP megacluster, which also gave rise to Hox and ParaHox genes, and at least some NK genes remain tightly linked in all animals examined so far. The protostome–deuterostome ancestor probably contained a cluster of nine Nk genes: (Msx)–(Nk4/tinman)–(Nk3/bagpipe)–(Lbx/ladybird)–(Tlx/c15)–(Nk7)–(Nk6/hgtx)–(Nk1/slouch)–(Nk5/Hmx). Of these genes, only NKX2.6–NKX3.1, LBX1–TLX1 and LBX2–TLX2 remain tightly linked in humans. However, it is currently unclear whether this is unique to the human genome as we do not know which of these Nk genes are clustered in other vertebrates. This makes it difficult to assess whether the remaining linkages are due to selective pressures or because chance rearrangements have “missed” certain genes. In this paper, we identify all of the paralogs of these ancestrally clustered NK genes in several distinct vertebrates. We demonstrate that tight linkages of Lbx1–Tlx1, Lbx2–Tlx2 and Nkx3.1–Nkx2.6 have been widely maintained in both the ray-finned and lobe-finned fish lineages. Moreover, the recently duplicated Hmx2–Hmx3 genes are also tightly linked. Finally, we show that Lbx1–Tlx1 and Hmx2–Hmx3 are flanked by highly conserved noncoding elements, suggesting that shared regulatory regions may have resulted in evolutionary pressure to maintain these linkages. Consistent with this, these pairs of genes have overlapping expression domains. In contrast, Lbx2–Tlx2 and Nkx3.1–Nkx2.6, which do not seem to be coexpressed, are also not associated with conserved noncoding sequences, suggesting that an alternative mechanism may be responsible for the continued clustering of these genes.     

Evidence for archaic Asian ancestry on the human X chromosome

The human RRM2P4 pseudogene has a pattern of nucleotide polymorphism that is unlike any locus published to date. A gene tree constructed from a 2.4-kb fragment of the RRM2P4 locus sequenced in a sample of 41 worldwide humans clearly roots in East Asia and has a most-recent common ancestor approximately 2 Myr before present. The presence of this basal lineage exclusively in Asia results in higher nucleotide diversity among non-Africans than among Africans. A global survey of a single-nucleotide polymorphism that is diagnostic for the basal, Asian lineage in 570 individuals shows that it occurs at frequencies up to 53% in south China, whereas only one of 177 surveyed Africans carries this archaic lineage. We suggest that this ancient lineage is a remnant of introgressive hybridization between expanding anatomically modern humans emerging from Africa and archaic populations in Eurasia.     

Evidence for extensive transmission distortion in the human genome

It is a basic principle of genetics that each chromosome is transmitted from parent to offspring with a probability that is given by Mendel's laws. However, several known biological processes lead to skewed transmission probabilities among surviving offspring and, therefore, to excess genetic sharing among relatives. Examples include in utero selection against deleterious mutations, meiotic drive, and maternal-fetal incompatibility. Although these processes affect our basic understanding of inheritance, little is known about their overall impact in humans or other mammals. In this study, we examined genome screen data from 148 nuclear families, collected without reference to phenotype, to look for departures from Mendelian transmission proportions. Using single-point and multipoint linkage analysis, we detected a modest but significant genomewide shift towards excess genetic sharing among siblings (average sharing of 50.43% for the autosomes; P=.009). Our calculations indicate that many loci with skewed transmission are required to produce a genomewide shift of this magnitude. Since transmission distortion loci are subject to strong selection, this raises interesting questions about the evolutionary forces that keep them polymorphic. Finally, our results also have implications for mapping disease genes and for the genetics of fertility.     

Widespread allele-specific topological domains in the human genome are not confined to imprinted gene clusters

Neural networks such as variational autoencoders (VAE) perform dimensionality reduction for the visualization and analysis of genomic data, but are limited in their interpretability: it is unknown which data features are represented by each embedding dimension. We present siVAE, a VAE that is interpretable by design, thereby enhancing downstream analysis tasks. Through interpretation, siVAE also identifies gene modules and hubs without explicit gene network inference. We use siVAE to identify gene modules whose connectivity is associated with diverse phenotypes such as iPSC neuronal differentiation efficiency and dementia, showcasing the wide applicability of interpretable generative models for genomic data analysis.

     

Molecular cloning of a human Vent-like homeobox gene

We have isolated a previously unknown human homeobox-containing cDNA, VENT-like homeobox-2 (VENTX2), using PCR with a bone marrow cDNA library and primers designed from the VENTX1 (alias HPX42) homeobox sequence. Here we describe the molecular cloning, chromosomal localization to 10q26.3, and functional analysis of this gene. The 2.4-kb human VENTX2 cDNA encoded a protein with a predicted molecular weight of 28 kDa containing a homeodomain with 65% identity to the Xenopus laevis ventralizing gene Xvent2B. VENTX2 antisera detected a 28-kDa protein in cells transfected with a VENTX2 expression construct, in a human erythroleukemic cell line and in bone marrow samples obtained from patients in recovery phase after chemotherapy. The similarity of the homeodomains from VENTX2 and the X. laevis Vent gene family places them in the same homeodomain class. Consistent with this structural classification, overexpression of VENTX2 in zebrafish embryos led to anterior truncations and failure to form a notochord, which are characteristics of ventralization.     

Evolutionary genomics of the Fox genes: Origin of gene families and the ancestry of gene clusters

Over the past decade genomic approaches have begun to revolutionise the study of animal diversity. In particular, genome sequencing programmes have spread beyond the traditional model species to encompass an increasing diversity of animals from many different phyla, as well as unicellular eukaryotes that are closely related to the animals. Whole genome sequences allow researchers to establish, with reasonable confidence, the full complement of any particular family of genes in a genome. Comparison of gene complements from appropriate genomes can reveal the evolutionary history of gene families, indicating when both gene diversification and gene loss have occurred. More than that, however, assembled genomes allow the genomic environment in which individual genes are found to be analysed and compared between species. This can reveal how gene diversification occurred. Here, we focus on the Fox genes, drawing from multiple animal genomes to develop an evolutionary framework explaining the timing and mechanism of origin of the diversity of animal Fox genes. Ancient linkages between genes are a prominent feature of the Fox genes, depicting a history of gene clusters, some of which may be relevant to understanding Fox gene function.     

Evidence for hitchhiking of deleterious mutations within the human genome

Deleterious mutations present a significant obstacle to adaptive evolution. Deleterious mutations can inhibit the spread of linked adaptive mutations through a population; conversely, adaptive substitutions can increase the frequency of linked deleterious mutations and even result in their fixation. To assess the impact of adaptive mutations on linked deleterious mutations, we examined the distribution of deleterious and neutral amino acid polymorphism in the human genome. Within genomic regions that show evidence of recent hitchhiking, we find fewer neutral but a similar number of deleterious SNPs compared to other genomic regions. The higher ratio of deleterious to neutral SNPs is consistent with simulated hitchhiking events and implies that positive selection eliminates some deleterious alleles and increases the frequency of others. The distribution of disease-associated alleles is also altered in hitchhiking regions. Disease alleles within hitchhiking regions have been associated with auto-immune disorders, metabolic diseases, cancers, and mental disorders. Our results suggest that positive selection has had a significant impact on deleterious polymorphism and may be partly responsible for the high frequency of certain human disease alleles.     

Proteorhodopsin photosystem gene clusters exhibit co-evolutionary trends and shared ancestry among diverse marine microbial phyla

Since the recent discovery of retinylidene proteins in marine bacteria (proteorhodopsins), the estimated abundance and diversity of this gene family has expanded rapidly. To explore proteorhodopsin photosystem evolutionary and distributional trends, we identified and compared 16 different proteorhodopsin-containing genome fragments recovered from naturally occurring bacterioplankton populations. In addition to finding several deep-branching proteorhodopsin sequences, proteorhodopsins were found in novel taxonomic contexts, including a betaproteobacterium and a planctomycete. Approximately one-third of the proteorhodopsin-containing genome fragments analysed, as well as a number of recently reported marine bacterial whole genome sequences, contained a linked set of genes required for biosynthesis of the rhodopsin chromophore, retinal. Phylogenetic analyses of the retinal biosynthetic genes suggested their co-evolution and probable coordinated lateral gene transfer into disparate lineages, including Euryarchaeota, Planctomycetales, and three different proteobacterial lineages. The lateral transfer and retention of genes required to assemble a functional proteorhodopsin photosystem appears to be a coordinated and relatively frequent evolutionary event. Strong selection pressure apparently acts to preserve these light-dependent photosystems in diverse marine microbial lineages.     

Study of polymorphism in human rRNA gene clusters

Human ribosomal DNA (rDNA) probe specific for the 3' end of the 28S rRNA gene was used for detecting standard restriction fragments' length polymorphism (RFLPs) in the non-transcribed spacer. The conditions for hybridization of rDNA probe which eliminate cross hybridization of parts of 28S rRNA gene were developed. A test for detecting incompletely restricted DNA was also developed which may be used in experiments for detecting new RFLPs. It was found that a set of standard RFLPs was identical in various human tissues for one individual. Frequency of standard RFLPs in the non-transcribed spacer of human rRNA gene clusters was calculated.     

Finding keywords for intergenic and gene regions for human genome

The analysis of functionally related sequences for conserved patterns is important for further research of different functional regions. This paper presents an analysis of genes and intergenic sequences from the point of view of linguistics analysis, where gene and intergenic regions are regarded as two different subjects written in the four-letter alphabet [A, C, G, T] and high-frequency simple sequences are taken as keywords. A measurement alpha[l(tau)] was introduced to describe the relative repeat ratio of simple sequences. Cutoff values were found for keywords selection. After eliminating "noise," 87 short sequences were selected as keywords for intergenic regions and 76 for gene regions.     

Comparative genome sequencing reveals chemotype-specific gene clusters in the toxigenic black mold Stachybotrys

Background

The fungal genus Stachybotrys produces several diverse toxins that affect human health. Its strains comprise two mutually-exclusive toxin chemotypes, one producing satratoxins, which are a subclass of trichothecenes, and the other producing the less-toxic atranones. To determine the genetic basis for chemotype-specific differences in toxin production, the genomes of four Stachybotrys strains were sequenced and assembled de novo. Two of these strains produce atranones and two produce satratoxins.

Results

Comparative analysis of these four 35-Mbp genomes revealed several chemotype-specific gene clusters that are predicted to make secondary metabolites. The largest, which was named the core atranone cluster, encodes 14 proteins that may suffice to produce all observed atranone compounds via reactions that include an unusual Baeyer-Villiger oxidation. Satratoxins are suggested to be made by products of multiple gene clusters that encode 21 proteins in all, including polyketide synthases, acetyltransferases, and other enzymes expected to modify the trichothecene skeleton. One such satratoxin chemotype-specific cluster is adjacent to the core trichothecene cluster, which has diverged from those of other trichothecene producers to contain a unique polyketide synthase.

Conclusions

The results suggest that chemotype-specific gene clusters are likely the genetic basis for the mutually-exclusive toxin chemotypes of Stachybotrys. A unified biochemical model for Stachybotrys toxin production is presented. Overall, the four genomes described here will be useful for ongoing studies of this mold’s diverse toxicity mechanisms.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-590) contains supplementary material, which is available to authorized users.     

A Gbx homeobox gene in amphioxus: insights into ancestry of the ANTP class and evolution of the midbrain/hindbrain boundary

In the vertebrate central nervous system (CNS), mutual antagonism between posteriorly expressed Gbx2 and anteriorly expressed Otx2 positions the midbrain/hindbrain boundary (MHB), but does not induce MHB organizer genes such as En, Pax2/5/8 and Wnt1. In the CNS of the cephalochordate amphioxus, Otx is also expressed anteriorly, but En, Pax2/5/8 and Wnt1 are not expressed near the caudal limit of Otx, raising questions about the existence of an MHB organizer in amphioxus. To investigate the evolutionary origins of the MHB, we cloned the single amphioxus Gbx gene. Fluorescence in situ hybridization showed that, as in vertebrates, amphioxus Gbx and the Hox cluster are on the same chromosome. From analysis of linked genes, we argue that during evolution a single ancestral Gbx gene duplicated fourfold in vertebrates, with subsequent loss of two duplicates. Amphioxus Gbx is expressed in all germ layers in the posterior 75% of the embryo, and in the CNS, the Gbx and Otx domains abut at the boundary between the cerebral vesicle (forebrain/midbrain) and the hindbrain. Thus, the genetic machinery to position the MHB was present in the protochordate ancestors of the vertebrates, but is insufficient for induction of organizer genes. Comparison with hemichordates suggests that anterior Otx and posterior Gbx domains were probably overlapping in the ancestral deuterostome and came to abut at the MHB early in the chordate lineage before MHB organizer properties evolved.     

Hox clusters as models for vertebrate genome evolution

The surprising variation in the number of Hox clusters and the genomic architecture within vertebrate lineages, especially within the ray-finned fish, reflects a history of duplications and subsequent lineage-specific gene loss. Recent research on the evolution of conserved non-coding sequences (CNS) in Hox clusters promises to reveal interesting results for functional and phenotypic diversification.     

The LIM homeobox gene ceh-14 confers thermosensory function to the AFD neurons in Caenorhabditis elegans

In Caenorhabditis elegans three pairs of neurons, AFD, AIY, and AIZ, play a key role in thermosensation. The LIM homeobox gene ceh-14 is expressed in the AFD thermosensory neurons. ceh-14 mutant animals display athermotactic behaviors, although the neurons are still present and differentiated. Two other LIM homeobox genes, ttx-3 and lin-11, function in the two interneurons AIY and AIZ, respectively. Thus, the three key thermosensory neurons are specified by three different LIM homeobox genes. ceh-14 ttx-3 lin-11 triple mutant animals have a basic cryophilic thermotaxis behavior indicative of a second thermotaxis pathway. Misexpression of ceh-14 in chemosensory neurons can restore thermotactic behavior without impairing the chemosensory function. Thus, ceh-14 confers thermosensory function to neurons.     

Organization and evolution of mitochondrial gene clusters in human

Currently, the spatial patterns of mitochondrial genes and how the genomic localization of (pseudo)genes originated from mitochondrial DNA remain largely unexplained. The aim of this study was to elucidate the organization of mitochondrial (pseudo)genes given their evolutionary origin. We used a keyword finding method and a bootstrapping method to estimate parameter values that represent the distribution pattern of mitochondrial genes in the nuclear genome. Almost half of mitochondrial genes showing physical clusters were located in the pericentromeric and subtelomeric regions of the chromosome. Most interestingly, the size of these clusters ranged from 0.085 to 3.2 Mb (average ± SD 1.3 ± 0.73 Mb), which coincides with the size of the evolutionary pocket, or the average size of evolutionary breakpoint regions. Our findings imply that the localization of mitochondrial genes in the human genome is determined independent of adaptation.     

Structure and sequence of the human homeobox gene HOX7.

A cosmid containing the human sequence HOX7, homologous to the murine Hox-7 gene, was isolated from a genomic library, and the positions of the coding sequences were determined by hybridization. DNA sequence analysis demonstrated two exons that code for a homeodomain-containing protein of 297 amino acids. The open reading frame is interrupted by a single intron of approximately 1.6 kb, the splice donor and acceptor sites of which conform to known consensus sequences. The human HOX7 coding sequence has a very high degree of identity with the murine Hox-7 cDNA. Within the homeobox, the two sequences share 94% identity at the DNA level, all substitutions being silent. This high level of sequence similarity is not confined to the homeodomain; overall the human and murine HOX7 gene products show 80% identity at the amino acid level. Both the 5' and 3' untranslated regions also show significant similarity to the murine gene, with 79 and 70% sequence identity, respectively. The sequence upstream of the coding sequence of exon 1 contains a GC-rich putative promoter region. There is no TATA box, but a CCAAT and numerous GC boxes are present. The region encompassing the promoter region, exon 1, and the 5' region of exon 2 have a higher than expected frequency of CpG dinucleotides; numerous sites for rare-cutter restriction enzymes are present, a characteristic of HTF islands.