Search

The adoption of agriculture triggered a rapid shift towards starch-rich diets in human populations1. Amylase genes facilitate starch digestion, and increased amylase copy number has been observed in some modern human populations with high-starch intake2, although evidence of recent selection is lacking3,4. Here, using 94 long-read haplotype-resolved assemblies and short-read data from approximately 5,600 contemporary and ancient humans, we resolve the diversity and evolutionary history of structural variation at the amylase locus. We find that amylase genes have higher copy numbers in agricultural populations than in fishing, hunting and pastoral populations. We identify 28 distinct amylase structural architectures and demonstrate that nearly identical structures have arisen recurrently on different haplotype backgrounds throughout recent human history. AMY1 and AMY2A genes each underwent multiple duplication/deletion events with mutation rates up to more than 10,000-fold the single-nucleotide polymorphism mutation rate, whereas AMY2B gene duplications share a single origin. Using a pangenome-based approach, we infer structural haplotypes across thousands of humans identifying extensively duplicated haplotypes at higher frequency in modern agricultural populations. Leveraging 533 ancient human genomes, we find that duplication-containing haplotypes (with more gene copies than the ancestral haplotype) have rapidly increased in frequency over the past 12,000 years in West Eurasians, suggestive of positive selection. Together, our study highlights the potential effects of the agricultural revolution on human genomes and the importance of structural variation in human adaptation.

Statin drugs lower blood cholesterol levels for cardiovascular disease prevention. Women are more likely than men to experience adverse statin effects, particularly new-onset diabetes (NOD) and muscle weakness. Here we find that impaired glucose homeostasis and muscle weakness in statin-treated female mice are associated with reduced levels of the omega-3 fatty acid, docosahexaenoic acid (DHA), impaired redox tone, and reduced mitochondrial respiration. Statin adverse effects are prevented in females by administering fish oil as a source of DHA, by reducing dosage of the X chromosome or the Kdm5c gene, which escapes X chromosome inactivation and is normally expressed at higher levels in females than males. As seen in female mice, we find that women experience more severe reductions than men in DHA levels after statin administration, and that DHA levels are inversely correlated with glucose levels. Furthermore, induced pluripotent stem cells from women who developed NOD exhibit impaired mitochondrial function when treated with statin, whereas cells from men do not. These studies identify X chromosome dosage as a genetic risk factor for statin adverse effects and suggest DHA supplementation as a preventive co-therapy.

Although the effects of genetic and environmental perturbations on multicellular organisms are rarely restricted to single phenotypic layers, our current understanding of how developmental programs react to these challenges remains limited. Here, we have examined the phenotypic consequences of disturbing the bicoid regulatory network in early Drosophila embryos. We generated flies with two extra copies of bicoid, which causes a posterior shift of the networks regulatory outputs and a decrease in fitness. We subjected these flies to EMS mutagenesis, followed by experimental evolution. After only 8-15 generations, experimental populations have normalized patterns of gene expression and increased survival. Using a phenomics approach, we find that populations were normalized through rapid increases in embryo size driven by maternal changes in metabolism and ovariole development. We extend our results to additional populations of flies, demonstrating predictability. Together, our results necessitate a broader view of regulatory network evolution at the systems level.

The 22q11.2 locus contains genes critical for brain development. Reciprocal Copy Number Variations (CNVs) at this locus impact risk for neurodevelopmental and psychiatric disorders. Both 22q11.2 deletions (22qDel) and duplications (22qDup) are associated with autism, but 22qDel uniquely elevates schizophrenia risk. Understanding brain phenotypes associated with these highly penetrant CNVs can provide insights into genetic pathways underlying neuropsychiatric disorders. Human neuroimaging and animal models indicate subcortical brain alterations in 22qDel, yet little is known about developmental differences across specific nuclei between reciprocal 22q11.2 CNV carriers and typically developing (TD) controls. We conducted a longitudinal MRI study in a total of 385 scans from 22qDel (n = 96, scans = 191, 53.1% female), 22qDup (n = 37, scans = 64, 45.9% female), and TD controls (n = 80, scans = 130, 51.2% female), across a wide age range (5.5-49.5 years). Volumes of the thalamus, hippocampus, amygdala, and anatomical subregions were estimated using FreeSurfer, and the linear effects of 22q11.2 gene dosage and non-linear effects of age were characterized with generalized additive mixed models (GAMMs). Positive gene dosage effects (volume increasing with copy number) were observed for total intracranial and whole hippocampus volumes, but not whole thalamus or amygdala volumes. Several amygdala subregions exhibited similar positive effects, with bi-directional effects found across thalamic nuclei. Distinct age-related trajectories were observed across the three groups. Notably, both 22qDel and 22qDup carriers exhibited flattened development of hippocampal CA2/3 subfields relative to TD controls. This study provides novel insights into the impact of 22q11.2 CNVs on subcortical brain structures and their developmental trajectories.

Loss of Cdx2 in vivo leads to stunted development of the allantois, an extraembryonic mesoderm-derived structure critical for nutrient delivery and waste removal in the early embryo. Here, we investigate how CDX2 dose-dependently influences the gene regulatory network underlying extraembryonic mesoderm development. By engineering human induced pluripotent stem cells (hiPSCs) consisting of wild-type (WT), heterozygous (CDX2-Het), and homozygous null CDX2 (CDX2-KO) genotypes, differentiating these cells in a 2D gastruloid model, and subjecting these cells to single-nucleus RNA and ATAC sequencing, we identify several pathways that are dose-dependently regulated by CDX2 including VEGF and non-canonical WNT. snATAC-seq reveals that CDX2-Het cells retain a WT-like chromatin accessibility profile, suggesting accessibility alone is not sufficient to drive this variability in gene expression. Because the loss of CDX2 or TBXT phenocopy one another in vivo, we compared differentially expressed genes in our CDX2-KO to those from TBXT-KO hiPSCs differentiated in an analogous experiment. This comparison identifies several communally misregulated genes that are critical for cytoskeletal integrity and tissue permeability. Together, these results clarify how CDX2 dose-dependently regulates gene expression in the extraembryonic mesoderm and reveal pathways that may underlie the defects in vascular development and allantoic elongation seen in vivo.

DYRK1A syndrome results from a reduction in copy number of the DYRK1A gene, which resides on human chromosome 21 (Hsa21). DYRK1A has been implicated in the development of cognitive phenotypes associated with many genetic disorders, including Down syndrome (DS) and Alzheimer's disease (AD). Additionally, overexpression of DYRK1A in DS has been implicated in the development of abnormal skeletal phenotypes in these individuals. Analyses of mouse models with Dyrk1a dosage imbalance (overexpression and underexpression) show skeletal deficits and abnormalities. Normalization of Dyrk1a copy number in an otherwise trisomic animal rescues some skeletal health parameters, and reduction of Dyrk1a copy number in an otherwise euploid (control) animal results in altered skeletal health measurements, including reduced bone mineral density (BMD) in the femur, mandible, and skull. However, little research has been conducted thus far on the implications of DYRK1A reduction on human skeletal health, specifically in individuals with DYRK1A syndrome. This review highlights the skeletal phenotypes of individuals with DYRK1A syndrome, as well as in murine models with reduced Dyrk1a copy number, and provides potential pathways altered by a reduction of DYRK1A copy number, which may impact skeletal health and phenotypes in these individuals. Understanding how decreased expression of DYRK1A in individuals with DYRK1A syndrome impacts bone health may increase awareness of skeletal traits and assist in the development of therapies to improve quality of life for these individuals. [ABSTRACT FROM AUTHOR]

The role of structurally dynamic genomic regions in speciation is poorly understood due to challenges inherent in diploid genome assembly. Here we reconstructed the evolutionary dynamics of structural variation in five cat species by phasing the genomes of three interspecies F1 hybrids to generate near-gapless single-haplotype assemblies. We discerned that cat genomes have a paucity of segmental duplications relative to great apes, explaining their remarkable karyotypic stability. X chromosomes were hotspots of structural variation, including enrichment with inversions in a large recombination desert with characteristics of a supergene. The X-linked macrosatellite DXZ4 evolves more rapidly than 99.5% of the genome clarifying its role in felid hybrid incompatibility. Resolved sensory gene repertoires revealed functional copy number changes associated with ecomorphological adaptations, sociality and domestication. This study highlights the value of gapless genomes to reveal structural mechanisms underpinning karyotypic evolution, reproductive isolation and ecological niche adaptation.

Loss-of-function variants in the triggering receptor expressed on myeloid cells 2 (TREM2) are responsible for a spectrum of neurodegenerative disorders. In the homozygous state, they cause severe pathologies with early onset dementia, such as Nasu-Hakola disease and behavioural variants of frontotemporal dementia (FTD), whereas heterozygous variants increase the risk of late-onset Alzheimer's disease (AD) and FTD. For over half of TREM2 variants found in families with recessive early onset dementia, the defect occurs at the transcript level via premature termination codons or aberrant splicing. The remaining variants are missense alterations thought to affect the protein; however, the underlying pathogenic mechanism is less clear. In this work, we tested whether these disease-associated TREM2 variants contribute to the pathology via altered splicing. Variants scored by SpliceAI algorithm were tested by a full-size TREM2 splicing reporter assay in different cell lines. The effect of variants was quantified by qRT-/RT-PCR and western blots. Nanostring nCounter was used to measure TREM2 RNA in the brains of NHD patients who carried spliceogenic variants. Exon skipping events were analysed from brain RNA-Seq datasets available through the Accelerating Medicines Partnership for Alzheimer's Disease Consortium. We found that for some Nasu-Hakola disease and early onset FTD-causing variants, splicing defects were the primary cause (D134G) or likely contributor to pathogenicity (V126G and K186N). Similar but milder effects on splicing of exons 2 and 3 were demonstrated for A130V, L133L and R136W enriched in patients with dementia. Moreover, the two most frequent missense variants associated with AD/FTD risk in European and African ancestries (R62H, 1% in Caucasians and T96K, 12% in Africans) had splicing defects via excessive skipping of exon 2 and overproduction of a potentially antagonistic TREM2 protein isoform. The effect of R62H on exon 2 skipping was confirmed in three independent brain RNA-Seq datasets. Our findings revealed an unanticipated complexity of pathogenic variation in TREM2 , in which effects on post-transcriptional gene regulation and protein function often coexist. This necessitates the inclusion of computational and experimental analyses of splicing and mRNA processing for a better understanding of genetic variation in disease. [ABSTRACT FROM AUTHOR]

The Gene Balance Hypothesis (GBH) proposes that selection acts on the dosage (i.e. copy number) of genes within dosage-sensitive portions of networks, pathways, and protein complexes to maintain balanced stoichiometry of interacting proteins, because perturbations to stoichiometric balance can result in reduced fitness. This selection has been called dosage balance selection. Dosage balance selection is also hypothesized to constrain expression responses to dosage changes, making dosage-sensitive genes (those encoding members of interacting proteins) experience more similar expression changes. In allopolyploids, where whole-genome duplication involves hybridization of diverged lineages, organisms often experience homoeologous exchanges (HEs) that recombine, duplicate, and delete homoeologous regions of the genome and alter the expression of homoeologous gene pairs. Although the GBH makes predictions about the expression response to HEs, they have not been empirically tested. We used genomic and transcriptomic data from six resynthesized, isogenic Brassica napus lines over ten generations to identify HEs, analyzed expression responses, and tested for patterns of genomic imbalance. Groups of dosage-sensitive genes had less variable expression responses to HEs than dosage-insensitive genes, a sign that their relative dosage is constrained. This difference was absent for homoeologous pairs whose expression was biased toward the BnA subgenome. Finally, the expression response to HEs was more variable than the response to WGD, suggesting HEs create genomic imbalance. These findings expand our knowledge of the impact of dosage balance selection on genome evolution and potentially connect patterns in polyploid genomes over time; from homoeolog expression bias to duplicate gene retention.

The loss of a single chromosome in a diploid organism halves the dosage of many genes and is usually accompanied by a substantial decrease in fitness. We asked whether this decrease simply reflects the joint damage caused by individual gene dosage deficiencies. We measured the fitness effects of single heterozygous gene deletions in yeast and combined them for each chromosome. This predicted a negative growth rate, that is, lethality, for multiple monosomies. However, monosomic strains remained alive and grew as if much (often most) of the damage caused by single mutations had disappeared, revealing an exceptionally large and positive epistatic component of fitness. We looked for functional explanations by analyzing the transcriptomes. There was no evidence of increased (compensatory) gene expression on the monosomic chromosomes. Nor were there signs of the cellular stress response that would be expected if monosomy led to protein destabilization and thus cytotoxicity. Instead, all monosomic strains showed extensive upregulation of genes encoding ribosomal proteins, but in an indiscriminate manner that did not correspond to their altered dosage. This response did not restore the stoichiometry required for efficient biosynthesis, which probably became growth limiting, making all other mutation-induced metabolic defects much less important. In general, the modular structure of the cell leads to an effective fragmentation of the total mutational load. Defects outside the module(s) currently defining fitness lose at least some of their relevance, producing the epiphenomenon of positive interactions between individually negative effects.

DYRK1A syndrome results from a reduction in copy number of the DYRK1A gene, which resides on human chromosome 21 (Hsa21). DYRK1A has been implicated in the development of cognitive phenotypes associated with many genetic disorders, including Down syndrome (DS) and Alzheimer’s disease (AD). Additionally, overexpression of DYRK1A in DS has been implicated in the development of abnormal skeletal phenotypes in these individuals. Analyses of mouse models with Dyrk1a dosage imbalance (overexpression and underexpression) show skeletal deficits and abnormalities. Normalization of Dyrk1a copy number in an otherwise trisomic animal rescues some skeletal health parameters, and reduction of Dyrk1a copy number in an otherwise euploid (control) animal results in altered skeletal health measurements, including reduced bone mineral density (BMD) in the femur, mandible, and skull. However, little research has been conducted thus far on the implications of DYRK1A reduction on human skeletal health, specifically in individuals with DYRK1A syndrome. This review highlights the skeletal phenotypes of individuals with DYRK1A syndrome, as well as in murine models with reduced Dyrk1a copy number, and provides potential pathways altered by a reduction of DYRK1A copy number, which may impact skeletal health and phenotypes in these individuals. Understanding how decreased expression of DYRK1A in individuals with DYRK1A syndrome impacts bone health may increase awareness of skeletal traits and assist in the development of therapies to improve quality of life for these individuals.

IntroductionPeople with Down syndrome (DS) are predisposed to Alzheimer's disease (AD). The amyloid hypothesis informs studies of AD. In AD-DS, but not sporadic AD, increased APP copy number is necessary, defining the APP gene dose hypothesis. Which amyloid precursor protein (APP) products contribute needs to be determined.MethodsBrain levels of full-length protein (fl-hAPP), C-terminal fragments (hCTFs), and amyloid beta (Aβ) peptides were measured in DS, AD-DS, non-demented controls (ND), and sporadic AD cases. The APP gene-dose hypothesis was evaluated in the Dp16 model.ResultsDS and AD-DS differed from ND and AD for all APP products. In AD-DS, Aβ42 and Aβ40 levels exceeded AD. APP products were increased in the Dp16 model; increased APP gene dose was necessary for loss of vulnerable neurons, tau pathology, and activation of astrocytes and microglia.DiscussionIncreases in APP products other than Aβ distinguished AD-DS from AD. Deciphering AD-DS pathogenesis necessitates deciphering which APP products contribute and how.

A first step in studying the molecular nature of dominant mutations is distinguishing between loss- and gain-of-function mutations. Loss-of-function mutations in haploinsufficient loci are dominant to their wild-type alleles. Haploinsufficient genes typically encode components of multimeric complexes or proteins that participate in highly connected nodes in regulatory or signaling pathways. The expression levels of haploinsufficient genes in yeast and plants are tightly controlled, and their reduced expression or overexpression affects fitness. Haploinsufficiency in flowering plants is not always easy to notice, because these plants contain large gene families whose members are functionally redundant. Raising awareness of this phenomenon among plant researchers will increase the detection of haploinsufficiency in plants. In diploid organisms, haploinsufficiency can be defined as the requirement for more than one fully functional copy of a gene. In contrast to most genes, whose loss-of-function alleles are recessive, loss-of-function alleles of haploinsufficient genes are dominant. However, forward and reverse genetic screens are biased toward obtaining recessive, loss-of-function mutations, and therefore, dominant mutations of all types are underrepresented in mutant collections. Despite this underrepresentation, haploinsufficient loci have intriguing implications for studies of genome evolution, gene dosage, stability of protein complexes, genetic redundancy, and gene expression. Here we review examples of haploinsufficiency in flowering plants and describe the underlying molecular mechanisms and evolutionary forces driving haploinsufficiency. Finally, we discuss the masking of haploinsufficiency by genetic redundancy, a widespread phenomenon among angiosperms. [ABSTRACT FROM AUTHOR]

In recent years, nanobodies have played an increasingly crucial role in virus neutralization, ELISA detection, and medical imaging. This study aimed to explore a universal expression strategy in Pichia pastoris using three nanobodies, denoted Va, Vb, and Vc, as model proteins. Initially, plasmids pLD-AOXα and pLD-AOX were engineered to minimize the risk of antibiotic resistance gene drift. Optimization of promoters and signal peptides resulted in a 1.38-fold and 1.89-fold increase in Va production. Further optimization of gene dosage led to an additional 1.39-fold enhancement in Va yield. Subsequently, 25 molecular chaperones were co-expressed with Va under the control of the wild-type AOX1 promoter, with HAC1 further increasing Va yield by 1.5-fold. By fine-tuning the promoter strength for HAC1, Va production was increased by 2.41-fold under the control of the 55p promoter. Finally, through high-density fermentation, the Va yield reached 2.13 g/L, representing a 49.8-fold increase compared to the initial strain 1-AOXα-Va in shake-flask culture. Integration of pLD-55p-HAC1 into the GS115 genome resulted in the H55 host, and the transformation of multicopy plasmids into this host led to a 1.98-fold increase in Vb yield and a 2.34-fold increase in Vc yield, respectively. The engineering of antibiotic-free parental plasmids, modification of expression components, gene dosage optimization, and the H55 host are regarded as a composite strategy which will pave the way for efficient expression of nanobodies in the future. [ABSTRACT FROM AUTHOR]

Nutrient scarcity is pervasive for natural microbial communities, affecting species reproduction and co-existence. However, it remains unclear whether there are general rules of how microbial species abundances are shaped by biotic and abiotic factors. Here we show that the ribosomal RNA gene operon (rrn) copy number, a genomic trait related to bacterial growth rate and nutrient demand, decreases from the abundant to the rare biosphere in the nutrient-rich coastal sediment but exhibits the opposite pattern in the nutrient-scarce pelagic zone of the global ocean. Both patterns are underlain by positive correlations between community-level rrn copy number and nutrients. Furthermore, inter-species co-exclusion inferred by negative network associations is observed more in coastal sediment than in ocean water samples. Nutrient manipulation experiments yield effects of nutrient availability on rrn copy numbers and network associations that are consistent with our field observations. Based on these results, we propose a "hunger games" hypothesis to define microbial species abundance rules using the rrn copy number, ecological interaction, and nutrient availability.

The development of multiple chromosome-scale reference genome sequences in many taxonomic groups has yielded a high-resolution view of the patterns and processes of molecular evolution. Nonetheless, leveraging information across multiple genomes remains a significant challenge in nearly all eukaryotic systems. These challenges range from studying the evolution of chromosome structure, to finding candidate genes for quantitative trait loci, to testing hypotheses about speciation and adaptation. Here, we present GENESPACE, which addresses these challenges by integrating conserved gene order and orthology to define the expected physical position of all genes across multiple genomes. We demonstrate this utility by dissecting presence-absence, copy-number, and structural variation at three levels of biological organization: spanning 300 million years of vertebrate sex chromosome evolution, across the diversity of the Poaceae (grass) plant family, and among 26 maize cultivars. The methods to build and visualize syntenic orthology in the GENESPACE R package offer a significant addition to existing gene family and synteny programs, especially in polyploid, outbred, and other complex genomes.

The last edition of the X-chromosome inactivation (XCI) meeting was held as an EMBO workshop in Berlin on 19-22 June 2023. The conference took place at the Harnack-haus in the Dahlem district, birthplace of the first modern research campus, where notable scientists such as Lise Meitner, Hans Krebs and, briefly, Albert Einstein conducted their research. This special edition, also accessible online, was organized by Rafael Galupa (Centre for Integrative Biology of Toulouse, France), Joost Gribnau (Erasmus MC Rotterdam, The Netherlands), Claire Rougeulle (Universite'Paris Cite'/CNRS, Epigenetics and Cell Fate Center, Paris, France), Edda Schulz (Max Planck Institute for Molecular Genetics, Berlin, Germany) and James Turner (The Francis Crick Institute, London, UK). Originally scheduled for 2021, to commemorate the 60th anniversary of Mary Lyon's hypothesis on X-chromosome inactivation in mammals and the 30th anniversary of XIST/Xist discovery, the meeting had to be postponed because of the COVID-19 pandemic. Seven years after the latest XCI meeting in London, the enthusiasm and expectations of the community were at their highest, bringing together over 160 scientists from around the world to share and discuss their research. Eighty posters and more than 40 talks were presented at this event, in a collegial and collaborative atmosphere.A historical session and several breakout discussions were also organized, as well as the now traditional boat trip, all thanks to great organization. Here, we debrief readers on this fantastic conference. [ABSTRACT FROM AUTHOR]

The gene balance hypothesis proposes that selection acts on the dosage (i.e. copy number) of genes within dosage-sensitive portions of networks, pathways, and protein complexes to maintain balanced stoichiometry of interacting proteins, because perturbations to stoichiometric balance can result in reduced fitness. This selection has been called dosage balance selection. Dosage balance selection is also hypothesized to constrain expression responses to dosage changes, making dosage-sensitive genes (those encoding members of interacting proteins) experience more similar expression changes. In allopolyploids, where whole-genome duplication involves hybridization of diverged lineages, organisms often experience homoeologous exchanges that recombine, duplicate, and delete homoeologous regions of the genome and alter the expression of homoeologous gene pairs. Although the gene balance hypothesis makes predictions about the expression response to homoeologous exchanges, they have not been empirically tested. We used genomic and transcriptomic data from 6 resynthesized, isogenic Brassica napus lines over 10 generations to identify homoeologous exchanges, analyzed expression responses, and tested for patterns of genomic imbalance. Groups of dosage-sensitive genes had less variable expression responses to homoeologous exchanges than dosage-insensitive genes, a sign that their relative dosage is constrained. This difference was absent for homoeologous pairs whose expression was biased toward the B. napus A subgenome. Finally, the expression response to homoeologous exchanges was more variable than the response to whole-genome duplication, suggesting homoeologous exchanges create genomic imbalance. These findings expand our knowledge of the impact of dosage balance selection on genome evolution and potentially connect patterns in polyploid genomes over time, from homoeolog expression bias to duplicate gene retention. [ABSTRACT FROM AUTHOR]

Sequence variants in gene regulatory regions alter gene expression and contribute to phenotypes of individual cells and the whole organism, including disease susceptibility and progression. Single-nucleotide variants in enhancers or promoters may affect gene transcription by altering transcription factor binding sites. Differential transcription factor binding in heterozygous genomic loci provides a natural source of information on such regulatory variants. We present a novel approach to call the allele-specific transcription factor binding events at single-nucleotide variants in ChIP-Seq data, taking into account the joint contribution of aneuploidy and local copy number variation, that is estimated directly from variant calls. We have conducted a meta-analysis of more than 7 thousand ChIP-Seq experiments and assembled the database of allele-specific binding events listing more than half a million entries at nearly 270 thousand single-nucleotide polymorphisms for several hundred human transcription factors and cell types. These polymorphisms are enriched for associations with phenotypes of medical relevance and often overlap eQTLs, making candidates for causality by linking variants with molecular mechanisms. Specifically, there is a special class of switching sites, where different transcription factors preferably bind alternative alleles, thus revealing allele-specific rewiring of molecular circuitry.

Haploinsufficiency of transcriptional regulators causes human congenital heart disease (CHD); however, the underlying CHD gene regulatory network (GRN) imbalances are unknown. Here, we define transcriptional consequences of reduced dosage of the CHD transcription factor, TBX5, in individual cells during cardiomyocyte differentiation from human induced pluripotent stem cells (iPSCs). We discovered highly sensitive dysregulation of TBX5-dependent pathways-including lineage decisions and genes associated with heart development, cardiomyocyte function, and CHD genetics-in discrete subpopulations of cardiomyocytes. Spatial transcriptomic mapping revealed chamber-restricted expression for many TBX5-sensitive transcripts. GRN analysis indicated that cardiac network stability, including vulnerable CHD-linked nodes, is sensitive to TBX5 dosage. A GRN-predicted genetic interaction between Tbx5 and Mef2c, manifesting as ventricular septation defects, was validated in mice. These results demonstrate exquisite and diverse sensitivity to TBX5 dosage in heterogeneous subsets of iPSC-derived cardiomyocytes and predicts candidate GRNs for human CHDs, with implications for quantitative transcriptional regulation in disease.

The DEFECTIVE EMBRYO AND MERISTEMS 1 (DEM1) gene encodes a protein of unknown biochemical function required for meristem formation and seedling development in tomato, but it was unclear whether DEM1's primary role was in cell division or alternatively, in defining the identity of meristematic cells. Genome sequence analysis indicates that flowering plants possess at least two DEM genes. Arabidopsis has two DEM genes, DEM1 and DEM2, which we show are expressed in developing embryos and meristems in a punctate pattern that is typical of genes involved in cell division. Homozygous dem1 dem2 double mutants were not recovered, and plants carrying a single functional DEM1 allele and no functional copies of DEM2, i.e. DEM1/dem1 dem2/dem2 plants, exhibit normal development through to the time of flowering but during male reproductive development, chromosomes fail to align on the metaphase plate at meiosis II and result in abnormal numbers of daughter cells following meiosis. Additionally, these plants show defects in both pollen and embryo sac development, and produce defective male and female gametes. In contrast, dem1/dem1 DEM2/dem2 plants showed normal levels of fertility, indicating that DEM2 plays a more important role than DEM1 in gamete viability. The increased importance of DEM2 in gamete viability correlated with higher mRNA levels of DEM2 compared to DEM1 in most tissues examined and particularly in the vegetative shoot apex, developing siliques, pollen and sperm. We also demonstrate that gamete viability depends not only on the number of functional DEM alleles inherited following meiosis, but also on the number of functional DEM alleles in the parent plant that undergoes meiosis. Furthermore, DEM1 interacts with RAS-RELATED NUCLEAR PROTEIN 1 (RAN1) in yeast two-hybrid and pull-down binding assays, and we show that fluorescent proteins fused to DEM1 and RAN1 co-localize transiently during male meiosis and pollen development. In eukaryotes, RAN is a highly conserved GTPase that plays key roles in cell cycle progression, spindle assembly during cell division, reformation of the nuclear envelope following cell division, and nucleocytoplasmic transport. Our results demonstrate that DEM proteins play an essential role in cell division in plants, most likely through an interaction with RAN1.

Tuberculosis (TB) is the largest infectious disease with 10 million new active-TB patients and1.7 million deaths per year. Active-TB is an inflammatory disease and is increasingly viewed as an imbalance of immune responses to M. tb. infection. The mechanisms of a switch from latent infection to active disease is not well worked out but a shift in the immune responses is thought to be responsible. Increasingly, the role of gut microbiota has been described as a major influencer of the immune system. And because the gut is the largest immune organ, we aimed to analyze the gut microbiome in active-TB patients in a TB-endemic country, Pakistan. The study revealed that Ruminococcacea, Enetrobactericeae, Erysipelotrichaceae, Bifidobacterium, etc. were the major genera associated with active-TB, also associated with chronic inflammatory disease. Plasma antibody profiles against several M. tb. antigens, as specific biomarkers for active-TB, correlated closely with the patient gut microbial profiles. Besides, bcoA gene copy number, indicative of the level of butyrate production by the gut microbiome was five-fold lower in TB patients compared to healthy individuals. These findings suggest that gut health in TB patients is compromised, with implications for disease morbidity (e.g., severe weight loss) as well as immune impairment.

BackgroundThe circadian clock drives endogenous 24-h rhythms that allow organisms to adapt and prepare for predictable and repeated changes in their environment throughout the day-night (diurnal) cycle. Many components of the circadian clock in Arabidopsis thaliana have been functionally characterized, but comparatively little is known about circadian clocks in grass species including major crops like maize and sorghum.ResultsComparative research based on protein homology and diurnal gene expression patterns suggests the function of some predicted clock components in grasses is conserved with their Arabidopsis counterparts, while others have diverged in function. Our analysis of diurnal gene expression in three panicoid grasses sorghum, maize, and foxtail millet revealed conserved and divergent evolution of expression for core circadian clock genes and for the overall transcriptome. We find that several classes of core circadian clock genes in these grasses differ in copy number compared to Arabidopsis, but mostly exhibit conservation of both protein sequence and diurnal expression pattern with the notable exception of maize paralogous genes. We predict conserved cis-regulatory motifs shared between maize, sorghum, and foxtail millet through identification of diurnal co-expression clusters for a subset of 27,196 orthologous syntenic genes. In this analysis, a Cochran-Mantel-Haenszel based method to control for background variation identified significant enrichment for both expected and novel 6-8 nucleotide motifs in the promoter regions of genes with shared diurnal regulation predicted to function in common physiological activities.ConclusionsThis study illustrates the divergence and conservation of circadian clocks and diurnal regulatory networks across syntenic orthologous genes in panacoid grass species. Further, conserved local regulatory sequences contribute to the architecture of these diurnal regulatory networks that produce conserved patterns of diurnal gene expression.

BackgroundIn plants, each ribosomal protein (RP) is encoded by a small gene family but it is largely unknown whether the family members are functionally diversified. There are two RPL23a paralogous genes (RPL23aA and RPL23aB) encoding cytoplasmic ribosomal proteins in Arabidopsis thaliana. Knock-down of RPL23aA using RNAi impeded growth and led to morphological abnormalities, whereas knock-out of RPL23aB had no observable phenotype, thus these two RPL23a paralogous proteins have been used as examples of ribosomal protein paralogues with functional divergence in many published papers.ResultsIn this study, we characterized T-DNA insertion mutants of RPL23aA and RPL23aB. A rare non-allelic non-complementation phenomenon was found in the F1 progeny of the rpl23aa X rpl23ab cross, which revealed a dosage effect of these two genes. Both RPL23aA and RPL23aB were found to be expressed almost in all examined tissues as revealed by GUS reporter analysis. Expression of RPL23aB driven by the RPL23aA promoter can rescue the phenotype of rpl23aa, indicating these two proteins are actually equivalent in function. Interestingly, based on the publicly available RNA-seq data, we found that these two RPL23a paralogues were expressed in a concerted manner and the expression level of RPL23aA was much higher than that of RPL23aB at different developmental stages and in different tissues.ConclusionsOur findings suggest that the two RPL23a paralogous proteins are functionally equivalent but the two genes are not. RPL23aA plays a predominant role due to its higher expression levels. RPL23aB plays a lesser role due to its lower expression. The presence of paralogous genes for the RPL23a protein in plants might be necessary to maintain its adequate dosage.

Bacteria that form long-term intracellular associations with host cells lose many genes, a process that often results in tiny, gene-dense, and stable genomes. Paradoxically, the some of the same evolutionary processes that drive genome reduction and simplification may also cause genome expansion and complexification. A bacterial endosymbiont of cicadas, Hodgkinia cicadicola , exemplifies this paradox. In many cicada species, a single Hodgkinia lineage with a tiny, gene-dense genome has split into several interdependent cell and genome lineages. Each new Hodgkinia lineage encodes a unique subset of the ancestral unsplit genome in a complementary way, such that the collective gene contents of all lineages match the total found in the ancestral single genome. This splitting creates genetically distinct Hodgkinia cells that must function together to carry out basic cellular processes. It also creates a gene dosage problem where some genes are encoded by only a small fraction of cells while others are much more abundant. Here, by sequencing DNA and RNA of Hodgkinia from different cicada species with different amounts of splitting—along with its structurally stable, unsplit partner endosymbiont Sulcia muelleri— we show that Hodgkinia does not transcriptionally compensate to rescue the wildly unbalanced gene and genome ratios that result from lineage splitting. We also find that Hodgkinia has a reduced capacity for basic transcriptional control independent of the splitting process. Our findings reveal another layer of degeneration further pushing the limits of canonical molecular and cell biology in Hodgkinia and may partially explain its propensity to go extinct through symbiont replacement. [ABSTRACT FROM AUTHOR]

Protein production in any expression system can be enhanced either by optimizing the culturing parameters or targeting the genetic factors that enhance protein production. One of those genetic factors is gene dosage, which can be increased by increasing the vector copy number. However, increased number of expression vector poses some metabolic burdens on bacterial cells. The current study provides gene multimerization strategy as an alternate way to amplify the gene dosage and explores its effect on the expression of human epidermal growth factor (hEGF) by constructing three types of expression vectors i.e., pET28-EGF-C1, pET28-EGF-C2 and pET28-EGF-C3, each containing single, double and triple expression cassettes, respectively. These vectors were transformed in E. coli strain Rosetta-gami 2 (DE3) to develop three different bacterial populations. Expression of human epidermal growth factor was analyzed by SDS-PAGE and subsequent densitometric analysis was done using ImageJ software. The difference in protein expression of all three populations was analyzed by using one-way ANOVA. Clones with double and triple copies of genes showed an increased expression up to 3.8 and 3.1 folds as compare to the clones having single copy of gene. Gene multimerization strategy can be used to enhance production of recombinant proteins in bacterial expression system. [ABSTRACT FROM AUTHOR]

Males and females differ in body composition and fat distribution. Using a mouse model that segregates gonadal sex (ovaries and testes) from chromosomal sex (XX and XY), we showed that XX chromosome complement in combination with a high-fat diet led to enhanced weight gain in the presence of male or female gonads. We identified the genomic dosage of Kdm5c, an X chromosome gene that escapes X chromosome inactivation, as a determinant of the X chromosome effect on adiposity. Modulating Kdm5c gene dosage in XX female mice to levels that are normally present in males resulted in reduced body weight, fat content, and food intake to a degree similar to that seen with altering the entire X chromosome dosage. In cultured preadipocytes, the levels of KDM5C histone demethylase influenced chromatin accessibility (ATAC-Seq), gene expression (RNA-Seq), and adipocyte differentiation. Both in vitro and in vivo, Kdm5c dosage influenced gene expression involved in extracellular matrix remodeling, which is critical for adipocyte differentiation and adipose tissue expansion. In humans, adipose tissue KDM5C mRNA levels and KDM5C genetic variants were associated with body mass. These studies demonstrate that the sex-dependent dosage of Kdm5c contributes to male/female differences in adipocyte biology and highlight X-escape genes as a critical component of female physiology.

Non-coding mutations at the far end of a large gene desert surrounding the SOX9 gene result in a human craniofacial disorder called Pierre Robin sequence (PRS). Leveraging a human stem cell differentiation model, we identify two clusters of enhancers within the PRS-associated region that regulate SOX9 expression during a restricted window of facial progenitor development at distances up to 1.45 Mb. Enhancers within the 1.45 Mb cluster exhibit highly synergistic activity that is dependent on the Coordinator motif. Using mouse models, we demonstrate that PRS phenotypic specificity arises from the convergence of two mechanisms: confinement of Sox9 dosage perturbation to developing facial structures through context-specific enhancer activity and heightened sensitivity of the lower jaw to Sox9 expression reduction. Overall, we characterize the longest-range human enhancers involved in congenital malformations, directly demonstrate that PRS is an enhanceropathy, and illustrate how small changes in gene expression can lead to morphological variation.

Marek's disease (MD) is a neoplastic disease of chickens caused by Marek's disease virus (MDV), a member of the subfamily Alphaherpesvirinae Like other alphaherpesviruses, MDV encodes a serine/threonine protein kinase, US3. The functions of US3 have been extensively studied in other alphaherpesviruses; however, the biological functions of MDV US3 and its substrates have not been studied in detail. In this study, we investigated potential cellular pathways that are regulated by MDV US3 and identified chicken CREB (chCREB) as a substrate of MDV US3. We show that wild-type MDV US3, but not kinase-dead US3 (US3-K220A), increases CREB phosphorylation, leading to recruitment of phospho-CREB (pCREB) to the promoter of the CREB-responsive gene and activation of CREB target gene expression. Using US3 deletion and US3 kinase-dead recombinant MDV, we identified US3-responsive MDV genes during infection and found that the majority of US3-responsive genes were located in the MDV repeat regions. Chromatin immunoprecipitation sequencing (ChIP-seq) studies determined that some US3-regulated genes colocalized with Meq (an MDV-encoded oncoprotein) recruitment sites. Chromatin immunoprecipitation-PCR (ChIP-PCR) further confirmed Meq binding to the ICP4/LAT region, which is also regulated by US3. Furthermore, biochemical studies demonstrated that MDV US3 interacts with Meq in transfected cells and MDV-infected chicken embryonic fibroblasts in a phosphorylation-dependent manner. Finally, in vitro kinase studies revealed that Meq is a US3 substrate. MDV US3 thus acts as an upstream kinase of the CREB signaling pathway to regulate the transcription function of the CREB/Meq heterodimer, which targets cellular and viral gene expression.IMPORTANCE MDV is a potent oncogenic herpesvirus that induces T-cell lymphoma in infected chickens. Marek's disease continues to have a significant economic impact on the poultry industry worldwide. US3 encoded by alphaherpesviruses is a multifunctional kinase involved in the regulation of various cellular pathways. Using an MDV genome quantitative reverse transcriptase PCR (qRT-PCR) array and chromatin immunoprecipitation, we elucidated the role of MDV US3 in viral and cellular gene regulation. Our results provide insights into how viral kinase regulates host cell signaling pathways to activate both viral and host gene expression. This is an important step toward understanding host-pathogen interaction through activation of signaling cascades.

Genomic instability contributes to tumorigenesis through the amplification and deletion of cancer driver genes. DNA copy number (CN) profiling of ensembles of tumors allows a thermodynamic analysis of the profile for each tumor. The free energy of the distribution of CNs is found to be a monotonically increasing function of the average chromosomal ploidy. The dependence is universal across several cancer types. Surprisal analysis distinguishes two main known subgroups: tumors with cells that have or have not undergone whole-genome duplication (WGD). The analysis uncovers that CN states having a narrower distribution are energetically more favorable toward the WGD transition. Surprisal analysis also determines the deviations from a fully stable-state distribution. These deviations reflect constraints imposed by tumor fitness selection pressures. The results point to CN changes that are more common in high-ploidy tumors and thus support altered selection pressures upon WGD.

A key goal of whole-genome sequencing for studies of human genetics is to interrogate all forms of variation, including single-nucleotide variants, small insertion or deletion (indel) variants and structural variants. However, tools and resources for the study of structural variants have lagged behind those for smaller variants. Here we used a scalable pipeline1 to map and characterize structural variants in 17,795 deeply sequenced human genomes. We publicly release site-frequency data to create the largest, to our knowledge, whole-genome-sequencing-based structural variant resource so far. On average, individuals carry 2.9 rare structural variants that alter coding regions; these variants affect the dosage or structure of 4.2 genes and account for 4.0-11.2% of rare high-impact coding alleles. Using a computational model, we estimate that structural variants account for 17.2% of rare alleles genome-wide, with predicted deleterious effects that are equivalent to loss-of-function coding alleles; approximately 90% of such structural variants are noncoding deletions (mean 19.1 per genome). We report 158,991 ultra-rare structural variants and show that 2% of individuals carry ultra-rare megabase-scale structural variants, nearly half of which are balanced or complex rearrangements. Finally, we infer the dosage sensitivity of genes and noncoding elements, and reveal trends that relate to element class and conservation. This work will help to guide the analysis and interpretation of structural variants in the era of whole-genome sequencing.

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO2 detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

Humanin is a member of a new family of peptides that are encoded by short open reading frames within the mitochondrial genome. It is conserved in animals and is both neuroprotective and cytoprotective. Here we report that in C. elegans the overexpression of humanin is sufficient to increase lifespan, dependent on daf-16/Foxo. Humanin transgenic mice have many phenotypes that overlap with the worm phenotypes and, similar to exogenous humanin treatment, have increased protection against toxic insults. Treating middle-aged mice twice weekly with the potent humanin analogue HNG, humanin improves metabolic healthspan parameters and reduces inflammatory markers. In multiple species, humanin levels generally decline with age, but here we show that levels are surprisingly stable in the naked mole-rat, a model of negligible senescence. Furthermore, in children of centenarians, who are more likely to become centenarians themselves, circulating humanin levels are much greater than age-matched control subjects. Further linking humanin to healthspan, we observe that humanin levels are decreased in human diseases such as Alzheimer's disease and MELAS (Mitochondrial Encephalopathy, Lactic Acidosis, and Stroke-like episodes). Together, these studies are the first to demonstrate that humanin is linked to improved healthspan and increased lifespan.

The invasive behavior of glioblastoma is essential to its aggressive potential. Here, we show that pleckstrin homology domain interacting protein (PHIP), acting through effects on the force transduction layer of the focal adhesion complex, drives glioblastoma motility and invasion. Immunofluorescence analysis localized PHIP to the leading edge of glioblastoma cells, together with several focal adhesion proteins: vinculin (VCL), talin 1 (TLN1), integrin beta 1 (ITGB1), as well as phosphorylated forms of paxillin (pPXN) and focal adhesion kinase (pFAK). Confocal microscopy specifically localized PHIP to the force transduction layer, together with TLN1 and VCL. Immunoprecipitation revealed a physical interaction between PHIP and VCL. Targeted suppression of PHIP resulted in significant down-regulation of these focal adhesion proteins, along with zyxin (ZYX), and produced profoundly disorganized stress fibers. Live-cell imaging of glioblastoma cells overexpressing a ZYX-GFP construct demonstrated a role for PHIP in regulating focal adhesion dynamics. PHIP silencing significantly suppressed the migratory and invasive capacity of glioblastoma cells, partially restored following TLN1 or ZYX cDNA overexpression. PHIP knockdown produced substantial suppression of tumor growth upon intracranial implantation, as well as significantly reduced microvessel density and secreted VEGF levels. PHIP copy number was elevated in the classical glioblastoma subtype and correlated with elevated EGFR levels. These results demonstrate PHIP's role in regulating the actin cytoskeleton, focal adhesion dynamics, and tumor cell motility, and identify PHIP as a key driver of glioblastoma migration and invasion.

Sex chromosome dosage (SCD) variation increases risk for neuropsychiatric impairment, which may reflect direct SCD effects on brain organization. Here, we 1) map cumulative X- and Y-chromosome dosage effects on regional cortical thickness (CT) and investigate potential functional implications of these effects using Neurosynth, 2) test if this map is organized by patterns of CT covariance that are evident in health, and 3) characterize SCD effects on CT covariance itself. We modeled SCD effects on CT and CT covariance for 308 equally sized regions of the cortical sheet using structural neuroimaging data from 301 individuals with varying numbers of sex chromosomes (169 euploid, 132 aneuploid). Mounting SCD increased CT in the rostral frontal cortex and decreased CT in the lateral temporal cortex, bilaterally. Regions targeted by SCD were associated with social functioning, language processing, and comprehension. Cortical regions with a similar degree of SCD-sensitivity showed heightened CT covariance in health. Finally, greater SCD also increased covariance among regions similarly affected by SCD. Our study both 1) develops novel methods for comparing typical and disease-related structural covariance networks in the brain and 2) uses these techniques to resolve and identify organizing principles for SCD effects on regional cortical anatomy and anatomical covariance.

Cancer develops through a process of somatic evolution1,2. Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes3. Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)4, we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

To date, over 1000 melanocytic neoplasms, spanning all stages of tumorigenesis, have been sequenced, offering detailed views into their -omic landscapes. This has coincided with advances in genetic engineering technologies that allow molecular biologists to edit the human genome with extreme precision and new mouse models to simulate disease progression. In this review, we describe how these technologies are being harnessed to provide insights into the evolution of melanoma at an unprecedented resolution, revealing that prior models of melanoma evolution, in which pathways are turned 'on' or 'off' in a binary fashion during the run-up to melanoma, are oversimplified.

Humans and orcas are among the very rare species that have a prolonged post-reproductive lifespan (PRLS), during which the aging process continues. Reactive oxygen species (ROS) derived from mitochondria and from the NADPH oxidase (NOX) enzymes of innate immune cells are known to contribute to aging, with the former thought to be dominant. CD33-related-Siglecs are immune receptors that recognize self-associated-molecular-patterns and modulate NOX-derived-ROS. We herewith demonstrate a strong correlation of lifespan with CD33rSIGLEC gene number in 26 species, independent of body weight or phylogeny. The correlation is stronger when considering total CD33rSIGLEC gene number rather than those encoding inhibitory and activating subsets, suggesting that lifetime balancing of ROS is important. Combining independent lines of evidence including the short half-life and spontaneous activation of neutrophils, we calculate that even without inter-current inflammation, a major source of lifetime ROS exposure may actually be neutrophil NOX-derived. However, genomes of human supercentenarians (>110 years) do not harbor a significantly higher number of functional CD33rSIGLEC genes. Instead, lifespan correlation with CD33rSIGLEC gene number was markedly strengthened by excluding the post-reproductive lifespan of humans and orcas (R2  = 0.83; P

Multi-omics datasets represent distinct aspects of the central dogma of molecular biology. Such high-dimensional molecular profiles pose challenges to data interpretation and hypothesis generation. ActivePathways is an integrative method that discovers significantly enriched pathways across multiple datasets using statistical data fusion, rationalizes contributing evidence and highlights associated genes. As part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2658 cancers across 38 tumor types, we integrated genes with coding and non-coding mutations and revealed frequently mutated pathways and additional cancer genes with infrequent mutations. We also analyzed prognostic molecular pathways by integrating genomic and transcriptomic features of 1780 breast cancers and highlighted associations with immune response and anti-apoptotic signaling. Integration of ChIP-seq and RNA-seq data for master regulators of the Hippo pathway across normal human tissues identified processes of tissue regeneration and stem cell regulation. ActivePathways is a versatile method that improves systems-level understanding of cellular organization in health and disease through integration of multiple molecular datasets and pathway annotations.

Mitochondrial DNA copy number (mtDNA-CN), a measure of the number of mitochondrial genomes per cell, is a minimally invasive proxy measure for mitochondrial function and has been associated with several aging-related diseases. Although quantitative real-time PCR (qPCR) is the current gold standard method for measuring mtDNA-CN, mtDNA-CN can also be measured from genotyping microarray probe intensities and DNA sequencing read counts. To conduct a comprehensive examination on the performance of these methods, we use known mtDNA-CN correlates (age, sex, white blood cell count, Duffy locus genotype, incident cardiovascular disease) to evaluate mtDNA-CN calculated from qPCR, two microarray platforms, as well as whole genome (WGS) and whole exome sequence (WES) data across 1,085 participants from the Atherosclerosis Risk in Communities (ARIC) study and 3,489 participants from the Multi-Ethnic Study of Atherosclerosis (MESA). We observe mtDNA-CN derived from WGS data is significantly more associated with known correlates compared to all other methods (p < 0.001). Additionally, mtDNA-CN measured from WGS is on average more significantly associated with traits by 5.6 orders of magnitude and has effect size estimates 5.8 times more extreme than the current gold standard of qPCR. We further investigated the role of DNA extraction method on mtDNA-CN estimate reproducibility and found mtDNA-CN estimated from cell lysate is significantly less variable than traditional phenol-chloroform-isoamyl alcohol (p = 5.44x10-4) and silica-based column selection (p = 2.82x10-7). In conclusion, we recommend the field moves towards more accurate methods for mtDNA-CN, as well as re-analyze trait associations as more WGS data becomes available from larger initiatives such as TOPMed.

Tumor DNA sequencing data can be interpreted by computational methods that analyze genomic heterogeneity to infer evolutionary dynamics. A growing number of studies have used these approaches to link cancer evolution with clinical progression and response to therapy. Although the inference of tumor phylogenies is rapidly becoming standard practice in cancer genome analyses, standards for evaluating them are lacking. To address this need, we systematically assess methods for reconstructing tumor subclonality. First, we elucidate the main algorithmic problems in subclonal reconstruction and develop quantitative metrics for evaluating them. Then we simulate realistic tumor genomes that harbor all known clonal and subclonal mutation types and processes. Finally, we benchmark 580 tumor reconstructions, varying tumor read depth, tumor type and somatic variant detection. Our analysis provides a baseline for the establishment of gold-standard methods to analyze tumor heterogeneity.

Sequences encoding Olduvai protein domains (formerly DUF1220) show the greatest human lineage-specific increase in copy number of any coding region in the genome and have been associated, in a dosage-dependent manner, with brain size, cognitive aptitude, autism, and schizophrenia. Tandem intragenic duplications of a three-domain block, termed the Olduvai triplet, in four NBPF genes in the chromosomal 1q21.1-0.2 region, are primarily responsible for the striking human-specific copy number increase. Interestingly, most of the Olduvai triplets are adjacent to, and transcriptionally coregulated with, three human-specific NOTCH2NL genes that have been shown to promote cortical neurogenesis. Until now, the underlying genomic events that drove the Olduvai hyperamplification in humans have remained unexplained. Here, we show that the presence or absence of an alternative first exon of the Olduvai triplet perfectly discriminates between amplified (58/58) and unamplified (0/12) triplets. We provide sequence and breakpoint analyses that suggest the alternative exon was produced by an nonallelic homologous recombination-based mechanism involving the duplicative transposition of an existing Olduvai exon found in the CON3 domain, which typically occurs at the C-terminal end of NBPF genes. We also provide suggestive in vitro evidence that the alternative exon may promote instability through a putative G-quadraplex (pG4)-based mechanism. Lastly, we use single-molecule optical mapping to characterize the intragenic structural variation observed in NBPF genes in 154 unrelated individuals and 52 related individuals from 16 families and show that the presence of pG4-containing Olduvai triplets is strongly correlated with high levels of Olduvai copy number variation. These results suggest that the same driver of genomic instability that allowed the evolutionarily recent, rapid, and extreme human-specific Olduvai expansion remains highly active in the human genome.

To examine the established link between DNA replication and telomere length, we tested whether firing of telomeric origins would cause telomere lengthening. We found that RIF1 mutants that block Protein Phosphatase 1 (PP1) binding activated telomeric origins but did not elongate telomeres. In a second approach, we found overexpression of ∆N-Dbf4 and Cdc7 increased DDK activity and activated telomeric origins, yet telomere length was unchanged. We tested a third mechanism to activate origins using the sld3-A mcm5-bob1 mutant that de-regulates the pre-replication complex, and again saw no change in telomere length. Finally, we tested whether mutations in RIF1 that cause telomere elongation would affect origin firing. We found that neither rif1-∆1322 nor rif1HOOK affected firing of telomeric origins. We conclude that telomeric origin firing does not cause telomere elongation, and the role of Rif1 in regulating origin firing is separable from its role in regulating telomere length.

BackgroundThe Hemiptera (aphids, cicadas, and true bugs) are a key insect order, with high diversity for feeding ecology and excellent experimental tractability for molecular genetics. Building upon recent sequencing of hemipteran pests such as phloem-feeding aphids and blood-feeding bed bugs, we present the genome sequence and comparative analyses centered on the milkweed bug Oncopeltus fasciatus, a seed feeder of the family Lygaeidae.ResultsThe 926-Mb Oncopeltus genome is well represented by the current assembly and official gene set. We use our genomic and RNA-seq data not only to characterize the protein-coding gene repertoire and perform isoform-specific RNAi, but also to elucidate patterns of molecular evolution and physiology. We find ongoing, lineage-specific expansion and diversification of repressive C2H2 zinc finger proteins. The discovery of intron gain and turnover specific to the Hemiptera also prompted the evaluation of lineage and genome size as predictors of gene structure evolution. Furthermore, we identify enzymatic gains and losses that correlate with feeding biology, particularly for reductions associated with derived, fluid nutrition feeding.ConclusionsWith the milkweed bug, we now have a critical mass of sequenced species for a hemimetabolous insect order and close outgroup to the Holometabola, substantially improving the diversity of insect genomics. We thereby define commonalities among the Hemiptera and delve into how hemipteran genomes reflect distinct feeding ecologies. Given Oncopeltus's strength as an experimental model, these new sequence resources bolster the foundation for molecular research and highlight technical considerations for the analysis of medium-sized invertebrate genomes.

SHOX deficiency is a common genetic cause of short stature of variable degree. SHOX haploinsufficiency causes Leri–Weill dyschondrosteosis (LWD) as well as nonspecific short stature. SHOX haploinsufficiency is known to result from heterozygous loss-of-function variants with pseudo-autosomal dominant inheritance, while biallelic SHOX loss-of-function variants cause the more severe skeletal dysplasia, Langer mesomelic dyschondrosteosis (LMD). Here we report for the first time the pseudo-autosomal recessive inheritance of LWD in two siblings caused by a novel homozygous non-canonical, leaky splice-site variant in intron 3 of SHOX: c.544+5G>C. Transcript analyses in patient-derived fibroblasts showed homozygous patients to produce approximately equal amounts of normally spliced mRNA and mRNA with the abnormal retention of intron 3 and containing a premature stop codon (p.Val183Glyfs*31). The aberrant transcript was shown to undergo nonsense-mediated mRNA decay, and thus resulting in SHOX haploinsufficiency in the homozygous patient. Six healthy relatives who are of normal height are heterozygous for this variant and fibroblasts from a heterozygote for the c.544+5G>C variant produced wild-type transcript amounts comparable to healthy control. The unique situation reported here highlights the fact that the dosage of SHOX determines the clinical phenotype rather than the Mendelian inheritance pattern of SHOX variants. This study extends the molecular and inheritance spectrum of SHOX deficiency disorder and highlights the importance of functional testing of SHOX variants of unknown significance in order to allow appropriate counseling and precision medicine for each family individual. [ABSTRACT FROM AUTHOR]

In recent years, nanobodies have played an increasingly crucial role in virus neutralization, ELISA detection, and medical imaging. This study aimed to explore a universal expression strategy in Pichia pastoris using three nanobodies, denoted Va, Vb, and Vc, as model proteins. Initially, plasmids pLD-AOXα and pLD-AOX were engineered to minimize the risk of antibiotic resistance gene drift. Optimization of promoters and signal peptides resulted in a 1.38-fold and 1.89-fold increase in Va production. Further optimization of gene dosage led to an additional 1.39-fold enhancement in Va yield. Subsequently, 25 molecular chaperones were co-expressed with Va under the control of the wild-type AOX1 promoter, with HAC1 further increasing Va yield by 1.5-fold. By fine-tuning the promoter strength for HAC1, Va production was increased by 2.41-fold under the control of the 55p promoter. Finally, through high-density fermentation, the Va yield reached 2.13 g/L, representing a 49.8-fold increase compared to the initial strain 1-AOXα-Va in shake-flask culture. Integration of pLD-55p-HAC1 into the GS115 genome resulted in the H55 host, and the transformation of multicopy plasmids into this host led to a 1.98-fold increase in Vb yield and a 2.34-fold increase in Vc yield, respectively. The engineering of antibiotic-free parental plasmids, modification of expression components, gene dosage optimization, and the H55 host are regarded as a composite strategy which will pave the way for efficient expression of nanobodies in the future.

Notholaenids are an unusual group of ferns that have adapted to, and diversified within, the deserts of Mexico and the southwestern United States. With approximately 40 species, this group is noted for being desiccation-tolerant and having "farina"-powdery exudates of lipophilic flavonoid aglycones-that occur on both the gametophytic and sporophytic phases of their life cycle. The most recent circumscription of notholaenids based on plastid markers surprisingly suggests that several morphological characters, including the expression of farina, are homoplasious. In a striking case of convergence, Notholaena standleyi appears to be distantly related to core Notholaena, with several taxa not before associated with Notholaena nested between them. Such conflicts can be due to morphological homoplasy resulting from adaptive convergence or, alternatively, the plastid phylogeny itself might be misleading, diverging from the true species tree due to incomplete lineage sorting, hybridization, or other factors. In this study, we present a species phylogeny for notholaenid ferns, using four low-copy nuclear loci and concatenated data from three plastid loci. A total of 61 individuals (49 notholaenids and 12 outgroup taxa) were sampled, including 31 out of 37 recognized notholaenid species. The homeologous/allelic nuclear sequences were retrieved using PacBio sequencing and the PURC bioinformatics pipeline. Each dataset was first analyzed individually using maximum likelihood and Bayesian inference, and the species phylogeny was inferred using *BEAST. Although we observed several incongruences between the nuclear and plastid phylogenies, our principal results are broadly congruent with previous inferences based on plastid data. By mapping the presence of farina and their biochemical constitutions on our consensus phylogenetic tree, we confirmed that the characters are indeed homoplastic and have complex evolutionary histories. Hybridization among recognized species of the notholaenid clade appears to be relatively rare compared to that observed in other well-studied fern genera.

Gene dosage variation and the associated changes in gene expression influence a wide variety of traits, ranging from cancer in humans to yield in plants. It is also expected to affect important traits of ecological and agronomic importance in forest trees, but this variation has not been systematically characterized or exploited. Here we performed a comprehensive scan of the Populus genome for dosage-sensitive loci affecting quantitative trait variation for spring and fall phenology and biomass production. The study population was a large collection of clonally propagated F1 hybrid lines of Populus that saturate the genome 10-fold with deletions and insertions (indels) of known sizes and positions. As a group, the phenotypic means of the indel lines consistently differed from control nonindel lines, with an overall negative effect of both insertions and deletions on all biomass-related traits but more diverse effects and an overall wider phenotypic distribution of the indel lines for the phenology-related traits. We also investigated the correlation between gene dosage at specific chromosomal locations and phenotype, to identify dosage quantitative trait loci (dQTL). Such dQTL were detected for most phenotypes examined, but stronger effect dQTL were identified for the phenology-related traits than for the biomass traits. Our genome-wide screen for dosage sensitivity in a higher eukaryote demonstrates the importance of global genomic balance and the impact of dosage on life history traits.