Phred-Phrap package to analyses tools: a pipeline to facilitate population genetics re-sequencing studies

BACKGROUND: Targeted re-sequencing is one of the most powerful and widely used strategies for population genetics studies because it allows an unbiased screening for variation that is suitable for a wide variety of organisms. Examples of studies that require re-sequencing data are evolutionary inferences, epidemiological studies designed to capture rare polymorphisms responsible for complex traits and screenings for mutations in families and small populations with high incidences of specific genetic diseases. Despite the advent of next-generation sequencing technologies, Sanger sequencing is still the most popular approach in population genetics studies because of the widespread availability of automatic sequencers based on capillary electrophoresis and because it is still less prone to sequencing errors, which is critical in population genetics studies. Two popular software applications for re-sequencing studies are Phred-Phrap-Consed-Polyphred, which performs base calling, alignment, graphical edition and genotype calling and DNAsp, which performs a set of population genetics analyses. These independent tools are the start and end points of basic analyses. In between the use of these tools, there is a set of basic but error-prone tasks to be performed with re-sequencing data. RESULTS: In order to assist with these intermediate tasks, we developed a pipeline that facilitates data handling typical of re-sequencing studies. Our pipeline: (1) consolidates different outputs produced by distinct Phred-Phrap-Consed contigs sharing a reference sequence; (2) checks for genotyping inconsistencies; (3) reformats genotyping data produced by Polyphred into a matrix of genotypes with individuals as rows and segregating sites as columns; (4) prepares input files for haplotype inferences using the popular software PHASE; and (5) handles PHASE output files that contain only polymorphic sites to reconstruct the inferred haplotypes including polymorphic and monomorphic sites as required by population genetics software for re-sequencing data such as DNAsp. CONCLUSION: We tested the pipeline in re-sequencing studies of haploid and diploid data in humans, plants, animals and microorganisms and observed that it allowed a substantial decrease in the time required for sequencing analyses, as well as being a more controlled process that eliminates several classes of error that may occur when handling datasets. The pipeline is also useful for investigators using other tools for sequencing and population genetics analyses

Measurement of the cross-section and charge asymmetry of WW bosons produced in proton-proton collisions at s=8\sqrt{s}=8 TeV with the ATLAS detector

This paper presents measurements of the $W^+ \rightarrow \mu^+\nu$ and $W^- \rightarrow \mu^-\nu$ cross-sections and the associated charge asymmetry as a function of the absolute pseudorapidity of the decay muon. The data were collected in proton--proton collisions at a centre-of-mass energy of 8 TeV with the ATLAS experiment at the LHC and correspond to a total integrated luminosity of 20.2~\mbox{fb^{-1}}. The precision of the cross-section measurements varies between 0.8% to 1.5% as a function of the pseudorapidity, excluding the 1.9% uncertainty on the integrated luminosity. The charge asymmetry is measured with an uncertainty between 0.002 and 0.003. The results are compared with predictions based on next-to-next-to-leading-order calculations with various parton distribution functions and have the sensitivity to discriminate between them.Comment: 38 pages in total, author list starting page 22, 5 figures, 4 tables, submitted to EPJC. All figures including auxiliary figures are available at https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/PAPERS/STDM-2017-13

Search for chargino-neutralino production with mass splittings near the electroweak scale in three-lepton final states in √s=13 TeV pp collisions with the ATLAS detector

A search for supersymmetry through the pair production of electroweakinos with mass splittings near the electroweak scale and decaying via on-shell W and Z bosons is presented for a three-lepton final state. The analyzed proton-proton collision data taken at a center-of-mass energy of √s=13  TeV were collected between 2015 and 2018 by the ATLAS experiment at the Large Hadron Collider, corresponding to an integrated luminosity of 139  fb−1. A search, emulating the recursive jigsaw reconstruction technique with easily reproducible laboratory-frame variables, is performed. The two excesses observed in the 2015–2016 data recursive jigsaw analysis in the low-mass three-lepton phase space are reproduced. Results with the full data set are in agreement with the Standard Model expectations. They are interpreted to set exclusion limits at the 95% confidence level on simplified models of chargino-neutralino pair production for masses up to 345 GeV

Avaliação da medida translucência nucal em gestantes de risco

27

Twin Town in South Brazil: A Nazi's Experiment or a Genetic Founder Effect?

Cândido Godói (CG) is a small municipality in South Brazil with approximately 6,000 inhabitants. It is known as the “Twins' Town” due to its high rate of twin births. Recently it was claimed that such high frequency of twinning would be connected to experiments performed by the German Nazi doctor Joseph Mengele. It is known, however, that this town was founded by a small number of families and therefore a genetic founder effect may represent an alternatively explanation for the high twinning prevalence in CG. In this study, we tested specific predictions of the “Nazi's experiment” and of the “founder effect” hypotheses. We surveyed a total of 6,262 baptism records from 1959–2008 in CG catholic churches, and identified 91 twin pairs and one triplet. Contrary to the “Nazi's experiment hypothesis”, there is no spurt in twinning between the years (1964–1968) when Mengele allegedly was in CG (P = 0.482). Moreover, there is no temporal trend for a declining rate of twinning since the 1960s (P = 0.351), and no difference in twinning among CG districts considering two different periods: 1927–1958 and 1959–2008 (P = 0.638). On the other hand, the “founder effect hypothesis” is supported by an isonymy analysis that shows that women who gave birth to twins have a higher inbreeding coefficient when compared to women who never had twins (0.0148, 0.0081, respectively, P = 0.019). In summary, our results show no evidence for the “Nazi's experiment hypothesis” and strongly suggest that the “founder effect hypothesis” is a much more likely alternative for explaining the high prevalence of twinning in CG. If this hypothesis is correct, then this community represents a valuable population where genetic factors linked to twinning may be identified

A genome-wide association study of asthma symptoms in Latin American children

Abstract\ud \ud Background\ud Asthma is a chronic disease of the airways and, despite the advances in the knowledge of associated genetic regions in recent years, their mechanisms have yet to be explored. Several genome-wide association studies have been carried out in recent years, but none of these have involved Latin American populations with a high level of miscegenation, as is seen in the Brazilian population.\ud \ud \ud Methods\ud 1246 children were recruited from a longitudinal cohort study in Salvador, Brazil. Asthma symptoms were identified in accordance with an International Study of Asthma and Allergies in Childhood (ISAAC) questionnaire. Following quality control, 1 877 526 autosomal SNPs were tested for association with childhood asthma symptoms by logistic regression using an additive genetic model. We complemented the analysis with an estimate of the phenotypic variance explained by common genetic variants. Replications were investigated in independent Mexican and US Latino samples.\ud \ud \ud Results\ud Two chromosomal regions reached genome-wide significance level for childhood asthma symptoms: the 14q11 region flanking the DAD1 and OXA1L genes (rs1999071, MAF 0.32, OR 1.78, 95 % CI 1.45–2.18, p-value 2.83 × 10−8) and 15q22 region flanking the FOXB1 gene (rs10519031, MAF 0.04, OR 3.0, 95 % CI 2.02–4.49, p-value 6.68 × 10−8 and rs8029377, MAF 0.03, OR 2.49, 95 % CI 1.76–3.53, p-value 2.45 × 10−7). eQTL analysis suggests that rs1999071 regulates the expression of OXA1L gene. However, the original findings were not replicated in the Mexican or US Latino samples.\ud \ud \ud Conclusions\ud We conclude that the 14q11 and 15q22 regions may be associated with asthma symptoms in childhood.This work was supported by the Department of Science and Technology\ud (DECIT, Ministry of Health), National Fund for Scientific and Technological\ud Development (FNDCT, Ministry of Science and Technology), Funding of\ud Studies and Projects (FINEP, Ministry of Science and Technology, Brazil), the\ud Brazilian National Research Council (CNPq) and the Wellcome Trust UK, Ref\ud 072405/Z/03/Z.\ud E.G.B. was funded by grants from National Institutes of Health (HL088133,\ud HL078885, HL004464, HL104608, HL117004, ES015794 and MD006902) and\ud by the American Asthma Foundation, the Sandler Foundation and the RWJF\ud Amos Medical Faculty Development Award.\ud Supported in part by the Intramural Research Program of the NIH, National\ud Institute of Environmental Health Sciences, USA

Basin-wide variation in tree hydraulic safety margins predicts the carbon balance of Amazon forests

Funding: Data collection was largely funded by the UK Natural Environment Research Council (NERC) project TREMOR (NE/N004655/1) to D.G., E.G. and O.P., with further funds from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES, finance code 001) to J.V.T. and a University of Leeds Climate Research Bursary Fund to J.V.T. D.G., E.G. and O.P. acknowledge further support from a NERC-funded consortium award (ARBOLES, NE/S011811/1). This paper is an outcome of J.V.T.’s doctoral thesis, which was sponsored by CAPES (GDE 99999.001293/2015-00). J.V.T. was previously supported by the NERC-funded ARBOLES project (NE/S011811/1) and is supported at present by the Swedish Research Council Vetenskapsrådet (grant no. 2019-03758 to R.M.). E.G., O.P. and D.G. acknowledge support from NERC-funded BIORED grant (NE/N012542/1). O.P. acknowledges support from an ERC Advanced Grant and a Royal Society Wolfson Research Merit Award. R.S.O. was supported by a CNPq productivity scholarship, the São Paulo Research Foundation (FAPESP-Microsoft 11/52072-0) and the US Department of Energy, project GoAmazon (FAPESP 2013/50531-2). M.M. acknowledges support from MINECO FUN2FUN (CGL2013-46808-R) and DRESS (CGL2017-89149-C2-1-R). C.S.-M., F.B.V. and P.R.L.B. were financed by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES, finance code 001). C.S.-M. received a scholarship from the Brazilian National Council for Scientific and Technological Development (CNPq 140353/2017-8) and CAPES (science without borders 88881.135316/2016-01). Y.M. acknowledges the Gordon and Betty Moore Foundation and ERC Advanced Investigator Grant (GEM-TRAITS, 321131) for supporting the Global Ecosystems Monitoring (GEM) network (gem.tropicalforests.ox.ac.uk), within which some of the field sites (KEN, TAM and ALP) are nested. The authors thank Brazil–USA Collaborative Research GoAmazon DOE-FAPESP-FAPEAM (FAPESP 2013/50533-5 to L.A.) and National Science Foundation (award DEB-1753973 to L. Alves). They thank Serrapilheira Serra-1709-18983 (to M.H.) and CNPq-PELD/POPA-441443/2016-8 (to L.G.) (P.I. Albertina Lima). They thank all the colleagues and grants mentioned elsewhere [8,36] that established, identified and measured the Amazon forest plots in the RAINFOR network analysed here. The authors particularly thank J. Lyod, S. Almeida, F. Brown, B. Vicenti, N. Silva and L. Alves. This work is an outcome approved Research Project no. 19 from ForestPlots.net, a collaborative initiative developed at the University of Leeds that unites researchers and the monitoring of their permanent plots from the world’s tropical forests [61]. The authros thank A. Levesley, K. Melgaço Ladvocat and G. Pickavance for ForestPlots.net management. They thank Y. Wang and J. Baker, respectively, for their help with the map and with the climatic data. The authors acknowledge the invaluable help of M. Brum for kindly providing the comparison of vulnerability curves based on PAD and on PLC shown in this manuscript. They thank J. Martinez-Vilalta for his comments on an early version of this manuscript. The authors also thank V. Hilares and the Asociación para la Investigación y Desarrollo Integral (AIDER, Puerto Maldonado, Peru); V. Saldaña and Instituto de Investigaciones de la Amazonía Peruana (IIAP) for local field campaign support in Peru; E. Chavez and Noel Kempff Natural History Museum for local field campaign support in Bolivia; ICMBio, INPA/NAPPA/LBA COOMFLONA (Cooperativa mista da Flona Tapajós) and T. I. Bragança-Marituba for the research support.Tropical forests face increasing climate risk1,2, yet our ability to predict their response to climate change is limited by poor understanding of their resistance to water stress. Although xylem embolism resistance thresholds (for example, Ψ50) and hydraulic safety margins (for example, HSM50) are important predictors of drought-induced mortality risk3-5, little is known about how these vary across Earth's largest tropical forest. Here, we present a pan-Amazon, fully standardized hydraulic traits dataset and use it to assess regional variation in drought sensitivity and hydraulic trait ability to predict species distributions and long-term forest biomass accumulation. Parameters Ψ50 and HSM50 vary markedly across the Amazon and are related to average long-term rainfall characteristics. Both Ψ50 and HSM50 influence the biogeographical distribution of Amazon tree species. However, HSM50 was the only significant predictor of observed decadal-scale changes in forest biomass. Old-growth forests with wide HSM50 are gaining more biomass than are low HSM50 forests. We propose that this may be associated with a growth-mortality trade-off whereby trees in forests consisting of fast-growing species take greater hydraulic risks and face greater mortality risk. Moreover, in regions of more pronounced climatic change, we find evidence that forests are losing biomass, suggesting that species in these regions may be operating beyond their hydraulic limits. Continued climate change is likely to further reduce HSM50 in the Amazon6,7, with strong implications for the Amazon carbon sink.Publisher PDFPeer reviewe

Geographic patterns of tree dispersal modes in Amazonia and their ecological correlates

Aim: To investigate the geographic patterns and ecological correlates in the geographic distribution of the most common tree dispersal modes in Amazonia (endozoochory, synzoochory, anemochory and hydrochory). We examined if the proportional abundance of these dispersal modes could be explained by the availability of dispersal agents (disperser-availability hypothesis) and/or the availability of resources for constructing zoochorous fruits (resource-availability hypothesis). Time period: Tree-inventory plots established between 1934 and 2019. Major taxa studied: Trees with a diameter at breast height (DBH) ≥ 9.55 cm. Location: Amazonia, here defined as the lowland rain forests of the Amazon River basin and the Guiana Shield. Methods: We assigned dispersal modes to a total of 5433 species and morphospecies within 1877 tree-inventory plots across terra-firme, seasonally flooded, and permanently flooded forests. We investigated geographic patterns in the proportional abundance of dispersal modes. We performed an abundance-weighted mean pairwise distance (MPD) test and fit generalized linear models (GLMs) to explain the geographic distribution of dispersal modes. Results: Anemochory was significantly, positively associated with mean annual wind speed, and hydrochory was significantly higher in flooded forests. Dispersal modes did not consistently show significant associations with the availability of resources for constructing zoochorous fruits. A lower dissimilarity in dispersal modes, resulting from a higher dominance of endozoochory, occurred in terra-firme forests (excluding podzols) compared to flooded forests. Main conclusions: The disperser-availability hypothesis was well supported for abiotic dispersal modes (anemochory and hydrochory). The availability of resources for constructing zoochorous fruits seems an unlikely explanation for the distribution of dispersal modes in Amazonia. The association between frugivores and the proportional abundance of zoochory requires further research, as tree recruitment not only depends on dispersal vectors but also on conditions that favour or limit seedling recruitment across forest types