Use of Bulk Segregant Analysis for Determining the Genetic Basis of Azole Resistance in the Opportunistic Pathogen Aspergillus fumigatus.

A sexual cycle was described in 2009 for the opportunistic fungal pathogen Aspergillus fumigatus, opening up for the first time the possibility of using techniques reliant on sexual crossing for genetic analysis. The present study was undertaken to evaluate whether the technique 'bulk segregant analysis' (BSA), which involves detection of differences between pools of progeny varying in a particular trait, could be applied in conjunction with next-generation sequencing to investigate the underlying basis of monogenic traits in A. fumigatus. Resistance to the azole antifungal itraconazole was chosen as a model, with a dedicated bioinformatic pipeline developed to allow identification of SNPs that differed between the resistant progeny pool and resistant parent compared to the sensitive progeny pool and parent. A clinical isolate exhibiting monogenic resistance to itraconazole of unknown basis was crossed to a sensitive parent and F1 progeny used in BSA. In addition, the use of backcrossing and increasing the number in progeny pools was evaluated as ways to enhance the efficiency of BSA. Use of F1 pools of 40 progeny led to the identification of 123 candidate genes with SNPs distributed over several contigs when aligned to an A1163 reference genome. Successive rounds of backcrossing enhanced the ability to identify specific genes and a genomic region, with BSA of progeny (using 40 per pool) from a third backcross identifying 46 genes with SNPs, and BSA of progeny from a sixth backcross identifying 20 genes with SNPs in a single 292 kb region of the genome. The use of an increased number of 80 progeny per pool also increased the resolution of BSA, with 29 genes demonstrating SNPs between the different sensitive and resistant groupings detected using progeny from just the second backcross with the majority of variants located on the same 292 kb region. Further bioinformatic analysis of the 292 kb region identified the presence of a cyp51A gene variant resulting in a methionine to lysine (M220K) change in the CYP51A protein, which was concluded to be the causal basis of the observed resistance to itraconazole. The future use of BSA in genetic analysis of A. fumigatus is discussed.

Subdividing Y-chromosome haplogroup R1a1 reveals Norse Viking dispersal lineages in Britain.

The influence of Viking-Age migrants to the British Isles is obvious in archaeological and place-names evidence, but their demographic impact has been unclear. Autosomal genetic analyses support Norse Viking contributions to parts of Britain, but show no signal corresponding to the Danelaw, the region under Scandinavian administrative control from the ninth to eleventh centuries. Y-chromosome haplogroup R1a1 has been considered as a possible marker for Viking migrations because of its high frequency in peninsular Scandinavia (Norway and Sweden). Here we select ten Y-SNPs to discriminate informatively among hg R1a1 sub-haplogroups in Europe, analyse these in 619 hg R1a1 Y chromosomes including 163 from the British Isles, and also type 23 short-tandem repeats (Y-STRs) to assess internal diversity. We find three specifically Western-European sub-haplogroups, two of which predominate in Norway and Sweden, and are also found in Britain; star-like features in the STR networks of these lineages indicate histories of expansion. We ask whether geographical distributions of hg R1a1 overall, and of the two sub-lineages in particular, correlate with regions of Scandinavian influence within Britain. Neither shows any frequency difference between regions that have higher (≥10%) or lower autosomal contributions from Norway and Sweden, but both are significantly overrepresented in the region corresponding to the Danelaw. These differences between autosomal and Y-chromosomal histories suggest either male-specific contribution, or the influence of patrilocality. Comparison of modern DNA with recently available ancient DNA data supports the interpretation that two sub-lineages of hg R1a1 spread with the Vikings from peninsular Scandinavia.

Human Y Chromosome Exerts Pleiotropic Effects on Susceptibility to Atherosclerosis.

OBJECTIVE: The male-specific region of the Y chromosome (MSY) remains one of the most unexplored regions of the genome. We sought to examine how the genetic variants of the MSY influence male susceptibility to coronary artery disease (CAD) and atherosclerosis. Approach and Results: Analysis of 129 133 men from UK Biobank revealed that only one of 7 common MSY haplogroups (haplogroup I1) was associated with CAD-carriers of haplogroup I1 had ≈11% increase in risk of CAD when compared with all other haplogroups combined (odds ratio, 1.11; 95% CI, 1.04-1.18; P=6.8×10-4). Targeted MSY sequencing uncovered 235 variants exclusive to this haplogroup. The haplogroup I1-specific variants showed 2.45- and 1.56-fold respective enrichment for promoter and enhancer chromatin states, in cells/tissues relevant to atherosclerosis, when compared with other MSY variants. Gene set enrichment analysis in CAD-relevant tissues showed that haplogroup I1 was associated with changes in pathways responsible for early and late stages of atherosclerosis development including defence against pathogens, immunity, oxidative phosphorylation, mitochondrial respiration, lipids, coagulation, and extracellular matrix remodeling. UTY was the only Y chromosome gene whose blood expression was associated with haplogroup I1. Experimental reduction of UTY expression in macrophages led to changes in expression of 59 pathways (28 of which overlapped with those associated with haplogroup I1) and a significant reduction in the immune costimulatory signal. CONCLUSIONS: Haplogroup I1 is enriched for regulatory chromatin variants in numerous cells of relevance to CAD and increases cardiovascular risk through proatherosclerotic reprogramming of the transcriptome, partly through UTY.

Population resequencing of European mitochondrial genomes highlights sex-bias in Bronze Age demographic expansions.

Interpretations of genetic data concerning the prehistory of Europe have long been a subject of great debate, but increasing amounts of ancient and modern DNA data are now providing new and more informative evidence. Y-chromosome resequencing studies in Europe have highlighted the prevalence of recent expansions of male lineages, and focused interest on the Bronze Age as a period of cultural and demographic change. These findings contrast with phylogeographic studies based on mitochondrial DNA (mtDNA), which have been interpreted as supporting expansions from glacial refugia. Here we have undertaken a population-based resequencing of complete mitochondrial genomes in Europe and the Middle East, in 340 samples from 17 populations for which Y-chromosome sequence data are also available. Demographic reconstructions show no signal of Bronze Age expansion, but evidence of Paleolithic expansions in all populations except the Saami, and with an absence of detectable geographical pattern. In agreement with previous inference from modern and ancient DNA data, the unbiased comparison between the mtDNA and Y-chromosome population datasets emphasizes the sex-biased nature of recent demographic transitions in Europe.

Signatures of human European Palaeolithic expansion shown by resequencing of non-recombining X-chromosome segments.

Human genetic diversity in Europe has been extensively studied using uniparentally inherited sequences (mitochondrial DNA (mtDNA) and the Y chromosome), which reveal very different patterns indicating sex-specific demographic histories. The X chromosome, haploid in males and inherited twice as often from mothers as from fathers, could provide insights into past female behaviours, but has not been extensively investigated. Here, we use HapMap single-nucleotide polymorphism data to identify genome-wide segments of the X chromosome in which recombination is historically absent and mutations are likely to be the only source of genetic variation, referring to these as phylogeographically informative haplotypes on autosomes and X chromosome (PHAXs). Three such sequences on the X chromosome spanning a total of ~49 kb were resequenced in 240 males from Europe, the Middle East and Africa at an average coverage of 181 ×. These PHAXs were confirmed to be essentially non-recombining across European samples. All three loci show highly homogeneous patterns across Europe and are highly differentiated from the African sample. Star-like structures of European-specific haplotypes in median-joining networks indicate past population expansions. Bayesian skyline plots and time-to-most-recent-common-ancestor estimates suggest expansions pre-dating the Neolithic transition, a finding that is more compatible with data on mtDNA than the Y chromosome, and with the female bias of X-chromosomal inheritance. This study demonstrates the potential of the use of X-chromosomal haplotype blocks, and the utility of the accurate ascertainment of rare variants for inferring human demographic history.

Great ape Y Chromosome and mitochondrial DNA phylogenies reflect subspecies structure and patterns of mating and dispersal.

The distribution of genetic diversity in great ape species is likely to have been affected by patterns of dispersal and mating. This has previously been investigated by sequencing autosomal and mitochondrial DNA (mtDNA), but large-scale sequence analysis of the male-specific region of the Y Chromosome (MSY) has not yet been undertaken. Here, we use the human MSY reference sequence as a basis for sequence capture and read mapping in 19 great ape males, combining the data with sequences extracted from the published whole genomes of 24 additional males to yield a total sample of 19 chimpanzees, four bonobos, 14 gorillas, and six orangutans, in which interpretable MSY sequence ranges from 2.61 to 3.80 Mb. This analysis reveals thousands of novel MSY variants and defines unbiased phylogenies. We compare these with mtDNA-based trees in the same individuals, estimating time-to-most-recent common ancestor (TMRCA) for key nodes in both cases. The two loci show high topological concordance and are consistent with accepted (sub)species definitions, but time depths differ enormously between loci and (sub)species, likely reflecting different dispersal and mating patterns. Gorillas and chimpanzees/bonobos present generally low and high MSY diversity, respectively, reflecting polygyny versus multimale-multifemale mating. However, particularly marked differences exist among chimpanzee subspecies: The western chimpanzee MSY phylogeny has a TMRCA of only 13.2 (10.8-15.8) thousand years, but that for central chimpanzees exceeds 1 million years. Cross-species comparison within a single MSY phylogeny emphasizes the low human diversity, and reveals species-specific branch length variation that may reflect differences in long-term generation times.

Large-scale recent expansion of European patrilineages shown by population resequencing.

The proportion of Europeans descending from Neolithic farmers ∼ 10 thousand years ago (KYA) or Palaeolithic hunter-gatherers has been much debated. The male-specific region of the Y chromosome (MSY) has been widely applied to this question, but unbiased estimates of diversity and time depth have been lacking. Here we show that European patrilineages underwent a recent continent-wide expansion. Resequencing of 3.7 Mb of MSY DNA in 334 males, comprising 17 European and Middle Eastern populations, defines a phylogeny containing 5,996 single-nucleotide polymorphisms. Dating indicates that three major lineages (I1, R1a and R1b), accounting for 64% of our sample, have very recent coalescent times, ranging between 3.5 and 7.3 KYA. A continuous swathe of 13/17 populations share similar histories featuring a demographic expansion starting ∼ 2.1-4.2 KYA. Our results are compatible with ancient MSY DNA data, and contrast with data on mitochondrial DNA, indicating a widespread male-specific phenomenon that focuses interest on the social structure of Bronze Age Europe.

The Y-chromosome tree bursts into leaf: 13,000 high-confidence SNPs covering the majority of known clades.

Many studies of human populations have used the male-specific region of the Y chromosome (MSY) as a marker, but MSY sequence variants have traditionally been subject to ascertainment bias. Also, dating of haplogroups has relied on Y-specific short tandem repeats (STRs), involving problems of mutation rate choice, and possible long-term mutation saturation. Next-generation sequencing can ascertain single nucleotide polymorphisms (SNPs) in an unbiased way, leading to phylogenies in which branch-lengths are proportional to time, and allowing the times-to-most-recent-common-ancestor (TMRCAs) of nodes to be estimated directly. Here we describe the sequencing of 3.7 Mb of MSY in each of 448 human males at a mean coverage of 51×, yielding 13,261 high-confidence SNPs, 65.9% of which are previously unreported. The resulting phylogeny covers the majority of the known clades, provides date estimates of nodes, and constitutes a robust evolutionary framework for analyzing the history of other classes of mutation. Different clades within the tree show subtle but significant differences in branch lengths to the root. We also apply a set of 23 Y-STRs to the same samples, allowing SNP- and STR-based diversity and TMRCA estimates to be systematically compared. Ongoing purifying selection is suggested by our analysis of the phylogenetic distribution of nonsynonymous variants in 15 MSY single-copy genes.

Structure-function analysis of the presumptive Arabidopsis auxin permease AUX1.

We have investigated the subcellular localization, the domain topology, and the amino acid residues that are critical for the function of the presumptive Arabidopsis thaliana auxin influx carrier AUX1. Biochemical fractionation experiments and confocal studies using an N-terminal yellow fluorescent protein (YFP) fusion observed that AUX1 colocalized with plasma membrane (PM) markers. Because of its PM localization, we were able to take advantage of the steep pH gradient that exists across the plant cell PM to investigate AUX1 topology using YFP as a pH-sensitive probe. The YFP-coding sequence was inserted in selected AUX1 hydrophilic loops to orient surface domains on either apoplastic or cytoplasmic faces of the PM based on the absence or presence of YFP fluorescence, respectively. We were able to demonstrate in conjunction with helix prediction programs that AUX1 represents a polytopic membrane protein composed of 11 transmembrane spanning domains. In parallel, a large aux1 allelic series containing null, partial-loss-of-function, and conditional mutations was characterized to identify the functionally important domains and amino acid residues within the AUX1 polypeptide. Whereas almost all partial-loss-of-function and null alleles cluster in the core permease region, the sole conditional allele aux1-7 modifies the function of the external C-terminal domain.