Sequencing human-gibbon breakpoints of synteny reveals mosaic new insertions at rearrangement sites

The gibbon genome exhibits extensive karyotypic diversity with an increased rate of chromosomal rearrangements during evolution. In an effort to understand the mechanistic origin and implications of these rearrangement events, we sequenced 24 synteny breakpoint regions in the white-cheeked gibbon (Nomascus leucogenys, NLE) in the form of high-quality BAC insert sequences (4.2 Mbp). While there is a significant deficit of breakpoints in genes, we identified seven human gene structures involved in signaling pathways (DEPDC4, GNG10), phospholipid metabolism (ENPP5, PLSCR2), beta-oxidation (ECH1), cellular structure and transport (HEATR4), and transcription (ZNF461), that have been disrupted in the NLE gibbon lineage. Notably, only three of these genes show the expected evolutionary signatures of pseudogenization. Sequence analysis of the breakpoints suggested both nonclassical nonhomologous end-joining (NHEJ) and replication-based mechanisms of rearrangement. A substantial number (11/24) of human-NLE gibbon breakpoints showed new insertions of gibbon-specific repeats and mosaic structures formed from disparate sequences including segmental duplications, LINE, SINE, and LTR elements. Analysis of these sites provides a model for a replication-dependent repair mechanism for double-strand breaks (DSBs) at rearrangement sites and insights into the structure and formation of primate segmental duplications at sites of genomic rearrangements during evolution.

Sensitive gene expression profiling of human T cell subsets reveals parallel post-thymic differentiation for CD4+ and CD8+ lineages.

The differentiation of CD4(+) or CD8(+) T cells following priming of naive cells is central in the establishment of the immune response against pathogens or tumors. However, our understanding of this complex process and the significance of the multiple subsets of differentiation remains controversial. Gene expression profiling has opened new directions of investigation in immunobiology. Nonetheless, the need for substantial amount of biological material often limits its application range. In this study, we have developed procedures to perform microarray analysis on amplified cDNA from low numbers of cells, including primary T lymphocytes, and applied this technology to the study of CD4 and CD8 lineage differentiation. Gene expression profiling was performed on samples of 1000 cells from 10 different subpopulations, defining the major stages of post-thymic CD4(+) or CD8(+) T cell differentiation. Surprisingly, our data revealed that while CD4(+) and CD8(+) T cell gene expression programs diverge at early stages of differentiation, they become increasingly similar as cells reach a late differentiation stage. This suggests that functional heterogeneity between Ag experienced CD4(+) and CD8(+) T cells is more likely to be located early during post-thymic differentiation, and that late stages of differentiation may represent a common end in the development of T-lymphocytes.

Numerous potentially functional but non-genic conserved sequences on human chromosome 21.

The use of comparative genomics to infer genome function relies on the understanding of how different components of the genome change over evolutionary time. The aim of such comparative analysis is to identify conserved, functionally transcribed sequences such as protein-coding genes and non-coding RNA genes, and other functional sequences such as regulatory regions, as well as other genomic features. Here, we have compared the entire human chromosome 21 with syntenic regions of the mouse genome, and have identified a large number of conserved blocks of unknown function. Although previous studies have made similar observations, it is unknown whether these conserved sequences are genes or not. Here we present an extensive experimental and computational analysis of human chromosome 21 in an effort to assign function to sequences conserved between human chromosome 21 (ref. 8) and the syntenic mouse regions. Our data support the presence of a large number of potentially functional non-genic sequences, probably regulatory and structural. The integration of the properties of the conserved components of human chromosome 21 to the rapidly accumulating functional data for this chromosome will improve considerably our understanding of the role of sequence conservation in mammalian genomes.

Detection of three Allexivirus species infecting garlic in Brazil

Garlic viruses often occur in mixed infections under field conditions. In this study, garlic samples collected in three geographical areas of Brazil were tested by Dot-ELISA for the detection of allexiviruses using monoclonal specific antibodies to detect Garlic virus A (GarV-A), Garlic virus B (GarV-B), Garlic virus C (GarV-C) and a polyclonal antiserum able to detect the three virus species mentioned plus Garlic virus D (GarV-D). The detected viruses were biologically isolated by successive passages through Chenopodium quinoa. Reverse Transcriptase Polimerase Chain Reaction (RT-PCR) was performed using primers designed from specific regions of the coat protein genes of Japanese allexiviruses available in the Genetic Bank of National Center of Biotechnology Information (NCBI). By these procedures, individual garlic virus genomes were isolated and sequenced. The nucleotide and amino acid sequence analysis and the one with serological data revealed the presence of three distinct allexiviruses GarV-C, GarV-D and a recently described allexivirus, named Garlic mite-borne filamentous virus (GarMbFV), in Brazil.

Adaptation of aerobically growing Pseudomonas aeruginosa to copper starvation.

Restricted bioavailability of copper in certain environments can interfere with cellular respiration because copper is an essential cofactor of most terminal oxidases. The global response of the metabolically versatile bacterium and opportunistic pathogen Pseudomonas aeruginosa to copper limitation was assessed under aerobic conditions. Expression of cioAB (encoding an alternative, copper-independent, cyanide-resistant ubiquinol oxidase) was upregulated, whereas numerous iron uptake functions (including the siderophores pyoverdine and pyochelin) were expressed at reduced levels, presumably reflecting a lower demand for iron by respiratory enzymes. Wild-type P. aeruginosa was able to grow aerobically in a defined glucose medium depleted of copper, whereas a cioAB mutant did not grow. Thus, P. aeruginosa relies on the CioAB enzyme to cope with severe copper deprivation. A quadruple cyo cco1 cco2 cox mutant, which was deleted for all known heme-copper terminal oxidases of P. aeruginosa, grew aerobically, albeit more slowly than did the wild type, indicating that the CioAB enzyme is capable of energy conservation. However, the expression of a cioA'-'lacZ fusion was less dependent on the copper status in the quadruple mutant than in the wild type, suggesting that copper availability might affect cioAB expression indirectly, via the function of the heme-copper oxidases.

Auxin response factors ARF6 and ARF8 promote jasmonic acid production and flower maturation.

Pollination in flowering plants requires that anthers release pollen when the gynoecium is competent to support fertilization. We show that in Arabidopsis thaliana, two paralogous auxin response transcription factors, ARF6 and ARF8, regulate both stamen and gynoecium maturation. arf6 arf8 double-null mutant flowers arrested as infertile closed buds with short petals, short stamen filaments, undehisced anthers that did not release pollen and immature gynoecia. Numerous developmentally regulated genes failed to be induced. ARF6 and ARF8 thus coordinate the transition from immature to mature fertile flowers. Jasmonic acid (JA) measurements and JA feeding experiments showed that decreased jasmonate production caused the block in pollen release, but not the gynoecium arrest. The double mutant had altered auxin responsive gene expression. However, whole flower auxin levels did not change during flower maturation, suggesting that auxin might regulate flower maturation only under specific environmental conditions, or in localized organs or tissues of flowers. arf6 and arf8 single mutants and sesquimutants (homozygous for one mutation and heterozygous for the other) had delayed stamen development and decreased fecundity, indicating that ARF6 and ARF8 gene dosage affects timing of flower maturation quantitatively.

The evolution of the thyroid hormone distributor protein transthyretin in the order insectivora, class mammalia.

Thyroid hormones are involved in the regulation of growth and metabolism in all vertebrates. Transthyretin is one of the extracellular proteins with high affinity for thyroid hormones which determine the partitioning of these hormones between extracellular compartments and intracellular lipids. During vertebrate evolution, both the tissue pattern of expression and the structure of the gene for transthyretin underwent characteristic changes. The purpose of this study was to characterize the position of Insectivora in the evolution of transthyretin in eutherians, a subclass of Mammalia. Transthyretin was identified by thyroxine binding and Western analysis in the blood of adult shrews, hedgehogs, and moles. Transthyretin is synthesized in the liver and secreted into the bloodstream, similar to the situation for other adult eutherians, birds, and diprotodont marsupials, but different from that for adult fish, amphibians, reptiles, monotremes, and Australian polyprotodont marsupials. For the characterization of the structure of the gene and the processing of mRNA for transthyretin, cDNA libraries were prepared from RNA from hedgehog and shrew livers, and full-length cDNA clones were isolated and sequenced. Sections of genomic DNA in the regions coding for the splice sites between exons 1 and 2 were synthesized by polymerase chain reaction and sequenced. The location of splicing was deduced from comparison of genomic with cDNA nucleotide sequences. Changes in the nucleotide sequence of the transthyretin gene during evolution are most pronounced in the region coding for the N-terminal region of the protein. Both the derived overall amino sequences and the N-terminal regions of the transthyretins in Insectivora were found to be very similar to those in other eutherians but differed from those found in marsupials, birds, reptiles, amphibians, and fish. Also, the pattern of transthyretin precursor mRNA splicing in Insectivora was more similar to that in other eutherians than to that in marsupials, reptiles, and birds. Thus, in contrast to the marsupials, with a different pattern of transthyretin gene expression in the evolutionarily "older" polyprotodonts compared with the evolutionarily "younger" diprotodonts, no separate lineages of transthyretin evolution could be identified in eutherians. We conclude that transthyretin gene expression in the liver of adult eutherians probably appeared before the branching of the lineages leading to modern eutherian species.

Leishmania major: histone H1 gene expression from the sw3 locus.

Histone H1 in the parasitic protozoan Leishmania is a developmentally regulated protein encoded by the sw3 gene. Here we report that histone H1 variants exist in different Leishmania species and strains of L. major and that they are encoded by polymorphic genes. Amplification of the sw3 gene from the genome of three strains of L. major gave rise to different products in each strain, suggesting the presence of a multicopy gene family. In L. major, these genes were all restricted to a 50-kb Bg/II fragment found on a chromosomal band of 1.3 Mb (chromosome 27). The detection of RFLPs in this locus demonstrated its heterogeneity within several species and strains of Leishmania. Two different copies of sw3 (sw3.0 and sw3.1) were identified after screening a cosmid library containing L. major strain Friedlin genomic DNA. They were identical in their 5' UTRs and open reading frames, but differed in their 3' UTRs. With respect to the originally cloned copy of sw3 from L. major strain LV39, their open reading frames lacked a repeat unit of 9 amino acids. Immunoblots of L. guyanensis parasites transfected with these cosmids revealed that both copies could give rise to the histone H1 protein. The characterization of this locus will now make possible a detailed analysis of the function of histone H1 in Leishmania, as well as permit the dissection of the molecular mechanisms governing the developmental regulation of the sw3 gene.

Nephrocalcinosis (enamel renal syndrome) caused by autosomal recessive FAM20A mutations.

BACKGROUND/AIMS: Calcium homeostasis requires regulated cellular and interstitial systems interacting to modulate the activity and movement of this ion. Disruption of these systems in the kidney results in nephrocalcinosis and nephrolithiasis, important medical problems whose pathogenesis is incompletely understood. METHODS: We investigated 25 patients from 16 families with unexplained nephrocalcinosis and characteristic dental defects (amelogenesis imperfecta, gingival hyperplasia, impaired tooth eruption). To identify the causative gene, we performed genome-wide linkage analysis, exome capture, next-generation sequencing, and Sanger sequencing. RESULTS: All patients had bi-allelic FAM20A mutations segregating with the disease; 20 different mutations were identified. CONCLUSIONS: This autosomal recessive disorder, also known as enamel renal syndrome, of FAM20A causes nephrocalcinosis and amelogenesis imperfecta. We speculate that all individuals with biallelic FAM20A mutations will eventually show nephrocalcinosis.

Indexing strategies for rapid searches of short words in genome sequences.

Searching for matches between large collections of short (14-30 nucleotides) words and sequence databases comprising full genomes or transcriptomes is a common task in biological sequence analysis. We investigated the performance of simple indexing strategies for handling such tasks and developed two programs, fetchGWI and tagger, that index either the database or the query set. Either strategy outperforms megablast for searches with more than 10,000 probes. FetchGWI is shown to be a versatile tool for rapidly searching multiple genomes, whose performance is limited in most cases by the speed of access to the filesystem. We have made publicly available a Web interface for searching the human, mouse, and several other genomes and transcriptomes with oligonucleotide queries.

Identification of HMX1 target genes: a predictive promoter model approach.

PURPOSE: A homozygous mutation in the H6 family homeobox 1 (HMX1) gene is responsible for a new oculoauricular defect leading to eye and auricular developmental abnormalities as well as early retinal degeneration (MIM 612109). However, the HMX1 pathway remains poorly understood, and in the first approach to better understand the pathway's function, we sought to identify the target genes. METHODS: We developed a predictive promoter model (PPM) approach using a comparative transcriptomic analysis in the retina at P15 of a mouse model lacking functional Hmx1 (dmbo mouse) and its respective wild-type. This PPM was based on the hypothesis that HMX1 binding site (HMX1-BS) clusters should be more represented in promoters of HMX1 target genes. The most differentially expressed genes in the microarray experiment that contained HMX1-BS clusters were used to generate the PPM, which was then statistically validated. Finally, we developed two genome-wide target prediction methods: one that focused on conserving PPM features in human and mouse and one that was based on the co-occurrence of HMX1-BS pairs fitting the PPM, in human or in mouse, independently. RESULTS: The PPM construction revealed that sarcoglycan, gamma (35kDa dystrophin-associated glycoprotein) (Sgcg), teashirt zinc finger homeobox 2 (Tshz2), and solute carrier family 6 (neurotransmitter transporter, glycine) (Slc6a9) genes represented Hmx1 targets in the mouse retina at P15. Moreover, the genome-wide target prediction revealed that mouse genes belonging to the retinal axon guidance pathway were targeted by Hmx1. Expression of these three genes was experimentally validated using a quantitative reverse transcription PCR approach. The inhibitory activity of Hmx1 on Sgcg, as well as protein tyrosine phosphatase, receptor type, O (Ptpro) and Sema3f, two targets identified by the PPM, were validated with luciferase assay. CONCLUSIONS: Gene expression analysis between wild-type and dmbo mice allowed us to develop a PPM that identified the first target genes of Hmx1.

Ecology and life history affect different aspects of the population structure of 27 high-alpine plants.

A plant species' genetic population structure is the result of a complex combination of its life history, ecological preferences, position in the ecosystem and historical factors. As a result, many different statistical methods exist that measure different aspects of species' genetic structure. However, little is known about how these methods are interrelated and how they are related to a species' ecology and life history. In this study, we used the IntraBioDiv amplified fragment length polymorphisms data set from 27 high-alpine species to calculate eight genetic summary statistics that we jointly correlate to a set of six ecological and life-history traits. We found that there is a large amount of redundancy among the calculated summary statistics and that there is a significant association with the matrix of species traits. In a multivariate analysis, two main aspects of population structure were visible among the 27 species. The first aspect is related to the species' dispersal capacities and the second is most likely related to the species' postglacial recolonization of the Alps. Furthermore, we found that some summary statistics, most importantly Mantel's r and Jost's D, show different behaviour than expected based on theory. We therefore advise caution in drawing too strong conclusions from these statistics.

Global invasion history of the fire ant Solenopsis invicta.

The fire ant Solenopsis invicta is a significant pest that was inadvertently introduced into the southern United States almost a century ago and more recently into California and other regions of the world. An assessment of genetic variation at a diverse set of molecular markers in 2144 fire ant colonies from 75 geographic sites worldwide revealed that at least nine separate introductions of S. invicta have occurred into newly invaded areas and that the main southern U.S. population is probably the source of all but one of these introductions. The sole exception involves a putative serial invasion from the southern United States to California to Taiwan. These results illustrate in stark fashion a severe negative consequence of an increasingly massive and interconnected global trade and travel system.

Multicentric carpotarsal osteolysis is caused by mutations clustering in the amino-terminal transcriptional activation domain of MAFB.

Multicentric carpotarsal osteolysis (MCTO) is a rare skeletal dysplasia characterized by aggressive osteolysis, particularly affecting the carpal and tarsal bones, and is frequently associated with progressive renal failure. Using exome capture and next-generation sequencing in five unrelated simplex cases of MCTO, we identified previously unreported missense mutations clustering within a 51 base pair region of the single exon of MAFB, validated by Sanger sequencing. A further six unrelated simplex cases with MCTO were also heterozygous for previously unreported mutations within this same region, as were affected members of two families with autosomal-dominant MCTO. MAFB encodes a transcription factor that negatively regulates RANKL-induced osteoclastogenesis and is essential for normal renal development. Identification of this gene paves the way for development of novel therapeutic approaches for this crippling disease and provides insight into normal bone and kidney development.

Signaling pathways controlling induced resistance to insect herbivores in Arabidopsis.

Insect attack triggers changes in transcript level in plants that are mediated predominantly by jasmonic acid (JA). The implication of ethylene (ET), salicylic acid (SA), and other signals in this response is less understood and was monitored with a microarray containing insect- and defense-regulated genes. Arabidopsis thaliana mutants coi1-1, ein2-1, and sid2-1 impaired in JA, ET, and SA signaling pathways were challenged with the specialist small cabbage white (Pieris rapae) and the generalist Egyptian cotton worm (Spodoptera littoralis). JA was shown to be a major signal controlling the upregulation of defense genes in response to either insect but was found to suppress changes in transcript level only in response to P. rapae. Larval growth was affected by the JA-dependent defenses, but S. littoralis gained much more weight on coi1-1 than P. rapae. ET and SA mutants had an altered transcript profile after S. littoralis herbivory but not after P. rapae herbivory. In contrast, both insects yielded similar transcript signatures in the abscisic acid (ABA)-biosynthetic mutants aba2-1 and aba3-1, and ABA controlled transcript levels both negatively and positively in insect-attacked plants. In accordance with the transcript signature, S. littoralis larvae performed better on aba2-1 mutants. This study reveals a new role for ABA in defense against insects in Arabidopsis and identifies some components important for plant resistance to herbivory.