Whole genome duplications and expansion of the vertebrate GATA transcription factor gene family

Genome-Wide Identification and Expression Analysis of GATA Family Genes in Dimocarpus longan Lour

GATA transcription factors, which are DNA-binding proteins with type IV zinc finger binding domains, have a role in transcriptional regulation in biological organisms. They have an indispensable role in the growth and development of plants, as well as in improvements in their ability to face various environmental stresses. To date, GATAs have been identified in many gene families, but the GATA gene in longan (Dimocarpus longan Lour) has not been studied in previous explorations. Various aspects of genes in the longan GATA family, including their identification and classification, the distribution of their positions on chromosomes, their exon/intron structures, a synteny analysis, their expression at different temperatures, concentration of PEG, early developmental stages of somatic embryos and their expression levels in different tissues, and concentrations of exogenous hormones, were investigated in this study. This study showed that the 22 DlGATAs could be divided into four subfamilies. There were 10 pairs of homologous GATA genes in the synteny analysis of DlGATA and AtGATA. Four segmental replication motifs and one pair of tandem duplication events were present among the DlGATA family members. The cis-acting elements located in promoter regions were also found to be enriched with light-responsive elements, which contained related hormone-responsive elements. In somatic embryos, DlGATA4 is upregulated for expression at the globular embryo (GE) stage. We also found that DlGATA expression was strongly up-regulated in roots and stems. The study demonstrated the expression of DlGATA under hormone (ABA and IAA) treatments in embryogenic callus of longan. Under ABA treatment, DlGATA4 was up-regulated and the other DlGATA genes did not respond significantly. Moreover, as demonstrated with qRT-PCR, the expression of DlGATA genes showed strong up-regulated expression levels under 100 μmol·L-1 concentration IAA treatment. This experiment further studied these and simulated their possible connections with a drought response mechanism, while correlating them with their expression under PEG treatment. Overall, this experiment explored the GATA genes and dug into their evolution, structure, function, and expression profile, thus providing more information for a more in-depth study of the characteristics of the GATA family of genes.

GATA1 in Normal and Pathologic Megakaryopoiesis and Platelet Development

GATA1 is a highly conserved hematopoietic transcription factor (TF), essential for normal erythropoiesis and megakaryopoiesis, that encodes a full-length, predominant isoform and an amino (N) terminus-truncated isoform GATA1s. It is consistently expressed throughout megakaryocyte development and interacts with its target genes either independently or in association with binding partners such as FOG1 (friend of GATA1). While the N-terminus and zinc finger have classically been demonstrated to be necessary for the normal regulation of platelet-specific genes, murine models, cell-line studies, and human case reports indicate that the carboxy-terminal activation domain and zinc finger also play key roles in precisely controlling megakaryocyte growth, proliferation, and maturation. Murine models have shown that disruptions to GATA1 increase the proliferation of immature megakaryocytes with abnormal architecture and impaired terminal differentiation into platelets. In humans, germline GATA1 mutations result in variable cytopenias, including macrothrombocytopenia with abnormal platelet aggregation and excessive bleeding tendencies, while acquired GATA1s mutations in individuals with trisomy 21 (T21) result in transient abnormal myelopoiesis (TAM) and myeloid leukemia of Down syndrome (ML-DS) arising from a megakaryocyte-erythroid progenitor (MEP). Taken together, GATA1 plays a key role in regulating megakaryocyte differentiation, maturation, and proliferative capacity. As sequencing and proteomic technologies expand, additional GATA1 mutations and regulatory mechanisms contributing to human diseases of megakaryocytes and platelets are likely to be revealed.

Genome-wide identification of CaARR-Bs transcription factor gene family in pepper and their expression patterns under salinity stress

In plants, ARRs-B transcription factors play a crucial role in regulating cytokinin signal transduction, abiotic stress resistance, and plant development. A number of adverse environmental conditions have caused severe losses for the pepper (Capsicum annuum L.)-a significant and economically important vegetable. Among the transcription factors of the type B-ARRs family, multiple members have different functions. In pepper, only a few members of the ARRs-B family have been reported and characterized. The current study aimed to characterize ARRs-B transcription factors in C. annuum, including phylogenetic relationships, gene structures, protein motif arrangement, and RT-qPCR expression analyses and their role in salinity stress. In total, ten genes encode CaARRs-B transcription factors (CaARR1 to CaARR10) from the largest subfamily of type-B ARRs were identified in C. annum. The genome-wide analyses of the CaARRs-B family in C. annuum were performed based on the reported ARRs-B genes in Arabidopsis. An analysis of homologous alignments of candidate genes, including their phylogenetic relationships, gene structures, conserved domains, and qPCR expression profiles, was conducted. In comparison with other plant ARRs-B proteins, CaARRs-B proteins showed gene conservation and potentially specialized functions. In addition, tissue-specific expression profiles showed that CaARRs-B genes were differentially expressed, suggesting functionally divergent. CaARRs-B proteins had a typical conserved domain, including AAR-like (pfam: PF00072) and Myb DNA binding (pfam: PF00249) domains. Ten of the CaARRs-B genes were asymmetrically mapped on seven chromosomes in Pepper. Additionally, the phylogenetic tree of CaARRs-B genes from C. annuum and other plant species revealed that CaARRs-B genes were classified into four clusters, which may have evolved conservatively. Further, using quantitative real-time qRT-PCR, the study assessed the expression patterns of CaARRs-B genes in Capsicum annuum seedlings subjected to salt stress. The study used quantitative real-time qRT-PCR to examine CaARRs-B gene expression in Capsicum annuum seedlings under salt stress. Roots exhibited elevated expression of CaARR2 and CaARR9, while leaves showed decreased expression for CaARR3, CaARR4, CaARR7, and CaARR8. Notably, no amplification was observed for CaARR10. This research sheds light on the roles of CaARRs-B genes in pepper's response to salinity stress. These findings enrich our comprehension of the functional implications of CaARRs-B genes in pepper, especially in responding to salinity stress, laying a solid groundwork for subsequent in-depth studies and applications in the growth and development of Capsicum annuum.

Genome-Wide Identification of GATA Family Genes in Phoebe bournei and Their Transcriptional Analysis under Abiotic Stresses

GATA transcription factors are crucial proteins in regulating transcription and are characterized by a type-IV zinc finger DNA-binding domain. They play a significant role in the growth and development of plants. While the GATA family gene has been identified in several plant species, it has not yet been reported in Phoebe bournei. In this study, 22 GATA family genes were identified from the P. bournei genome, and their physicochemical properties, chromosomal distribution, subcellular localization, phylogenetic tree, conserved motif, gene structure, cis-regulatory elements in promoters, and expression in plant tissues were analyzed. Phylogenetic analysis showed that the PbGATAs were clearly divided into four subfamilies. They are unequally distributed across 11 out of 12 chromosomes, except chromosome 9. Promoter cis-elements are mostly involved in environmental stress and hormonal regulation. Further studies showed that PbGATA11 was localized to chloroplasts and expressed in five tissues, including the root bark, root xylem, stem bark, stem xylem, and leaf, which means that PbGATA11 may have a potential role in the regulation of chlorophyll synthesis. Finally, the expression profiles of four representative genes, PbGATA5, PbGATA12, PbGATA16, and PbGATA22, under drought, salinity, and temperature stress, were detected by qRT-PCR. The results showed that PbGATA5, PbGATA22, and PbGATA16 were significantly expressed under drought stress. PbGATA12 and PbGATA22 were significantly expressed after 8 h of low-temperature stress at 10 °C. This study concludes that the growth and development of the PbGATA family gene in P. bournei in coping with adversity stress are crucial. This study provides new ideas for studying the evolution of GATAs, provides useful information for future functional analysis of PbGATA genes, and helps better understand the abiotic stress response of P. bournei.

Long noncoding RNA GATA2-AS1 augments endothelial Hypoxia Inducible Factor 1-α induction and regulates hypoxic signaling

Vascular endothelial cells form the inner cellular lining of blood vessels and have myriad physiologic functions including angiogenesis and response to hypoxia. We recently identified a set of endothelial cell (EC)-enriched long noncoding RNAs (lncRNAs) in differentiated human primary cell types and described the role of the STEEL lncRNA in angiogenic patterning. We sought to further understand the role of EC-enriched lncRNAs in physiologic adaptation of the vascular endothelium. In this work, we describe an abundant, cytoplasmic, and EC-enriched lncRNA, GATA2-AS1, that is divergently transcribed from the EC-enriched developmental regulator, GATA2. While GATA2-AS1 is largely co-expressed with GATA2 in ECs, GATA2-AS1 and GATA2 appear to be complementary rather than synergistic as they have mostly distinct target genes. Common single nucleotide variants (SNVs) in GATA2-AS1 exons are associated with early onset coronary artery disease (CAD) and decreased expression of GATA2-AS1 in endothelial cell lines. In most cells, HIF1-α is central to the transcriptional response to hypoxia, while in ECs, both HIF1-α and HIF2-α are required to coordinate an acute and chronic response respectively. In this setting, GATA2-AS1 contributes to the "HIF switch" and augments HIF1-α induction in acute hypoxia to regulate HIF1-α/ HIF2-α balance. In hypoxia, GATA2-AS1 orchestrates HIF1-α-dependent induction of the glycolytic pathway, and HIF1-α-independent maintenance of mitochondrial biogenesis. Similarly, GATA2-AS1 coordinates both metabolism and "tip/stalk" cell signaling to regulate angiogenesis in hypoxic ECs. Furthermore, we find that GATA2-AS1 expression patterns are perturbed in atherosclerotic disease. Together, these results define a role for GATA2-AS1 in the EC-specific response to hypoxia.

GATA2 deficiency and MDS/AML: Experimental strategies for disease modelling and future therapeutic prospects

The importance of predisposition to leukaemia in clinical practice is being increasingly recognized. This is emphasized by the establishment of a novel WHO disease category in 2016 called "myeloid neoplasms with germline predisposition". A major syndrome within this group is GATA2 deficiency, a heterogeneous immunodeficiency syndrome with a very high lifetime risk to develop myelodysplastic syndrome (MDS) and acute myeloid leukaemia (AML). GATA2 deficiency has been identified as the most common hereditary cause of MDS in adolescents with monosomy 7. Allogenic haematopoietic stem cell transplantation is the only curative option; however, chances of survival decrease with progression of immunodeficiency and MDS evolution. Penetrance and expressivity within families carrying GATA2 mutations is often variable, suggesting that co-operating extrinsic events are required to trigger the disease. Predictive tools are lacking, and intrafamilial heterogeneity is poorly understood; hence there is a clear unmet medical need. On behalf of the ERAPerMed GATA2 HuMo consortium, in this review we describe the genetic, clinical, and biological aspects of familial GATA2-related MDS, highlighting the importance of developing robust disease preclinical models to improve early detection and clinical decision-making of GATA2 carriers.

Genome-Wide Identification and Characterization of the BRD Family in Nile Tilapia (Oreochromis niloticus)

Simple Summary

Nile tilapia is a good model for genome-wide identification and examination of the expression and role of gene families. In this study, we identified 54 bromodomain genes (BRDs) divided into eight subfamilies in Nile tilapia. Phylogenetic analysis revealed a high conservation of the BRDs family in vertebrates, with BRDs expansion due to fish-specific duplications. Most of the BRDs displayed sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and 9 ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testes of treated fish.

Abstract

The bromodomain (BRD) proteins specifically recognize the N-acetyllysine motifs, which is a key event in the reading process of epigenetic marks. BRDs are evolutionarily highly conserved. Over recent years, BRDs attracted great interest because of their important roles in biological processes. However, the genome-wide identification of this family was not carried out in many animal groups, in particular, in teleosts. Moreover, the expression patterns were not reported for any of the members in this family, and the role of the BRD family was not extensively studied in fish reproduction. In this study, we identified 16 to 120 BRD genes in 24 representative species. BRDs expanded significantly in vertebrates. Phylogenetic analysis showed that the BRD family was divided into eight subfamilies (I–VIII). Transcriptome analysis showed that BRDs in Nile tilapia (Oreochromis niloticus) exhibited different expression patterns in different tissues, suggesting that these genes may play different roles in growth and development. Gonadal transcriptome analysis showed that most of the BRDs display sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and nine ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Consistent with transcriptomic data, the results of qRT-PCR and fluorescence in situ hybridization showed that brdt expression was higher in the testis than in the ovary, suggesting its critical role in the spermatogenesis of the tilapia. Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced, and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testis tissue of treated fish. Our data provide insights into the evolution and expression of BRD genes, which is helpful for understanding their critical roles in sex differentiation and gonadal development in teleosts.

Two Isoforms of serpent Containing Either One or Two GATA Zinc Fingers Provide Functional Diversity During Drosophila Development

GATA transcription factors play crucial roles in various developmental processes in organisms ranging from flies to humans. In mammals, GATA factors are characterized by the presence of two highly conserved domains, the N-terminal (N-ZnF) and the C-terminal (C-ZnF) zinc fingers. The Drosophila GATA factor Serpent (Srp) is produced in different isoforms that contains either both N-ZnF and C-ZnF (SrpNC) or only the C-ZnF (SrpC). Here, we investigated the functional roles ensured by each of these isoforms during Drosophila development. Using the CRISPR/Cas9 technique, we generated new mutant fly lines deleted for one (ΔsrpNC) or the other (ΔsrpC) encoded isoform, and a third one with a single point mutation in the N-ZnF that alters its interaction with its cofactor, the Drosophila FOG homolog U-shaped (Ush). Analysis of these mutants revealed that the Srp zinc fingers are differentially required for Srp to fulfill its functions. While SrpC is essential for embryo to adult viability, SrpNC, which is the closest conserved isoform to that of vertebrates, is not. However, to ensure its specific functions in larval hematopoiesis and fertility, Srp requires the presence of both N- and C-ZnF (SrpNC) and interaction with its cofactor Ush. Our results also reveal that in vivo the presence of N-ZnF restricts rather than extends the ability of GATA factors to regulate the repertoire of C-ZnF bound target genes.

Differential Requirement of Gata2a and Gata2b for Primitive and Definitive Myeloid Development in Zebrafish

Germline loss or mutation of one copy of the transcription factor GATA2 in humans leads to a range of clinical phenotypes affecting hematopoietic, lymphatic and vascular systems. GATA2 heterozygous mice show only a limited repertoire of the features observed in humans. Zebrafish have two copies of the Gata2 gene as a result of an additional round of ancestral whole genome duplication. These genes, Gata2a and Gata2b, show distinct but overlapping expression patterns, and between them, highlight a significantly broader range of the phenotypes observed in GATA2 deficient syndromes, than each one alone. In this manuscript, we use mutants for Gata2a and Gata2b to interrogate the effects on hematopoiesis of these two ohnologs, alone and in combination, during development in order to further define the role of GATA2 in developmental hematopoiesis. We define unique roles for each ohnolog at different stages of developmental myelopoiesis and for the emergence of hematopoietic stem and progenitor cells. These effects are not additive in the haploinsufficient state suggesting a redundancy between these two genes in hematopoietic stem and progenitor cells. Rescue studies additionally support that Gata2b can compensate for the effects of Gata2a loss. Finally we show that adults with loss of combined heterozygosity show defects in the myeloid compartment consistent with GATA2 loss in humans. These results build on existing knowledge from other models of GATA2 deficiency and refine our understanding of the early developmental effects of GATA2. In addition, these studies shed light on the complexity and potential structure-function relationships as well as sub-functionalization of Gata2 genes in the zebrafish model.

Cdan1 Is Essential for Primitive Erythropoiesis

Congenital dyserythropoietic anemia type I (CDA I) is an autosomal recessive disease characterized by moderate to severe macrocytic anemia and pathognomonic morphologic abnormalities of the erythroid precursors, including spongy heterochromatin. The disease is mainly caused by mutations in CDAN1 (encoding for Codanin-1). No patients with homozygous null type mutations have been described, and mouse null mutants die during early embryogenesis prior to the initiation of erythropoiesis. The cellular functions of Codanin-1 and the erythroid specificity of the phenotype remain elusive. To investigate the role of Codanin-1 in erythropoiesis, we crossed mice carrying the Cdan1 floxed allele (Cdan^fl/fl) with mice expressing Cre-recombinase under regulation of the erythropoietin receptor promoter (ErGFPcre). The resulting Cdan^ΔEry transgenic embryos died at mid-gestation (E12.5–E13.5) from severe anemia, with very low numbers of circulating erythroblast. Transmission electron microscopy studies of primitive erythroblasts (E9.5) revealed the pathognomonic spongy heterochromatin. The morphology of Cdan^ΔEry primitive erythroblasts demonstrated progressive development of dyserythropoiesis. Annexin V staining showed increases in both early and late-apoptotic erythroblasts compared to controls. Flow cytometry studies using the erythroid-specific cell-surface markers CD71 and Ter119 demonstrated that Cdan^ΔEry erythroid progenitors do not undergo the semi-synchronous maturation characteristic of primitive erythroblasts. Gene expression studies aimed to evaluate the effect of Cdan1 depletion on erythropoiesis revealed a delay of ζ to α globin switch compared to controls. We also found increased expression of Gata2, Pu.1, and Runx1, which are known to inhibit terminal erythroid differentiation. Consistent with this data, our zebrafish model showed increased gata2 expression upon cdan1 knockdown. In summary, we demonstrated for the first time that Cdan1 is required for primitive erythropoiesis, while providing two experimental models for studying the role of Codanin-1 in erythropoiesis and in the pathogenesis of CDA type I.

Genome-wide identification and characterization of 14-3-3 genes in fishes

The 14-3-3 family genes are highly conserved regulatory factors in eukaryotes with involvement in multiple important cellular processes. However, detailed investigations of this family in fishes are very limited. Here, a comparative genomic and transcriptomic survey were performed to investigate the 14-3-3 family in fishes. We confirmed that the numbers of 14-3-3 genes ranged from 5 to 7 in non-teleost fishes, as well as additional 14-3-3 genes (9 to 11) in teleost fishes. In addition, some special teleost fishes possess 17 to 25 14-3-3s, which undergone the fourth whole-genome duplication (WGD). We also found that six pairs of fish 14-3-3 genes were clustered with mammalian ε, γ, ς, η, τand β isotypes, respectively, while σ was absent with a potential specificity within mammals, on the basis of their phylogenetic and synteny analyses. According to our results, we inferred that the diversity of 14-3-3 genes in fishes seems to be generated from a combination of WGD and gene loss. Comparative transcriptomic analysis revealed that there are differences in tissue distribution, and we speculated that 14-3-3 genes may contribute to terrestrial adaptations in mudskippers. In addition, protein sequence alignments of 14-3-3s supported their differential roles in fishes. In summary, our present comparative genomic and transcriptomic survey will benefit for further functional investigations of these fish genes.

Zebrafish Cancer Predisposition Models

Cancer predisposition syndromes are rare, typically monogenic disorders that result from germline mutations that increase the likelihood of developing cancer. Although these disorders are individually rare, resulting cancers collectively represent 5-10% of all malignancies. In addition to a greater incidence of cancer, affected individuals have an earlier tumor onset and are frequently subjected to long-term multi-modal cancer screening protocols for earlier detection and initiation of treatment. In vivo models are needed to better understand tumor-driving mechanisms, tailor patient screening approaches and develop targeted therapies to improve patient care and disease prognosis. The zebrafish (Danio rerio) has emerged as a robust model for cancer research due to its high fecundity, time- and cost-efficient genetic manipulation and real-time high-resolution imaging. Tumors developing in zebrafish cancer models are histologically and molecularly similar to their human counterparts, confirming the validity of these models. The zebrafish platform supports both large-scale random mutagenesis screens to identify potential candidate/modifier genes and recently optimized genome editing strategies. These techniques have greatly increased our ability to investigate the impact of certain mutations and how these lesions impact tumorigenesis and disease phenotype. These unique characteristics position the zebrafish as a powerful in vivo tool to model cancer predisposition syndromes and as such, several have already been created, including those recapitulating Li-Fraumeni syndrome, familial adenomatous polyposis, RASopathies, inherited bone marrow failure syndromes, and several other pathogenic mutations in cancer predisposition genes. In addition, the zebrafish platform supports medium- to high-throughput preclinical drug screening to identify compounds that may represent novel treatment paradigms or even prevent cancer evolution. This review will highlight and synthesize the findings from zebrafish cancer predisposition models created to date. We will discuss emerging trends in how these zebrafish cancer models can improve our understanding of the genetic mechanisms driving cancer predisposition and their potential to discover therapeutic and/or preventative compounds that change the natural history of disease for these vulnerable children, youth and adults.

Genome-wide identification and function characterization of GATA transcription factors during development and in response to abiotic stresses and hormone treatments in pepper

Pepper (Capsicum annuum L.) is an economically important vegetable crop whose production and quality are severely reduced under adverse environmental stress conditions. The GATA transcription factors belonging to type IV zinc-finger proteins, play a significant role in regulating light morphogenesis, nitrate assimilation, and organ development in plants. However, the functional characteristics of GATA gene family during development and in response to environmental stresses have not yet been investigated in pepper. In this study, a total of 28 pepper GATA (CaGATA) genes were identified. To gain an overview of the CaGATAs, we analyzed their chromosomal distribution, gene structure, conservative domains, cis-elements, phylogeny, and evolutionary relationship. We divided 28 CaGATAs into four groups distributed on 10 chromosomes, and identified 7 paralogs in CaGATA family of pepper and 35 orthologous gene pairs between CaGATAs and Arabidopsis GATAs (AtGATAs). The results of promoter cis-element analysis and the quantitative real-time PCR (qRT-PCR) analysis revealed that CaGATA genes were involved in regulating the plant growth and development and the responses to various abiotic stresses and hormone treatments in pepper. Tissue-specific expression analysis showed that most CaGATA genes were preferentially expressed in flower buds, flowers, and leaves. Several CaGATA genes, especially CaGATA14, were significantly regulated under multiple abiotic stresses, and CaGATA21 and CaGATA27 were highly responsive to phytohormone treatments. Taken together, our results lay a foundation for the biological function analysis of GATA gene family in pepper.

Evolution of Developmental GATA Factors in Nematodes

GATA transcription factors are found in animals, plants, and fungi. In animals, they have important developmental roles in controlling specification of cell identities and executing tissue-specific differentiation. The Phylum Nematoda is a diverse group of vermiform animals that inhabit ecological niches all over the world. Both free-living and parasitic species are known, including those that cause human infectious disease. To date, GATA factors in nematodes have been studied almost exclusively in the model system C. elegans and its close relatives. In this study, we use newly available sequences to identify GATA factors across the nematode phylum. We find that most species have fewer than six GATA factors, but some species have 10 or more. Comparisons of gene and protein structure suggest that there were at most two GATA factors at the base of the phylum, which expanded by duplication and modification to result in a core set of four factors. The high degree of structural similarity with the corresponding orthologues in C. elegans suggests that the nematode GATA factors share similar functions in development.

Deletion of a conserved Gata2 enhancer impairs haemogenic endothelium programming and adult Zebrafish haematopoiesis

Gata2 is a key transcription factor required to generate Haematopoietic Stem and Progenitor Cells (HSPCs) from haemogenic endothelium (HE); misexpression of Gata2 leads to haematopoietic disorders. Here we deleted a conserved enhancer (i4 enhancer) driving pan-endothelial expression of the zebrafish gata2a and showed that Gata2a is required for HE programming by regulating expression of runx1 and of the second Gata2 orthologue, gata2b. By 5 days, homozygous gata2aΔi4/Δi4 larvae showed normal numbers of HSPCs, a recovery mediated by Notch signalling driving gata2b and runx1 expression in HE. However, gata2aΔi4/Δi4 adults showed oedema, susceptibility to infections and marrow hypo-cellularity, consistent with bone marrow failure found in GATA2 deficiency syndromes. Thus, gata2a expression driven by the i4 enhancer is required for correct HE programming in embryos and maintenance of steady-state haematopoietic stem cell output in the adult. These enhancer mutants will be useful in exploring further the pathophysiology of GATA2-related deficiencies in vivo.

Evolutionary Dynamics of the SKN-1 → MED → END-1,3 Regulatory Gene Cascade in Caenorhabditis Endoderm Specification

Gene regulatory networks and their evolution are important in the study of animal development. In the nematode, Caenorhabditis elegans, the endoderm (gut) is generated from a single embryonic precursor, E. Gut is specified by the maternal factor SKN-1, which activates the MED → END-1,3 → ELT-2,7 cascade of GATA transcription factors. In this work, genome sequences from over two dozen species within the Caenorhabditis genus are used to identify MED and END-1,3 orthologs. Predictions are validated by comparison of gene structure, protein conservation, and putative cis-regulatory sites. All three factors occur together, but only within the Elegans supergroup, suggesting they originated at its base. The MED factors are the most diverse and exhibit an unexpectedly extensive gene amplification. In contrast, the highly conserved END-1 orthologs are unique in nearly all species and share extended regions of conservation. The END-1,3 proteins share a region upstream of their zinc finger and an unusual amino-terminal poly-serine domain exhibiting high codon bias. Compared with END-1, the END-3 proteins are otherwise less conserved as a group and are typically found as paralogous duplicates. Hence, all three factors are under different evolutionary constraints. Promoter comparisons identify motifs that suggest the SKN-1, MED, and END factors function in a similar gut specification network across the Elegans supergroup that has been conserved for tens of millions of years. A model is proposed to account for the rapid origin of this essential kernel in the gut specification network, by the upstream intercalation of duplicate genes into a simpler ancestral network.

Cluster expansion of apolipoprotein D (ApoD) genes in teleost fishes

BACKGROUND

Gene and genome duplication play important roles in the evolution of gene function. Compared to individual duplicated genes, gene clusters attract particular attention considering their frequent associations with innovation and adaptation. Here, we report for the first time the expansion of the apolipoprotein D (ApoD) ligand-transporter genes in a cluster manner specific to teleost fishes.

RESULTS

Based on comparative genomic and transcriptomic analyses, protein 3D structure comparison, positive selection detection and breakpoints detection, the single ApoD gene in the ancestor expanded into two clusters following a dynamic evolutionary pattern in teleost fishes. Orthologous genes show conserved expression patterns, whereas lineage-specific duplicated genes show tissue-specific expression patterns and even evolve new gene expression profiles. Positive selection occurred in branches before and after gene duplication, especially for lineage-specific duplicated genes. Cluster analyses based on protein 3D structure comparisons, especially comparisons of the four loops at the opening side, show gene duplication-segregating patterns. Duplicated ApoD genes are predicted to be associated with forkhead transcription factors and MAPK genes. ApoD clusters are located next to the breakpoints of genome rearrangements.

CONCLUSIONS

Here, we report the expansion of ApoD genes specific to teleost fishes in a cluster manner for the first time. Neofunctionalization and subfunctionalization were observed at both the protein and expression levels after duplication. Evidence from different aspects-i.e., abnormal expression-induced disease in humans, fish-specific expansion, predicted associations with forkhead transcription factors and MAPK genes, specific expression patterns in tissues related to sexual selection and adaptation, duplicated genes under positive selection and their location next to the breakpoints of genome rearrangements-suggests the potentially advantageous roles of ApoD genes in teleost fishes. The cluster expansion of ApoD genes specific to teleost fishes provides thus an ideal evo-devo model for studying gene duplication, cluster maintenance and new gene function emergence.

An Anopheles aquasalis GATA factor Serpent is required for immunity against Plasmodium and bacteria

Innate immunity is an ancient and conserved defense system that provides an early effective response against invaders. Many immune genes of Anopheles mosquitoes have been implicated in defense against a variety of pathogens, including plasmodia. Nevertheless, only recent work identified some immune genes of Anopheles aquasalis mosquitoes upon P. vivax infection. Among these was a GATA transcription factor gene, which is described here. This is an ortholog of GATA factor Serpent genes described in Drosophila melanogaster and Anopheles gambiae. Gene expression analyses showed an increase of GATA-Serpent mRNA in P. vivax-infected A. aquasalis and functional RNAi experiments identified this transcription factor as an important immune gene of A. aquasalis against both bacteria and P. vivax. Besides, we were able to identify an effect of GATA-Serpent knockdown on A. aquasalis hemocyte proliferation and differentiation. These findings expand our understanding of the poorly studied A. aquasalis-P. vivax interactions and uncover GATA-Serpent as a key player of the mosquito innate immune response.

The ancient role for GATA2/3 transcription factor homolog in the hemocyte production of oyster

Hemocytes, the cellular component of invertebrate hemolymph, are essential for invertebrate immunity, but the hematopoiesis and regulation mechanism are still largely unknown. In the present study, a conserved hematopoietic transcription factor Cg-GATA2/3 was identified in Pacific oyster Crassotrea gigas, which was evolutionarily close to the vertebrate GATA1/2/3. Cg-GATA2/3 was mainly distributed in the immune organs, such as gill, hemocytes, and mantle. After Cg-GATA2/3 was interferenced by dsRNA, the mRNA expressions of hemocytes specific gene (EcSOD) and hematopoietic transcription factor (C-Myb) were all significant down-regulated, and the hemocyte renewal rates also decreased both in hemolymph and gill. During the larval developmental stages, the mRNA transcripts of Cg-GATA2/3 increased immediately after fertilization and kept a high level during blastula and early trochophore larvae stage (4-10 hpf, hours post fertilization), then decreased sharply in early D-veliger larvae stage (15 hpf). Whole-mount immunofluorescence assay further revealed that the abundant immunoreactivity of Cg-GATA2/3 was distributed in the whole body of blastula and gastrula embryos, while specialized gradually to a ring structure around the dorsal region in trochophore larvae. In the D-veliger and umbo larvae, scattered positive signals appeared in the specific sinus structure on the dorsal side and velum region. These results demonstrated that Cg-GATA2/3 was a hematopoietic lineage-specific transcription factor to regulate the hemocyte production, and it could also be used as hematopoietic specific marker to trace potential developmental events of hematopoiesis during ontogenesis of oyster.

Genome-wide transcriptomics analysis identifies sox7 and sox18 as specifically regulated by gata4 in cardiomyogenesis

The transcription factors GATA4, GATA5 and GATA6 are important regulators of heart muscle differentiation (cardiomyogenesis), which function in a partially redundant manner. We identified genes specifically regulated by individual cardiogenic GATA factors in a genome-wide transcriptomics analysis. The genes regulated by gata4 are particularly interesting because GATA4 is able to induce differentiation of beating cardiomyocytes in Xenopus and in mammalian systems. Among the specifically gata4-regulated transcripts we identified two SoxF family members, sox7 and sox18. Experimental reinstatement of gata4 restores sox7 and sox18 expression, and loss of cardiomyocyte differentiation due to gata4 knockdown is partially restored by reinstating sox7 or sox18 expression, while (as previously reported) knockdown of sox7 or sox18 interferes with heart muscle formation. In order to test for conservation in mammalian cardiomyogenesis, we confirmed in mouse embryonic stem cells (ESCs) undergoing cardiomyogenesis that knockdown of Gata4 leads to reduced Sox7 (and Sox18) expression and that Gata4 is also uniquely capable of promptly inducing Sox7 expression. Taken together, we identify an important and conserved gene regulatory axis from gata4 to the SoxF paralogs sox7 and sox18 and further to heart muscle cell differentiation.

•

Gata 4, 5 and 6 have redundant and non-redundant functions in heart development.

•

RNA-seq analysis of Gata4, 5 and 6 knockdown experiments was carried out.

•

Genes specifically regulated by Gata4, 5 and 6 were identified.

•

The SoxF genes sox7 and sox18 were identified as specifically regulated by Gata4.

•

Epistasis demonstrates a regulatory axis from Gata4 to Sox7/18 to cardiomyogenesis.

An evolutionary perspective on the systems of adaptive immunity

We propose an evolutionary perspective to classify and characterize the diverse systems of adaptive immunity that have been discovered across all major domains of life. We put forward a new function-based classification according to the way information is acquired by the immune systems: Darwinian immunity (currently known from, but not necessarily limited to, vertebrates) relies on the Darwinian process of clonal selection to 'learn' by cumulative trial-and-error feedback; Lamarckian immunity uses templated targeting (guided adaptation) to internalize heritable information on potential threats; finally, shotgun immunity operates through somatic mechanisms of variable targeting without feedback. We argue that the origin of Darwinian (but not Lamarckian or shotgun) immunity represents a radical innovation in the evolution of individuality and complexity, and propose to add it to the list of major evolutionary transitions. While transitions to higher-level units entail the suppression of selection at lower levels, Darwinian immunity re-opens cell-level selection within the multicellular organism, under the control of mechanisms that direct, rather than suppress, cell-level evolution for the benefit of the individual. From a conceptual point of view, the origin of Darwinian immunity can be regarded as the most radical transition in the history of life, in which evolution by natural selection has literally re-invented itself. Furthermore, the combination of clonal selection and somatic receptor diversity enabled a transition from limited to practically unlimited capacity to store information about the antigenic environment. The origin of Darwinian immunity therefore comprises both a transition in individuality and the emergence of a new information system - the two hallmarks of major evolutionary transitions. Finally, we present an evolutionary scenario for the origin of Darwinian immunity in vertebrates. We propose a revival of the concept of the 'Big Bang' of vertebrate immunity, arguing that its origin involved a 'difficult' (i.e. low-probability) evolutionary transition that might have occurred only once, in a common ancestor of all vertebrates. In contrast to the original concept, we argue that the limiting innovation was not the generation of somatic diversity, but the regulatory circuitry needed for the safe operation of amplifiable immune responses with somatically acquired targeting. Regulatory complexity increased abruptly by genomic duplications at the root of the vertebrate lineage, creating a rare opportunity to establish such circuitry. We discuss the selection forces that might have acted at the origin of the transition, and in the subsequent stepwise evolution leading to the modern immune systems of extant vertebrates.

Roles of GATA6 during Gonadal Development in Japanese Flounder: Gonadogenesis, Regulation of Gender-Related Genes, Estrogen Formation and Gonadal Function Maintenance

GATA-binding protein 6 (GATA6), a highly-conserved transcription factor of the GATA family plays an important role in gonadal cell proliferation, differentiation and endoderm development. In this study, the full-length cDNA of GATA6 of Paralichthys olivaceus (Japanese flounder) was obtained. Phylogenetic, gene structure and synteny analyses demonstrated that GATA6 of P. olivaceus is homologous to that of teleosts and tetrapods. The P. olivaceus GATA6 transcript showed higher expression in testis than in ovary, demonstrating a sexually dimorphic gene expression. During embryonic development, the expression of P. olivaceus GATA6 increased at the blastula stage, demonstrating that GATA6 is involved in morphogenesis. Results of in situ hybridization showed that GATA6 signals were detected in Sertoli cells, oogonia and oocytes. Moreover, 17α methyl testosterone, a male hormone, could moderately upregulate P. olivaceus GATA6 and downregulate P. olivaceus aromatase CYP19A1 in testis cells. These results suggest that GATA6 may play an important role in gonadal development in P. olivaceus. This study provides valuable information on the function of P. olivaceus GATA6, laying the foundation for further development of breeding techniques in this species.

Origin and evolution of GATA2a and GATA2b in teleosts: insights from tongue sole, Cynoglossus semilaevis

Background. Following the two rounds of whole-genome duplication that occurred during deuterostome evolution, a third genome duplication occurred in the lineage of teleost fish and is considered to be responsible for much of the biological diversification within the lineage. GATA2, a member of GATA family of transcription factors, is an important regulator of gene expression in hematopoietic cell in mammals, yet the role of this gene or its putative paralogs in ray-finned fishes remains relatively unknown.

Methods. In this study, we attempted to identify GATA2 sequences from the transcriptomes and genomes of multiple teleosts using the bioinformatic tools MrBayes, MEME, and PAML. Following identification, comparative analysis of genome structure, molecular evolution rate, and expression by real-time qPCR were used to predict functional divergence of GATA2 paralogs and their relative transcription in organs of female and male tongue soles (Cynoglossus semilaevis).

Results. Two teleost GATA2 genes were identified in the transcriptomes of tongue sole and Japanese flounder (Paralichthysolivaceus). Synteny and phylogenetic analysis confirmed that the two genes likely originated from the teleost-specific genome duplication . Additionally, selection pressure analysis predicted these gene duplicates to have undergone purifying selection and possible divergent new functions. This was supported by differential expression pattern of GATA2a and GATA2b observed in organs of female and male tongue soles.

Discussion. Our results indicate that two GATA2 genes originating from the first teleost-specific genome duplication have remained transcriptionally active in some fish species and have likely undergone neofunctionalization. This knowledge provides novel insights into the evolution of the teleost GATA2 genes and constituted important groundwork for further research on the GATA gene family.

Gata2b is a restricted early regulator of hemogenic endothelium in the zebrafish embryo

The adult blood system is established by hematopoietic stem cells (HSCs), which arise during development from an endothelial-to-hematopoietic transition of cells comprising the floor of the dorsal aorta. Expression of aortic runx1 has served as an early marker of HSC commitment in the zebrafish embryo, but recent studies have suggested that HSC specification begins during the convergence of posterior lateral plate mesoderm (PLM), well before aorta formation and runx1 transcription. Further understanding of the earliest stages of HSC specification necessitates an earlier marker of hemogenic endothelium. Studies in mice have suggested that GATA2 might function at early stages within hemogenic endothelium. Two orthologs of Gata2 exist in zebrafish: gata2a and gata2b. Here, we report that gata2b expression initiates during the convergence of PLM, becoming restricted to emerging HSCs. We observe Notch-dependent gata2b expression within the hemogenic subcompartment of the dorsal aorta that is in turn required to initiate runx1 expression. Our results indicate that Gata2b functions within hemogenic endothelium from an early stage, whereas Gata2a functions more broadly throughout the vascular system.

Regulation of Cardiac Cell Fate by microRNAs: Implications for Heart Regeneration

microRNAs are post-transcriptional regulators of gene expression that have been shown to be central players in the establishment of cellular programs, often acting as switches that control the choice between proliferation and differentiation during development and in adult tissues. The heart develops from two small patches of cells in the mesoderm, the heart fields, which originate the different cardiac cell types, including cardiomyocytes, vascular smooth muscle and endothelial cells. These progenitors proliferate and differentiate to establish a highly connected three-dimensional structure, involving a robust succession of gene expression programs strongly influenced by microRNAs. Although the mammalian heart has conventionally been viewed as a post-mitotic organ, cardiac cells have recently been shown to display some regenerative potential, which is nonetheless insufficient to regenerate heart lesions, in contrast with other vertebrates like the zebrafish. Both the proliferation of adult cardiac stem cells and the ability of cardiomyocytes to re-enter the cell cycle have been proposed to sustain these regenerative processes. Here we review the role of microRNAs in the control of stem cell and cardiomyocyte dependent cardiac regeneration processes, and discuss potential applications for the treatment of cardiac injury.

AnimalTFDB 2.0: a resource for expression, prediction and functional study of animal transcription factors

Transcription factors (TFs) are key regulators for gene expression. Here we updated the animal TF database AnimalTFDB to version 2.0 (http://bioinfo.life.hust.edu.cn/AnimalTFDB/). Using the improved prediction pipeline, we identified 72 336 TF genes, 21 053 transcription co-factor genes and 6502 chromatin remodeling factor genes from 65 species covering main animal lineages. Besides the abundant annotations (basic information, gene model, protein functional domain, gene ontology, pathway, protein interaction, ortholog and paralog, etc.) in the previous version, we made several new features and functions in the updated version. These new features are: (i) gene expression from RNA-Seq for nine model species, (ii) gene phenotype information, (iii) multiple sequence alignment of TF DNA-binding domains, and the weblogo and phylogenetic tree based on the alignment, (iv) a TF prediction server to identify new TFs from input sequences and (v) a BLAST server to search against TFs in AnimalTFDB. A new nice web interface was designed for AnimalTFDB 2.0 allowing users to browse and search all data in the database. We aim to maintain the AnimalTFDB as a solid resource for TF identification and studies of transcription regulation and comparative genomics.

A synergism between adaptive effects and evolvability drives whole genome duplication to fixation

Whole genome duplication has shaped eukaryotic evolutionary history and has been associated with drastic environmental change and species radiation. While the most common fate of WGD duplicates is a return to single copy, retained duplicates have been found enriched for highly interacting genes. This pattern has been explained by a neutral process of subfunctionalization and more recently, dosage balance selection. However, much about the relationship between environmental change, WGD and adaptation remains unknown. Here, we study the duplicate retention pattern postWGD, by letting virtual cells adapt to environmental changes. The virtual cells have structured genomes that encode a regulatory network and simple metabolism. Populations are under selection for homeostasis and evolve by point mutations, small indels and WGD. After populations had initially adapted fully to fluctuating resource conditions re-adaptation to a broad range of novel environments was studied by tracking mutations in the line of descent. WGD was established in a minority (≈30%) of lineages, yet, these were significantly more successful at re-adaptation. Unexpectedly, WGD lineages conserved more seemingly redundant genes, yet had higher per gene mutation rates. While WGD duplicates of all functional classes were significantly over-retained compared to a model of neutral losses, duplicate retention was clearly biased towards highly connected TFs. Importantly, no subfunctionalization occurred in conserved pairs, strongly suggesting that dosage balance shaped retention. Meanwhile, singles diverged significantly. WGD, therefore, is a powerful mechanism to cope with environmental change, allowing conservation of a core machinery, while adapting the peripheral network to accommodate change.

The evolution of eukaryotes is characterized by drastic changes in their genome content. Genome expansions have often occurred by duplication of the entire genome. It is generally not know whether organisms gain any adaptive advantage from these mutations. However, they appear to become fixed in response to environmental change. Many interesting whole genome duplications happened long ago in eukaryotic evolutionary history during periods of turbulent genome and species evolution. Genomic data analysis alone cannot resolve the evolutionary mechanisms and consequences of whole genome duplication. Here, we modeled evolution with whole genome duplications in a Virtual Cell model. Simulating populations that undergo a range of different environmental changes we found that next to often increasing fitness directly, whole genome duplications made lineages more evolvable and hence more able to adapt to harsh new environments. Although most duplicates are deleted in subsequent evolution, genes with many interaction partners were retained preferentially, increasing regulatory complexity. Interestingly however, we found that innovation happened most likely in the more loosely connected and less essential genes.

GATA transcription factors in vertebrates: evolutionary, structural and functional interplay

GATA transcription factors perform conserved and essential roles during animal development, including germ-layer specification, hematopoiesis, and cardiogenesis. The evolutionary history and the changes in selection pressures following duplication of the six GATA family members in vertebrates have not been completely understood. Recently, we explored multiple databases to find GATAs in different vertebrate species. Using these sequences, we have performed molecular phylogenetic analyses using Maximum Likelihood and Bayesian methods, and statistical tests of tree topologies, to ascertain the phylogenetic relationship and selection pressures among GATA proteins. Seventy-one full-length cDNA sequences from 24 vertebrate species were extracted from multiple databases. By phylogenetic analyses, we investigated the origin, conservation, and evolution of the GATAs. Six GATA genes in vertebrates might be formed by gene duplication. The inferred evolutionary transitions that separate members which belong to different gene clusters correlated with changes in functional properties. Selection analysis and protein structure analysis were combined to explain Darwinian selection in GATA sequences and these changes brought putative biological significance. 26 positive selection sites were detected in this process. This study reveals the evolutionary history of vertebrate GATA paralogous and positively selected sites likely relevant for the distinct functional properties of the paralogs. It provides a new perspective for understanding the origin and evolution and biological functions of GATAs, which will help to uncover the GATAs' biological roles, evolution and their relationship with associated diseases; in addition, other complex multidomain families and also larger superfamilies can be investigated in a similar way.

Age distribution patterns of human gene families: divergent for Gene Ontology categories and concordant between different subcellular localizations

The age distribution of gene duplication events within the human genome exhibits two waves of duplications along with an ancient component. However, because of functional constraint differences, genes in different functional categories might show dissimilar retention patterns after duplication. It is known that genes in some functional categories are highly duplicated in the early stage of vertebrate evolution. However, the correlations of the age distribution pattern of gene duplication between the different functional categories are still unknown. To investigate this issue, we developed a robust pipeline to date the gene duplication events in the human genome. We successfully estimated about three-quarters of the duplication events within the human genome, along with the age distribution pattern in each Gene Ontology (GO) slim category. We found that some GO slim categories show different distribution patterns when compared to the whole genome. Further hierarchical clustering of the GO slim functional categories enabled grouping into two main clusters. We found that human genes located in the duplicated copy number variant regions, whose duplicate genes have not been fixed in the human population, were mainly enriched in the groups with a high proportion of recently duplicated genes. Moreover, we used a phylogenetic tree-based method to date the age of duplications in three signaling-related gene superfamilies: transcription factors, protein kinases and G-protein coupled receptors. These superfamilies were expressed in different subcellular localizations. They showed a similar age distribution as the signaling-related GO slim categories. We also compared the differences between the age distributions of gene duplications in multiple subcellular localizations. We found that the distribution patterns of the major subcellular localizations were similar to that of the whole genome. This study revealed the whole picture of the evolution patterns of gene functional categories in the human genome.

How life changes itself: the Read-Write (RW) genome

The genome has traditionally been treated as a Read-Only Memory (ROM) subject to change by copying errors and accidents. In this review, I propose that we need to change that perspective and understand the genome as an intricately formatted Read-Write (RW) data storage system constantly subject to cellular modifications and inscriptions. Cells operate under changing conditions and are continually modifying themselves by genome inscriptions. These inscriptions occur over three distinct time-scales (cell reproduction, multicellular development and evolutionary change) and involve a variety of different processes at each time scale (forming nucleoprotein complexes, epigenetic formatting and changes in DNA sequence structure). Research dating back to the 1930s has shown that genetic change is the result of cell-mediated processes, not simply accidents or damage to the DNA. This cell-active view of genome change applies to all scales of DNA sequence variation, from point mutations to large-scale genome rearrangements and whole genome duplications (WGDs). This conceptual change to active cell inscriptions controlling RW genome functions has profound implications for all areas of the life sciences.

LGI proteins in the nervous system

The development and function of the vertebrate nervous system depend on specific interactions between different cell types. Two examples of such interactions are synaptic transmission and myelination. LGI1-4 (leucine-rich glioma inactivated proteins) play important roles in these processes. They are secreted proteins consisting of an LRR (leucine-rich repeat) domain and a so-called epilepsy-associated or EPTP (epitempin) domain. Both domains are thought to function in protein–protein interactions. The first LGI gene to be identified, LGI1, was found at a chromosomal translocation breakpoint in a glioma cell line. It was subsequently found mutated in ADLTE (autosomal dominant lateral temporal (lobe) epilepsy) also referred to as ADPEAF (autosomal dominant partial epilepsy with auditory features). LGI1 protein appears to act at synapses and antibodies against LGI1 may cause the autoimmune disorder limbic encephalitis. A similar function in synaptic remodelling has been suggested for LGI2, which is mutated in canine Benign Familial Juvenile Epilepsy. LGI4 is required for proliferation of glia in the peripheral nervous system and binds to a neuronal receptor, ADAM22, to foster ensheathment and myelination of axons by Schwann cells. Thus, LGI proteins play crucial roles in nervous system development and function and their study is highly important, both to understand their biological functions and for their therapeutic potential. Here, we review our current knowledge about this important family of proteins, and the progress made towards understanding their functions.

The evolutionary origins of chordate hematopoiesis and vertebrate endothelia

The vertebrate circulatory system is the most complex vascular system among those of metazoans, with key innovations including a multi-chambered heart and highly specialized blood cells. Invertebrate vessels, on the other hand, consist of hemal spaces between the basal laminae of epithelia. How the evolutionary transition from an invertebrate-type system to the complex vertebrate one occurred is, however, poorly understood. We investigate here the development of the cardiovascular system of the cephalochordate amphioxus Branchiostoma lanceolatum in order to gain insight into the origin of the vertebrate cardiovascular system. The cardiac markers Hand, Csx (Nkx2-5) and Tbx4/5 reveal a broad cardiac-like domain in amphioxus; such a decentralized organization during development parallels that seen in the adult anatomy. Our data therefore support the hypothesis that amphioxus never possessed a proper heart, even transiently during development. We also define a putative hematopoietic domain, supported by the expression of the hematopoietic markers Scl and Pdvegfr. We show that this area is closed to the dorsal aorta anlages, partially linked to excretory tissues, and that its development is regulated by retinoic acid, thus recalling the aorta-gonads-mesonephros (AGM) area of vertebrates. This region probably produces Pdvegfr+ hemal cells, with an important role in amphioxus vessel formation, since treatments with an inhibitor of PDGFR/VEGFR lead to a decrease of Laminin in the basal laminae of developing vessels. Our results point to a chordate origin of hematopoiesis in an AGM-like area from where hemal Pdvegfr+ cells are produced. These Pdvegfr+ cells probably resemble the ancestral chordate blood cells from which the vertebrate endothelium later originated.

Evolution of the tbx6/16 subfamily genes in vertebrates: insights from zebrafish

In any comparative studies striving to understand the similarities and differences of the living organisms at the molecular genetic level, the crucial first step is to establish the homology (orthology and paralogy) of genes between different organisms. Determination of the homology of genes becomes complicated when the genes have undergone a rapid divergence in sequence or when the involved genes are members of a gene family that has experienced a differential gain or loss of its constituents in different taxonomic groups. Organisms with duplicated genomes such as teleost fishes might have been especially prone to these problems because the functional redundancies provided by the duplicate copies of genes would have allowed a rapid divergence or loss of genes during evolution. In this study, we will demonstrate that much of the ambiguities in the determination of the homology between fish and tetrapod genes resulting from the problems like these can be eliminated by complementing the sequence-based phylogenies with nonsequence information, such as the exon-intron structure of a gene or the composition of a gene's genomic neighbors. We will use the Tbx6/16 subfamily genes of zebrafish (tbx6, tbx16, tbx24, and mga genes), which have been well known for the ambiguity of their evolutionary relationships to the Tbx6/16 subfamily genes of tetrapods, as an illustrative example. We will show that, despite the similarity of sequence and expression to the tetrapod Tbx6 genes, zebrafish tbx6 gene is actually a novel T-box gene more closely related to the tetrapod Tbx16 genes, whereas the zebrafish tbx24 gene, hitherto considered to be a novel gene due to the high level of sequence divergence, is actually an ortholog of tetrapod Tbx6 genes. We will also show that, after their initial appearance by the multiplication of a common ancestral gene at the beginning of vertebrate evolution, the Tbx6/16 subfamily of vertebrate T-box genes might have experienced differential losses of member genes in different vertebrate groups and gradual pooling of member gene's functions in surviving members, which might have prevented the revelation of the true identity of member genes by way of the comparison of sequence and function.

An evolutionarily conserved kernel of gata5, gata6, otx2 and prdm1a operates in the formation of endoderm in zebrafish

An evolutionarily conserved subcircuit (kernel) dedicated to a specific developmental function is found at the top of the gene regulatory networks (GRNs) hierarchy. Here we comprehensively demonstrate that a pan-deuterostome endoderm specification kernel exists in zebrafish. We analyzed interactions among gata5, gata6, otx2 and prdm1a using specific morpholino knockdowns and measured the gene expression profiles by quantitative real-time RT-PCR and in situ hybridization. The mRNA rescue experiment validated the specificity of the morpholino knockdown. We found that the interactions among gata5, gata6, otx2 and prdm1a determine the initial specification of the zebrafish endoderm. Although otx2 can activate both gata5 and gata6, and the prdm1a/krox homologue also activates some endoderm transcription factors, a feedback loop from Gata to otx2 and prdm1a is missing. Furthermore, we found the positive regulation between gata5 and gata6 to further lock-on the mesendoderm specification by the Gata family. Chromatin immunoprecipitation was used to further validate the recruitment of Otx2 to the gata5 and gata6 loci. Functional assays revealed that module B of gata6 and the basal promoter of gata5 drive the gene at the mesendoderm, and mutational analysis demonstrated that Otx2 and Gata5/6 contribute to reporter gene activation. This is the first direct evidence for evolutionarily conserved endoderm specification across echinoderms and vertebrates.

Evolution of the regulatory control of vertebrate striated muscle: the roles of troponin I and myosin binding protein-C

Troponin I (TnI) and myosin binding protein-C (MyBP-C) are key regulatory proteins of contractile function in vertebrate muscle. TnI modulates the Ca2+ activation signal, while MyBP-C regulates cross-bridge cycling kinetics. In vertebrates, each protein is distributed as tissue-specific paralogs in fast skeletal (fs), slow skeletal (ss), and cardiac (c) muscles. The purpose of this study is to characterize how TnI and MyBP-C have changed during the evolution of vertebrate striated muscle and how tissue-specific paralogs have adapted to different physiological conditions. To accomplish this we have completed phylogenetic analyses using the amino acid sequences of all known TnI and MyBP-C isoforms. This includes 99 TnI sequences (fs, ss, and c) from 51 different species and 62 MyBP-C sequences from 26 species, with representatives from each vertebrate group. Results indicate that the role of protein kinase A (PKA) and protein kinase C (PKC) in regulating contractile function has changed during the evolution of vertebrate striated muscle. This is reflected in an increased number of phosphorylatable sites in cTnI and cMyBP-C in endothermic vertebrates and the loss of two PKC sites in fsTnI in a common ancestor of mammals, birds, and reptiles. In addition, we find that His132, Val134, and Asn141 in human ssTnI, previously identified as enabling contractile function during cellular acidosis, are present in all vertebrate cTnI isoforms except those from monotremes, marsupials, and eutherian mammals. This suggests that the replacement of these residues with alternative residues coincides with the evolution of endothermy in the mammalian lineage.

Development of the C. elegans digestive tract

The C. elegans digestive tract (pharynx, intestine, and rectum) contains only approximately 100 cells but develops under the control of the same types of transcription factors (e.g. FoxA and GATA factors) that control digestive tract development in far more complex animals. The GATA-factor dominated core regulatory hierarchy directing development of the homogenous clonal intestine from oocyte to mature organ is now known with some degree of certainty, setting the stage for more biochemical experiments to understand developmental mechanisms. The FoxA-factor dominated development of the pharynx (and rectum) is less well understood but is beginning to reveal how transcription factor combinations produce unique cell types within organs.

Genome duplication and T cell immunity

The adaptive immune system (AIS) mediated by T cells and B cells arose ~450 million years ago in a common ancestor of jawed vertebrates. This system was so successful that, once established, it has been maintained in all classes of jawed vertebrates with only minor modifications. One event thought to have contributed to the emergence of this form of AIS is two rounds of whole-genome duplication. This event enabled jawed vertebrate ancestors to acquire many paralogous genes, known as ohnologs, with essential roles in T cell and B cell immunity. Ohnologs encode the key components of the antigen presentation machinery and signal transduction pathway for lymphocyte activation as well as numerous transcription factors important for lymphocyte development. Recently, it has been discovered that jawless vertebrates have developed an AIS employing antigen receptors unrelated to T/B cell receptors, but with marked overall similarities to the AIS of jawed vertebrates. Emerging evidence suggests that a common ancestor of all vertebrates was equipped with T-lymphoid and B-lymphoid lineages.