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 multi-copy Mega-NUMTs in the human genome

The maternal mode of mitochondrial DNA (mtDNA) inheritance is central to human genetics. Recently, evidence for bi-parental inheritance of mtDNA was claimed for individuals of three pedigrees that suffered mitochondrial disorders. We sequenced mtDNA using both direct Sanger and Massively Parallel Sequencing in several tissues of eleven maternally related and other affiliated healthy individuals of a family pedigree and observed mixed mitotypes in eight individuals. Cells without nuclear DNA, i.e. thrombocytes and hair shafts, only showed the mitotype of haplogroup (hg) V. Skin biopsies were prepared to generate ρ° cells void of mtDNA, sequencing of which resulted in a hg U4c1 mitotype. The position of the Mega-NUMT sequence was determined by fluorescence in situ hybridization and two different quantitative PCR assays were used to determine the number of contributing mtDNA copies. Thus, evidence for the presence of repetitive, full mitogenome Mega-NUMTs matching haplogroup U4c1 in various tissues of eight maternally related individuals was provided. Multi-copy Mega-NUMTs mimic mixtures of mtDNA that cannot be experimentally avoided and thus may appear in diverse fields of mtDNA research and diagnostics. We demonstrate that hair shaft mtDNA sequencing provides a simple but reliable approach to exclude NUMTs as source of misleading results.     

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.     

Evidence for 26 distinct acyl-coenzyme A synthetase genes in the human genome

Acyl-coenzyme A synthetases (ACSs) catalyze the fundamental, initial reaction in fatty acid metabolism. "Activation" of fatty acids by thioesterification to CoA allows their participation in both anabolic and catabolic pathways. The availability of the sequenced human genome has facilitated the investigation of the number of ACS genes present. Using two conserved amino acid sequence motifs to probe human DNA databases, 26 ACS family genes/proteins were identified. ACS activity in either humans or rodents was demonstrated previously for 20 proteins, but 6 remain candidate ACSs. For two candidates, cDNA was cloned, protein was expressed in COS-1 cells, and ACS activity was detected. Amino acid sequence similarities were used to assign enzymes into subfamilies, and subfamily assignments were consistent with acyl chain length preference. Four of the 26 proteins did not fit into a subfamily, and bootstrap analysis of phylograms was consistent with evolutionary divergence. Three additional conserved amino acid sequence motifs were identified that likely have functional or structural roles. The existence of many ACSs suggests that each plays a unique role, directing the acyl-CoA product to a specific metabolic fate. Knowing the full complement of ACS genes in the human genome will facilitate future studies to characterize their specific biological functions.     

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.     

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.     

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.     

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.