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.     

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.     

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.     

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.     

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.     

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.     

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.     

Copy Number Loss of the Interferon Gene Cluster in Melanomas Is Linked to Reduced T Cell Infiltrate and Poor Patient Prognosis

While immunotherapies are rapidly becoming mainstays of cancer treatment, significant gaps remain in our understanding of how to optimally target them, alone or in combination. Here we describe a novel method to monitor levels of immune cells and pathways in expression data from solid tumors using pre-defined groups or modules of co-regulated immune genes. We show that expression of an interconnected sub-network of type I interferon-stimulated genes (ISGs) in melanomas at the time of diagnosis significantly predicted patient survival, as did, to a lesser extent, sub-networks of T helper/T regulatory and NK/T Cytotoxic cell genes. As a group, poor prognosis tumors with reduced ISG and immune gene levels exhibited significant copy number loss of the interferon gene cluster located at chromosome 9p21.3. Our studies demonstrate a link between type I interferon action and immune cell levels in melanomas, and suggest that therapeutic approaches augmenting both activities may be most beneficial.     

水稻逆转录转座子RIRE10拷贝数与LTR区转录活性的鉴定

在水稻第四号染色体的长臂上鉴定了一个结构完整的Ty3型逆转录转座子RIRE10。RIRE10两LTR间的中间区域在gag pol的上游还包含另一个开放阅读框。通过RT PCR与Northern印迹杂交检测到来自LTR区的转录产物 ;根据点杂交结果 ,鉴定出包含中间区域的RIRE10成员的个数以及LTR区的拷贝数。除了 6 5个完整的逆转录转座子所具备的两个LTR外 ,水稻基因组还含有近 90 0个RIRE10的solo LTR。LTR区的转录以及导致solo LTR产生的同源重组可能影响了RIRE10成员在水稻基因组中的转座活性