Genomic evidence for the Pleistocene and recent population history of Native Americans

How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative “Paleoamerican” relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.

Integration of steroid hormone initiated membrane action to genomic function in the brain

Estrogen is a ligand for the estrogen receptor (ER), which on binding 17beta-estradiol, functions as a ligand-activated transcription factor and regulates the transcription of target genes. This is the slow genomic mode of action. However, rapid non-genomic actions of estrogen also exist at the cell membrane. Using a novel two-pulse paradigm in which the first pulse rapidly initiates non-genomic actions using a membrane-limited estrogen conjugate (E-BSA), while the second pulse promotes genomic transcription from a consensus estrogen response element (ERE), we have demonstrated that rapid actions of estrogen potentiate the slower transcriptional response from an ERE-reporter in neuroblastoma cells. Since rapid actions of estrogen activate kinases, we used selective inhibitors in the two-pulse paradigm to determine the intracellular signaling cascades important in such potentiation. Inhibition of protein kinase A (PKA), PKC, mitogen activated protein kinase (MAPK) or phosphatidylinositol 3-OH kinase (PI-3K) in the first pulse decreases potentiation of transcription. Also, our data with both dominant negative and constitutive mutants of Galpha subunits show that Galpha(q) initiates the rapid signaling cascade at the membrane in SK-N-BE(2)C neuroblastoma cells. We discuss two models of multiple kinase activation at the membrane Pulses of estrogen induce lordosis behavior in female rats. Infusion of E-BSA into the ventromedial hypothalamus followed by 17beta-estradiol in the second pulse could induce lordosis behavior, demonstrating the applicability of this paradigm in vivo. A model where non-genomic actions of estrogen couple to genomic actions unites both aspects of hormone action.

Proteomic analysis of the medicinal plant Artemisia annua: data from leaf and trichome extracts

This article contains raw and processed data related to research published by Bryant et al. [1]. Data was obtained by MS-based proteomics, analysing trichome-enriched, trichome-depleted and whole leaf samples taken from the medicinal plant Artemisia annua and searching the acquired MS/MS data against a recently published contig database [2] and other genomic and proteomic sequence databases for comparison. The processed data shows that an order-of-magnitude more proteins have been identified from trichome-enriched Artemisia annua samples in comparison to previously published data. Proteins known to have a role in the biosynthesis of artemisinin and other highly abundant proteins were found which imply additional enzymatically driven processes occurring within the trichomes that are significant for the biosynthesis of artemisinin.

Genomic and nongenomic stimulatory effect of aldosterone on H(+)-ATPase in proximal S3 segments

Leite-Dellova DC, Malnic G, Mello-Aires M. Genomic and non-genomic stimulatory effect of aldosterone on H(+)-ATPase in proximal S3 segments. Am J Physiol Renal Physiol 300: F682-F691, 2011. First published December 29, 2010; doi:10.1152/ajprenal.00172.2010.-The genomic and nongenomic effects of aldosterone on the intracellular pH recovery rate (pHirr) via H(-)(+)ATPase and on cytosolic free calcium concentration ([Ca(2+)](i)) were investigated in isolated proximal S3 segments of rats during superfusion with an Na(+)-free solution, by using the fluorescent probes BCECF-AM and FLUO-4-AM, respectively. The pHirr, after cellular acidification with a NH(4)Cl pulse, was 0.064 +/- 0.003 pH units/min (n = 17/74) and was abolished with concanamycin. Aldosterone (10(-12), 10(-10),10(-8), or 10(-6) M with 1-h or 15- or 2-min preincubation) increased the pHirr. The baseline [Ca(2+)](i) was 103 +/- 2 nM (n = 58). After 1 min of aldosterone preincubation, there was a transient and dose-dependent increase in [Ca(2+)](i) and after 6-min preincubation there was a new increase in [Ca(2+)](i) that persisted after 1 h. Spironolactone [mineralocorticoid (MR) antagonist], actinomycin D, or cycloheximide did not affect the effects of aldosterone (15- or 2-min preincubation) on pHirr and on [Ca(2+)](i) but inhibited the effects of aldosterone (1-h preincubation) on these parameters. RU 486 [glucocorticoid (GR) antagonist] and dimethyl-BAPTA (Ca(2+) chelator) prevented the effect of aldosterone on both parameters. The data indicate a genomic (1 h, via MR) and a nongenomic action (15 or 2 min, probably via GR) on the H(+)-ATPase and on [Ca(2+)](i). The results are compatible with stimulation of the H(+)-ATPase by increases in [Ca(2+)](i) (at 10(-12)-10(-6) M aldosterone) and inhibition of the H(+)-ATPase by decreases in [Ca(2+)](i) (at 10(-12) or 10(-6) M aldosterone plus RU 486).

Nonrecurrent MECP2 duplications mediated by genomic architecture-driven DNA breaks and break-induced replication repair

Recurrent submicroscopic genomic copy number changes are the result of nonallelic homologous recombination (NAHR). Nonrecurrent aberrations, however, can result from different nonexclusive recombination-repair mechanisms. We previously described small microduplications at Xq28 containing MECP2 in four male patients with a severe neurological phenotype. Here, we report on the fine-mapping and breakpoint analysis of 16 unique microduplications. The size of the overlapping copy number changes varies between 0.3 and 2.3 Mb, and FISH analysis on three patients demonstrated a tandem orientation. Although eight of the 32 breakpoint regions coincide with low-copy repeats, none of the duplications are the result of NAHR. Bioinformatics analysis of the breakpoint regions demonstrated a 2.5-fold higher frequency of Alu interspersed repeats as compared with control regions, as well as a very high GC content (53%). Unexpectedly, we obtained the junction in only one patient by long-range PCR, which revealed nonhomologous end joining as the mechanism. Breakpoint analysis in two other patients by inverse PCR and subsequent array comparative genomic hybridization analysis demonstrated the presence of a second duplicated region more telomeric at Xq28, of which one copy was inserted in between the duplicated MECP2 regions. These data suggest a two-step mechanism in which part of Xq28 is first inserted near the MECP2 locus, followed by breakage-induced replication with strand invasion of the normal sister chromatid. Our results indicate that the mechanism by which copy number changes occur in regions with a complex genomic architecture can yield complex rearrangements.

Genomic Analysis of Wild Tomato Introgressions Determining Metabolism- and Yield-Associated Traits

With the aim of determining the genetic basis of metabolic regulation in tomato fruit, we constructed a detailed physical map of genomic regions spanning previously described metabolic quantitative trait loci of a Solanum pennellii introgression line population. Two genomic libraries from S. pennellii were screened with 104 colocated markers from five selected genomic regions, and a total of 614 bacterial artificial chromosome (BAC)/cosmids were identified as seed clones. Integration of sequence data with the genetic and physical maps of Solanum lycopersicum facilitated the anchoring of 374 of these BAC/cosmid clones. The analysis of this information resulted in a genome-wide map of a nondomesticated plant species and covers 10% of the physical distance of the selected regions corresponding to approximately 1% of the wild tomato genome. Comparative analyses revealed that S. pennellii and domesticated tomato genomes can be considered as largely colinear. A total of 1,238,705 bp from both BAC/cosmid ends and nine large insert clones were sequenced, annotated, and functionally categorized. The sequence data allowed the evaluation of the level of polymorphism between the wild and cultivated tomato species. An exhaustive microsynteny analysis allowed us to estimate the divergence date of S. pennellii and S. lycopersicum at 2.7 million years ago. The combined results serve as a reference for comparative studies both at the macrosyntenic and microsyntenic levels. They also provide a valuable tool for fine-mapping of quantitative trait loci in tomato. Furthermore, they will contribute to a deeper understanding of the regulatory factors underpinning metabolism and hence defining crop chemical composition.

Genomic and nongenomic dose-dependent biphasic effect of aldosterone on Na(+)/H(+) exchanger in proximal S3 segment: role of cytosolic calcium

Leite-Dellova DC, Oliveira-Souza M, Malnic G, Mello-Aires M. Genomic and nongenomic dose-dependent biphasic effect of aldosterone on Na(+)/H(+) exchanger in proximal S3 segment: role of cytosolic calcium. Am J Physiol Renal Physiol 295: F1342-F1352, 2008. First published August 20, 2008; doi:10.1152/ajprenal.00048.2008.-The effects of aldosterone on the intracellular pH recovery rate (pHirr) via Na(+)/H(+) exchanger and on the [Ca(2+)](i) were investigated in isolated rat S3 segment. Aldosterone [10(-12), 10(-10), or 10(-8) M with 1-h, 15- or 2-min preincubation (pi)] caused a dose-dependent increase in the pHirr, but aldosterone (10(-6) M with 1-h, 15- or 2-min pi) decreased it (these effects were prevented by HOE694 but not by S3226). After 1 min of aldosterone pi, there was a transient and dose-dependent increase of the [Ca(2+)](i) and after 6-min pi there was a new increase of [Ca(2+)](i) that persisted after 1 h. Spironolactone, actinomycin D, or cycloheximide did not affect the effects of aldosterone (15 -or 2-min pi) but inhibited the effects of aldosterone (1-h pi) on pHirr and on [Ca(2+)](i). RU 486 prevented the stimulatory effect of aldosterone (10(-12) M, 15 -or 2-min pi) on both parameters and maintained the inhibitory effect of aldosterone (10(-6) M, 15- or 2-min pi) on the pHirr but reversed its stimulatory effect on the [Ca(2+)](i) to an inhibitory effect. The data indicate a genomic (1 h, via MR) and a nongenomic action (15 or 2 min, probably via GR) on [Ca(2+)](i) and on the basolateral NHE1 and are compatible with stimulation of the NHE1 by increases in [Ca(2+)](i) in the lower range (at 10(-12) M aldosterone) and inhibition by increases at high levels (at 10(-6) M aldosterone) or decreases in [Ca(2+)](i) (at 10(-6) M aldosterone plus RU 486).

DGGE with genomic DNA: Suitable for detection of numerically important organisms but not for identification of the most abundant organisms

Identification of all important community members as well as of the numerically dominant members of a community are key aspects of microbial community analysis of bioreactor samples. A systematic study was conducted with artificial consortia to test whether denaturing gradient gel electrophoresis (DGCE) is a reliable technique to obtain such community data under conditions where results would not be affected by differences in DNA extraction efficiency from cells. A total of 27 consortia were established by mixing DNA extracted from Escherichia coli K12, Burkholderia cepacia and Stenotrophomonas maltophilia in different proportions. Concentrations of DNA of single organisms in the consortia were either 0.04, 0.4 or 4 ng/mu l. DGGE-PCR of genomic DNA with primer sets targeted at the V3 and V6-V8 regions of the 16S rDNA failed to detect the three community members in only 7% of consortia, but provided incorrect information about dominance or co-dominance for 85% and 89% of consortia with the primer sets for the V6-V8 and V3 regions, respectively. The high failure rate in detection of dominant B. cepacia with the primers for the V6-V8 region was attributable to a single nucleoticle primer mismatch in the target sequences of both, the forward and reverse primer. Amplification bias in PCR of E. coli and S. maltophilia for the V6-V8 region and for all three organisms for the V3 region occurred due to interference of genomic DNA in PCR-DGGE, since a nested PCR approach, where PCR-DGGE was started from mixtures of 16S rRNA genes of the organisms, provided correct information about the relative abundance of original DNA in the sample. Multiple bands were not observed in pure culture amplicons produced with the V6-V8 primer pair, but pure culture V3 DGGE profiles of E. coli, S. maltophilia and B. cepacia contained 5, 3 and 3 bands, respectively. These results demonstrate DGGE was suitable for identification of all important community members in the three-membered artificial consortium, but not for identification of the dominant organisms in this small community. Multiple DGGE bands obtained for single organisms with the V3 primer pair could greatly confound interpretation of DGGE profiles. (C) 2008 Elsevier Ltd. All rights reserved.

Genomic organization and transcription analysis of the 195-bp satellite DNA in Trypanosoma cruzi

The 195-bp satellite DNA is the most abundant Trypanosoma cruzi repetitive sequence. Here we show by RNA blotting and RT-PCR that 195 SAT is intensely transcribed. We observed a positive correlation between the level of satellite RNA and the abundance of the satellite copies in the genome of T cruzi strains and that the satellite expression is not developmentally regulated. By analyzing CL Brener individual reads, we estimated that 195 SAT corresponds to approximately 5% of the CL Brener genome. 195 SAT elements were found in only 37 annotated contigs, indicating that a large number of satellite copies were not incorporated into the assembled data. The assembled satellite units are distributed in non-syntenic regions with Trypanosoma brucei and Leishmania major genomes, enriched with surface proteins, retroelements, RHS and hypothetical proteins. Satellite repeats were not observed in annotated subtelomeric regions. We report that 12 satellite sequences are truncated by the retroelement VIPER. (C) 2008 Elsevier B.V. All rights reserved.

Defining the role of common variation in the genomic and biological architecture of adult human height

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated ∼2,000, ∼3,700 and ∼9,500 SNPs explained ∼21%, ∼24% and ∼29% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/β-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.2014

A comparison of statistical methods for genomic selection in a mice population

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genomic Strategy Identifies a Missense Mutation in WD-Repeat Domain 65 (WDR65) in an Individual With Van der Woude Syndrome

Genetic variation in the transcription factor interferon regulatory factor 6 (IRF6) causes and contributes risk for oral clefting disorders. We hypothesized that genes regulated by IRF6 are also involved in oral clefting disorders. We used five criteria to identify potential IRF6 target genes; differential gene expression in skin taken from wild-type and Irf6-deficient murine embryos, localization to the Van der Woude syndrome 2 (VWS2) locus at 1p36-1p32, overlapping expression with Irf6, presence of a conserved predicted-binding site in the promoter region, and a mutant murine phenotype that was similar to the Irf6 mutant mouse. Previously, we observed altered expression for 573 genes; 13 were located in the murine region syntenic to the VWS2 locus. Two of these genes, Wdr65 and Stratifin, met 4 of 5 criteria. Wdr65 was a novel gene that encoded a predicted protein of 1,250 amino acids with two WD domains. As potential targets for Irf6 regulation, we hypothesized that disease-causing mutations will be found in WDR65 and Stratifin in individuals with VWS or VWS-like syndromes. We identified a potentially etiologic missense mutation in WDR65 in a person with VWS who does not have an exonic mutation in IRF6. The expression and mutation data were consistent with the hypothesis that WDR65 was a novel gene involved in oral clefting. (C) 2011 Wiley-Liss, Inc.

Genomic content and new insights on the origin of the B chromosome of the cichlid fish Astatotilapia latifasciata

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Comparative chromosome mapping of repetitive sequences. Implications for genomic evolution in the fish, Hoplias malabaricus

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Surprising genomic diversity in the Neotropical fish Synbranchus marmoratus (Teleostei: Synbranchidae): how many species?

Os cromossomos e a variação de conteúdo de DNA nuclear foram estudados em sete amostras de Synbranchus marmoratus das bacias dos rios Paraguai e Paraná com o objetivo de avaliar se diferenças no cariótipo e no conteúdo de DNA nuclear poderiam fornecer informações para a caracterização de linhagens evolutivas independentes no gênero e para a elaboração de hipóteses evolutivas e biogeográficas. A ocorrência de diferentes cariótipos já foi descrita para essa espécie, entretanto, um novo citótipo foi encontrado no rio Miranda. O conteúdo de DNA nuclear mostrou uma ampla variação entre as amostras e indivíduos, com valores entre 5,2 a 9,1 pg de DNA/núcleo. Uma análise de variância confirmou a ocorrência de diferenças significativas entre as amostras. em uma série de análises, um padrão multimodal foi encontrado na distribuição do conteúdo de DNA nuclear, pelo qual várias unidades, mais ou menos discretas, foram identificadas. Combinando a fórmula cariotípica com o conteúdo de DNA nuclear, uma relação complexa entre os rios foi observada. Com base nos dados disponíveis, sugerimos que existem várias linhagens evolutivas independentes de Synbranchus marmoratus nos rios amostrados. Hipóteses biogeográficas são propostas e discutidas.