Gene copy number variations in adaptive evolution: The genomic distribution of gene copy number variations revealed by genetic mapping and their adaptive role in an undomesticated species,white spruce (Picea glauca)

Gene copy number variation (CNV) has been associated with phenotypic variability in animals and plants, but a genomewide understanding of their impacts on phenotypes is largely restricted to human and agricultural systems. As such, CNVs have rarely been considered in investigations of the genomic architecture of adaptation in wild species. Here, we report on the genetic mapping of gene CNVs in white spruce, which lacks a contiguous assembly of its large genome (~20 Gb), and their relationships with adaptive phenotypic variation. We detected 3,911 gene CNVs including de novo structural variations using comparative genome hybridization on arrays (aCGH) in a large progeny set. We inferred the heterozygosity at CNV loci within parents by comparing haploid and diploid tissues and genetically mapped 82 gene CNVs. Our analysis showed that CNVs were distributed over 10 linkage groups and identified four CNV hotspots that we predict to occur in other species of the Pinaceae. Significant relationships were found between 29 of the gene CNVs and adaptive traits based on regression analyses with timings of bud set and bud flush, and height growth, suggesting a role for CNVs in climate adaptation. The importance of CNVs in adaptive evolution of white spruce was also indicated by functional gene annotations and the clustering of 31% of the mapped adaptive gene CNVs in CNV hotspots. Taken together, these results illustrate the feasibility of studying CNVs in undomesticated species and represent a major step towards a better understanding of the roles of CNVs in adaptive evolution.     

Gene copy number variations as signatures of adaptive evolution in the parthenogenetic,plant‐parasitic nematode Meloidogyne incognita

Adaptation to changing environmental conditions represents a challenge to parthenogenetic organisms, and until now, how phenotypic variants are generated in clones in response to the selection pressure of their environment remains poorly known. The obligatory parthenogenetic root‐knot nematode species Meloidogyne incognita has a worldwide distribution and is the most devastating plant‐parasitic nematode. Despite its asexual reproduction, this species exhibits an unexpected capacity of adaptation to environmental constraints, for example, resistant hosts. Here, we used a genomewide comparative hybridization strategy to evaluate variations in gene copy numbers between genotypes of M. incognita resulting from two parallel experimental evolution assays on a susceptible vs. resistant host plant. We detected gene copy number variations (CNVs) associated with the ability of the nematodes to overcome resistance of the host plant, and this genetic variation may reflect an adaptive response to host resistance in this parthenogenetic species. The CNV distribution throughout the nematode genome is not random and suggests the occurrence of genomic regions more prone to undergo duplications and losses in response to the selection pressure of the host resistance. Furthermore, our analysis revealed an outstanding level of gene loss events in nematode genotypes that have overcome the resistance. Overall, our results support the view that gene loss could be a common class of adaptive genetic mechanism in response to a challenging new biotic environment in clonal animals.     

Liebenberg syndrome is caused by a deletion upstream to the PITX1 gene resulting in transformation of the upper limbs to reflect lower limb characteristics

Liebenberg syndrome (MIM 186550) is a very rare autosomal dominant condition characterized by three main features: dysplasia of all of the bony components of the elbow joint, abnormalities in the shape of carpal bones, and brachydactyly. In this paper, we report a Saudi Arabian family with Liebenberg syndrome. Comparative genomic hybridization (CGH) revealed a 275-kb deletion within the cytogenetic band 5q31.1 which contains the H2AFY gene and 190,428 bp of its downstream region. The deleted region is upstream to the PITX1 gene. The radiological features in the upper limbs of all affected members of the family were almost identical to the phenotype in the mouse model with ectopic expression of Pitx1 in the forelimbs. We therefore re-define the phenotype of Liebenberg syndrome as a transformation of the upper limbs to reflect lower limb characteristics and speculate that the area of deletion contains a regulatory sequence that suppresses the expression of PITX1 in the upper limb buds.     

Compound heterozygosity for TNXB genetic variants in a mixed‐breed dog with Ehlers‐Danlos syndrome

The Ehlers‐Danlos syndromes (EDSs) are a heterogeneous group of inherited connective tissue disorders characterized by skin hyperextensibility, joint hypermobility and tissue fragility. Inherited disorders similar to human EDS have been reported in different mammalian species. In the present study, we investigated a female mixed‐breed dog with clinical signs of EDS. Whole‐genome sequencing of the affected dog revealed two missense variants in the TNXB gene, encoding the extracellular matrix protein tenascin XB. In humans, TNXB genetic variants cause classical‐like EDS or the milder hypermobile EDS. The affected dog was heterozygous at both identified variants. Each variant allele was transmitted from one of the case's parents, consistent with compound heterozygosity. Although one of the variant alleles, XM_003431680.3:c.2012G>A, p.(Ser671Asn), was private to the family of the affected dog and absent from whole‐genome sequencing data of 599 control dogs, the second variant allele, XM_003431680.3:c.2900G>A, p.(Gly967Asp), is present at a low frequency in the Chihuahua and Poodle population. Given that TNXB is a functional candidate gene for EDS, we suggest that compound heterozygosity for the identified TNXB variants may have caused the EDS‐like phenotype in the affected dog. Chihuahuas and Poodles should be monitored for EDS cases, which might confirm the hypothesized pathogenic effect of the segregating TNXB variant.     

Genome‐scale models of metabolism and gene expression extend and refine growth phenotype prediction

Growth is a fundamental process of life. Growth requirements are well‐characterized experimentally for many microbes; however, we lack a unified model for cellular growth. Such a model must be predictive of events at the molecular scale and capable of explaining the high‐level behavior of the cell as a whole. Here, we construct an ME‐Model for Escherichia coli—a genome‐scale model that seamlessly integrates metabolic and gene product expression pathways. The model computes ～80% of the functional proteome (by mass), which is used by the cell to support growth under a given condition. Metabolism and gene expression are interdependent processes that affect and constrain each other. We formalize these constraints and apply the principle of growth optimization to enable the accurate prediction of multi‐scale phenotypes, ranging from coarse‐grained (growth rate, nutrient uptake, by‐product secretion) to fine‐grained (metabolic fluxes, gene expression levels). Our results unify many existing principles developed to describe bacterial growth.     

Testing genotyping strategies for ultra‐deep sequencing of a co‐amplifying gene family: MHC class I in a passerine bird

Characterization of highly duplicated genes, such as genes of the major histocompatibility complex (MHC), where multiple loci often co‐amplify, has until recently been hindered by insufficient read depths per amplicon. Here, we used ultra‐deep Illumina sequencing to resolve genotypes at exon 3 of MHC class I genes in the sedge warbler (Acrocephalus schoenobaenus). We sequenced 24 individuals in two replicates and used this data, as well as a simulated data set, to test the effect of amplicon coverage (range: 500–20 000 reads per amplicon) on the repeatability of genotyping using four different genotyping approaches. A third replicate employed unique barcoding to assess the extent of tag jumping, that is swapping of individual tag identifiers, which may confound genotyping. The reliability of MHC genotyping increased with coverage and approached or exceeded 90% within‐method repeatability of allele calling at coverages of >5000 reads per amplicon. We found generally high agreement between genotyping methods, especially at high coverages. High reliability of the tested genotyping approaches was further supported by our analysis of the simulated data set, although the genotyping approach relying primarily on replication of variants in independent amplicons proved sensitive to repeatable errors. According to the most repeatable genotyping method, the number of co‐amplifying variants per individual ranged from 19 to 42. Tag jumping was detectable, but at such low frequencies that it did not affect the reliability of genotyping. We thus demonstrate that gene families with many co‐amplifying genes can be reliably genotyped using HTS, provided that there is sufficient per amplicon coverage.     

Extreme‐phenotype genome‐wide association study (XP‐GWAS): a method for identifying trait‐associated variants by sequencing pools of individuals selected from a diversity panel

Although approaches for performing genome‐wide association studies (GWAS) are well developed, conventional GWAS requires high‐density genotyping of large numbers of individuals from a diversity panel. Here we report a method for performing GWAS that does not require genotyping of large numbers of individuals. Instead XP‐GWAS (extreme‐phenotype GWAS) relies on genotyping pools of individuals from a diversity panel that have extreme phenotypes. This analysis measures allele frequencies in the extreme pools, enabling discovery of associations between genetic variants and traits of interest. This method was evaluated in maize (Zea mays) using the well‐characterized kernel row number trait, which was selected to enable comparisons between the results of XP‐GWAS and conventional GWAS. An exome‐sequencing strategy was used to focus sequencing resources on genes and their flanking regions. A total of 0.94 million variants were identified and served as evaluation markers; comparisons among pools showed that 145 of these variants were statistically associated with the kernel row number phenotype. These trait‐associated variants were significantly enriched in regions identified by conventional GWAS. XP‐GWAS was able to resolve several linked QTL and detect trait‐associated variants within a single gene under a QTL peak. XP‐GWAS is expected to be particularly valuable for detecting genes or alleles responsible for quantitative variation in species for which extensive genotyping resources are not available, such as wild progenitors of crops, orphan crops, and other poorly characterized species such as those of ecological interest.     

Chromosome 19p13.3 deletion in a child with Peutz-Jeghers syndrome, congenital heart defect, high myopia, learning difficulties and dysmorphic features: Clinical and molecular characterization of a new contiguous gene syndrome

The Peutz-Jeghers syndrome (PJS) is an autosomal-dominant hamartomatous polyposis syndrome characterized by mucocutaneous pigmentation, gastrointestinal polyps and the increased risk of multiple cancers. The causative point mutation in the STK11 gene of most patients accounts for about 30% of the cases of partial and complete gene deletion. This is a report on a girl with PJS features, learning difficulties, dysmorphic features and cardiac malformation, bearing a de novo 1.1 Mb deletion at 19p13.3. This deletion encompasses at least 47 genes, including STK11. This is the first report on 19p13.3 deletion associated with a PJS phenotype, as well as other atypical manifestations, thereby implying a new contiguous gene syndrome.     

A submicroscopic deletion involving part of the CREBBP gene detected by array-CGH in a patient with Rubinstein-Taybi syndrome

We report a girl with Rubinstein-Taybi syndrome (RSTS) who was found to have copy number loss on 16p13.3 by array-CGH. She has developmental delay and other features of RSTS including downslanting palpebral fissures, a prominent nose with the nasal septum extending below the alae nasi, broad thumbs and big toes, postaxial polydactyly of the right foot and constipation from birth. We report the junction sequence across the breakpoint region for a microdeletion in RSTS. The sequencing results also showed that the deletion was 81.4kb involving three genes DNASE 1, TRAP 1, and CREBBP.     

A 1.1 Mb deletion in distal 13q deletion syndrome region with congenital heart defect and postaxial polydactyly: Additional support for a CHD locus at distal 13q34 region

13q deletion syndrome is a rare genetic disorder, especially for group 3 deletion (13q33–q34 deletion). Previously we described a patient with congenital heart defect and mental retardation and proposed that a distal 6 Mb region might contain the causative gene of congenital heart defect. Here we present a new patient with congenital heart defects (CHD), hand and foot anomalies and mild mental retardation. We identified a 1.1 Mb deletion at chromosome 13q34 with high resolution SNP-array BeadChips (HumanOmni1-Quad, Illumina, USA). This chromosome region contains ten annotated genes, including GRK1, TFDP1, RASA3 and GAS6. To our knowledge, this represents the smallest 13q34 deletion identified to date. Our study provides additional support that distal 13q34 deletion region might contain key gene(s) responsible for cardiac development.     

Evolutionary conservation,diversity and specificity of LTR‐retrotransposons in flowering plants: insights from genome‐wide analysis and multi‐specific comparison

The availability of complete or nearly complete genome sequences from several plant species permits detailed discovery and cross‐species comparison of transposable elements (TEs) at the whole genome level. We initially investigated 510 long terminal repeat‐retrotransposon (LTR‐RT) families comprising 32 370 elements in soybean (Glycine max (L.) Merr.). Approximately 87% of these elements were located in recombination‐suppressed pericentromeric regions, where the ratio (1.26) of solo LTRs to intact elements (S/I) is significantly lower than that of chromosome arms (1.62). Further analysis revealed a significant positive correlation between S/I and LTR sizes, indicating that larger LTRs facilitate solo LTR formation. Phylogenetic analysis revealed seven Copia and five Gypsy evolutionary lineages that were present before the divergence of eudicot and monocot species, but the scales and timeframes within which they proliferated vary dramatically across families, lineages and species, and notably, a Copia lineage has been lost in soybean. Analysis of the physical association of LTR‐RTs with centromere satellite repeats identified two putative centromere retrotransposon (CR) families of soybean, which were grouped into the CR (e.g. CRR and CRM) lineage found in grasses, indicating that the ‘functional specification’ of CR pre‐dates the bifurcation of eudicots and monocots. However, a number of families of the CR lineage are not concentrated in centromeres, suggesting that their CR roles may now be defunct. Our data also suggest that the envelope‐like genes in the putative Copia retrovirus‐like family are probably derived from the Gypsy retrovirus‐like lineage, and thus we propose the hypothesis of a single ancient origin of envelope‐like genes in flowering plants.     

Tuning a ménage à trois: Co‐evolution and co‐adaptation of nuclear and organellar genomes in plants

Plastids and mitochondria arose through endosymbiotic acquisition of formerly free‐living bacteria. During more than a billion years of subsequent concerted evolution, the three genomes of plant cells have undergone dramatic structural changes to optimize the expression of the compartmentalized genetic material and to fine‐tune the communication between the nucleus and the organelles. The chimeric composition of many multiprotein complexes in plastids and mitochondria (one part of the subunits being nuclear encoded and another one being encoded in the organellar genome) provides a paradigm for co‐evolution at the cellular level. In this paper, we discuss the co‐evolution of nuclear and organellar genomes in the context of environmental adaptation in species and populations. We highlight emerging genetic model systems and new experimental approaches that are particularly suitable to elucidate the molecular basis of co‐adaptation processes and describe how nuclear‐cytoplasmic co‐evolution can cause genetic incompatibilities that contribute to the establishment of hybridization barriers, ultimately leading to the formation of new species.     

Selection and sex‐biased dispersal in a coastal shark: the influence of philopatry on adaptive variation

Sex‐biased dispersal is expected to homogenize nuclear genetic variation relative to variation in genetic material inherited through the philopatric sex. When site fidelity occurs across a heterogeneous environment, local selective regimes may alter this pattern. We assessed spatial patterns of variation in nuclear‐encoded, single nucleotide polymorphisms (SNPs) and sequences of the mitochondrial control region in bonnethead sharks (Sphyrna tiburo), a species thought to exhibit female philopatry, collected from summer habitats used for gestation. Geographic patterns of mtDNA haplotypes and putatively neutral SNPs confirmed female philopatry and male‐mediated gene flow along the northeastern coast of the Gulf of Mexico. A total of 30 outlier SNP loci were identified; alleles at over half of these loci exhibited signatures of latitude‐associated selection. Our results indicate that in species with sex‐biased dispersal, philopatry can facilitate sorting of locally adaptive variation, with the dispersing sex facilitating movement of potentially adaptive variation among locations and environments.     

A complex microcephaly syndrome in a Pakistani family associated with a novel missense mutation in RBBP8 and a heterozygous deletion in NRXN1

We report on a consanguineous Pakistani family with a severe congenital microcephaly syndrome resembling the Seckel syndrome and Jawad syndrome. The affected individuals in this family were born to consanguineous parents of whom the mother presented with mild intellectual disability (ID), epilepsy and diabetes mellitus. The two living affected brothers presented with microcephaly, white matter disease of the brain, hyponychia, dysmorphic facial features with synophrys, epilepsy, diabetes mellitus and ID. Genotyping with a 250K SNP array in both affected brothers revealed an 18 MB homozygous region on chromosome 18p11.21-q12.1 encompassing the SCKL2 locus of the Seckel and Jawad syndromes. Sequencing of the RBBP8 gene, underlying the Seckel and Jawad syndromes, identified the novel mutation c.919A > G, p.Arg307Gly, segregating in a recessive manner in the family. In addition, in the two affected brothers and their mother we have also found a heterozygous 607 kb deletion, encompassing exons 13–19 of NRXN1. Bidirectional sequencing of the coding exons of NRXN1 did not reveal any other mutation on the other allele. It thus appears that the phenotype of the mildly affected mother can be explained by the NRXN1 deletion, whereas the more severe and complex microcephalic phenotype of the two affected brothers is due to the simultaneous deletion in NRXN1 and the homozygous missense mutation affecting RBBP8.     

Expanding the mutation spectrum for Fraser syndrome: Identification of a novel heterozygous deletion in FRAS1

Fraser syndrome (FS) is a rare autosomal recessive inherited disorder characterized by cryptophthalmos, laryngeal defects and oral clefting, mental retardation, syndactyly, and urogenital defects. To date, 250 patients have been described in the literature. Mutations in the FRAS1 gene on chromosome 4 have been identified in patients with Fraser syndrome. So far, 26 mutations have been identified, most of them are truncating mutations. The mutational spectrum includes nucleotide substitutions, splicing defects, a large insertion, and small deletions/insertions. Moreover, single heterozygous missense mutations in FRAS1 seem to be responsible for non-syndromic unilateral renal agenesis.