Chromosomal Copy Number Variation in Saccharomyces pastorianus Is Evidence for Extensive Genome Dynamics in Industrial Lager Brewing Strains

Lager brewing strains of Saccharomyces pastorianus are natural interspecific hybrids originating from the spontaneous hybridization of Saccharomyces cerevisiae and Saccharomyces eubayanus. Over the past 500 years, S. pastorianus has been domesticated to become one of the most important industrial microorganisms. Production of lager-type beers requires a set of essential phenotypes, including the ability to ferment maltose and maltotriose at low temperature, the production of flavors and aromas, and the ability to flocculate. Understanding of the molecular basis of complex brewing-related phenotypic traits is a prerequisite for rational strain improvement. While genome sequences have been reported, the variability and dynamics of S. pastorianus genomes have not been investigated in detail. Here, using deep sequencing and chromosome copy number analysis, we showed that S. pastorianus strain CBS1483 exhibited extensive aneuploidy. This was confirmed by quantitative PCR and by flow cytometry. As a direct consequence of this aneuploidy, a massive number of sequence variants was identified, leading to at least 1,800 additional protein variants in S. pastorianus CBS1483. Analysis of eight additional S. pastorianus strains revealed that the previously defined group I strains showed comparable karyotypes, while group II strains showed large interstrain karyotypic variability. Comparison of three strains with nearly identical genome sequences revealed substantial chromosome copy number variation, which may contribute to strain-specific phenotypic traits. The observed variability of lager yeast genomes demonstrates that systematic linking of genotype to phenotype requires a three-dimensional genome analysis encompassing physical chromosomal structures, the copy number of individual chromosomes or chromosomal regions, and the allelic variation of copies of individual genes.     

Copy Number Variation Is a Fundamental Aspect of the Placental Genome

Discovery of lineage-specific somatic copy number variation (CNV) in mammals has led to debate over whether CNVs are mutations that propagate disease or whether they are a normal, and even essential, aspect of cell biology. We show that 1,000N polyploid trophoblast giant cells (TGCs) of the mouse placenta contain 47 regions, totaling 138 Megabases, where genomic copies are underrepresented (UR). UR domains originate from a subset of late-replicating heterochromatic regions containing gene deserts and genes involved in cell adhesion and neurogenesis. While lineage-specific CNVs have been identified in mammalian cells, classically in the immune system where V(D)J recombination occurs, we demonstrate that CNVs form during gestation in the placenta by an underreplication mechanism, not by recombination nor deletion. Our results reveal that large scale CNVs are a normal feature of the mammalian placental genome, which are regulated systematically during embryogenesis and are propagated by a mechanism of underreplication.     

Whole Genome Distribution and Ethnic Differentiation of Copy Number Variation in Caucasian and Asian Populations

Although copy number variation (CNV) has recently received much attention as a form of structure variation within the human genome, knowledge is still inadequate on fundamental CNV characteristics such as occurrence rate, genomic distribution and ethnic differentiation. In the present study, we used the Affymetrix GeneChip® Mapping 500K Array to discover and characterize CNVs in the human genome and to study ethnic differences of CNVs between Caucasians and Asians. Three thousand and nineteen CNVs, including 2381 CNVs in autosomes and 638 CNVs in X chromosome, from 985 Caucasian and 692 Asian individuals were identified, with a mean length of 296 kb. Among these CNVs, 190 had frequencies greater than 1% in at least one ethnic group, and 109 showed significant ethnic differences in frequencies (p<0.01). After merging overlapping CNVs, 1135 copy number variation regions (CNVRs), covering approximately 439 Mb (14.3%) of the human genome, were obtained. Our findings of ethnic differentiation of CNVs, along with the newly constructed CNV genomic map, extend our knowledge on the structural variation in the human genome and may furnish a basis for understanding the genomic differentiation of complex traits across ethnic groups.     

Maize Inbreds Exhibit High Levels of Copy Number Variation (CNV) and Presence/Absence Variation (PAV) in Genome Content

Following the domestication of maize over the past ∼10,000 years, breeders have exploited the extensive genetic diversity of this species to mold its phenotype to meet human needs. The extent of structural variation, including copy number variation (CNV) and presence/absence variation (PAV), which are thought to contribute to the extraordinary phenotypic diversity and plasticity of this important crop, have not been elucidated. Whole-genome, array-based, comparative genomic hybridization (CGH) revealed a level of structural diversity between the inbred lines B73 and Mo17 that is unprecedented among higher eukaryotes. A detailed analysis of altered segments of DNA conservatively estimates that there are several hundred CNV sequences among the two genotypes, as well as several thousand PAV sequences that are present in B73 but not Mo17. Haplotype-specific PAVs contain hundreds of single-copy, expressed genes that may contribute to heterosis and to the extraordinary phenotypic diversity of this important crop.     

The Necessity of DNA Taxonomy to Reveal Cryptic Diversity and Spatial Distribution of Meiofauna,with a Focus on Nemertea

Meiofauna represent one of the most abundant and diverse communities in marine benthic ecosystems. However, an accurate assessment of diversity at the level of species has been and remains challenging for these microscopic organisms. Therefore, for many taxa, especially the soft body forms such as nemerteans, which often lack clear diagnostic morphological traits, DNA taxonomy is an effective means to assess species diversity. Morphological taxonomy of Nemertea is well documented as complicated by scarcity of unambiguous character states and compromised by diagnoses of a majority of species (and higher clades) being inadequate or based on ambiguous characters and character states. Therefore, recent studies have advocated for the primacy of molecular tools to solve the taxonomy of this group. DNA taxonomy uncovers possible hidden cryptic species, provides a coherent means to systematize taxa in definite clades, and also reveals possible biogeographic patterns. Here, we analyze diversity of nemertean species by considering the barcode region of the mitochondrial gene Cytochrome Oxidase subunit I (COI) and different species delineation approaches in order to infer evolutionarily significant units. In the aim to uncover actual diversity of meiofaunal nemerteans across different sites in Central America, COI sequences were obtained for specimens assigned here to the genera Cephalothrix, Ototyphlonemertes, and Tetrastemma-like worms, each commonly encountered in our sampling. Additional genetic, taxonomic, and geographic data of other specimens belonging to these genera were added from GenBank. Results are consistent across different DNA taxonomy approaches, and revealed (i) the presence of several hidden cryptic species and (ii) numerous potential misidentifications due to traditional taxonomy. (iii) We additionally test a possible biogeographic pattern of taxonomic units revealed by this study, and, except for a few cases, the putative species seem not to be widely distributed, in contrast to what traditional taxonomy would suggest for the recognized morphotypes.     

Algorithms to Model Single Gene,Single Chromosome,and Whole Genome Copy Number Changes Jointly in Tumor Phylogenetics

We present methods to construct phylogenetic models of tumor progression at the cellular level that include copy number changes at the scale of single genes, entire chromosomes, and the whole genome. The methods are designed for data collected by fluorescence in situ hybridization (FISH), an experimental technique especially well suited to characterizing intratumor heterogeneity using counts of probes to genetic regions frequently gained or lost in tumor development. Here, we develop new provably optimal methods for computing an edit distance between the copy number states of two cells given evolution by copy number changes of single probes, all probes on a chromosome, or all probes in the genome. We then apply this theory to develop a practical heuristic algorithm, implemented in publicly available software, for inferring tumor phylogenies on data from potentially hundreds of single cells by this evolutionary model. We demonstrate and validate the methods on simulated data and published FISH data from cervical cancers and breast cancers. Our computational experiments show that the new model and algorithm lead to more parsimonious trees than prior methods for single-tumor phylogenetics and to improved performance on various classification tasks, such as distinguishing primary tumors from metastases obtained from the same patient population.     

Selection,Recombination, and Virulence Gene Diversity among Group B Streptococcal Genotypes

Transmission of group B Streptococcus (GBS) from mothers to neonates during childbirth is a leading cause of neonatal sepsis and meningitis. Although subtyping tools have identified specific GBS phylogenetic lineages that are important in neonatal disease, little is known about the genetic diversity of these lineages or the roles that recombination and selection play in the generation of emergent genotypes. Here, we examined genetic variation, selection, and recombination in seven multilocus sequence typing (MLST) loci from 94 invasive, colonizing, and bovine strains representing 38 GBS sequence types and performed DNA sequencing and PCR-based restriction fragment length polymorphism analysis of several putative virulence genes to identify gene content differences between genotypes. Despite the low level of diversity in the MLST loci, a neighbor net analysis revealed a variable range of genetic exchange among the seven clonal complexes (CCs) identified, suggesting that recombination is partly responsible for the diversity observed between genotypes. Recombination is also important for several virulence genes, as some gene alleles had evidence for lateral gene exchange across divergent genotypes. The CC-17 lineage, which is associated with neonatal disease, is relatively homogeneous and therefore appears to have diverged independently with an exclusive set of virulence characteristics. These data suggest that different GBS genetic backgrounds have distinct virulence gene profiles that may be important for disease pathogenesis. Such profiles could be used as markers for the rapid detection of strains with an increased propensity to cause neonatal disease and may be considered useful vaccine targets.Group B Streptococcus (GBS) is a leading cause of neonatal sepsis, pneumonia, and meningitis (51) and causes infections in pregnant women, nonpregnant adults, and the elderly with underlying medical conditions. Maternal GBS colonization is a main risk factor for neonatal disease, and roughly 20 to 40% of pregnant women are colonized (14, 23). Colonization rates of up to 31% and 34% have been documented in young men (4) and nonpregnant women (4, 42), respectively, whereas a rate of 22% has been observed in individuals over 65 years of age (18). GBS has also been identified as the cause of bovine mastitis in up to 45% of symptomatic bovines (30). Nine distinct polysaccharide capsule types (serotypes) are known, and the serotype distribution varies by population.The genetic diversity of GBS populations has been studied using a variety of different methods, including restriction fragment length polymorphism (RFLP) (24), ribotyping (5, 25), pulsed-field gel electrophoresis (49), multilocus enzyme electrophoresis (MLEE) (45), random amplification of polymorphic DNA (36), restriction digestion pattern (RDP) typing (53), and multilocus sequence typing (MLST) (28). By utilizing methods that focus on conserved genetic changes within GBS strains, virulent GBS clones that have diversified genetically can be identified. Both MLEE and MLST can distinguish the major GBS serotype III clones associated with neonatal invasive disease as sequence type 17 (ST-17) in the MLST system (28, 29, 38) or electrophoretic type 1 in the MLEE system (45). This clone is also evident in the RDP system as RDP-III (53).A recent study of 75 GBS strains representing different sources and STs reported that the ST-17 lineage is relatively homogeneous and contains a unique set of surface proteins (9). Homogeneity within a GBS lineage that is significantly associated with neonatal disease is likely important for disease pathogenesis, though few studies have been conducted to identify specific differences in virulence characteristics between lineages. Similarly, the roles of selection and recombination in the generation of STs, as well as known virulence genes, have only recently been explored and require further investigation (9a). Here, we assess the genomic diversity of GBS strains representing a variety of common clonal genotypes, examine evidence for selection and recombination, and evaluate the extent of DNA polymorphism and allelic variation in several putative virulence genes.     

Dietary Variation and Evolution of Gene Copy Number among Dog Breeds

Prolonged human interactions and artificial selection have influenced the genotypic and phenotypic diversity among dog breeds. Because humans and dogs occupy diverse habitats, ecological contexts have likely contributed to breed-specific positive selection. Prior to the advent of modern dog-feeding practices, there was likely substantial variation in dietary landscapes among disparate dog breeds. As such, we investigated one type of genetic variant, copy number variation, in three metabolic genes: glucokinase regulatory protein (GCKR), phytanol-CoA 2-hydroxylase (PHYH), and pancreatic α-amylase 2B (AMY2B). These genes code for proteins that are responsible for metabolizing dietary products that originate from distinctly different food types: sugar, meat, and starch, respectively. After surveying copy number variation among dogs with diverse dietary histories, we found no correlation between diet and positive selection in either GCKR or PHYH. Although it has been previously demonstrated that dogs experienced a copy number increase in AMY2B relative to wolves during or after the dog domestication process, we demonstrate that positive selection continued to act on amylase copy number in dog breeds that consumed starch-rich diets in time periods after domestication. Furthermore, we found that introgression with wolves is not responsible for deterioration of positive selection on AMY2B among diverse dog breeds. Together, this supports the hypothesis that the amylase copy number expansion is found universally in dogs.     

SomatiCA: Identifying,Characterizing and Quantifying Somatic Copy Number Aberrations from Cancer Genome Sequencing Data

Whole genome sequencing of matched tumor-normal sample pairs is becoming routine in cancer research. However, analysis of somatic copy-number changes from sequencing data is still challenging because of insufficient sequencing coverage, unknown tumor sample purity and subclonal heterogeneity. Here we describe a computational framework, named SomatiCA, which explicitly accounts for tumor purity and subclonality in the analysis of somatic copy-number profiles. Taking read depths (RD) and lesser allele frequencies (LAF) as input, SomatiCA will output 1) admixture rate for each tumor sample, 2) somatic allelic copy-number for each genomic segment, 3) fraction of tumor cells with subclonal change in each somatic copy number aberration (SCNA), and 4) a list of substantial genomic aberration events including gain, loss and LOH. SomatiCA is available as a Bioconductor R package at http://www.bioconductor.org/packages/2.13/bioc/html/SomatiCA.html.     

MCM8- and MCM9-Deficient Mice Reveal Gametogenesis Defects and Genome Instability Due to Impaired Homologous Recombination

We generated knockout mice for MCM8 and MCM9 and show that deficiency for these genes impairs homologous recombination (HR)-mediated DNA repair during gametogenesis and somatic cells cycles. MCM8(-/-) mice are sterile because spermatocytes are blocked in meiotic prophase I, and females have only arrested primary follicles and frequently develop ovarian tumors. MCM9(-/-) females also are sterile as ovaries are completely devoid of oocytes. In contrast, MCM9(-/-) testes produce spermatozoa, albeit in much reduced quantity. Mcm8(-/-) and Mcm9(-/-) embryonic fibroblasts show growth defects and chromosomal damage and cannot overcome a transient inhibition of replication fork progression. In these cells, chromatin recruitment of HR factors like Rad51 and RPA is impaired and HR strongly reduced. We further demonstrate that MCM8 and MCM9 form a complex and that they coregulate their stability. Our work uncovers essential functions of MCM8 and MCM9 in HR-mediated DSB repair during gametogenesis, replication fork maintenance, and DNA repair.     

Selection and Characterization of Phage-Resistant Mutant Strains of Listeria monocytogenes Reveal Host Genes Linked to Phage Adsorption

Listeria-infecting phages are readily isolated from Listeria-containing environments, yet little is known about the selective forces they exert on their host. Here, we identified that two virulent phages, LP-048 and LP-125, adsorb to the surface of Listeria monocytogenes strain 10403S through different mechanisms. We isolated and sequenced, using whole-genome sequencing, 69 spontaneous mutant strains of 10403S that were resistant to either one or both phages. Mutations from 56 phage-resistant mutant strains with only a single mutation mapped to 10 genes representing five loci on the 10403S chromosome. An additional 12 mutant strains showed two mutations, and one mutant strain showed three mutations. Two of the loci, containing seven of the genes, accumulated the majority (n = 64) of the mutations. A representative mutant strain for each of the 10 genes was shown to resist phage infection through mechanisms of adsorption inhibition. Complementation of mutant strains with the associated wild-type allele was able to rescue phage susceptibility for 6 out of the 10 representative mutant strains. Wheat germ agglutinin, which specifically binds to N-acetylglucosamine, bound to 10403S and mutant strains resistant to LP-048 but did not bind to mutant strains resistant to only LP-125. We conclude that mutant strains resistant to only LP-125 lack terminal N-acetylglucosamine in their wall teichoic acid (WTA), whereas mutant strains resistant to both phages have disruptive mutations in their rhamnose biosynthesis operon but still possess N-acetylglucosamine in their WTA.     

Extreme Recombination Frequencies Shape Genome Variation and Evolution in the Honeybee,Apis mellifera

Meiotic recombination is a fundamental cellular process, with important consequences for evolution and genome integrity. However, we know little about how recombination rates vary across the genomes of most species and the molecular and evolutionary determinants of this variation. The honeybee, Apis mellifera, has extremely high rates of meiotic recombination, although the evolutionary causes and consequences of this are unclear. Here we use patterns of linkage disequilibrium in whole genome resequencing data from 30 diploid honeybees to construct a fine-scale map of rates of crossing over in the genome. We find that, in contrast to vertebrate genomes, the recombination landscape is not strongly punctate. Crossover rates strongly correlate with levels of genetic variation, but not divergence, which indicates a pervasive impact of selection on the genome. Germ-line methylated genes have reduced crossover rate, which could indicate a role of methylation in suppressing recombination. Controlling for the effects of methylation, we do not infer a strong association between gene expression patterns and recombination. The site frequency spectrum is strongly skewed from neutral expectations in honeybees: rare variants are dominated by AT-biased mutations, whereas GC-biased mutations are found at higher frequencies, indicative of a major influence of GC-biased gene conversion (gBGC), which we infer to generate an allele fixation bias 5 – 50 times the genomic average estimated in humans. We uncover further evidence that this repair bias specifically affects transitions and favours fixation of CpG sites. Recombination, via gBGC, therefore appears to have profound consequences on genome evolution in honeybees and interferes with the process of natural selection. These findings have important implications for our understanding of the forces driving molecular evolution.     

Morphology and Molecules Reveal Unexpected Cryptic Diversity in the Enigmatic Genus Sinobirma Bryk, 1944 (Lepidoptera: Saturniidae)

The wild silkmoth genus Sinobirma Bryk, 1944 is a poorly known monotypic taxon from the eastern end of the Himalaya Range. It was convincingly proposed to be closely related to some members of an exclusively Afro-tropical group of Saturniidae, but its biogeographical and evolutionary history remains enigmatic. After examining recently collected material from Tibet, northern India, and northeastern Myanmar, we realized that this unique species, S. malaisei Bryk, 1944 only known so far from a few specimens and from a very restricted area near the border between north-eastern Myanmar and the Yunnan province of China, may in fact belong to a group of closely related cryptic species. In this work, we combined morphological comparative study, DNA barcoding, and the sequences of a nuclear marker (D2 expansion segment of the 28S rRNA gene) to unequivocally delimit three distinct species in the genus Sinobirma, of which two are described as new to science: S. myanmarensis sp. n. and S. bouyeri sp. n. An informative DNA barcode sequence was obtained from the female holotype of S. malaisei—collected in 1934—ensuring the proper assignation of this name to the newly collected and studied specimens. Our findings represent another example of the potential of coupling traditional taxonomy and DNA barcoding for revealing and solving difficult cases of cryptic diversity. This approach is now being generalized to the world fauna of Saturniidae, with the participation of most of the taxonomists studying these moths.     

Global Spectrum of Copy Number Variations Reveals Genome Organizational Plasticity and Proposes New Migration Routes

Global spectrum of CNVs is required to catalog variations to provide a high-resolution on the dynamics of genome-organization and human migration. In this study, we performed genome-wide genotyping using high-resolution arrays and identified 44,109 CNVs from 1,715 genomes across 12 populations. The study unraveled the force of independent evolutionary dynamics on genome-organizational plasticity across populations. We demonstrated the use of CNV tool to study human migration and identified a second major settlement establishing new migration routes in addition to existing ones.     

Incorporating 16S Gene Copy Number Information Improves Estimates of Microbial Diversity and Abundance

The abundance of different SSU rRNA (“16S”) gene sequences in environmental samples is widely used in studies of microbial ecology as a measure of microbial community structure and diversity. However, the genomic copy number of the 16S gene varies greatly – from one in many species to up to 15 in some bacteria and to hundreds in some microbial eukaryotes. As a result of this variation the relative abundance of 16S genes in environmental samples can be attributed both to variation in the relative abundance of different organisms, and to variation in genomic 16S copy number among those organisms. Despite this fact, many studies assume that the abundance of 16S gene sequences is a surrogate measure of the relative abundance of the organisms containing those sequences. Here we present a method that uses data on sequences and genomic copy number of 16S genes along with phylogenetic placement and ancestral state estimation to estimate organismal abundances from environmental DNA sequence data. We use theory and simulations to demonstrate that 16S genomic copy number can be accurately estimated from the short reads typically obtained from high-throughput environmental sequencing of the 16S gene, and that organismal abundances in microbial communities are more strongly correlated with estimated abundances obtained from our method than with gene abundances. We re-analyze several published empirical data sets and demonstrate that the use of gene abundance versus estimated organismal abundance can lead to different inferences about community diversity and structure and the identity of the dominant taxa in microbial communities. Our approach will allow microbial ecologists to make more accurate inferences about microbial diversity and abundance based on 16S sequence data.